4 void ITSPlotAllTracks(){
6 ifstream in ("AliITSTra.out");
10 ///////////////////////////////// Histograms definition ///////////////////////////////////////////
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);
16 // TH1F *hp=new TH1F("hp","PHI resolution",50,-5.,5.); hp->SetFillColor(4);
17 //TH1F *hl=new TH1F("hl","LAMBDA resolution",50,-5.,5.); hl->SetFillColor(4);
19 // TH1F *hpt=new TH1F("hpt","Relative Pt resolution",50,-5.,5.);
20 TH1F *hpt=new TH1F("hpt","Relative Pt resolution",40,-10.,10.);
22 TH1F *hd=new TH1F("hd","Impact parameter distribution ",100,0,2000);
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);
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);
37 TH1F *hgood=new TH1F("hgood","Good tracks",10,0,2);
38 hgood->Sumw2(); // aggiunto il 15-01-01
39 TH1F *hfound=new TH1F("hfound","Found tracks",10,0,2);
40 hfound->Sumw2(); // aggiunto il 15-01-01
41 TH1F *hfake=new TH1F("hfake","Fake tracks",10,0,2);
43 TH1F *hg=new TH1F("hg","",10,0,2); //efficiency for good tracks
44 hg->SetLineColor(4); hg->SetLineWidth(2);
45 TH1F *hf=new TH1F("hf","Efficiency for fake tracks",10,0,2);
46 /*hf->SetFillColor(1); hf->SetFillStyle(3013);*/ hf->SetLineColor(4); hf->SetLineWidth(2);
48 ///////////////////////////////////////////////////////////////////////////////////////////////////
50 ifstream in1 ("AliITSTrag.out");
54 if( in1.eof() ) break;
59 Int_t neglabs=0; // provvisoria
61 for (int r=0; r<9; r++) in>>DataOut(r);
64 Double_t ptg=DataOut(0); Double_t labITS=DataOut(1); Double_t labTPC=DataOut(2); Double_t ptperc=DataOut(3);
65 Double_t deltalam=DataOut(4); Double_t deltaphi=DataOut(5);
66 Double_t Dtot=DataOut(6); Double_t Dr=DataOut(7); Double_t Dz=DataOut(8);
68 if(labITS<0) neglabs++; // provvisoria
69 if(labITS>=0) hfound->Fill(ptg); else { hfake->Fill(ptg);}
83 Stat_t ngood=hgood->GetEntries(); cerr<<"Good tracks "<<ngood<<endl;
84 Stat_t nfound=hfound->GetEntries(); cerr<<"Found tracks "<<nfound<<endl;
85 Stat_t nfake=hfake->GetEntries(); cerr<<"Fake tracks "<<nfake<<endl;
86 gStyle->SetOptStat(111110);
88 TCanvas *c1=new TCanvas("c1","",0,0,700,700);
89 TPad *p1=new TPad("p1","",0,0.5,0.5,1); p1->Draw(); hp->SetXTitle("(mrad)");
90 p1->cd(); hp->Draw(); hp->Fit("gaus"); c1->cd();
91 TPad *p2=new TPad("p2","",0.5,0.5,1,1); p2->Draw(); hl->SetXTitle("(mrad)");
92 p2->cd(); hl->Draw(); hl->Fit("gaus"); c1->cd();
93 TPad *p3=new TPad("p3","",0,0,0.5,0.5); p3->Draw(); hpt->SetXTitle("(%)");
94 p3->cd(); hpt->Draw(); hpt->Fit("gaus"); c1->cd();
95 TPad *p4=new TPad("p4","",0.5,0,1,0.5); p4->Draw(); hd->SetXTitle("(micron)");
96 p4->cd(); hd->Draw(); c1->cd();
98 TCanvas *c3=new TCanvas("c3","",200,200,800,500);
99 hfound->Print("all"); // aggiunto il 16-01-01
100 hgood->Print("all"); // aggiunto il 16-01-01
101 TPad *p7=new TPad("p7","",0,0,0.333,1); p7->Draw(); p7->cd(); hfound->Draw(); c3->cd();
102 TPad *p8=new TPad("p8","",0.333,0,0.666,1); p8->Draw(); p8->cd(); hfake->Draw(); c3->cd();
103 TPad *p9=new TPad("p9","",0.666,0,1,1); p9->Draw(); p9->cd(); hgood->Draw(); c3->cd();
105 TCanvas *c4=new TCanvas("c4","",300,300,800,500);
106 hg->Divide(hfound,hgood,1.,1.); //,"b");
107 hf->Divide(hfake,hgood,1.,1.); //,"b");
109 hg->SetYTitle("Tracking efficiency");
110 hg->SetXTitle("Pt (GeV/c)");
111 hg->Print("all"); // aggiunto il 16-01-01
112 hg->Draw(); // to not plot the erro bars hg->Draw("histo");
113 TLine *line1 = new TLine(0,1.0,2,1.0); line1->SetLineStyle(4);
115 TLine *line2 = new TLine(0,0.9,2,0.9); line2->SetLineStyle(4);
116 line2->Draw("histosame");
120 hf->SetFillStyle(3013);
123 hf->Draw("same"); // to not plot the error bars hf->Draw("histosame");
125 TText *text = new TText(0.461176,0.248448,"Fake tracks");
126 text->SetTextSize(0.05);
128 text = new TText(0.453919,1.11408,"Good tracks");
129 text->SetTextSize(0.05);
132 TCanvas *c2=new TCanvas("c2","",100,100,700,400);
133 TPad *p5=new TPad("p5","",0,0,0.5,1); p5->Draw(); hdr->SetXTitle("(micron)");
134 p5->cd(); hdr->Draw(); hdr->Fit("gaus"); c2->cd();
135 TPad *p6=new TPad("p6","",0.5,0,1,1); p6->Draw(); hdz->SetXTitle("(micron)");
136 p6->cd(); hdz->Draw(); hdz->Fit("gaus"); c2->cd();
138 cout<<"neglabs = "<<neglabs<<"\n"; // provvisoria