#include <math.h>

void pos_plot(){
   gROOT->SetStyle("Plain");
   gStyle->SetOptStat(kFALSE);

   const int num = 1000000;
   double x[num], y[num], z[num], dummy;
   int n = 0;
   int n2 = 0;

   TCanvas *c1 = new TCanvas("c1", "dd", 3, 300, 750, 600);    

//   c1->SetGridy();

   TGraph2D *dt1 = new TGraph2D();
   TGraph2D *dt2 = new TGraph2D();

   ifstream widthfile("output_file1.txt");
   while((widthfile >> x[n] >> y[n] >> z[n] >> dummy >> dummy >> dummy >> dummy >> dummy)){
      if(dummy == 1){		// particle type 1
         dt1->SetPoint(n,x[n],y[n],z[n]);
      }
      if(dummy == 2){		// particle type 2
         dt2->SetPoint(n2,x[n],y[n],z[n]);
         n2++;
      }
      n++;
   }  
   widthfile.close();

   dt1->SetMarkerSize(1);
   dt1->SetMarkerColor(1);
   dt1->SetMarkerColor(4);
   dt1->SetMarkerStyle(1);
   //dt1->GetXaxis()->SetRange(-200,200);
   dt1->SetTitle("");

   dt1->Draw("P");
   dt2->Draw("P,same");
}