ITS/ITSPlotTracks.C

   1 #include <iostream.h>
   2 #include <fstream.h>
   3
   4 void ITSPlotTracks(){
   5
   6   ifstream in ("AliITSTra.out");
   7
   8   TVector DataOut(10);
   9
  10   ///////////////////////////////// Histograms definition ///////////////////////////////////////////
  11
  12
  13   TH1F *hp=new TH1F("hp","PHI resolution",50,-15.,15.); hp->SetFillColor(4);
  14   TH1F *hl=new TH1F("hl","LAMBDA resolution",50,-15.,15.); hl->SetFillColor(4);
  15
  16  //TH1F *hp=new TH1F("hp","PHI resolution",100,-50.,50.); hp->SetFillColor(4);
  17  // TH1F *hl=new TH1F("hl","LAMBDA resolution",100,-50.,50.); hl->SetFillColor(4);
  18
  19  // TH1F *hpt=new TH1F("hpt","Relative Pt resolution",40,-20.,20.);
  20  TH1F *hpt=new TH1F("hpt","Relative Pt resolution",40,-10.,10.);
  21   hpt->SetFillColor(2);
  22   TH1F *hd=new TH1F("hd","Impact parameter distribution ",100,0,2000);
  23   hd->SetFillColor(6);
  24
  25   TH1F *hdr=new TH1F("hdr","Dr ",50,-1000,1000);
  26   hdr->SetFillColor(kGreen);
  27   TH1F *hdz=new TH1F("hdz","Dz ",50,-1000,1000);
  28   hdz->SetFillColor(kBlue);
  29
  30   /*
  31   TH1F *hdr=new TH1F("hdr","Dr ",50,-100,100);
  32   hdr->SetFillColor(kGreen);
  33   TH1F *hdz=new TH1F("hdz","Dz ",50,-1000,1000);
  34   hdz->SetFillColor(kBlue);
  35   */
  36
  37   TH1F *hgood=new TH1F("hgood","Good tracks",10,0,2);
  38   TH1F *hfound=new TH1F("hfound","Found tracks",10,0,2);
  39   TH1F *hfake=new TH1F("hfake","Fake tracks",10,0,2);
  40   TH1F *hg=new TH1F("hg","",10,0,2); //efficiency for good tracks
  41   hg->SetLineColor(4); hg->SetLineWidth(2);
  42   TH1F *hf=new TH1F("hf","Efficiency for fake tracks",10,0,2);
  43   hf->SetFillColor(1); hf->SetFillStyle(3013); hf->SetLineWidth(2);
  44
  45   ///////////////////////////////////////////////////////////////////////////////////////////////////
  46
  47   ifstream in1 ("AliITSTrag.out");
  48   Double_t ptg;
  49   for(;;) {
  50     in1 >> ptg;
  51     if( in1.eof() ) break;
  52     hgood->Fill(ptg);
  53   }
  54   in1.close();
  55
  56   for (;;){
  57     for (int r=0; r<9; r++) in>>DataOut(r);
  58     if( in.eof() ) break;
  59
  60     Double_t ptg=DataOut(0); Double_t labITS=DataOut(1); Double_t labTPC=DataOut(2); Double_t ptperc=DataOut(3);
  61     Double_t deltalam=DataOut(4); Double_t deltaphi=DataOut(5);
  62     Double_t Dtot=DataOut(6); Double_t Dr=DataOut(7); Double_t Dz=DataOut(8);
  63
  64           if(labITS>=0) hfound->Fill(ptg); else         {   hfake->Fill(ptg);}
  65
  66           if(labITS>=0 ) {
  67       hpt->Fill(ptperc);
  68       hl->Fill(deltalam);
  69       hp->Fill(deltaphi);
  70       hd->Fill(Dtot);
  71       hdr->Fill(Dr);
  72       hdz->Fill(Dz);
  73    }
  74   }
  75
  76
  77   in.close();
  78   Stat_t ngood=hgood->GetEntries(); cerr<<"Good tracks "<<ngood<<endl;
  79   Stat_t nfound=hfound->GetEntries(); cerr<<"Found tracks "<<nfound<<endl;
  80   Stat_t nfake=hfake->GetEntries(); cerr<<"Fake tracks "<<nfake<<endl;
  81   gStyle->SetOptStat(111110);
  82   gStyle->SetOptFit(1);
  83   TCanvas *c1=new TCanvas("c1","",0,0,700,700);
  84   TPad *p1=new TPad("p1","",0,0.5,0.5,1); p1->Draw(); hp->SetXTitle("(mrad)");
  85   p1->cd(); hp->Draw();  hp->Fit("gaus"); c1->cd();
  86   TPad *p2=new TPad("p2","",0.5,0.5,1,1); p2->Draw(); hl->SetXTitle("(mrad)");
  87   p2->cd(); hl->Draw(); hl->Fit("gaus"); c1->cd();
  88   TPad *p3=new TPad("p3","",0,0,0.5,0.5); p3->Draw(); hpt->SetXTitle("(%)");
  89   p3->cd(); hpt->Draw(); hpt->Fit("gaus"); c1->cd();
  90   TPad *p4=new TPad("p4","",0.5,0,1,0.5); p4->Draw(); hd->SetXTitle("(micron)");
  91   p4->cd(); hd->Draw(); c1->cd();
  92
  93   TCanvas *c3=new TCanvas("c3","",200,200,800,500);
  94   TPad *p7=new TPad("p7","",0,0,0.333,1); p7->Draw(); p7->cd(); hfound->Draw(); c3->cd();
  95   TPad *p8=new TPad("p8","",0.333,0,0.666,1); p8->Draw(); p8->cd(); hfake->Draw(); c3->cd();
  96   TPad *p9=new TPad("p9","",0.666,0,1,1); p9->Draw(); p9->cd(); hgood->Draw(); c3->cd();
  97
  98   TCanvas *c4=new TCanvas("c4","",300,300,800,500);
  99   hg->Divide(hfound,hgood,1,1.,"b");
 100   hf->Divide(hfake,hgood,1,1.,"b");
 101   hg->SetMaximum(1.4);
 102   hg->SetYTitle("Tracking efficiency");
 103   hg->SetXTitle("Pt (GeV/c)");
 104   hg->Draw();
 105   TLine *line1 = new TLine(0,1.0,2,1.0); line1->SetLineStyle(4);
 106   line1->Draw("same");
 107   TLine *line2 = new TLine(0,0.9,2,0.9); line2->SetLineStyle(4);
 108   line2->Draw("same");
 109
 110   hf->SetFillColor(1);
 111   hf->SetFillStyle(3013);
 112   hf->SetLineColor(2);
 113   hf->SetLineWidth(2);
 114   hf->Draw("histsame");
 115   TText *text = new TText(0.461176,0.248448,"Fake tracks");
 116   text->SetTextSize(0.05);
 117   text->Draw();
 118   text = new TText(0.453919,1.11408,"Good tracks");
 119   text->SetTextSize(0.05);
 120   text->Draw();
 121
 122   TCanvas *c2=new TCanvas("c2","",100,100,700,400);
 123   TPad *p5=new TPad("p5","",0,0,0.5,1); p5->Draw(); hdr->SetXTitle("(micron)");
 124   p5->cd(); hdr->Draw();  hdr->Fit("gaus"); c2->cd();
 125   TPad *p6=new TPad("p6","",0.5,0,1,1); p6->Draw(); hdz->SetXTitle("(micron)");
 126   p6->cd(); hdz->Draw(); hdz->Fit("gaus"); c2->cd();
 127
 128 }