]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGPP/macros/PlotITSTrackingHists.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Moving PWG1 to PWGPP
[u/mrichter/AliRoot.git] / PWGPP / macros / PlotITSTrackingHists.C
CommitLineData
e75c9f34 1void Test(TString fname="QAresults.root") {
2
3 TCanvas *clist;
4 Int_t cnum=1;
5 TFile *f=new TFile(fname.Data());
6 PlotITSTPCMatchingEff(f,clist,cnum);
7
8
9 return;
10}
11//--------------------------------------------------------------------------
12Bool_t PlotITSTPCMatchingEff(TFile *f, TCanvas*& clist,Int_t& cnum) {
13
14 Int_t cnum=1;
15 //clist = new TCanvas[2];
16
17
18 TCanvas* cITSTPCmatch = new TCanvas("cITSTPCmatch","cITSTPCmatch",10,10,1200,600);
19 cITSTPCmatch->Divide(2,1);
20 cITSTPCmatch_1->SetGridy();
21 cITSTPCmatch_2->SetGridy();
22 cITSTPCmatch_1->SetLogx();
23 cITSTPCmatch_2->SetLogx();
24
25 clist = cITSTPCmatch;
26
27 if(!f) return kFALSE;
28
29 TList *list=0;
30 TList *listSPD=0;
31 TDirectoryFile *dir=0;
32
33 // count active SPD HSs
34 dir=(TDirectoryFile*)f->GetDirectory("SPD_Performance");
35 if(dir) listSPD = (TList*)dir->Get("coutput1");
36 if(!dir) return kFALSE;
37
38 Float_t spdFrac[2];
39 TH1F *hnHSsSPD=new TH1F("hnHSsSPD","Active HSs in SPD layers 1 and 2; layer; HSs",2,0.5,2.5);
40 if(listSPD) {
41 //listSPD->Print();
42 TH1F *hFiredChip = (TH1F*)listSPD->FindObject("hFiredChip");
43 Int_t nHSsInner=0,nHSsOuter=0;
44 for(Int_t i=0;i<400;i++) if(hFiredChip->GetBinContent(i)>0) nHSsInner++;
45 for(Int_t i=400;i<1200;i++) if(hFiredChip->GetBinContent(i)>0) nHSsOuter++;
46 nHSsInner = (Int_t)(nHSsInner/10);
47 nHSsOuter = (Int_t)(nHSsOuter/10);
48 hnHSsSPD->SetBinContent(1,nHSsInner);
49 hnHSsSPD->SetBinContent(2,nHSsOuter);
50 spdFrac[0]=(Float_t)nHSsInner/40.;
51 spdFrac[1]=(Float_t)nHSsOuter/80.;
52 }
53 TGraph *spdFrac0=new TGraph(1);
54 spdFrac0->SetPoint(0,0.08,spdFrac[0]);
55 spdFrac0->SetMarkerColor(1); spdFrac0->SetMarkerStyle(20);
56 TGraph *spdFrac1=new TGraph(1);
57 spdFrac1->SetPoint(0,0.08,spdFrac[1]);
58 spdFrac1->SetMarkerColor(1); spdFrac1->SetMarkerStyle(24);
59 TLegend *l2=new TLegend(0.1,0.62,0.5,0.93);
60 l2->SetBorderSize(1);
61 l2->AddEntry(spdFrac0,"Frac. active SPD0","p");
62 l2->AddEntry(spdFrac1,"Frac. active SPD1","p");
63
64 //
65 // Efficiencies for CENTRAL
66 //
67 dir=(TDirectoryFile*)f->GetDirectory("ITS_Performance");
68 if(dir) list = (TList*)dir->Get("cOutputITS_3500_10000");
69 if(!list) return kFALSE;
70
71 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
72 TH1F *fHistPtITSMI6InAcc = (TH1F*)list->FindObject("fHistPtITSMI6InAcc");
73 TH1F *fHistPtITSMI5InAcc = (TH1F*)list->FindObject("fHistPtITSMI5InAcc");
74 TH1F *fHistPtITSMI4InAcc = (TH1F*)list->FindObject("fHistPtITSMI4InAcc");
75 TH1F *fHistPtITSMI3InAcc = (TH1F*)list->FindObject("fHistPtITSMI3InAcc");
76 TH1F *fHistPtITSMI2InAcc = (TH1F*)list->FindObject("fHistPtITSMI2InAcc");
77 TH1F *fHistPtITSMISPDInAcc = (TH1F*)list->FindObject("fHistPtITSMISPDInAcc");
78 TH1F *fHistPtITSMIoneSPDInAcc = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAcc");
79 TH1F *fHistPtITSTPCsel = (TH1F*)list->FindObject("fHistPtITSTPCsel");
80 TH1F *fHistPtITSMIge2InAcc = (TH1F*)fHistPtITSMI6InAcc->Clone("fHistPtITSMIge2InAcc");
81 fHistPtITSMIge2InAcc->Add(fHistPtITSMI5InAcc);
82 fHistPtITSMIge2InAcc->Add(fHistPtITSMI4InAcc);
83 fHistPtITSMIge2InAcc->Add(fHistPtITSMI3InAcc);
84 fHistPtITSMIge2InAcc->Add(fHistPtITSMI2InAcc);
85
86
87 TLegend *l3=new TLegend(0.5,0.62,0.95,0.93);
88 l3->SetBorderSize(1);
89 cITSTPCmatch->cd(1);
90 fHistPtITSMIge2InAcc->SetTitle("Fraction of prolonged tracks with N ITS points: central");
91 fHistPtITSMIge2InAcc->SetYTitle("ITS+TPC / TPC");
92 fHistPtITSMIge2InAcc->Divide(fHistPtITSMIge2InAcc,fHistPtTPCInAcc,1,1,"B");
93 fHistPtITSMIge2InAcc->SetMaximum(1.6);
94 fHistPtITSMIge2InAcc->SetMinimum(0);
95 fHistPtITSMIge2InAcc->GetXaxis()->SetRangeUser(0.1,30);
96 fHistPtITSMIge2InAcc->Draw();
97 l3->AddEntry(fHistPtITSMIge2InAcc,">=2 cls","l");
98 fHistPtITSMI6InAcc->Divide(fHistPtITSMI6InAcc,fHistPtTPCInAcc,1,1,"B");
99 fHistPtITSMI6InAcc->SetLineColor(2);
100 l3->AddEntry(fHistPtITSMI6InAcc,"6 cls","l");
101 fHistPtITSMI6InAcc->Draw("same");
102 fHistPtITSMI5InAcc->Divide(fHistPtITSMI5InAcc,fHistPtTPCInAcc,1,1,"B");
103 fHistPtITSMI5InAcc->SetLineColor(3);
104 l3->AddEntry(fHistPtITSMI5InAcc,"5 cls","l");
105 fHistPtITSMI5InAcc->Draw("same");
106 fHistPtITSMI4InAcc->Divide(fHistPtITSMI4InAcc,fHistPtTPCInAcc,1,1,"B");
107 fHistPtITSMI4InAcc->SetLineColor(4);
108 l3->AddEntry(fHistPtITSMI4InAcc,"4 cls","l");
109 fHistPtITSMI4InAcc->Draw("same");
110 fHistPtITSMI3InAcc->Divide(fHistPtITSMI3InAcc,fHistPtTPCInAcc,1,1,"B");
111 fHistPtITSMI3InAcc->SetLineColor(6);
112 l3->AddEntry(fHistPtITSMI3InAcc,"3 cls","l");
113 fHistPtITSMI3InAcc->Draw("same");
114 fHistPtITSMI2InAcc->Divide(fHistPtITSMI2InAcc,fHistPtTPCInAcc,1,1,"B");
115 fHistPtITSMI2InAcc->SetLineColor(7);
116 l3->AddEntry(fHistPtITSMI2InAcc,"2 cls","l");
117 fHistPtITSMI2InAcc->Draw("same");
118 fHistPtITSMISPDInAcc->Divide(fHistPtITSMISPDInAcc,fHistPtTPCInAcc,1,1,"B");
119 fHistPtITSMISPDInAcc->SetLineColor(9);
120 l3->AddEntry(fHistPtITSMISPDInAcc,"2SPD + any","l");
121 fHistPtITSMISPDInAcc->Draw("same");
122 fHistPtITSMIoneSPDInAcc->Divide(fHistPtITSMIoneSPDInAcc,fHistPtTPCInAcc,1,1,"B");
123 fHistPtITSMIoneSPDInAcc->SetLineColor(15);
124 l3->AddEntry(fHistPtITSMIoneSPDInAcc,">=1SPD + any","l");
125 fHistPtITSMIoneSPDInAcc->Draw("same");
126 fHistPtITSTPCsel->Divide(fHistPtITSTPCsel,fHistPtTPCInAcc,1,1,"B");
127 fHistPtITSTPCsel->SetLineColor(kAzure+1);
128 l3->AddEntry(fHistPtITSTPCsel,">=1SPD + any + d_{0} cut","l");
129 fHistPtITSTPCsel->Draw("same");
130 fHistPtITSMIge2InAcc->Draw("same");
131 l3->Draw();
132 l2->Draw();
133 spdFrac0->Draw("p");
134 spdFrac1->Draw("p");
135 /*
136 if(ioValues) {
137 Int_t index,ptbin;
138 ptbin=fHistPtITSMIge2InAcc->FindBin(0.201);
139 index=2;
140 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
141 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
142 ptbin=fHistPtITSMIge2InAcc->FindBin(1.001);
143 index=3;
144 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
145 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
146 ptbin=fHistPtITSMIge2InAcc->FindBin(10.001);
147 index=4;
148 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
149 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
150 ptbin=fHistPtITSMI6InAcc->FindBin(0.201);
151 index=5;
152 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
153 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
154 ptbin=fHistPtITSMI6InAcc->FindBin(1.001);
155 index=6;
156 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
157 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
158 ptbin=fHistPtITSMI6InAcc->FindBin(10.001);
159 index=7;
160 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
161 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
162 ptbin=fHistPtITSMISPDInAcc->FindBin(0.201);
163 index=8;
164 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
165 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
166 ptbin=fHistPtITSMISPDInAcc->FindBin(1.001);
167 index=9;
168 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
169 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
170 ptbin=fHistPtITSMISPDInAcc->FindBin(10.001);
171 index=10;
172 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
173 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
174 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(0.201);
175 index=11;
176 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
177 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
178 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(1.001);
179 index=12;
180 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
181 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
182 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(10.001);
183 index=13;
184 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
185 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
186 }
187 */
188
189
190
191 //
192 // Efficiencies for PERIPHERAL
193 //
194 dir=(TDirectoryFile*)f->GetDirectory("ITS_Performance");
195 if(dir) list = (TList*)dir->Get("cOutputITS_70_310");
196 if(!list) return kFALSE;
197
198 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
199 TH1F *fHistPtITSMI6InAcc = (TH1F*)list->FindObject("fHistPtITSMI6InAcc");
200 TH1F *fHistPtITSMI5InAcc = (TH1F*)list->FindObject("fHistPtITSMI5InAcc");
201 TH1F *fHistPtITSMI4InAcc = (TH1F*)list->FindObject("fHistPtITSMI4InAcc");
202 TH1F *fHistPtITSMI3InAcc = (TH1F*)list->FindObject("fHistPtITSMI3InAcc");
203 TH1F *fHistPtITSMI2InAcc = (TH1F*)list->FindObject("fHistPtITSMI2InAcc");
204 TH1F *fHistPtITSMISPDInAcc = (TH1F*)list->FindObject("fHistPtITSMISPDInAcc");
205 TH1F *fHistPtITSMIoneSPDInAcc = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAcc");
206 TH1F *fHistPtITSTPCsel = (TH1F*)list->FindObject("fHistPtITSTPCsel");
207 TH1F *fHistPtITSMIge2InAcc = (TH1F*)fHistPtITSMI6InAcc->Clone("fHistPtITSMIge2InAcc");
208 fHistPtITSMIge2InAcc->Add(fHistPtITSMI5InAcc);
209 fHistPtITSMIge2InAcc->Add(fHistPtITSMI4InAcc);
210 fHistPtITSMIge2InAcc->Add(fHistPtITSMI3InAcc);
211 fHistPtITSMIge2InAcc->Add(fHistPtITSMI2InAcc);
212
213
214 cITSTPCmatch->cd(2);
215 fHistPtITSMIge2InAcc->SetTitle("Fraction of prolonged tracks with N ITS points: peripheral");
216 fHistPtITSMIge2InAcc->SetYTitle("ITS+TPC / TPC");
217 fHistPtITSMIge2InAcc->Divide(fHistPtITSMIge2InAcc,fHistPtTPCInAcc,1,1,"B");
218 fHistPtITSMIge2InAcc->SetMaximum(1.6);
219 fHistPtITSMIge2InAcc->SetMinimum(0);
220 fHistPtITSMIge2InAcc->GetXaxis()->SetRangeUser(0.1,30);
221 fHistPtITSMIge2InAcc->Draw();
222 fHistPtITSMI6InAcc->Divide(fHistPtITSMI6InAcc,fHistPtTPCInAcc,1,1,"B");
223 fHistPtITSMI6InAcc->SetLineColor(2);
224 fHistPtITSMI6InAcc->Draw("same");
225 fHistPtITSMI5InAcc->Divide(fHistPtITSMI5InAcc,fHistPtTPCInAcc,1,1,"B");
226 fHistPtITSMI5InAcc->SetLineColor(3);
227 fHistPtITSMI5InAcc->Draw("same");
228 fHistPtITSMI4InAcc->Divide(fHistPtITSMI4InAcc,fHistPtTPCInAcc,1,1,"B");
229 fHistPtITSMI4InAcc->SetLineColor(4);
230 fHistPtITSMI4InAcc->Draw("same");
231 fHistPtITSMI3InAcc->Divide(fHistPtITSMI3InAcc,fHistPtTPCInAcc,1,1,"B");
232 fHistPtITSMI3InAcc->SetLineColor(6);
233 fHistPtITSMI3InAcc->Draw("same");
234 fHistPtITSMI2InAcc->Divide(fHistPtITSMI2InAcc,fHistPtTPCInAcc,1,1,"B");
235 fHistPtITSMI2InAcc->SetLineColor(7);
236 fHistPtITSMI2InAcc->Draw("same");
237 fHistPtITSMISPDInAcc->Divide(fHistPtITSMISPDInAcc,fHistPtTPCInAcc,1,1,"B");
238 fHistPtITSMISPDInAcc->SetLineColor(9);
239 fHistPtITSMISPDInAcc->Draw("same");
240 fHistPtITSMIoneSPDInAcc->Divide(fHistPtITSMIoneSPDInAcc,fHistPtTPCInAcc,1,1,"B");
241 fHistPtITSMIoneSPDInAcc->SetLineColor(15);
242 fHistPtITSMIoneSPDInAcc->Draw("same");
243 fHistPtITSTPCsel->Divide(fHistPtITSTPCsel,fHistPtTPCInAcc,1,1,"B");
244 fHistPtITSTPCsel->SetLineColor(kAzure+1);
245 fHistPtITSTPCsel->Draw("same");
246 fHistPtITSMIge2InAcc->Draw("same");
247 l3->Draw();
248 l2->Draw();
249 spdFrac0->Draw("p");
250 spdFrac1->Draw("p");
251 /*
252 if(ioValues) {
253 Int_t index,ptbin;
254 ptbin=fHistPtITSMIge2InAcc->FindBin(0.201);
255 index=2;
256 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
257 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
258 ptbin=fHistPtITSMIge2InAcc->FindBin(1.001);
259 index=3;
260 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
261 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
262 ptbin=fHistPtITSMIge2InAcc->FindBin(10.001);
263 index=4;
264 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
265 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
266 ptbin=fHistPtITSMI6InAcc->FindBin(0.201);
267 index=5;
268 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
269 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
270 ptbin=fHistPtITSMI6InAcc->FindBin(1.001);
271 index=6;
272 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
273 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
274 ptbin=fHistPtITSMI6InAcc->FindBin(10.001);
275 index=7;
276 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
277 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
278 ptbin=fHistPtITSMISPDInAcc->FindBin(0.201);
279 index=8;
280 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
281 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
282 ptbin=fHistPtITSMISPDInAcc->FindBin(1.001);
283 index=9;
284 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
285 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
286 ptbin=fHistPtITSMISPDInAcc->FindBin(10.001);
287 index=10;
288 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
289 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
290 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(0.201);
291 index=11;
292 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
293 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
294 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(1.001);
295 index=12;
296 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
297 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
298 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(10.001);
299 index=13;
300 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
301 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
302 }
303 */
304
305
306
307 return kTRUE;
308}
309
b51e971f 310Bool_t PlotITSTrackingHists(TString fname="ITS.Performance.root",
311 Float_t *ioValues=0,Float_t *ioErrors=0)
312{
2308c814 313 //
314 // Macro to plot the histos from the task AliAnalysisTaskITSTrackingCheck
315 // A. Dainese 28.11.09
316 //
317
b51e971f 318 Bool_t plotAlignmentChecks=kFALSE;
fa516524 319 //gStyle->SetOptStat(0);
2308c814 320
d388aa11 321 //if(fname.Contains("alien")) TGrid::Connect("alien://");
8fa56655 322
323 TFile *f= TFile::Open(fname.Data());
b51e971f 324 if(!f) return kFALSE;
2308c814 325
fa516524 326 TList *list=(TList*)f->Get("cOutputITS");
b51e971f 327 TList *listSPD=0;
af874c53 328 TDirectoryFile *dir=0;
8fa56655 329 if(!list) {
af874c53 330 dir=(TDirectoryFile*)f->GetDirectory("ITS_Performance");
e75c9f34 331 if(dir) list = (TList*)dir->Get("cOutputITS_70_310");
b51e971f 332 // count active SPD HSs
333 dir=(TDirectoryFile*)f->GetDirectory("SPD_Performance");
334 if(dir) listSPD = (TList*)dir->Get("coutput1");
8fa56655 335 }
b51e971f 336
af874c53 337 if(!list) return kFALSE;
b51e971f 338
339 TH1F *hnHSsSPD=new TH1F("hnHSsSPD","Active HSs in SPD layers 1 and 2; layer; HSs",2,0.5,2.5);
340 if(listSPD) {
af874c53 341 //listSPD->Print();
b51e971f 342 TH1F *hFiredChip = (TH1F*)listSPD->FindObject("hFiredChip");
343 Int_t nHSsInner=0,nHSsOuter=0;
344 for(Int_t i=0;i<400;i++) if(hFiredChip->GetBinContent(i)>0) nHSsInner++;
345 for(Int_t i=400;i<1200;i++) if(hFiredChip->GetBinContent(i)>0) nHSsOuter++;
346 nHSsInner = (Int_t)(nHSsInner/10);
347 nHSsOuter = (Int_t)(nHSsOuter/10);
348 hnHSsSPD->SetBinContent(1,nHSsInner);
349 hnHSsSPD->SetBinContent(2,nHSsOuter);
350 if(ioValues) ioValues[0]=(Float_t)nHSsInner/40.;
351 if(ioValues) ioValues[1]=(Float_t)nHSsOuter/80.;
352 }
353
354 TCanvas *cSPD=new TCanvas("cSPD","cSPD");
355 hnHSsSPD->SetMaximum(90);
356 hnHSsSPD->SetMinimum(0);
357 hnHSsSPD->Draw();
358
2308c814 359 TH1F *fHistNclsITSSA = (TH1F*)list->FindObject("fHistNclsITSSA");
360 TH1F *fHistClusterMapITSSA = (TH1F*)list->FindObject("fHistClusterMapITSSA");
361 TH1F *fHistClusterMapITSSAok = (TH1F*)list->FindObject("fHistClusterMapITSSAok");
362 TH1F *fHistClusterMapITSSAbad = (TH1F*)list->FindObject("fHistClusterMapITSSAbad");
363 TH1F *fHistClusterMapITSSAskipped = (TH1F*)list->FindObject("fHistClusterMapITSSAskipped");
364 TH1F *fHistClusterMapITSSAoutinz = (TH1F*)list->FindObject("fHistClusterMapITSSAoutinz");
365 TH1F *fHistClusterMapITSSAokoutinzbad = (TH1F*)list->FindObject("fHistClusterMapITSSAokoutinzbad");
366 TH1F *fHistClusterMapITSSAnorefit = (TH1F*)list->FindObject("fHistClusterMapITSSAnorefit");
367 TH1F *fHistClusterMapITSSAnocls = (TH1F*)list->FindObject("fHistClusterMapITSSAnocls");
368 TH1F *fHistNclsITSSAInAcc = (TH1F*)list->FindObject("fHistNclsITSSAInAcc");
369 TH1F *fHistClusterMapITSSAInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAInAcc");
370 TH1F *fHistClusterMapITSSAokInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAokInAcc");
371 TH1F *fHistClusterMapITSSAbadInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAbadInAcc");
fa516524 372 TH1F *fHistClusterMapModuleITSMIokInAcc = (TH1F*)list->FindObject("fHistClusterMapModuleITSMIokInAcc");
373 TH1F *fHistClusterMapModuleITSMIbadInAcc = (TH1F*)list->FindObject("fHistClusterMapModuleITSMIbadInAcc");
2308c814 374 TH1F *fHistClusterMapITSSAskippedInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAskippedInAcc");
375 TH1F *fHistClusterMapITSSAoutinzInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAoutinzInAcc");
376 TH1F *fHistClusterMapITSSAokoutinzbadInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAokoutinzbadInAcc");
377 TH1F *fHistClusterMapITSSAnorefitInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAnorefitInAcc");
378 TH1F *fHistClusterMapITSSAnoclsInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAnoclsInAcc");
fa516524 379 TH1F *fHistClusterMapModuleITSMInoclsInAcc = (TH1F*)list->FindObject("fHistClusterMapModuleITSMInoclsInAcc");
2308c814 380 TH1F *fHistNclsITSMI = (TH1F*)list->FindObject("fHistNclsITSMI");
381 TH1F *fHistClusterMapITSMI = (TH1F*)list->FindObject("fHistClusterMapITSMI");
382 TH1F *fHistClusterMapITSMIok = (TH1F*)list->FindObject("fHistClusterMapITSMIok");
383 TH1F *fHistClusterMapITSMIbad = (TH1F*)list->FindObject("fHistClusterMapITSMIbad");
384 TH1F *fHistClusterMapITSMIskipped = (TH1F*)list->FindObject("fHistClusterMapITSMIskipped");
385 TH1F *fHistClusterMapITSMIoutinz = (TH1F*)list->FindObject("fHistClusterMapITSMIoutinz");
386 TH1F *fHistClusterMapITSMIokoutinzbad = (TH1F*)list->FindObject("fHistClusterMapITSMIokoutinzbad");
387 TH1F *fHistClusterMapITSMInorefit = (TH1F*)list->FindObject("fHistClusterMapITSMInorefit");
388 TH1F *fHistClusterMapITSMInocls = (TH1F*)list->FindObject("fHistClusterMapITSMInocls");
389
c2893e82 390 TH1F *fHistPhiTPCInAcc = (TH1F*)list->FindObject("fHistPhiTPCInAcc");
391 TH1F *fHistPhiITSMIokbadoutinz6InAcc = (TH1F*)list->FindObject("fHistPhiITSMIokbadoutinz6InAcc");
392
393 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
8fa56655 394 TH1F *fHistPtTPCInAccSfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccSfromStrange");
395 TH1F *fHistPtTPCInAccPfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccPfromStrange");
fa516524 396 TH1F *fHistPtTPCInAccMCtwoSPD = (TH1F*)list->FindObject("fHistPtTPCInAccMCtwoSPD");
397 TH1F *fHistPtTPCInAccMConeSPD = (TH1F*)list->FindObject("fHistPtTPCInAccMConeSPD");
c2893e82 398 TH1F *fHistPtITSMI6InAcc = (TH1F*)list->FindObject("fHistPtITSMI6InAcc");
399 TH1F *fHistPtITSMI5InAcc = (TH1F*)list->FindObject("fHistPtITSMI5InAcc");
400 TH1F *fHistPtITSMI4InAcc = (TH1F*)list->FindObject("fHistPtITSMI4InAcc");
401 TH1F *fHistPtITSMI3InAcc = (TH1F*)list->FindObject("fHistPtITSMI3InAcc");
402 TH1F *fHistPtITSMI2InAcc = (TH1F*)list->FindObject("fHistPtITSMI2InAcc");
403 TH1F *fHistPtITSMISPDInAcc = (TH1F*)list->FindObject("fHistPtITSMISPDInAcc");
fa516524 404 TH1F *fHistPtITSMIoneSPDInAcc = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAcc");
af874c53 405 TH1F *fHistPtITSMI6InAccFake = (TH1F*)list->FindObject("fHistPtITSMI6InAccFake");
406 TH1F *fHistPtITSMI5InAccFake = (TH1F*)list->FindObject("fHistPtITSMI5InAccFake");
407 TH1F *fHistPtITSMI4InAccFake = (TH1F*)list->FindObject("fHistPtITSMI4InAccFake");
408 TH1F *fHistPtITSMI3InAccFake = (TH1F*)list->FindObject("fHistPtITSMI3InAccFake");
409 TH1F *fHistPtITSMI2InAccFake = (TH1F*)list->FindObject("fHistPtITSMI2InAccFake");
410 TH1F *fHistPtITSMISPDInAccFake = (TH1F*)list->FindObject("fHistPtITSMISPDInAccFake");
411 TH1F *fHistPtITSMIoneSPDInAccFake = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAccFake");
8fa56655 412 TH1F *fHistPtITSMIoneSPDthreeSDDSSDInAcc = (TH1F*)list->FindObject("fHistPtITSMIoneSPDthreeSDDSSDInAcc");
c2893e82 413 TH1F *fHistPtITSMIokbadoutinz6InAcc = (TH1F*)list->FindObject("fHistPtITSMIokbadoutinz6InAcc");
414 TH1F *fHistPtITSMIokbadoutinz5InAcc = (TH1F*)list->FindObject("fHistPtITSMIokbadoutinz5InAcc");
415 TH1F *fHistPtITSMIokbadoutinz4InAcc = (TH1F*)list->FindObject("fHistPtITSMIokbadoutinz4InAcc");
fa516524 416 TH1F *fHistPtTPCInAccP = (TH1F*)list->FindObject("fHistPtTPCInAccP");
417 TH1F *fHistPtITSMI6InAccP = (TH1F*)list->FindObject("fHistPtITSMI6InAccP");
418 TH1F *fHistPtITSMI5InAccP = (TH1F*)list->FindObject("fHistPtITSMI5InAccP");
419 TH1F *fHistPtITSMI4InAccP = (TH1F*)list->FindObject("fHistPtITSMI4InAccP");
420 TH1F *fHistPtITSMI3InAccP = (TH1F*)list->FindObject("fHistPtITSMI3InAccP");
421 TH1F *fHistPtITSMI2InAccP = (TH1F*)list->FindObject("fHistPtITSMI2InAccP");
422 TH1F *fHistPtITSMISPDInAccP = (TH1F*)list->FindObject("fHistPtITSMISPDInAccP");
423 TH1F *fHistPtITSMIoneSPDInAccP = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAccP");
424 TH1F *fHistPtTPCInAccS = (TH1F*)list->FindObject("fHistPtTPCInAccS");
425 TH1F *fHistPtITSMI6InAccS = (TH1F*)list->FindObject("fHistPtITSMI6InAccS");
426 TH1F *fHistPtITSMI5InAccS = (TH1F*)list->FindObject("fHistPtITSMI5InAccS");
427 TH1F *fHistPtITSMI4InAccS = (TH1F*)list->FindObject("fHistPtITSMI4InAccS");
428 TH1F *fHistPtITSMI3InAccS = (TH1F*)list->FindObject("fHistPtITSMI3InAccS");
429 TH1F *fHistPtITSMI2InAccS = (TH1F*)list->FindObject("fHistPtITSMI2InAccS");
430 TH1F *fHistPtITSMISPDInAccS = (TH1F*)list->FindObject("fHistPtITSMISPDInAccS");
431 TH1F *fHistPtITSMIoneSPDInAccS = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAccS");
432
433 TH1F *fHistRProdVtxInAccP = (TH1F*)list->FindObject("fHistRProdVtxInAccP");
434 TH1F *fHistRProdVtxInAccS = (TH1F*)list->FindObject("fHistRProdVtxInAccS");
435
436 TH1F *fHistxlocSDDok = (TH1F*)list->FindObject("fHistxlocSDDok");
437 TH1F *fHistxlocSDDall = (TH1F*)list->FindObject("fHistxlocSDDall");
438 TH1F *fHistzlocSDDok = (TH1F*)list->FindObject("fHistzlocSDDok");
439 TH1F *fHistzlocSDDall = (TH1F*)list->FindObject("fHistzlocSDDall");
440
441 TH1F *fHistPtITSTPCsel = (TH1F*)list->FindObject("fHistPtITSTPCsel");
442 TH1F *fHistPtITSTPCselP = (TH1F*)list->FindObject("fHistPtITSTPCselP");
443 TH1F *fHistPtITSTPCselS = (TH1F*)list->FindObject("fHistPtITSTPCselS");
af874c53 444 TH1F *fHistPtITSTPCselFake = (TH1F*)list->FindObject("fHistPtITSTPCselFake");
fa516524 445 TH1F *fHistPtITSTPCselPfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
446 TH1F *fHistPtITSTPCselSfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
447 TH1F *fHistPtITSTPCselSfromMat = (TH1F*)list->FindObject("fHistPtITSTPCselSfromMat");
448
fa516524 449 TH1F *fHistPtTPCInAccSfromMat = (TH1F*)list->FindObject("fHistPtTPCInAccSfromMat");
450
af874c53 451
8fa56655 452 //ReweightStrange(fHistPtTPCInAcc,fHistPtTPCInAccPfromStrange,fHistPtTPCInAccSfromStrange);
453 //ReweightStrange(fHistPtITSTPCsel,fHistPtITSTPCselPfromStrange,fHistPtITSTPCselSfromStrange);
c2893e82 454
455 //---------------------------------------------------------------
fa516524 456 TH1F *fHistPtITSMIge2InAcc = (TH1F*)fHistPtITSMI6InAcc->Clone("fHistPtITSMIge2InAcc");
457 fHistPtITSMIge2InAcc->Add(fHistPtITSMI5InAcc);
458 fHistPtITSMIge2InAcc->Add(fHistPtITSMI4InAcc);
459 fHistPtITSMIge2InAcc->Add(fHistPtITSMI3InAcc);
460 fHistPtITSMIge2InAcc->Add(fHistPtITSMI2InAcc);
af874c53 461 TH1F *fHistPtITSMIge2InAccFake = 0;
462 if(fHistPtITSMI6InAccFake) {
463 fHistPtITSMIge2InAccFake=(TH1F*)fHistPtITSMI6InAccFake->Clone("fHistPtITSMIge2InAccFake");
464 fHistPtITSMIge2InAccFake->Add(fHistPtITSMI5InAccFake);
465 fHistPtITSMIge2InAccFake->Add(fHistPtITSMI4InAccFake);
466 fHistPtITSMIge2InAccFake->Add(fHistPtITSMI3InAccFake);
467 fHistPtITSMIge2InAccFake->Add(fHistPtITSMI2InAccFake);
468 }
fa516524 469 TH1F *fHistPtITSMIge2InAccP = (TH1F*)fHistPtITSMI6InAccP->Clone("fHistPtITSMIge2InAccP");
470 fHistPtITSMIge2InAccP->Add(fHistPtITSMI5InAccP);
471 fHistPtITSMIge2InAccP->Add(fHistPtITSMI4InAccP);
472 fHistPtITSMIge2InAccP->Add(fHistPtITSMI3InAccP);
473 fHistPtITSMIge2InAccP->Add(fHistPtITSMI2InAccP);
474 TH1F *fHistPtITSMIge2InAccS = (TH1F*)fHistPtITSMI6InAccS->Clone("fHistPtITSMIge2InAccS");
475 fHistPtITSMIge2InAccS->Add(fHistPtITSMI5InAccS);
476 fHistPtITSMIge2InAccS->Add(fHistPtITSMI4InAccS);
477 fHistPtITSMIge2InAccS->Add(fHistPtITSMI3InAccS);
478 fHistPtITSMIge2InAccS->Add(fHistPtITSMI2InAccS);
479
af874c53 480 // fake fraction
481 if(fHistPtITSMIge2InAccFake) {
482 fHistPtITSMIge2InAccFake->Divide(fHistPtITSMIge2InAccFake,fHistPtITSMIge2InAcc,1,1,"B");
483 fHistPtITSMI2InAccFake->Divide(fHistPtITSMI2InAccFake,fHistPtITSMI2InAcc,1,1,"B");
484 fHistPtITSMI3InAccFake->Divide(fHistPtITSMI3InAccFake,fHistPtITSMI3InAcc,1,1,"B");
485 fHistPtITSMI4InAccFake->Divide(fHistPtITSMI4InAccFake,fHistPtITSMI4InAcc,1,1,"B");
486 fHistPtITSMI5InAccFake->Divide(fHistPtITSMI5InAccFake,fHistPtITSMI5InAcc,1,1,"B");
487 fHistPtITSMI6InAccFake->Divide(fHistPtITSMI6InAccFake,fHistPtITSMI6InAcc,1,1,"B");
488 fHistPtITSMISPDInAccFake->Divide(fHistPtITSMISPDInAccFake,fHistPtITSMISPDInAcc,1,1,"B");
489 fHistPtITSMIoneSPDInAccFake->Divide(fHistPtITSMIoneSPDInAccFake,fHistPtITSMIoneSPDInAcc,1,1,"B");
490 if (fHistPtITSTPCselFake)fHistPtITSTPCselFake->Divide(fHistPtITSTPCselFake,fHistPtITSTPCsel,1,1,"B");
491 }
492
fa516524 493 TH1F* fHistPtITSMISPDInAccMC=(TH1F*)fHistPtITSMISPDInAcc->Clone("fHistPtITSMISPDInAccMC");
494 TH1F* fHistPtITSMIoneSPDInAccMC=(TH1F*)fHistPtITSMIoneSPDInAcc->Clone("fHistPtITSMIoneSPDInAccMC");
495
c2893e82 496
2308c814 497 TLegend *l1=new TLegend(0.5,0.5,0.9,0.9);
498 TLegend *l2=new TLegend(0.5,0.5,0.9,0.9);
499
500 TCanvas *c1= new TCanvas("c1","c1",10,10,600,500);
501 fHistNclsITSSA->SetMinimum(0);
502 fHistNclsITSSA->SetLineColor(1);
503 l1->AddEntry(fHistNclsITSSA,"ITS-SA","l");
504 fHistNclsITSSA->Draw();
505 fHistNclsITSSAInAcc->SetLineColor(4);
506 l1->AddEntry(fHistNclsITSSAInAcc,"ITS-SA in acc.","l");
507 fHistNclsITSSAInAcc->Draw("same");
508 fHistNclsITSMI->SetLineColor(2);
509 l1->AddEntry(fHistNclsITSMI,"ITS from TPC","l");
510 fHistNclsITSMI->Draw("same");
511 l1->Draw();
512
513
514 TCanvas *c2 =new TCanvas("c2","c2",10,10,1200,800);
515 c2->Divide(3,2);
516 c2->cd(1);
517 //
518 fHistClusterMapITSSAokoutinzbad->SetLineColor(1);
519 fHistClusterMapITSSAokoutinzbad->SetMarkerColor(1);
520 fHistClusterMapITSSAokoutinzbad->SetMarkerStyle(20);
521 fHistClusterMapITSSAokoutinzbad->Draw();
c2893e82 522 fHistClusterMapITSMIokoutinzbad->SetLineColor(2);
523 fHistClusterMapITSMIokoutinzbad->SetMarkerColor(2);
2308c814 524 fHistClusterMapITSMIokoutinzbad->SetMarkerStyle(20);
525 fHistClusterMapITSMIokoutinzbad->Draw("same");
526 l1->Draw();
527 //
528 c2->cd(2);
529 fHistClusterMapITSSAok->SetLineColor(1);
530 fHistClusterMapITSSAok->SetMarkerColor(1);
531 fHistClusterMapITSSAok->SetMarkerStyle(21);
532 fHistClusterMapITSSAok->Draw();
533 fHistClusterMapITSMIok->SetLineColor(2);
534 fHistClusterMapITSMIok->SetMarkerColor(2);
535 fHistClusterMapITSMIok->SetMarkerStyle(21);
536 fHistClusterMapITSMIok->Draw("same");
537 //
538 c2->cd(3);
539 fHistClusterMapITSSAoutinz->SetLineColor(1);
540 fHistClusterMapITSSAoutinz->SetMarkerColor(1);
541 fHistClusterMapITSSAoutinz->SetMarkerStyle(22);
542 fHistClusterMapITSSAoutinz->Draw();
543 fHistClusterMapITSMIoutinz->SetLineColor(2);
544 fHistClusterMapITSMIoutinz->SetMarkerColor(2);
545 fHistClusterMapITSMIoutinz->SetMarkerStyle(22);
546 fHistClusterMapITSMIoutinz->Draw("same");
547 //
548 c2->cd(4);
549 fHistClusterMapITSSAbad->SetLineColor(1);
550 fHistClusterMapITSSAbad->SetMarkerColor(1);
551 fHistClusterMapITSSAbad->SetMarkerStyle(23);
552 fHistClusterMapITSSAbad->Draw();
553 fHistClusterMapITSMIbad->SetLineColor(2);
554 fHistClusterMapITSMIbad->SetMarkerColor(2);
555 fHistClusterMapITSMIbad->SetMarkerStyle(23);
556 fHistClusterMapITSMIbad->Draw("same");
557 //
558 c2->cd(5);
559 fHistClusterMapITSSAnocls->SetLineColor(1);
560 fHistClusterMapITSSAnocls->SetMarkerColor(1);
561 fHistClusterMapITSSAnocls->SetMarkerStyle(24);
562 fHistClusterMapITSSAnocls->Draw();
563 fHistClusterMapITSMInocls->SetLineColor(2);
564 fHistClusterMapITSMInocls->SetMarkerColor(2);
565 fHistClusterMapITSMInocls->SetMarkerStyle(24);
566 fHistClusterMapITSMInocls->Draw("same");
567
fa516524 568 TCanvas *c3 =new TCanvas("c3","c3",10,10,1200,800);
569 c3->Divide(1,2);
570 c3->cd(1);
571 fHistClusterMapModuleITSMIokInAcc->SetLineColor(1);
572 fHistClusterMapModuleITSMIokInAcc->Draw();
573 l2->AddEntry(fHistClusterMapModuleITSMIokInAcc,"ok","l");
574 fHistClusterMapModuleITSMIbadInAcc->SetLineColor(3);
575 fHistClusterMapModuleITSMIbadInAcc->Draw("same");
576 l2->AddEntry(fHistClusterMapModuleITSMIbadInAcc,"bad","l");
577 fHistClusterMapModuleITSMInoclsInAcc->SetLineColor(2);
578 fHistClusterMapModuleITSMInoclsInAcc->Draw("same");
579 l2->AddEntry(fHistClusterMapModuleITSMInoclsInAcc,"no cls","l");
2308c814 580 l2->Draw();
fa516524 581 c3->cd(2);
582 TH1F *fHistClusterMapModuleITSMIallInAcc=(TH1F*)fHistClusterMapModuleITSMIokInAcc->Clone("fHistClusterMapModuleITSMIallInAcc");
583 fHistClusterMapModuleITSMIallInAcc->Add(fHistClusterMapModuleITSMIbadInAcc);
584 fHistClusterMapModuleITSMIallInAcc->Add(fHistClusterMapModuleITSMInoclsInAcc);
585 TH1F *fHistClusterMapModuleITSMIokRatioInAcc=(TH1F*)fHistClusterMapModuleITSMIokInAcc->Clone("fHistClusterMapModuleITSMIokRatioInAcc");
586 fHistClusterMapModuleITSMIokRatioInAcc->Divide(fHistClusterMapModuleITSMIallInAcc);
587 for(Int_t ib=1;ib<=fHistClusterMapModuleITSMIokRatioInAcc->GetNbinsX();ib++) {
588 fHistClusterMapModuleITSMIokRatioInAcc->SetBinError(ib,0);
589 if(fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)) fHistClusterMapModuleITSMIokRatioInAcc->SetBinError(ib,TMath::Sqrt(fHistClusterMapModuleITSMIokRatioInAcc->GetBinContent(ib)*(1.-fHistClusterMapModuleITSMIokRatioInAcc->GetBinContent(ib))/fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)));
590 }
591 fHistClusterMapModuleITSMIokRatioInAcc->Draw();
592 TH1F *fHistClusterMapModuleITSMIbadRatioInAcc=(TH1F*)fHistClusterMapModuleITSMIbadInAcc->Clone("fHistClusterMapModuleITSMIbadRatioInAcc");
593 fHistClusterMapModuleITSMIbadRatioInAcc->Divide(fHistClusterMapModuleITSMIallInAcc);
594 for(Int_t ib=1;ib<=fHistClusterMapModuleITSMIokRatioInAcc->GetNbinsX();ib++) {
595 fHistClusterMapModuleITSMIbadRatioInAcc->SetBinError(ib,0);
596 if(fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)) fHistClusterMapModuleITSMIbadRatioInAcc->SetBinError(ib,TMath::Sqrt(fHistClusterMapModuleITSMIbadRatioInAcc->GetBinContent(ib)*(1.-fHistClusterMapModuleITSMIbadRatioInAcc->GetBinContent(ib))/fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)));
597 }
598 fHistClusterMapModuleITSMIbadRatioInAcc->Draw("same");
599 TH1F *fHistClusterMapModuleITSMInoclsRatioInAcc=(TH1F*)fHistClusterMapModuleITSMInoclsInAcc->Clone("fHistClusterMapModuleITSMInoclsRatioInAcc");
600 fHistClusterMapModuleITSMInoclsRatioInAcc->Divide(fHistClusterMapModuleITSMIallInAcc);
601 for(Int_t ib=1;ib<=fHistClusterMapModuleITSMIokRatioInAcc->GetNbinsX();ib++) {
602 fHistClusterMapModuleITSMInoclsRatioInAcc->SetBinError(ib,0);
603 if(fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)) fHistClusterMapModuleITSMInoclsRatioInAcc->SetBinError(ib,TMath::Sqrt(fHistClusterMapModuleITSMInoclsRatioInAcc->GetBinContent(ib)*(1.-fHistClusterMapModuleITSMInoclsRatioInAcc->GetBinContent(ib))/fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)));
604 }
605 fHistClusterMapModuleITSMInoclsRatioInAcc->Draw("same");
2308c814 606
fa516524 607 TCanvas *c3b = new TCanvas("c3b","c3b");
608 c3b->Divide(2,1);
609 c3b->cd(1);
610 fHistxlocSDDok->Divide(fHistxlocSDDall);
611 fHistxlocSDDok->Draw();
612 c3b->cd(2);
613 fHistzlocSDDok->Divide(fHistzlocSDDall);
614 fHistzlocSDDok->Draw();
2308c814 615
c2893e82 616 TCanvas *c4 =new TCanvas("c4","c4",10,10,500,500);
617 c4->SetGridy();
618 fHistPhiITSMIokbadoutinz6InAcc->Divide(fHistPhiTPCInAcc);
619 fHistPhiITSMIokbadoutinz6InAcc->SetMinimum(0);
620 fHistPhiITSMIokbadoutinz6InAcc->SetMaximum(1.5);
621 fHistPhiITSMIokbadoutinz6InAcc->SetYTitle("ITS+TPC / TPC");
622 fHistPhiITSMIokbadoutinz6InAcc->SetTitle("Fraction of tracks with 6 layers ok");
623 fHistPhiITSMIokbadoutinz6InAcc->Draw();
624
fa516524 625 TCanvas *c5c =new TCanvas("c5c","c5c",10,10,600,600);
626 c5c->SetGridy();
627 TH1F *fHistPtTPCInAccSecFrac = (TH1F*)fHistPtTPCInAccS->Clone("fHistPtTPCInAccSecFrac");
628 fHistPtTPCInAccSecFrac->Divide(fHistPtTPCInAcc);
629 fHistPtTPCInAccSecFrac->SetLineColor(kOrange+7);
630 fHistPtTPCInAccSecFrac->SetTitle("Fraction of secondaries");
631 fHistPtTPCInAccSecFrac->SetYTitle("sec/all");
632 fHistPtTPCInAccSecFrac->Draw();
633 TH1F *fHistPtITSMI2InAccSecFrac = (TH1F*)fHistPtITSMI2InAccS->Clone("fHistPtITSMI2InAccSecFrac");
634 fHistPtITSMI2InAccSecFrac->Divide(fHistPtITSMI2InAcc);
635 fHistPtITSMI2InAccSecFrac->SetLineColor(7);
636 fHistPtITSMI2InAccSecFrac->Draw("same");
637 TH1F *fHistPtITSMIge2InAccSecFrac = (TH1F*)fHistPtITSMIge2InAccS->Clone("fHistPtITSMIge2InAccSecFrac");
638 fHistPtITSMIge2InAccSecFrac->Divide(fHistPtITSMIge2InAcc);
639 fHistPtITSMIge2InAccSecFrac->Draw("same");
640 TH1F *fHistPtITSMI3InAccSecFrac = (TH1F*)fHistPtITSMI3InAccS->Clone("fHistPtITSMI3InAccSecFrac");
641 fHistPtITSMI3InAccSecFrac->Divide(fHistPtITSMI3InAcc);
642 fHistPtITSMI3InAccSecFrac->SetLineColor(6);
643 fHistPtITSMI3InAccSecFrac->Draw("same");
644 TH1F *fHistPtITSMI4InAccSecFrac = (TH1F*)fHistPtITSMI4InAccS->Clone("fHistPtITSMI4InAccSecFrac");
645 fHistPtITSMI4InAccSecFrac->Divide(fHistPtITSMI4InAcc);
646 fHistPtITSMI4InAccSecFrac->SetLineColor(4);
647 fHistPtITSMI4InAccSecFrac->Draw("same");
648 TH1F *fHistPtITSMI5InAccSecFrac = (TH1F*)fHistPtITSMI5InAccS->Clone("fHistPtITSMI5InAccSecFrac");
649 fHistPtITSMI5InAccSecFrac->Divide(fHistPtITSMI5InAcc);
650 fHistPtITSMI5InAccSecFrac->SetLineColor(3);
651 fHistPtITSMI5InAccSecFrac->Draw("same");
652 TH1F *fHistPtITSMI6InAccSecFrac = (TH1F*)fHistPtITSMI6InAccS->Clone("fHistPtITSMI6InAccSecFrac");
653 fHistPtITSMI6InAccSecFrac->Divide(fHistPtITSMI6InAcc);
654 fHistPtITSMI6InAccSecFrac->SetLineColor(2);
655 fHistPtITSMI6InAccSecFrac->SetLineColor(2);
656 fHistPtITSMI6InAccSecFrac->Draw("same");
657 TH1F *fHistPtITSMISPDInAccSecFrac = (TH1F*)fHistPtITSMISPDInAccS->Clone("fHistPtITSMISPDInAccSecFrac");
658 fHistPtITSMISPDInAccSecFrac->Divide(fHistPtITSMISPDInAcc);
659 fHistPtITSMISPDInAccSecFrac->SetLineColor(9);
660 fHistPtITSMISPDInAccSecFrac->Draw("same");
661 TH1F *fHistPtITSMIoneSPDInAccSecFrac = (TH1F*)fHistPtITSMIoneSPDInAccS->Clone("fHistPtITSMIoneSPDInAccSecFrac");
662 fHistPtITSMIoneSPDInAccSecFrac->Divide(fHistPtITSMIoneSPDInAcc);
663 fHistPtITSMIoneSPDInAccSecFrac->SetLineColor(15);
664 fHistPtITSMIoneSPDInAccSecFrac->Draw("same");
665 TH1F *fHistPtITSTPCselSecFrac = (TH1F*)fHistPtITSTPCselS->Clone("fHistPtITSTPCselSecFrac");
666 fHistPtITSTPCselSecFrac->Divide(fHistPtITSTPCsel);
667 fHistPtITSTPCselSecFrac->SetLineColor(15);
668 fHistPtITSTPCselSecFrac->Draw("same");
669
670
c2893e82 671 TLegend *l3=new TLegend(0.5,0.5,0.9,0.9);
672 TLegend *l4=new TLegend(0.5,0.5,0.9,0.9);
673 TCanvas *c5 =new TCanvas("c5","c5",10,10,1200,600);
674 c5->Divide(2,1);
675 c5_1->SetGridy();
676 c5_2->SetGridy();
8fa56655 677 c5_1->SetLogx();
678 c5_2->SetLogx();
c2893e82 679 c5->cd(1);
c2893e82 680 fHistPtITSMIge2InAcc->SetMaximum(1.5);
681 fHistPtITSMIge2InAcc->SetMinimum(0);
682 fHistPtITSMIge2InAcc->SetTitle("Fraction of prolonged tracks with N ITS points");
683 fHistPtITSMIge2InAcc->SetYTitle("ITS+TPC / TPC");
8fa56655 684 fHistPtITSMIge2InAcc->Divide(fHistPtITSMIge2InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 685 fHistPtITSMIge2InAcc->Draw();
686 l3->AddEntry(fHistPtITSMIge2InAcc,">=2 cls","l");
8fa56655 687 fHistPtITSMI6InAcc->Divide(fHistPtITSMI6InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 688 fHistPtITSMI6InAcc->SetLineColor(2);
689 l3->AddEntry(fHistPtITSMI6InAcc,"6 cls","l");
690 fHistPtITSMI6InAcc->Draw("same");
8fa56655 691 fHistPtITSMI5InAcc->Divide(fHistPtITSMI5InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 692 fHistPtITSMI5InAcc->SetLineColor(3);
693 l3->AddEntry(fHistPtITSMI5InAcc,"5 cls","l");
694 fHistPtITSMI5InAcc->Draw("same");
8fa56655 695 fHistPtITSMI4InAcc->Divide(fHistPtITSMI4InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 696 fHistPtITSMI4InAcc->SetLineColor(4);
697 l3->AddEntry(fHistPtITSMI4InAcc,"4 cls","l");
698 fHistPtITSMI4InAcc->Draw("same");
8fa56655 699 fHistPtITSMI3InAcc->Divide(fHistPtITSMI3InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 700 fHistPtITSMI3InAcc->SetLineColor(6);
701 l3->AddEntry(fHistPtITSMI3InAcc,"3 cls","l");
702 fHistPtITSMI3InAcc->Draw("same");
8fa56655 703 fHistPtITSMI2InAcc->Divide(fHistPtITSMI2InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 704 fHistPtITSMI2InAcc->SetLineColor(7);
705 l3->AddEntry(fHistPtITSMI2InAcc,"2 cls","l");
706 fHistPtITSMI2InAcc->Draw("same");
8fa56655 707 fHistPtITSMISPDInAcc->Divide(fHistPtITSMISPDInAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 708 fHistPtITSMISPDInAcc->SetLineColor(9);
709 l3->AddEntry(fHistPtITSMISPDInAcc,"2SPD + any","l");
710 fHistPtITSMISPDInAcc->Draw("same");
8fa56655 711 fHistPtITSMIoneSPDInAcc->Divide(fHistPtITSMIoneSPDInAcc,fHistPtTPCInAcc,1,1,"B");
fa516524 712 fHistPtITSMIoneSPDInAcc->SetLineColor(15);
713 l3->AddEntry(fHistPtITSMIoneSPDInAcc,">=1SPD + any","l");
714 fHistPtITSMIoneSPDInAcc->Draw("same");
8fa56655 715 //fHistPtITSMIoneSPDthreeSDDSSDInAcc->Divide(fHistPtITSMIoneSPDthreeSDDSSDInAcc,fHistPtTPCInAcc,1,1,"B");
716 //fHistPtITSMIoneSPDthreeSDDSSDInAcc->SetLineColor(kOrange);
717 //l3->AddEntry(fHistPtITSMIoneSPDthreeSDDSSDInAcc,">=1SPD + >=3SDDSSD","l");
718 //fHistPtITSMIoneSPDthreeSDDSSDInAcc->Draw("same");
719 fHistPtITSTPCsel->Divide(fHistPtITSTPCsel,fHistPtTPCInAcc,1,1,"B");
fa516524 720 fHistPtITSTPCsel->SetLineColor(kAzure+1);
721 l3->AddEntry(fHistPtITSTPCsel,">=1SPD + any + d_{0} cut","l");
722 fHistPtITSTPCsel->Draw("same");
c2893e82 723 fHistPtITSMIge2InAcc->Draw("same");
724 l3->Draw();
b51e971f 725 if(ioValues) {
726 Int_t index,ptbin;
727 ptbin=fHistPtITSMIge2InAcc->FindBin(0.201);
728 index=2;
729 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
730 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
731 ptbin=fHistPtITSMIge2InAcc->FindBin(1.001);
732 index=3;
733 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
734 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
735 ptbin=fHistPtITSMIge2InAcc->FindBin(10.001);
736 index=4;
737 ioValues[index]=fHistPtITSMIge2InAcc->GetBinContent(ptbin);
738 ioErrors[index]=fHistPtITSMIge2InAcc->GetBinError(ptbin);
739 ptbin=fHistPtITSMI6InAcc->FindBin(0.201);
740 index=5;
741 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
742 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
743 ptbin=fHistPtITSMI6InAcc->FindBin(1.001);
744 index=6;
745 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
746 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
747 ptbin=fHistPtITSMI6InAcc->FindBin(10.001);
748 index=7;
749 ioValues[index]=fHistPtITSMI6InAcc->GetBinContent(ptbin);
750 ioErrors[index]=fHistPtITSMI6InAcc->GetBinError(ptbin);
751 ptbin=fHistPtITSMISPDInAcc->FindBin(0.201);
752 index=8;
753 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
754 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
755 ptbin=fHistPtITSMISPDInAcc->FindBin(1.001);
756 index=9;
757 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
758 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
759 ptbin=fHistPtITSMISPDInAcc->FindBin(10.001);
760 index=10;
761 ioValues[index]=fHistPtITSMISPDInAcc->GetBinContent(ptbin);
762 ioErrors[index]=fHistPtITSMISPDInAcc->GetBinError(ptbin);
763 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(0.201);
764 index=11;
765 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
766 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
767 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(1.001);
768 index=12;
769 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
770 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
771 ptbin=fHistPtITSMIoneSPDInAcc->FindBin(10.001);
772 index=13;
773 ioValues[index]=fHistPtITSMIoneSPDInAcc->GetBinContent(ptbin);
774 ioErrors[index]=fHistPtITSMIoneSPDInAcc->GetBinError(ptbin);
775 }
c2893e82 776 c5->cd(2);
777 TH1F *fHistPtITSMIokbadoutinzge4InAcc = (TH1F*)fHistPtITSMIokbadoutinz6InAcc->Clone("fHistPtITSMIokbadoutinzge4InAcc");
778 fHistPtITSMIokbadoutinzge4InAcc->Add(fHistPtITSMIokbadoutinz5InAcc);
779 fHistPtITSMIokbadoutinzge4InAcc->Add(fHistPtITSMIokbadoutinz4InAcc);
780 fHistPtITSMIokbadoutinzge4InAcc->SetMaximum(1.5);
781 fHistPtITSMIokbadoutinzge4InAcc->SetMinimum(0);
782 fHistPtITSMIokbadoutinzge4InAcc->SetTitle("Fraction of prolonged tracks with N ITS layers \"ok\"");
783 fHistPtITSMIokbadoutinzge4InAcc->SetYTitle("ITS+TPC / TPC");
8fa56655 784 fHistPtITSMIokbadoutinzge4InAcc->Divide(fHistPtITSMIokbadoutinzge4InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 785 fHistPtITSMIokbadoutinzge4InAcc->SetLineColor(1);
786 fHistPtITSMIokbadoutinzge4InAcc->Draw();
8fa56655 787 fHistPtITSMIokbadoutinz6InAcc->Divide(fHistPtITSMIokbadoutinz6InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 788 fHistPtITSMIokbadoutinz6InAcc->SetLineColor(2);
789 fHistPtITSMIokbadoutinz6InAcc->Draw("same");
8fa56655 790 fHistPtITSMIokbadoutinz5InAcc->Divide(fHistPtITSMIokbadoutinz5InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 791 fHistPtITSMIokbadoutinz5InAcc->SetLineColor(3);
792 fHistPtITSMIokbadoutinz5InAcc->Draw("same");
8fa56655 793 fHistPtITSMIokbadoutinz4InAcc->Divide(fHistPtITSMIokbadoutinz4InAcc,fHistPtTPCInAcc,1,1,"B");
c2893e82 794 fHistPtITSMIokbadoutinz4InAcc->SetLineColor(4);
795 fHistPtITSMIokbadoutinz4InAcc->Draw("same");
796 fHistPtITSMIokbadoutinzge4InAcc->Draw("same");
797 l4->AddEntry(fHistPtITSMIokbadoutinzge4InAcc,">=4 layers","l");
798 l4->AddEntry(fHistPtITSMIokbadoutinz6InAcc,"6 layers","l");
799 l4->AddEntry(fHistPtITSMIokbadoutinz5InAcc,"5 layers","l");
800 l4->AddEntry(fHistPtITSMIokbadoutinz4InAcc,"4 layers","l");
801 l4->Draw();
802
af874c53 803 if(fHistPtITSMIge2InAccFake) {
804 TCanvas *c5f =new TCanvas("c5f","c5f",10,10,600,600);
805 c5f->SetGridy();
806 c5f->SetLogx();
807 fHistPtITSMIge2InAccFake->SetMaximum(1.5);
808 fHistPtITSMIge2InAccFake->SetMinimum(0);
809 fHistPtITSMIge2InAccFake->SetTitle("Fraction of fake tracks with N ITS points");
810 fHistPtITSMIge2InAccFake->SetYTitle("Fraction of fakes");
811 fHistPtITSMIge2InAccFake->Draw();
812 fHistPtITSMI6InAccFake->SetLineColor(2);
813 fHistPtITSMI6InAccFake->Draw("same");
814 fHistPtITSMI5InAccFake->SetLineColor(3);
815 //fHistPtITSMI5InAccFake->Draw("same");
816 fHistPtITSMI4InAccFake->SetLineColor(4);
817 //fHistPtITSMI4InAccFake->Draw("same");
818 fHistPtITSMI3InAccFake->SetLineColor(6);
819 //fHistPtITSMI3InAccFake->Draw("same");
820 fHistPtITSMI2InAccFake->SetLineColor(7);
821 //fHistPtITSMI2InAccFake->Draw("same");
822 fHistPtITSMISPDInAccFake->SetLineColor(9);
823 fHistPtITSMISPDInAccFake->Draw("same");
824 fHistPtITSMIoneSPDInAccFake->SetLineColor(15);
825 fHistPtITSMIoneSPDInAccFake->Draw("same");
826 if(fHistPtITSTPCselFake) {
827 fHistPtITSTPCselFake->SetLineColor(kAzure+1);
828 fHistPtITSTPCselFake->Draw("same");
829 }
830 fHistPtITSMIge2InAccFake->Draw("same");
831 l3->Draw();
832 }
833
c2893e82 834
fa516524 835 TLegend *l4d=new TLegend(0.5,0.5,0.9,0.9);
836 TCanvas *c5d =new TCanvas("c5d","c5d",10,10,600,600);
837 c5d->SetGridy();
838 c5d->SetGridy();
839 fHistPtITSMISPDInAccMC->SetMaximum(1.5);
840 fHistPtITSMISPDInAccMC->SetMinimum(0);
841 fHistPtITSMISPDInAccMC->SetTitle("Fraction of prolonged tracks with N ITS points");
842 fHistPtITSMISPDInAccMC->SetYTitle("ITS+TPC / TPC");
8fa56655 843 fHistPtITSMISPDInAccMC->Divide(fHistPtITSMISPDInAccMC,fHistPtTPCInAccMCtwoSPD,1,1,"B");
fa516524 844 fHistPtITSMISPDInAccMC->SetLineColor(9);
845 l4d->AddEntry(fHistPtITSMISPDInAccMC,"2SPD + any","l");
846 fHistPtITSMISPDInAccMC->Draw();
8fa56655 847 fHistPtITSMIoneSPDInAccMC->Divide(fHistPtITSMIoneSPDInAccMC,fHistPtTPCInAccMConeSPD,1,1,"B");
fa516524 848 fHistPtITSMIoneSPDInAccMC->SetLineColor(15);
849 l4d->AddEntry(fHistPtITSMIoneSPDInAccMC,">=1SPD + any","l");
850 fHistPtITSMIoneSPDInAccMC->Draw("same");
851 l4d->Draw();
852
853 TFile *eff=new TFile("eff.root","recreate");
854 fHistPtITSMIge2InAcc->Write();
855 fHistPtITSMI6InAcc->Write();
856 fHistPtITSMI5InAcc->Write();
857 fHistPtITSMI4InAcc->Write();
858 fHistPtITSMI3InAcc->Write();
859 fHistPtITSMI2InAcc->Write();
860 fHistPtITSMISPDInAcc->Write();
861 fHistPtITSMIoneSPDInAcc->Write();
8fa56655 862 //fHistPtITSMIoneSPDthreeSDDSSDInAcc->Write();
fa516524 863 fHistPtITSTPCsel->Write();
8fa56655 864 fHistNclsITSMI->Write();
fa516524 865 fHistClusterMapModuleITSMIokRatioInAcc->Write();
866 fHistClusterMapModuleITSMIbadRatioInAcc->Write();
867 fHistClusterMapModuleITSMInoclsRatioInAcc->Write();
868 eff->Close();
869
870 c5c->cd();
871 l3->Draw();
872
873
874 TCanvas *c5b =new TCanvas("c5b","c5b",10,10,600,600);
875 c5b->SetLogy();
876 fHistRProdVtxInAccP->SetLineColor(2);
877 fHistRProdVtxInAccP->SetMinimum(1);
878 fHistRProdVtxInAccS->SetMinimum(1);
879 fHistRProdVtxInAccP->Draw();
880 fHistRProdVtxInAccS->Draw("sames");
881
882 TCanvas *c5a =new TCanvas("c5a","c5a",10,10,1200,600);
883 c5a->Divide(2,1);
8fa56655 884 c5a_1->SetLogx();
885 c5a_2->SetLogx();
fa516524 886 c5a_1->SetGridy();
887 c5a_2->SetGridy();
888 c5a->cd(1);
889 fHistPtITSMIge2InAccP->SetMaximum(1.5);
890 fHistPtITSMIge2InAccP->SetMinimum(0);
891 fHistPtITSMIge2InAccP->SetTitle("Fraction of prolonged tracks with N ITS points");
892 fHistPtITSMIge2InAccP->SetYTitle("ITS+TPC / TPC");
893 fHistPtITSMIge2InAccP->Draw();
8fa56655 894 fHistPtITSMIge2InAccP->Divide(fHistPtITSMIge2InAccP,fHistPtTPCInAccP,1,1,"B");
895 fHistPtITSMI6InAccP->Divide(fHistPtITSMI6InAccP,fHistPtTPCInAccP,1,1,"B");
fa516524 896 fHistPtITSMI6InAccP->SetLineColor(2);
897 fHistPtITSMI6InAccP->Draw("same");
8fa56655 898 fHistPtITSMI5InAccP->Divide(fHistPtITSMI5InAccP,fHistPtTPCInAccP,1,1,"B");
fa516524 899 fHistPtITSMI5InAccP->SetLineColor(3);
900 fHistPtITSMI5InAccP->Draw("same");
8fa56655 901 fHistPtITSMI4InAccP->Divide(fHistPtITSMI4InAccP,fHistPtTPCInAccP,1,1,"B");
fa516524 902 fHistPtITSMI4InAccP->SetLineColor(4);
903 fHistPtITSMI4InAccP->Draw("same");
8fa56655 904 fHistPtITSMI3InAccP->Divide(fHistPtITSMI3InAccP,fHistPtTPCInAccP,1,1,"B");
fa516524 905 fHistPtITSMI3InAccP->SetLineColor(6);
906 fHistPtITSMI3InAccP->Draw("same");
8fa56655 907 fHistPtITSMI2InAccP->Divide(fHistPtITSMI2InAccP,fHistPtTPCInAccP,1,1,"B");
fa516524 908 fHistPtITSMI2InAccP->SetLineColor(7);
909 fHistPtITSMI2InAccP->Draw("same");
8fa56655 910 fHistPtITSMISPDInAccP->Divide(fHistPtITSMISPDInAccP,fHistPtTPCInAccP,1,1,"B");
fa516524 911 fHistPtITSMISPDInAccP->SetLineColor(9);
912 fHistPtITSMISPDInAccP->Draw("same");
8fa56655 913 fHistPtITSMIoneSPDInAccP->Divide(fHistPtITSMIoneSPDInAccP,fHistPtTPCInAccP,1,1,"B");
fa516524 914 fHistPtITSMIoneSPDInAccP->SetLineColor(15);
915 fHistPtITSMIoneSPDInAccP->Draw("same");
8fa56655 916 fHistPtITSTPCselP->Divide(fHistPtITSTPCselP,fHistPtTPCInAccP,1,1,"B");
fa516524 917 fHistPtITSTPCselP->SetLineColor(kAzure+1);
918 fHistPtITSTPCselP->Draw("same");
919 fHistPtITSMIge2InAccP->Draw("same");
920 l3->Draw();
921 c5a->cd(2);
922 fHistPtITSMIge2InAccS->SetMaximum(1.5);
923 fHistPtITSMIge2InAccS->SetMinimum(0);
924 fHistPtITSMIge2InAccS->SetTitle("Fraction of prolonged tracks with N ITS points");
925 fHistPtITSMIge2InAccS->SetYTitle("ITS+TPC / TPC");
8fa56655 926 fHistPtITSMIge2InAccS->Divide(fHistPtITSMIge2InAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 927 fHistPtITSMIge2InAccS->Draw();
8fa56655 928 fHistPtITSMI6InAccS->Divide(fHistPtITSMI6InAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 929 fHistPtITSMI6InAccS->SetLineColor(2);
930 fHistPtITSMI6InAccS->Draw("same");
8fa56655 931 fHistPtITSMI5InAccS->Divide(fHistPtITSMI5InAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 932 fHistPtITSMI5InAccS->SetLineColor(3);
933 fHistPtITSMI5InAccS->Draw("same");
8fa56655 934 fHistPtITSMI4InAccS->Divide(fHistPtITSMI4InAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 935 fHistPtITSMI4InAccS->SetLineColor(4);
936 fHistPtITSMI4InAccS->Draw("same");
8fa56655 937 fHistPtITSMI3InAccS->Divide(fHistPtITSMI3InAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 938 fHistPtITSMI3InAccS->SetLineColor(6);
939 fHistPtITSMI3InAccS->Draw("same");
8fa56655 940 fHistPtITSMI2InAccS->Divide(fHistPtITSMI2InAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 941 fHistPtITSMI2InAccS->SetLineColor(7);
942 fHistPtITSMI2InAccS->Draw("same");
8fa56655 943 fHistPtITSMISPDInAccS->Divide(fHistPtITSMISPDInAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 944 fHistPtITSMISPDInAccS->SetLineColor(9);
945 fHistPtITSMISPDInAccS->Draw("same");
8fa56655 946 fHistPtITSMIoneSPDInAccS->Divide(fHistPtITSMIoneSPDInAccS,fHistPtTPCInAccS,1,1,"B");
fa516524 947 fHistPtITSMIoneSPDInAccS->SetLineColor(15);
948 fHistPtITSMIoneSPDInAccS->Draw("same");
8fa56655 949 fHistPtITSTPCselS->Divide(fHistPtITSTPCselS,fHistPtTPCInAccS,1,1,"B");
fa516524 950 fHistPtITSTPCselS->SetLineColor(kAzure+1);
951 fHistPtITSTPCselS->Draw("same");
952 fHistPtITSMIge2InAccS->Draw("same");
953 l3->Draw();
954
955
b51e971f 956 if(!plotAlignmentChecks) return kTRUE;
fa516524 957
c2893e82 958 // PLOT ALIGNMENT CHECKS
959 //
960 TH1F *hSPDTrackletsTBdxy = new TH1F("hSPDTrackletsTBdxy","SPD tracklets; SPD tracklet to SPD vertex distance in (x,y) [cm]; entries",100,-0.1,0.1);
961 TH1F *hSPDTrackletsLRdxy = new TH1F("hSPDTrackletsLRdxy","SPD tracklets; SPD tracklet to SPD vertex distance in (x,y) [cm]; entries",100,-0.1,0.1);
962 TH1F *hSPDTrackletsTBdz = new TH1F("hSPDTrackletsTBdz","SPD tracklets; SPD tracklet to SPD vertex distance in z [cm]; entries",100,-0.1,0.1);
963 TH1F *hSPDTrackletsLRdz = new TH1F("hSPDTrackletsLRdz","SPD tracklets; SPD tracklet to SPD vertex distance in z [cm]; entries",100,-0.1,0.1);
964
965 TNtuple *fNtupleITSAlignSPDTracklets = (TNtuple*)list->FindObject("fNtupleITSAlignSPDTracklets");
966 Float_t dxy,dz,phi,pt;
967 fNtupleITSAlignSPDTracklets->SetBranchAddress("pt",&pt);
968 fNtupleITSAlignSPDTracklets->SetBranchAddress("phi",&phi);
969 fNtupleITSAlignSPDTracklets->SetBranchAddress("dxy",&dxy);
970 fNtupleITSAlignSPDTracklets->SetBranchAddress("dz",&dz);
971
972 for(Int_t i=0;i<fNtupleITSAlignSPDTracklets->GetEntries();i++) {
973 fNtupleITSAlignSPDTracklets->GetEvent(i);
fa516524 974 if(pt<.01) continue;
c2893e82 975 if(TMath::Abs(TMath::Abs(phi)-0.5*TMath::Pi())<0.25*TMath::Pi()) {
976 // top-bottom
977 hSPDTrackletsTBdxy->Fill(dxy);
978 hSPDTrackletsTBdz->Fill(dz);
979 } else {
980 // left-right
981 hSPDTrackletsLRdxy->Fill(dxy);
982 hSPDTrackletsLRdz->Fill(dz);
983 }
984 }
985
986 TLegend *l6=new TLegend(0.5,0.5,0.9,0.9);
987
988 TCanvas *c6 = new TCanvas("c6","c6",0,0,1000,500);
989 c6->Divide(2,1);
990 c6->cd(1);
991 hSPDTrackletsTBdxy->SetLineColor(4);
992 hSPDTrackletsTBdxy->Draw();
993 hSPDTrackletsLRdxy->SetLineColor(2);
fa516524 994 hSPDTrackletsLRdxy->Draw("sames");
c2893e82 995 l6->AddEntry(hSPDTrackletsTBdxy,"top-bottom","l");
996 l6->AddEntry(hSPDTrackletsLRdxy,"left-right","l");
997 l6->Draw();
998 c6->cd(2);
999 hSPDTrackletsTBdz->SetLineColor(4);
1000 hSPDTrackletsTBdz->Draw();
1001 hSPDTrackletsLRdz->SetLineColor(2);
fa516524 1002 hSPDTrackletsLRdz->Draw("sames");
c2893e82 1003 l6->Draw();
1004
1005
1006 TH1F *hSPDExtraClsTBdxy = new TH1F("hSPDExtraClsTBdxy","SPD extra clusters; track-to-point distance in (x,y) [cm]; entries",100,-0.1,0.1);
1007 TH1F *hSPDExtraClsLRdxy = new TH1F("hSPDExtraClsLRdxy","SPD extra clusters; track-to-point distance in (x,y) [cm]; entries",100,-0.1,0.1);
1008 TH1F *hSPDExtraClsTBdz = new TH1F("hSPDExtraClsTBdz","SPD extra clusters; track-to-point distance in z [cm]; entries",100,-0.1,0.1);
1009 TH1F *hSPDExtraClsLRdz = new TH1F("hSPDExtraClsLRdz","SPD extra clusters; track-to-point distance in z [cm]; entries",100,-0.1,0.1);
1010
1011 TNtuple *fNtupleITSAlignExtra = (TNtuple*)list->FindObject("fNtupleITSAlignExtra");
1012 Float_t layer,npoints,x,y,z;
1013 fNtupleITSAlignExtra->SetBranchAddress("layer",&layer);
1014 fNtupleITSAlignExtra->SetBranchAddress("npoints",&npoints);
1015 fNtupleITSAlignExtra->SetBranchAddress("x",&x);
1016 fNtupleITSAlignExtra->SetBranchAddress("y",&y);
1017 fNtupleITSAlignExtra->SetBranchAddress("z",&z);
1018 fNtupleITSAlignExtra->SetBranchAddress("dxy",&dxy);
1019 fNtupleITSAlignExtra->SetBranchAddress("dz",&dz);
1020 fNtupleITSAlignExtra->SetBranchAddress("pt",&pt);
1021
1022 for(Int_t i=0;i<fNtupleITSAlignExtra->GetEntries();i++) {
1023 fNtupleITSAlignExtra->GetEvent(i);
fa516524 1024 if(pt<0.1) continue;
c2893e82 1025 if(layer!=0 && layer!=1) continue;
fa516524 1026 if(npoints<3) continue;
c2893e82 1027 phi = TMath::ATan2(y,x);
1028 if(TMath::Abs(TMath::Abs(phi)-0.5*TMath::Pi())<0.25*TMath::Pi()) {
1029 // top-bottom
1030 hSPDExtraClsTBdxy->Fill(dxy);
1031 hSPDExtraClsTBdz->Fill(dz);
1032 } else {
1033 // left-right
1034 hSPDExtraClsLRdxy->Fill(dxy);
1035 hSPDExtraClsLRdz->Fill(dz);
1036 }
1037 }
1038
1039 TCanvas *c7 = new TCanvas("c7","c7",0,0,1000,500);
1040 c7->Divide(2,1);
1041 c7->cd(1);
1042 hSPDExtraClsTBdxy->SetLineColor(4);
1043 hSPDExtraClsTBdxy->Draw();
1044 hSPDExtraClsLRdxy->SetLineColor(2);
1045 hSPDExtraClsLRdxy->Draw("same");
1046 l6->Draw();
1047 c7->cd(2);
1048 hSPDExtraClsTBdz->SetLineColor(4);
1049 hSPDExtraClsTBdz->Draw();
1050 hSPDExtraClsLRdz->SetLineColor(2);
1051 hSPDExtraClsLRdz->Draw("same");
1052 l6->Draw();
1053
1054
b51e971f 1055 return kTRUE;
fa516524 1056}
1057//-----------------------------------------------------------------------------
1058void PlotEffRatio() {
1059
e75c9f34 1060 TFile *f= new TFile("eff.root");
fa516524 1061
1062 TH1F *fHistPtITSMI6InAcc = (TH1F*)f->Get("fHistPtITSMI6InAcc");
1063 TH1F *fHistPtITSMI5InAcc = (TH1F*)f->Get("fHistPtITSMI5InAcc");
1064 TH1F *fHistPtITSMI4InAcc = (TH1F*)f->Get("fHistPtITSMI4InAcc");
1065 TH1F *fHistPtITSMI3InAcc = (TH1F*)f->Get("fHistPtITSMI3InAcc");
1066 TH1F *fHistPtITSMI2InAcc = (TH1F*)f->Get("fHistPtITSMI2InAcc");
1067 TH1F *fHistPtITSMIge2InAcc = (TH1F*)f->Get("fHistPtITSMIge2InAcc");
1068 TH1F *fHistPtITSMISPDInAcc = (TH1F*)f->Get("fHistPtITSMISPDInAcc");
1069 TH1F *fHistPtITSMIoneSPDInAcc = (TH1F*)f->Get("fHistPtITSMIoneSPDInAcc");
8fa56655 1070 TH1F *fHistPtITSMIoneSPDthreeSDDSSDInAcc = (TH1F*)f->Get("fHistPtITSMIoneSPDthreeSDDSSDInAcc");
fa516524 1071 TH1F *fHistPtITSTPCsel = (TH1F*)f->Get("fHistPtITSTPCsel");
1072 TH1F *fHistClusterMapModuleITSMIokRatioInAcc = (TH1F*)f->Get("fHistClusterMapModuleITSMIokRatioInAcc");
8fa56655 1073 TH1F *fHistNclsITSMI = (TH1F*)f->Get("fHistNclsITSMI");
1074
fa516524 1075
e75c9f34 1076 TFile *fMC= new TFile("../MC/MCLHC11b10a/eff.root");
fa516524 1077
1078 TH1F *fHistPtITSMI6InAccMC = (TH1F*)fMC->Get("fHistPtITSMI6InAcc");
1079 TH1F *fHistPtITSMI5InAccMC = (TH1F*)fMC->Get("fHistPtITSMI5InAcc");
1080 TH1F *fHistPtITSMI4InAccMC = (TH1F*)fMC->Get("fHistPtITSMI4InAcc");
1081 TH1F *fHistPtITSMI3InAccMC = (TH1F*)fMC->Get("fHistPtITSMI3InAcc");
1082 TH1F *fHistPtITSMI2InAccMC = (TH1F*)fMC->Get("fHistPtITSMI2InAcc");
1083 TH1F *fHistPtITSMIge2InAccMC = (TH1F*)fMC->Get("fHistPtITSMIge2InAcc");
1084 TH1F *fHistPtITSMISPDInAccMC = (TH1F*)fMC->Get("fHistPtITSMISPDInAcc");
1085 TH1F *fHistPtITSMIoneSPDInAccMC = (TH1F*)fMC->Get("fHistPtITSMIoneSPDInAcc");
8fa56655 1086 TH1F *fHistPtITSMIoneSPDthreeSDDSSDInAccMC = (TH1F*)fMC->Get("fHistPtITSMIoneSPDthreeSDDSSDInAcc");
fa516524 1087 TH1F *fHistPtITSTPCselMC = (TH1F*)fMC->Get("fHistPtITSTPCsel");
1088 TH1F *fHistClusterMapModuleITSMIokRatioInAccMC = (TH1F*)fMC->Get("fHistClusterMapModuleITSMIokRatioInAcc");
8fa56655 1089 TH1F *fHistNclsITSMIMC = (TH1F*)fMC->Get("fHistNclsITSMIMC");
fa516524 1090
8fa56655 1091 /*
1092 TCanvas *c0 = new TCanvas("c0","c0",0,0,600,500);
1093 c0->SetGridy();
1094 fHistNclsITSMI->SetLineColor(4);
1095 fHistNclsITSMI->SetMarkerColor(4);
1096 fHistNclsITSMI->SetLineStyle(1);
1097 fHistNclsITSMI->SetMarkerStyle(20);
1098 fHistNclsITSMI->Draw();
1099 fHistNclsITSMIMC->SetLineColor(4);
1100 fHistNclsITSMIMC->SetMarkerColor(4);
1101 fHistNclsITSMIMC->SetLineStyle(2);
1102 fHistNclsITSMIMC->SetMarkerStyle(24);
1103 fHistNclsITSMIMC->Draw("same");
1104 */
fa516524 1105 TCanvas *c = new TCanvas("c","c",0,0,600,500);
1106 c->SetGridy();
8fa56655 1107 c->SetLogx();
1108
fa516524 1109 fHistPtITSMI6InAcc->Divide(fHistPtITSMI6InAccMC);
1110 fHistPtITSMI6InAcc->Draw();
1111 fHistPtITSMI5InAcc->Divide(fHistPtITSMI5InAccMC);
1112 fHistPtITSMI5InAcc->Draw("same");
1113 fHistPtITSMI4InAcc->Divide(fHistPtITSMI4InAccMC);
1114 fHistPtITSMI4InAcc->Draw("same");
1115 fHistPtITSMI3InAcc->Divide(fHistPtITSMI3InAccMC);
1116 fHistPtITSMI3InAcc->Draw("same");
1117 fHistPtITSMI2InAcc->Divide(fHistPtITSMI2InAccMC);
1118 fHistPtITSMI2InAcc->Draw("same");
1119 fHistPtITSMIge2InAcc->Divide(fHistPtITSMIge2InAccMC);
1120 fHistPtITSMIge2InAcc->Draw("same");
1121 fHistPtITSMISPDInAcc->Divide(fHistPtITSMISPDInAccMC);
1122 fHistPtITSMISPDInAcc->Draw("same");
1123 fHistPtITSMIoneSPDInAcc->Divide(fHistPtITSMIoneSPDInAccMC);
1124 fHistPtITSMIoneSPDInAcc->Draw("same");
8fa56655 1125 //fHistPtITSMIoneSPDthreeSDDSSDInAcc->Divide(fHistPtITSMIoneSPDthreeSDDSSDInAccMC);
1126 //fHistPtITSMIoneSPDthreeSDDSSDInAcc->Draw("same");
1127
fa516524 1128 fHistPtITSTPCsel->Divide(fHistPtITSTPCselMC);
1129 fHistPtITSTPCsel->Draw("same");
8fa56655 1130
fa516524 1131 TLegend *l3=new TLegend(0.5,0.5,0.9,0.9);
1132 l3->AddEntry(fHistPtITSMIge2InAcc,">=2 cls","l");
1133 l3->AddEntry(fHistPtITSMI6InAcc,"6 cls","l");
1134 l3->AddEntry(fHistPtITSMI5InAcc,"5 cls","l");
1135 l3->AddEntry(fHistPtITSMI4InAcc,"4 cls","l");
1136 l3->AddEntry(fHistPtITSMI3InAcc,"3 cls","l");
1137 l3->AddEntry(fHistPtITSMI2InAcc,"2 cls","l");
1138 l3->AddEntry(fHistPtITSMISPDInAcc,"2SPD + any","l");
1139 l3->AddEntry(fHistPtITSMIoneSPDInAcc,">=1SPD + any","l");
8fa56655 1140 //l3->AddEntry(fHistPtITSMIoneSPDthreeSDDSSDInAcc,">=1SPD + >=3SDDSSD","l");
fa516524 1141 l3->AddEntry(fHistPtITSTPCsel,">=1SPD + any + d_{0} cut","l");
1142 l3->Draw();
8fa56655 1143
fa516524 1144 TCanvas *cc = new TCanvas("cc","cc",0,0,600,500);
1145 cc->Divide(1,2);
1146 cc->cd(1);
1147 fHistClusterMapModuleITSMIokRatioInAccMC->SetLineColor(4);
1148 fHistClusterMapModuleITSMIokRatioInAccMC->Draw();
1149 fHistClusterMapModuleITSMIokRatioInAcc->Draw("same");
1150 cc->cd(2);
1151 for(Int_t i=1;i<=fHistClusterMapModuleITSMIokRatioInAccMC->GetNbinsX();i++) {
1152 if(fHistClusterMapModuleITSMIokRatioInAccMC->GetBinContent(i)<0.02) {
1153 fHistClusterMapModuleITSMIokRatioInAccMC->SetBinContent(i,1);
1154 if(fHistClusterMapModuleITSMIokRatioInAcc->GetBinContent(i)<0.02) {
1155 fHistClusterMapModuleITSMIokRatioInAcc->SetBinContent(i,1);
1156 }
1157 }
1158 }
1159 TH1F* fHistClusterMapDataOverMC=(TH1F*)fHistClusterMapModuleITSMIokRatioInAcc->Clone("fHistClusterMapDataOverMC");
1160 fHistClusterMapDataOverMC->Divide(fHistClusterMapModuleITSMIokRatioInAccMC);
1161 fHistClusterMapDataOverMC->Draw();
1162
1163 return;
1164}
1165
8fa56655 1166//-----------------------------------------------------------------------------
1167void PlotEffOfficial(Bool_t drawRatio=kTRUE) {
1168
1169 gStyle->SetOptStat(0);
1170
1171 gROOT->LoadMacro("/Users/dainesea/ALICEWorkInProgress.C");
1172
1173 TFile *f= new TFile("eff_115328_115393_minipass.root");
1174
1175 TH1F *fHistPtITSMI6InAcc = (TH1F*)f->Get("fHistPtITSMI6InAcc");
1176 TH1F *fHistPtITSMI5InAcc = (TH1F*)f->Get("fHistPtITSMI5InAcc");
1177 TH1F *fHistPtITSMI4InAcc = (TH1F*)f->Get("fHistPtITSMI4InAcc");
1178 TH1F *fHistPtITSMI3InAcc = (TH1F*)f->Get("fHistPtITSMI3InAcc");
1179 TH1F *fHistPtITSMI2InAcc = (TH1F*)f->Get("fHistPtITSMI2InAcc");
1180 TH1F *fHistPtITSMIge2InAcc = (TH1F*)f->Get("fHistPtITSMIge2InAcc");
1181 fHistPtITSMIge2InAcc->SetLineColor(1);
1182 fHistPtITSMIge2InAcc->SetLineStyle(1);
1183 fHistPtITSMIge2InAcc->SetMarkerColor(1);
1184 fHistPtITSMIge2InAcc->SetMarkerStyle(21);
1185 TH1F *fHistPtITSMISPDInAcc = (TH1F*)f->Get("fHistPtITSMISPDInAcc");
1186 TH1F *fHistPtITSMIoneSPDInAcc = (TH1F*)f->Get("fHistPtITSMIoneSPDInAcc");
1187 fHistPtITSMIoneSPDInAcc->SetLineColor(2);
1188 fHistPtITSMIoneSPDInAcc->SetLineStyle(1);
1189 fHistPtITSMIoneSPDInAcc->SetMarkerColor(2);
1190 fHistPtITSMIoneSPDInAcc->SetMarkerStyle(20);
1191 TH1F *fHistPtITSTPCsel = (TH1F*)f->Get("fHistPtITSTPCsel");
1192 TH1F *fHistClusterMapModuleITSMIokRatioInAcc = (TH1F*)f->Get("fHistClusterMapModuleITSMIokRatioInAcc");
1193 TH1F *fHistNclsITSMI = (TH1F*)f->Get("fHistNclsITSMI");
1194
1195
1196 TFile *fMC= new TFile("eff_lhc10b2.root");
1197
1198 TH1F *fHistPtITSMI6InAccMC = (TH1F*)fMC->Get("fHistPtITSMI6InAcc");
1199 TH1F *fHistPtITSMI5InAccMC = (TH1F*)fMC->Get("fHistPtITSMI5InAcc");
1200 TH1F *fHistPtITSMI4InAccMC = (TH1F*)fMC->Get("fHistPtITSMI4InAcc");
1201 TH1F *fHistPtITSMI3InAccMC = (TH1F*)fMC->Get("fHistPtITSMI3InAcc");
1202 TH1F *fHistPtITSMI2InAccMC = (TH1F*)fMC->Get("fHistPtITSMI2InAcc");
1203 TH1F *fHistPtITSMIge2InAccMC = (TH1F*)fMC->Get("fHistPtITSMIge2InAcc");
1204 fHistPtITSMIge2InAccMC->SetLineColor(1);
1205 fHistPtITSMIge2InAccMC->SetLineStyle(2);
1206 fHistPtITSMIge2InAccMC->SetMarkerColor(1);
1207 fHistPtITSMIge2InAccMC->SetMarkerStyle(25);
1208 TH1F *fHistPtITSMISPDInAccMC = (TH1F*)fMC->Get("fHistPtITSMISPDInAcc");
1209 TH1F *fHistPtITSMIoneSPDInAccMC = (TH1F*)fMC->Get("fHistPtITSMIoneSPDInAcc");
1210 fHistPtITSMIoneSPDInAccMC->SetLineColor(2);
1211 fHistPtITSMIoneSPDInAccMC->SetLineStyle(2);
1212 fHistPtITSMIoneSPDInAccMC->SetMarkerColor(2);
1213 fHistPtITSMIoneSPDInAccMC->SetMarkerStyle(24);
1214 TH1F *fHistPtITSTPCselMC = (TH1F*)fMC->Get("fHistPtITSTPCsel");
1215 TH1F *fHistClusterMapModuleITSMIokRatioInAccMC = (TH1F*)fMC->Get("fHistClusterMapModuleITSMIokRatioInAcc");
1216 TH1F *fHistNclsITSMIMC = (TH1F*)fMC->Get("fHistNclsITSMI");
1217
1218
1219 TCanvas *c0 = new TCanvas("c0","c0",0,0,600,500);
1220 c0->SetGridy();
1221 fHistNclsITSMI->SetLineColor(4);
1222 fHistNclsITSMI->SetMarkerColor(4);
1223 fHistNclsITSMI->SetLineStyle(1);
1224 fHistNclsITSMI->SetMarkerStyle(20);
1225 fHistNclsITSMI->Scale(1./fHistNclsITSMI->Integral());
1226 fHistNclsITSMI->Draw();
1227 fHistNclsITSMIMC->SetLineColor(4);
1228 fHistNclsITSMIMC->SetMarkerColor(4);
1229 fHistNclsITSMIMC->SetLineStyle(2);
1230 fHistNclsITSMIMC->SetMarkerStyle(24);
1231 fHistNclsITSMIMC->Scale(1./fHistNclsITSMIMC->Integral());
1232 fHistNclsITSMIMC->Draw("same");
1233
1234 TLegend *l0=new TLegend(0.5,0.5,0.9,0.9);
1235 l0->SetFillStyle(0);
1236 l0->SetBorderSize(0);
1237 l0->AddEntry(fHistNclsITSMI,"Data","lp");
1238 l0->AddEntry(fHistNclsITSMIMC,"MC","lp");
1239 l0->Draw();
1240
1241 ALICEWorkInProgress(c0,"20/05/2010");
1242
1243 TCanvas *c = new TCanvas("c","c",0,0,600,500);
1244 c->SetGridy();
1245 c->SetLogx();
1246
1247 TLegend *l3=new TLegend(0.5,0.5,0.9,0.9);
1248 l3->SetFillStyle(0);
1249 l3->SetBorderSize(0);
1250
1251 if(drawRatio) {
1252 //fHistPtITSMI6InAcc->Divide(fHistPtITSMI6InAccMC);
1253 //fHistPtITSMI6InAcc->Draw();
1254 //fHistPtITSMI5InAcc->Divide(fHistPtITSMI5InAccMC);
1255 //fHistPtITSMI5InAcc->Draw("same");
1256 //fHistPtITSMI4InAcc->Divide(fHistPtITSMI4InAccMC);
1257 //fHistPtITSMI4InAcc->Draw("same");
1258 //fHistPtITSMI3InAcc->Divide(fHistPtITSMI3InAccMC);
1259 //fHistPtITSMI3InAcc->Draw("same");
1260 //fHistPtITSMI2InAcc->Divide(fHistPtITSMI2InAccMC);
1261 //fHistPtITSMI2InAcc->Draw("same");
1262 fHistPtITSMIge2InAcc->SetYTitle("Data efficiency / MC efficiency");
1263 fHistPtITSMIge2InAcc->Divide(fHistPtITSMIge2InAccMC);
1264 fHistPtITSMIge2InAcc->Draw();
1265 //fHistPtITSMISPDInAcc->Divide(fHistPtITSMISPDInAccMC);
1266 //fHistPtITSMISPDInAcc->Draw("same");
1267 fHistPtITSMIoneSPDInAcc->Divide(fHistPtITSMIoneSPDInAccMC);
1268 fHistPtITSMIoneSPDInAcc->Draw("same");
1269
1270 //fHistPtITSTPCsel->Divide(fHistPtITSTPCselMC);
1271 //fHistPtITSTPCsel->Draw("same");
1272 l3->AddEntry(fHistPtITSMIge2InAcc,"at least 2 ITS hits","lp");
1273 //l3->AddEntry(fHistPtITSMI6InAcc,"6 cls","l");
1274 //l3->AddEntry(fHistPtITSMI5InAcc,"5 cls","l");
1275 //l3->AddEntry(fHistPtITSMI4InAcc,"4 cls","l");
1276 //l3->AddEntry(fHistPtITSMI3InAcc,"3 cls","l");
1277 //l3->AddEntry(fHistPtITSMI2InAcc,"2 cls","l");
1278 //l3->AddEntry(fHistPtITSMISPDInAcc,"2 SPD hits","lp");
1279 l3->AddEntry(fHistPtITSMIoneSPDInAcc,"at least 1 SPD hit","lp");
1280 //l3->AddEntry(fHistPtITSTPCsel,">=1SPD + any + d_{0} cut","l");
1281
1282 } else {
1283 fHistPtITSMIge2InAcc->SetYTitle("ITS prolongation efficiency");
1284 fHistPtITSMIge2InAcc->Draw();
1285 fHistPtITSMIge2InAccMC->Draw("same");
1286 fHistPtITSMIoneSPDInAcc->Draw("same");
1287 fHistPtITSMIoneSPDInAccMC->Draw("same");
1288 l3->AddEntry(fHistPtITSMIge2InAcc,"at least 2 ITS hits (Data)","lp");
1289 l3->AddEntry(fHistPtITSMIoneSPDInAcc,"at least 1 SPD hit (Data)","lp");
1290 l3->AddEntry(fHistPtITSMIge2InAccMC,"at least 2 ITS hits (MC)","lp");
1291 l3->AddEntry(fHistPtITSMIoneSPDInAccMC,"at least 1 SPD hit (MC)","lp");
1292 }
1293
1294 l3->Draw();
1295
1296 ALICEWorkInProgress(c,"20/05/2010");
1297
1298 return;
1299}
1300
fa516524 1301void PlotImpPar_rphi(Int_t rebin=1) {
1302
1303
71f9d353 1304 TFile *fMC= new TFile("AnalysisResults_onlynonfakes.root");
fa516524 1305
1306 TList *list=(TList*)fMC->Get("cOutputITS");
71f9d353 1307 TDirectoryFile *dir=0;
1308 if(!list) {
1309 dir=(TDirectoryFile*)fMC->GetDirectory("ITS_Performance");
1310 if(dir) list = (TList*)dir->Get("cOutputITS");
1311 }
1312
1313 if(!list) return kFALSE;
1314
fa516524 1315 TH1F *fHistd0rphiITSMIoneSPDInAccP150200MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
1316 TH1F *fHistd0rphiITSMIoneSPDInAccS150200MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
1317 TH1F *fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200fromStrange");
1318 TH1F *fHistd0rphiITSMIoneSPDInAccS150200fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200fromMat");
8fa56655 1319 TH1F *fHistd0rphiITSMIoneSPDInAccP350450MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP350450");
fa516524 1320 TH1F *fHistd0rphiITSMIoneSPDInAccS350450MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450");
1321 TH1F *fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450fromStrange");
8fa56655 1322 TH1F *fHistd0rphiITSMIoneSPDInAccS350450fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450fromMat");
fa516524 1323 TH1F *fHistd0rphiITSMIoneSPDInAccP500700MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP500700");
1324 TH1F *fHistd0rphiITSMIoneSPDInAccS500700MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700");
1325 TH1F *fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700fromStrange");
1326 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from211MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from211");
1327 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from22MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from22");
1328 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from310MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from310");
1329 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from321MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from321");
1330 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from3122MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from3122");
1331 TH1F *fHistd0rphiITSMIoneSPDInAccS500700fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700fromMat");
1332 TH1F *fHistd0rphiITSMIoneSPDInAccP10001500MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP10001500");
1333 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500");
1334 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500fromStrange");
1335 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500fromMat");
8fa56655 1336 TH1F *fHistd0rphiITSMIoneSPDInAccP25004000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP25004000");
fa516524 1337 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000");
1338 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000fromStrange");
1339 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000fromMat");
1340 TH1F *fHistd0rphiITSMIoneSPDInAccP40008000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP40008000");
1341 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000");
1342 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000fromStrange");
8fa56655 1343 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000fromMat");
fa516524 1344
1345 TH1F *fHistd0rphiITSMIoneSPDInAcc150200MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP150200MC->Clone("fHistd0rphiITSMIoneSPDInAcc150200MC");
1346 fHistd0rphiITSMIoneSPDInAcc150200MC->Add(fHistd0rphiITSMIoneSPDInAccS150200MC);
1347 fHistd0rphiITSMIoneSPDInAcc150200MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc150200MC->GetEntries());
8fa56655 1348 TH1F *fHistd0rphiITSMIoneSPDInAcc350450MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP350450MC->Clone("fHistd0rphiITSMIoneSPDInAcc350450MC");
1349 fHistd0rphiITSMIoneSPDInAcc350450MC->Add(fHistd0rphiITSMIoneSPDInAccS350450MC);
1350 fHistd0rphiITSMIoneSPDInAcc350450MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc350450MC->GetEntries());
fa516524 1351 TH1F *fHistd0rphiITSMIoneSPDInAcc500700MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP500700MC->Clone("fHistd0rphiITSMIoneSPDInAcc500700MC");
1352 fHistd0rphiITSMIoneSPDInAcc500700MC->Add(fHistd0rphiITSMIoneSPDInAccS500700MC);
1353 fHistd0rphiITSMIoneSPDInAcc500700MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc500700MC->GetEntries());
1354 TH1F *fHistd0rphiITSMIoneSPDInAcc10001500MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP10001500MC->Clone("fHistd0rphiITSMIoneSPDInAcc10001500MC");
1355 fHistd0rphiITSMIoneSPDInAcc10001500MC->Add(fHistd0rphiITSMIoneSPDInAccS10001500MC);
1356 fHistd0rphiITSMIoneSPDInAcc10001500MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc10001500MC->GetEntries());
8fa56655 1357 TH1F *fHistd0rphiITSMIoneSPDInAcc25004000MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP25004000MC->Clone("fHistd0rphiITSMIoneSPDInAcc25004000MC");
fa516524 1358 fHistd0rphiITSMIoneSPDInAcc25004000MC->Add(fHistd0rphiITSMIoneSPDInAccS25004000MC);
1359 fHistd0rphiITSMIoneSPDInAcc25004000MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc25004000MC->GetEntries());
1360 TH1F *fHistd0rphiITSMIoneSPDInAcc40008000MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP40008000MC->Clone("fHistd0rphiITSMIoneSPDInAcc40008000MC");
1361 fHistd0rphiITSMIoneSPDInAcc40008000MC->Add(fHistd0rphiITSMIoneSPDInAccS40008000MC);
1362 fHistd0rphiITSMIoneSPDInAcc40008000MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc40008000MC->GetEntries());
8fa56655 1363
fa516524 1364
71f9d353 1365 TFile *f= new TFile("AnalysisResults_onlynonfakes.root");
1366 list=(TList*)f->Get("cOutputITS");
1367 if(!list) {
1368 dir=(TDirectoryFile*)f->GetDirectory("ITS_Performance");
1369 if(dir) list = (TList*)dir->Get("cOutputITS");
1370 }
1371
1372 if(!list) return kFALSE;
fa516524 1373
fa516524 1374 TH1F *fHistd0rphiITSMIoneSPDInAccP150200 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
1375 TH1F *fHistd0rphiITSMIoneSPDInAccS150200 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
8fa56655 1376 TH1F *fHistd0rphiITSMIoneSPDInAccP350450 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP350450");
1377 TH1F *fHistd0rphiITSMIoneSPDInAccS350450 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450");
fa516524 1378 TH1F *fHistd0rphiITSMIoneSPDInAccP500700 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP500700");
1379 TH1F *fHistd0rphiITSMIoneSPDInAccS500700 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700");
1380 TH1F *fHistd0rphiITSMIoneSPDInAccP10001500 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP10001500");
1381 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500");
8fa56655 1382 TH1F *fHistd0rphiITSMIoneSPDInAccP25004000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP25004000");
fa516524 1383 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000");
1384 TH1F *fHistd0rphiITSMIoneSPDInAccP40008000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP40008000");
8fa56655 1385 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000");
fa516524 1386 TH1F *fHistd0rphiITSMIoneSPDInAcc150200=(TH1F*)fHistd0rphiITSMIoneSPDInAccP150200->Clone("fHistd0rphiITSMIoneSPDInAcc150200");
1387 fHistd0rphiITSMIoneSPDInAcc150200->Add(fHistd0rphiITSMIoneSPDInAccS150200);
1388 fHistd0rphiITSMIoneSPDInAcc150200->Scale(1./fHistd0rphiITSMIoneSPDInAcc150200->GetEntries());
8fa56655 1389 TH1F *fHistd0rphiITSMIoneSPDInAcc350450=(TH1F*)fHistd0rphiITSMIoneSPDInAccP350450->Clone("fHistd0rphiITSMIoneSPDInAcc350450");
fa516524 1390 fHistd0rphiITSMIoneSPDInAcc350450->Add(fHistd0rphiITSMIoneSPDInAccS350450);
8fa56655 1391 fHistd0rphiITSMIoneSPDInAcc350450->Scale(1./fHistd0rphiITSMIoneSPDInAcc350450->GetEntries());
fa516524 1392 TH1F *fHistd0rphiITSMIoneSPDInAcc500700=(TH1F*)fHistd0rphiITSMIoneSPDInAccP500700->Clone("fHistd0rphiITSMIoneSPDInAcc500700");
1393 fHistd0rphiITSMIoneSPDInAcc500700->Add(fHistd0rphiITSMIoneSPDInAccS500700);
1394 fHistd0rphiITSMIoneSPDInAcc500700->Scale(1./fHistd0rphiITSMIoneSPDInAcc500700->GetEntries());
1395 TH1F *fHistd0rphiITSMIoneSPDInAcc10001500=(TH1F*)fHistd0rphiITSMIoneSPDInAccP10001500->Clone("fHistd0rphiITSMIoneSPDInAcc10001500");
1396 fHistd0rphiITSMIoneSPDInAcc10001500->Add(fHistd0rphiITSMIoneSPDInAccS10001500);
1397 fHistd0rphiITSMIoneSPDInAcc10001500->Scale(1./fHistd0rphiITSMIoneSPDInAcc10001500->GetEntries());
8fa56655 1398 TH1F *fHistd0rphiITSMIoneSPDInAcc25004000=(TH1F*)fHistd0rphiITSMIoneSPDInAccP25004000->Clone("fHistd0rphiITSMIoneSPDInAcc25004000");
fa516524 1399 fHistd0rphiITSMIoneSPDInAcc25004000->Add(fHistd0rphiITSMIoneSPDInAccS25004000);
1400 fHistd0rphiITSMIoneSPDInAcc25004000->Scale(1./fHistd0rphiITSMIoneSPDInAcc25004000->GetEntries());
1401 TH1F *fHistd0rphiITSMIoneSPDInAcc40008000=(TH1F*)fHistd0rphiITSMIoneSPDInAccP40008000->Clone("fHistd0rphiITSMIoneSPDInAcc40008000");
1402 fHistd0rphiITSMIoneSPDInAcc40008000->Add(fHistd0rphiITSMIoneSPDInAccS40008000);
8fa56655 1403 fHistd0rphiITSMIoneSPDInAcc40008000->Scale(1./fHistd0rphiITSMIoneSPDInAcc40008000->GetEntries());
fa516524 1404
1405
1406 TCanvas *c1 = new TCanvas("c1","c1");
1407 c1->Divide(3,2);
1408 c1->cd(1);
1409 fHistd0rphiITSMIoneSPDInAcc150200MC->Rebin(rebin);
1410 fHistd0rphiITSMIoneSPDInAcc150200->Rebin(rebin);
1411 fHistd0rphiITSMIoneSPDInAcc150200MC->SetLineColor(2);
1412 fHistd0rphiITSMIoneSPDInAcc150200MC->Draw();
1413 fHistd0rphiITSMIoneSPDInAcc150200->SetLineColor(4);
1414 fHistd0rphiITSMIoneSPDInAcc150200->Draw("same");
8fa56655 1415 c1->cd(2);
fa516524 1416 fHistd0rphiITSMIoneSPDInAcc350450MC->Rebin(rebin);
1417 fHistd0rphiITSMIoneSPDInAcc350450->Rebin(rebin);
1418 fHistd0rphiITSMIoneSPDInAcc350450MC->SetLineColor(2);
1419 fHistd0rphiITSMIoneSPDInAcc350450MC->Draw();
1420 fHistd0rphiITSMIoneSPDInAcc350450->SetLineColor(4);
8fa56655 1421 fHistd0rphiITSMIoneSPDInAcc350450->Draw("same");
fa516524 1422 c1->cd(3);
1423 fHistd0rphiITSMIoneSPDInAcc500700MC->Rebin(rebin);
1424 fHistd0rphiITSMIoneSPDInAcc500700->Rebin(rebin);
1425 fHistd0rphiITSMIoneSPDInAcc500700MC->SetLineColor(2);
1426 fHistd0rphiITSMIoneSPDInAcc500700MC->Draw();
1427 fHistd0rphiITSMIoneSPDInAcc500700->SetLineColor(4);
1428 fHistd0rphiITSMIoneSPDInAcc500700->Draw("same");
1429 c1->cd(4);
1430 fHistd0rphiITSMIoneSPDInAcc10001500MC->Rebin(rebin);
1431 fHistd0rphiITSMIoneSPDInAcc10001500->Rebin(rebin);
1432 fHistd0rphiITSMIoneSPDInAcc10001500MC->SetLineColor(2);
1433 fHistd0rphiITSMIoneSPDInAcc10001500MC->Draw();
1434 fHistd0rphiITSMIoneSPDInAcc10001500->SetLineColor(4);
1435 fHistd0rphiITSMIoneSPDInAcc10001500->Draw("same");
8fa56655 1436 c1->cd(5);
fa516524 1437 fHistd0rphiITSMIoneSPDInAcc25004000MC->Rebin(rebin);
1438 fHistd0rphiITSMIoneSPDInAcc25004000->Rebin(rebin);
1439 fHistd0rphiITSMIoneSPDInAcc25004000MC->SetLineColor(2);
1440 fHistd0rphiITSMIoneSPDInAcc25004000MC->Draw();
1441 fHistd0rphiITSMIoneSPDInAcc25004000->SetLineColor(4);
1442 fHistd0rphiITSMIoneSPDInAcc25004000->Draw("same");
1443 c1->cd(6);
1444 fHistd0rphiITSMIoneSPDInAcc40008000MC->Rebin(rebin);
1445 fHistd0rphiITSMIoneSPDInAcc40008000->Rebin(rebin);
1446 fHistd0rphiITSMIoneSPDInAcc40008000MC->SetLineColor(2);
1447 fHistd0rphiITSMIoneSPDInAcc40008000MC->Draw();
1448 fHistd0rphiITSMIoneSPDInAcc40008000->SetLineColor(4);
1449 fHistd0rphiITSMIoneSPDInAcc40008000->Draw("same");
8fa56655 1450
1451
fa516524 1452 TCanvas *c1a = new TCanvas("c1a","c1a");
1453 c1a->Divide(3,1);
1454 c1a->cd(1);
1455 TH1F* fHistd0rphiITSMIoneSPDInAcc150200Ratio=(TH1F*)fHistd0rphiITSMIoneSPDInAcc150200->Clone("fHistd0rphiITSMIoneSPDInAcc150200Ratio");
1456 fHistd0rphiITSMIoneSPDInAcc150200Ratio->Divide(fHistd0rphiITSMIoneSPDInAcc150200MC);
1457 fHistd0rphiITSMIoneSPDInAcc150200Ratio->Draw();
1458 c1a->cd(2);
1459 TH1F* fHistd0rphiITSMIoneSPDInAcc500700Ratio=(TH1F*)fHistd0rphiITSMIoneSPDInAcc500700->Clone("fHistd0rphiITSMIoneSPDInAcc500700Ratio");
1460 fHistd0rphiITSMIoneSPDInAcc500700Ratio->Divide(fHistd0rphiITSMIoneSPDInAcc500700MC);
1461 fHistd0rphiITSMIoneSPDInAcc500700Ratio->Draw();
1462 c1a->cd(3);
1463 TH1F* fHistd0rphiITSMIoneSPDInAcc10001500Ratio=(TH1F*)fHistd0rphiITSMIoneSPDInAcc10001500->Clone("fHistd0rphiITSMIoneSPDInAcc10001500Ratio");
1464 fHistd0rphiITSMIoneSPDInAcc10001500Ratio->Divide(fHistd0rphiITSMIoneSPDInAcc10001500MC);
1465 fHistd0rphiITSMIoneSPDInAcc10001500Ratio->Draw();
1466
1467
1468
1469 TCanvas *c2 = new TCanvas("c2","c2");
1470 c2->Divide(3,2);
1471 c2_1->SetLogy();
1472 c2_2->SetLogy();
1473 c2_3->SetLogy();
1474 c2_4->SetLogy();
1475 c2_5->SetLogy();
1476 c2_6->SetLogy();
1477 c2->cd(1);
1478 fHistd0rphiITSMIoneSPDInAccP150200MC->SetMinimum(1);
1479 fHistd0rphiITSMIoneSPDInAccS150200MC->SetMinimum(1);
1480 fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->SetMinimum(1);
1481 fHistd0rphiITSMIoneSPDInAccS150200fromMatMC->SetMinimum(1);
1482 fHistd0rphiITSMIoneSPDInAccP150200MC->SetLineColor(1);
1483 fHistd0rphiITSMIoneSPDInAccP150200MC->Draw();
1484 fHistd0rphiITSMIoneSPDInAccS150200MC->SetLineColor(6);
1485 fHistd0rphiITSMIoneSPDInAccS150200MC->Draw("same");
1486 fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->SetLineColor(8);
1487 fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->Draw("same");
1488 fHistd0rphiITSMIoneSPDInAccS150200fromMatMC->SetLineColor(9);
1489 fHistd0rphiITSMIoneSPDInAccS150200fromMatMC->Draw("same");
8fa56655 1490 c2->cd(2);
fa516524 1491 fHistd0rphiITSMIoneSPDInAccP350450MC->SetMinimum(1);
1492 fHistd0rphiITSMIoneSPDInAccS350450MC->SetMinimum(1);
1493 fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->SetMinimum(1);
1494 fHistd0rphiITSMIoneSPDInAccS350450fromMatMC->SetMinimum(1);
1495 fHistd0rphiITSMIoneSPDInAccP350450MC->SetLineColor(1);
1496 fHistd0rphiITSMIoneSPDInAccP350450MC->Draw();
1497 fHistd0rphiITSMIoneSPDInAccS350450MC->SetLineColor(6);
1498 fHistd0rphiITSMIoneSPDInAccS350450MC->Draw("same");
1499 fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->SetLineColor(8);
1500 fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->Draw("same");
1501 fHistd0rphiITSMIoneSPDInAccS350450fromMatMC->SetLineColor(9);
8fa56655 1502 fHistd0rphiITSMIoneSPDInAccS350450fromMatMC->Draw("same");
fa516524 1503 c2->cd(3);
1504 fHistd0rphiITSMIoneSPDInAccP500700MC->SetMinimum(1);
1505 fHistd0rphiITSMIoneSPDInAccS500700MC->SetMinimum(1);
1506 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetMinimum(1);
1507 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->SetMinimum(1);
1508 fHistd0rphiITSMIoneSPDInAccP500700MC->SetLineColor(1);
1509 fHistd0rphiITSMIoneSPDInAccP500700MC->Draw();
1510 fHistd0rphiITSMIoneSPDInAccS500700MC->SetLineColor(6);
1511 fHistd0rphiITSMIoneSPDInAccS500700MC->Draw("same");
1512 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetLineColor(8);
1513 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->Draw("same");
1514 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->SetLineColor(9);
1515 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->Draw("same");
1516 c2->cd(4);
1517 fHistd0rphiITSMIoneSPDInAccP10001500MC->SetMinimum(1);
1518 fHistd0rphiITSMIoneSPDInAccS10001500MC->SetMinimum(1);
1519 fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->SetMinimum(1);
1520 fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC->SetMinimum(1);
1521 fHistd0rphiITSMIoneSPDInAccP10001500MC->SetLineColor(1);
1522 fHistd0rphiITSMIoneSPDInAccP10001500MC->Draw();
1523 fHistd0rphiITSMIoneSPDInAccS10001500MC->SetLineColor(6);
1524 fHistd0rphiITSMIoneSPDInAccS10001500MC->Draw("same");
1525 fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->SetLineColor(8);
1526 fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->Draw("same");
1527 fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC->SetLineColor(9);
1528 fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC->Draw("same");
fa516524 1529 c2->cd(5);
1530 fHistd0rphiITSMIoneSPDInAccP25004000MC->SetMinimum(1);
1531 fHistd0rphiITSMIoneSPDInAccS25004000MC->SetMinimum(1);
1532 fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->SetMinimum(1);
1533 fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC->SetMinimum(1);
1534 fHistd0rphiITSMIoneSPDInAccP25004000MC->SetLineColor(1);
1535 fHistd0rphiITSMIoneSPDInAccP25004000MC->Draw();
1536 fHistd0rphiITSMIoneSPDInAccS25004000MC->SetLineColor(6);
1537 fHistd0rphiITSMIoneSPDInAccS25004000MC->Draw("same");
1538 fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->SetLineColor(8);
1539 fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->Draw("same");
1540 fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC->SetLineColor(9);
1541 fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC->Draw("same");
1542 c2->cd(6);
1543 fHistd0rphiITSMIoneSPDInAccP40008000MC->SetMinimum(1);
1544 fHistd0rphiITSMIoneSPDInAccS40008000MC->SetMinimum(1);
1545 fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->SetMinimum(1);
1546 fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC->SetMinimum(1);
1547 fHistd0rphiITSMIoneSPDInAccP40008000MC->SetLineColor(1);
1548 fHistd0rphiITSMIoneSPDInAccP40008000MC->Draw();
1549 fHistd0rphiITSMIoneSPDInAccS40008000MC->SetLineColor(6);
1550 fHistd0rphiITSMIoneSPDInAccS40008000MC->Draw("same");
1551 fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->SetLineColor(8);
1552 fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->Draw("same");
1553 fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC->SetLineColor(9);
1554 fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC->Draw("same");
8fa56655 1555
fa516524 1556 TCanvas *c2b = new TCanvas("c2b","c2b");
1557 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetMinimum(1);
1558 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetLineColor(8);
1559 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->Draw();
1560 fHistd0rphiITSMIoneSPDInAccS500700from310MC->SetMinimum(1);
1561 fHistd0rphiITSMIoneSPDInAccS500700from310MC->SetLineColor(2);
1562 fHistd0rphiITSMIoneSPDInAccS500700from310MC->Draw("same");
1563 fHistd0rphiITSMIoneSPDInAccS500700from321MC->SetMinimum(1);
1564 fHistd0rphiITSMIoneSPDInAccS500700from321MC->SetLineColor(5);
1565 fHistd0rphiITSMIoneSPDInAccS500700from321MC->Draw("same");
1566 fHistd0rphiITSMIoneSPDInAccS500700from3122MC->SetMinimum(1);
1567 fHistd0rphiITSMIoneSPDInAccS500700from3122MC->SetLineColor(1);
1568 fHistd0rphiITSMIoneSPDInAccS500700from3122MC->Draw("same");
1569
1570 TCanvas *c2c = new TCanvas("c2c","c2c");
1571 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->SetMinimum(1);
1572 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->Draw();
1573 fHistd0rphiITSMIoneSPDInAccS500700from211MC->SetMinimum(1);
1574 fHistd0rphiITSMIoneSPDInAccS500700from211MC->SetLineColor(2);
1575 fHistd0rphiITSMIoneSPDInAccS500700from211MC->Draw("same");
1576 fHistd0rphiITSMIoneSPDInAccS500700from22MC->SetMinimum(1);
1577 fHistd0rphiITSMIoneSPDInAccS500700from22MC->SetLineColor(1);
1578 fHistd0rphiITSMIoneSPDInAccS500700from22MC->Draw("same");
1579
1580
1581
1582 TCanvas *c3 = new TCanvas("c3","c3");
1583 c3->Divide(3,2);
1584 c3_1->SetLogx();
1585 c3_2->SetLogx();
1586 c3_3->SetLogx();
1587 c3_4->SetLogx();
1588 c3_5->SetLogx();
1589 c3_6->SetLogx();
1590 Float_t d0cut[15]={0.0301,0.0401,0.0501,0.0601,0.0801,0.101,0.151,0.201,0.251,0.301,0.401,0.601,0.801,1.001,1.401};
8fa56655 1591
fa516524 1592 c3->cd(1);
1593 Float_t fracP150200[15],fracS150200[15],fracSfromStrange150200[15],fracS150200Strp30[15],fracS150200Strm30[15],fracS150200Matm10[15],fracS150200Matp10[15];
1594 Float_t intPtot150200=fHistd0rphiITSMIoneSPDInAccP150200MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP150200MC->GetNbinsX());
1595 for(Int_t i=0;i<15;i++) {
1596 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP150200MC->FindBin(-d0cut[i]);
1597 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP150200MC->FindBin(+d0cut[i]);
1598 Float_t intPcut150200=fHistd0rphiITSMIoneSPDInAccP150200MC->Integral(bin1,bin2);
1599 Float_t intScut150200=fHistd0rphiITSMIoneSPDInAccS150200MC->Integral(bin1,bin2);
1600 Float_t intSfromStrangecut150200=fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->Integral(bin1,bin2);
1601 Float_t intScut150200Strp30 = intScut150200 + 0.3*intSfromStrangecut150200;
1602 Float_t intScut150200Strm30 = intScut150200 - 0.3*intSfromStrangecut150200;
1603 Float_t intScut150200Matp10 = intScut150200 + 0.1*(intScut150200-intSfromStrangecut150200);
1604 Float_t intScut150200Matm10 = intScut150200 - 0.1*(intScut150200-intSfromStrangecut150200);
1605 fracP150200[i]=intPcut150200/intPtot150200;
1606 fracS150200[i]=1.-intScut150200/(intPcut150200+intScut150200);
1607 fracS150200Strp30[i]=1.-intScut150200Strp30/(intPcut150200+intScut150200Strp30);
1608 fracS150200Strm30[i]=1.-intScut150200Strm30/(intPcut150200+intScut150200Strm30);
1609 fracS150200Matp10[i]=1.-intScut150200Matp10/(intPcut150200+intScut150200Matp10);
1610 fracS150200Matm10[i]=1.-intScut150200Matm10/(intPcut150200+intScut150200Matm10);
1611 fracSfromStrange150200[i]=1.-intSfromStrangecut150200/(intPcut150200+intScut150200);
1612 }
1613 TGraph *gfracP150200=new TGraph(15,d0cut,fracP150200);
1614 gfracP150200->SetMarkerColor(2);
1615 gfracP150200->SetMarkerStyle(20);
1616 gfracP150200->Draw("ap");
1617 TGraph *gfracS150200=new TGraph(15,d0cut,fracS150200);
1618 gfracS150200->SetMarkerColor(4);
1619 gfracS150200->SetMarkerStyle(21);
1620 gfracS150200->Draw("p");
1621 TGraph *gfracSfromStrange150200=new TGraph(15,d0cut,fracSfromStrange150200);
1622 gfracSfromStrange150200->SetMarkerColor(8);
1623 gfracSfromStrange150200->SetMarkerStyle(22);
8fa56655 1624 gfracSfromStrange150200->Draw("p");
1625
fa516524 1626 c3->cd(2);
1627 Float_t fracP350450[15],fracS350450[15],fracSfromStrange350450[15],fracS350450Strp30[15],fracS350450Strm30[15],fracS350450Matm10[15],fracS350450Matp10[15];
1628 Float_t intPtot350450=fHistd0rphiITSMIoneSPDInAccP350450MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP350450MC->GetNbinsX());
1629 for(Int_t i=0;i<15;i++) {
1630 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP350450MC->FindBin(-d0cut[i]);
1631 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP350450MC->FindBin(+d0cut[i]);
1632 Float_t intPcut350450=fHistd0rphiITSMIoneSPDInAccP350450MC->Integral(bin1,bin2);
1633 Float_t intScut350450=fHistd0rphiITSMIoneSPDInAccS350450MC->Integral(bin1,bin2);
1634 Float_t intSfromStrangecut350450=fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->Integral(bin1,bin2);
1635 Float_t intScut350450Strp30 = intScut350450 + 0.3*intSfromStrangecut350450;
1636 Float_t intScut350450Strm30 = intScut350450 - 0.3*intSfromStrangecut350450;
1637 Float_t intScut350450Matp10 = intScut350450 + 0.1*(intScut350450-intSfromStrangecut350450);
1638 Float_t intScut350450Matm10 = intScut350450 - 0.1*(intScut350450-intSfromStrangecut350450);
1639 fracP350450[i]=intPcut350450/intPtot350450;
1640 fracS350450[i]=1.-intScut350450/(intPcut350450+intScut350450);
1641 fracS350450Strp30[i]=1.-intScut350450Strp30/(intPcut350450+intScut350450Strp30);
1642 fracS350450Strm30[i]=1.-intScut350450Strm30/(intPcut350450+intScut350450Strm30);
1643 fracS350450Matp10[i]=1.-intScut350450Matp10/(intPcut350450+intScut350450Matp10);
1644 fracS350450Matm10[i]=1.-intScut350450Matm10/(intPcut350450+intScut350450Matm10);
1645 fracSfromStrange350450[i]=1.-intSfromStrangecut350450/(intPcut350450+intScut350450);
1646 }
1647 TGraph *gfracP350450=new TGraph(15,d0cut,fracP350450);
1648 gfracP350450->SetMarkerColor(2);
1649 gfracP350450->SetMarkerStyle(20);
1650 gfracP350450->Draw("ap");
1651 TGraph *gfracS350450=new TGraph(15,d0cut,fracS350450);
1652 gfracS350450->SetMarkerColor(4);
1653 gfracS350450->SetMarkerStyle(21);
1654 gfracS350450->Draw("p");
1655 TGraph *gfracSfromStrange350450=new TGraph(15,d0cut,fracSfromStrange350450);
1656 gfracSfromStrange350450->SetMarkerColor(8);
1657 gfracSfromStrange350450->SetMarkerStyle(22);
1658 gfracSfromStrange350450->Draw("p");
8fa56655 1659
fa516524 1660 c3->cd(3);
1661 Float_t fracP500700[15],fracS500700[15],fracSfromStrange500700[15],fracS500700Strp30[15],fracS500700Strm30[15],fracS500700Matm10[15],fracS500700Matp10[15];
1662 Float_t intPtot500700=fHistd0rphiITSMIoneSPDInAccP500700MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP500700MC->GetNbinsX());
1663 for(Int_t i=0;i<15;i++) {
1664 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP500700MC->FindBin(-d0cut[i]);
1665 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP500700MC->FindBin(+d0cut[i]);
1666 Float_t intPcut500700=fHistd0rphiITSMIoneSPDInAccP500700MC->Integral(bin1,bin2);
1667 Float_t intScut500700=fHistd0rphiITSMIoneSPDInAccS500700MC->Integral(bin1,bin2);
1668 Float_t intSfromStrangecut500700=fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->Integral(bin1,bin2);
1669 Float_t intScut500700Strp30 = intScut500700 + 0.3*intSfromStrangecut500700;
1670 Float_t intScut500700Strm30 = intScut500700 - 0.3*intSfromStrangecut500700;
1671 Float_t intScut500700Matp10 = intScut500700 + 0.1*(intScut500700-intSfromStrangecut500700);
1672 Float_t intScut500700Matm10 = intScut500700 - 0.1*(intScut500700-intSfromStrangecut500700);
1673 fracP500700[i]=intPcut500700/intPtot500700;
1674 fracS500700[i]=1.-intScut500700/(intPcut500700+intScut500700);
1675 fracS500700Strp30[i]=1.-intScut500700Strp30/(intPcut500700+intScut500700Strp30);
1676 fracS500700Strm30[i]=1.-intScut500700Strm30/(intPcut500700+intScut500700Strm30);
1677 fracS500700Matp10[i]=1.-intScut500700Matp10/(intPcut500700+intScut500700Matp10);
1678 fracS500700Matm10[i]=1.-intScut500700Matm10/(intPcut500700+intScut500700Matm10);
1679 fracSfromStrange500700[i]=1.-intSfromStrangecut500700/(intPcut500700+intScut500700);
1680 }
1681 TGraph *gfracP500700=new TGraph(15,d0cut,fracP500700);
1682 gfracP500700->SetMarkerColor(2);
1683 gfracP500700->SetMarkerStyle(20);
1684 gfracP500700->Draw("ap");
1685 TGraph *gfracS500700=new TGraph(15,d0cut,fracS500700);
1686 gfracS500700->SetMarkerColor(4);
1687 gfracS500700->SetMarkerStyle(21);
1688 gfracS500700->Draw("p");
1689 TGraph *gfracSfromStrange500700=new TGraph(15,d0cut,fracSfromStrange500700);
1690 gfracSfromStrange500700->SetMarkerColor(8);
1691 gfracSfromStrange500700->SetMarkerStyle(22);
1692 gfracSfromStrange500700->Draw("p");
1693
1694 c3->cd(4);
1695 Float_t fracP10001500[15],fracS10001500[15],fracSfromStrange10001500[15],fracS10001500Strp30[15],fracS10001500Strm30[15],fracS10001500Matm10[15],fracS10001500Matp10[15];
1696 Float_t intPtot10001500=fHistd0rphiITSMIoneSPDInAccP10001500MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP10001500MC->GetNbinsX());
1697 for(Int_t i=0;i<15;i++) {
1698 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP10001500MC->FindBin(-d0cut[i]);
1699 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP10001500MC->FindBin(+d0cut[i]);
1700 Float_t intPcut10001500=fHistd0rphiITSMIoneSPDInAccP10001500MC->Integral(bin1,bin2);
1701 Float_t intScut10001500=fHistd0rphiITSMIoneSPDInAccS10001500MC->Integral(bin1,bin2);
1702 Float_t intSfromStrangecut10001500=fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->Integral(bin1,bin2);
1703 Float_t intScut10001500Strp30 = intScut10001500 + 0.3*intSfromStrangecut10001500;
1704 Float_t intScut10001500Strm30 = intScut10001500 - 0.3*intSfromStrangecut10001500;
1705 Float_t intScut10001500Matp10 = intScut10001500 + 0.1*(intScut10001500-intSfromStrangecut10001500);
1706 Float_t intScut10001500Matm10 = intScut10001500 - 0.1*(intScut10001500-intSfromStrangecut10001500);
1707 fracP10001500[i]=intPcut10001500/intPtot10001500;
1708 fracS10001500[i]=1.-intScut10001500/(intPcut10001500+intScut10001500);
1709 fracS10001500Strp30[i]=1.-intScut10001500Strp30/(intPcut10001500+intScut10001500Strp30);
1710 fracS10001500Strm30[i]=1.-intScut10001500Strm30/(intPcut10001500+intScut10001500Strm30);
1711 fracS10001500Matp10[i]=1.-intScut10001500Matp10/(intPcut10001500+intScut10001500Matp10);
1712 fracS10001500Matm10[i]=1.-intScut10001500Matm10/(intPcut10001500+intScut10001500Matm10);
1713 fracSfromStrange10001500[i]=1.-intSfromStrangecut10001500/(intPcut10001500+intScut10001500);
1714 }
1715 TGraph *gfracP10001500=new TGraph(15,d0cut,fracP10001500);
1716 gfracP10001500->SetMarkerColor(2);
1717 gfracP10001500->SetMarkerStyle(20);
1718 gfracP10001500->Draw("ap");
1719 TGraph *gfracS10001500=new TGraph(15,d0cut,fracS10001500);
1720 gfracS10001500->SetMarkerColor(4);
1721 gfracS10001500->SetMarkerStyle(21);
1722 gfracS10001500->Draw("p");
1723 TGraph *gfracSfromStrange10001500=new TGraph(15,d0cut,fracSfromStrange10001500);
1724 gfracSfromStrange10001500->SetMarkerColor(8);
1725 gfracSfromStrange10001500->SetMarkerStyle(22);
1726 gfracSfromStrange10001500->Draw("p");
8fa56655 1727
fa516524 1728 c3->cd(5);
1729 Float_t fracP25004000[15],fracS25004000[15],fracSfromStrange25004000[15],fracS25004000Strp30[15],fracS25004000Strm30[15],fracS25004000Matm10[15],fracS25004000Matp10[15];
1730 Float_t intPtot25004000=fHistd0rphiITSMIoneSPDInAccP25004000MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP25004000MC->GetNbinsX());
1731 for(Int_t i=0;i<15;i++) {
1732 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP25004000MC->FindBin(-d0cut[i]);
1733 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP25004000MC->FindBin(+d0cut[i]);
1734 Float_t intPcut25004000=fHistd0rphiITSMIoneSPDInAccP25004000MC->Integral(bin1,bin2);
1735 Float_t intScut25004000=fHistd0rphiITSMIoneSPDInAccS25004000MC->Integral(bin1,bin2);
1736 Float_t intSfromStrangecut25004000=fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->Integral(bin1,bin2);
1737 Float_t intScut25004000Strp30 = intScut25004000 + 0.3*intSfromStrangecut25004000;
1738 Float_t intScut25004000Strm30 = intScut25004000 - 0.3*intSfromStrangecut25004000;
1739 Float_t intScut25004000Matp10 = intScut25004000 + 0.1*(intScut25004000-intSfromStrangecut25004000);
1740 Float_t intScut25004000Matm10 = intScut25004000 - 0.1*(intScut25004000-intSfromStrangecut25004000);
1741 fracP25004000[i]=intPcut25004000/intPtot25004000;
1742 fracS25004000[i]=1.-intScut25004000/(intPcut25004000+intScut25004000);
1743 fracS25004000Strp30[i]=1.-intScut25004000Strp30/(intPcut25004000+intScut25004000Strp30);
1744 fracS25004000Strm30[i]=1.-intScut25004000Strm30/(intPcut25004000+intScut25004000Strm30);
1745 fracS25004000Matp10[i]=1.-intScut25004000Matp10/(intPcut25004000+intScut25004000Matp10);
1746 fracS25004000Matm10[i]=1.-intScut25004000Matm10/(intPcut25004000+intScut25004000Matm10);
1747 fracSfromStrange25004000[i]=1.-intSfromStrangecut25004000/(intPcut25004000+intScut25004000);
1748 }
1749 TGraph *gfracP25004000=new TGraph(15,d0cut,fracP25004000);
1750 gfracP25004000->SetMarkerColor(2);
1751 gfracP25004000->SetMarkerStyle(20);
1752 gfracP25004000->Draw("ap");
1753 TGraph *gfracS25004000=new TGraph(15,d0cut,fracS25004000);
1754 gfracS25004000->SetMarkerColor(4);
1755 gfracS25004000->SetMarkerStyle(21);
1756 gfracS25004000->Draw("p");
1757 TGraph *gfracSfromStrange25004000=new TGraph(15,d0cut,fracSfromStrange25004000);
1758 gfracSfromStrange25004000->SetMarkerColor(8);
1759 gfracSfromStrange25004000->SetMarkerStyle(22);
1760 gfracSfromStrange25004000->Draw("p");
1761
1762 c3->cd(6);
1763 Float_t fracP40008000[15],fracS40008000[15],fracSfromStrange40008000[15],fracS40008000Strp30[15],fracS40008000Strm30[15],fracS40008000Matm10[15],fracS40008000Matp10[15];
1764 Float_t intPtot40008000=fHistd0rphiITSMIoneSPDInAccP40008000MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP40008000MC->GetNbinsX());
1765 for(Int_t i=0;i<15;i++) {
1766 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP40008000MC->FindBin(-d0cut[i]);
1767 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP40008000MC->FindBin(+d0cut[i]);
1768 Float_t intPcut40008000=fHistd0rphiITSMIoneSPDInAccP40008000MC->Integral(bin1,bin2);
1769 Float_t intScut40008000=fHistd0rphiITSMIoneSPDInAccS40008000MC->Integral(bin1,bin2);
1770 Float_t intSfromStrangecut40008000=fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->Integral(bin1,bin2);
1771 Float_t intScut40008000Strp30 = intScut40008000 + 0.3*intSfromStrangecut40008000;
1772 Float_t intScut40008000Strm30 = intScut40008000 - 0.3*intSfromStrangecut40008000;
1773 Float_t intScut40008000Matp10 = intScut40008000 + 0.1*(intScut40008000-intSfromStrangecut40008000);
1774 Float_t intScut40008000Matm10 = intScut40008000 - 0.1*(intScut40008000-intSfromStrangecut40008000);
1775 fracP40008000[i]=intPcut40008000/intPtot40008000;
1776 fracS40008000[i]=1.-intScut40008000/(intPcut40008000+intScut40008000);
1777 fracS40008000Strp30[i]=1.-intScut40008000Strp30/(intPcut40008000+intScut40008000Strp30);
1778 fracS40008000Strm30[i]=1.-intScut40008000Strm30/(intPcut40008000+intScut40008000Strm30);
1779 fracS40008000Matp10[i]=1.-intScut40008000Matp10/(intPcut40008000+intScut40008000Matp10);
1780 fracS40008000Matm10[i]=1.-intScut40008000Matm10/(intPcut40008000+intScut40008000Matm10);
1781 fracSfromStrange40008000[i]=1.-intSfromStrangecut40008000/(intPcut40008000+intScut40008000);
1782 }
1783 TGraph *gfracP40008000=new TGraph(15,d0cut,fracP40008000);
1784 gfracP40008000->SetMarkerColor(2);
1785 gfracP40008000->SetMarkerStyle(20);
1786 gfracP40008000->Draw("ap");
1787 TGraph *gfracS40008000=new TGraph(15,d0cut,fracS40008000);
1788 gfracS40008000->SetMarkerColor(4);
1789 gfracS40008000->SetMarkerStyle(21);
1790 gfracS40008000->Draw("p");
1791 TGraph *gfracSfromStrange40008000=new TGraph(15,d0cut,fracSfromStrange40008000);
1792 gfracSfromStrange40008000->SetMarkerColor(8);
1793 gfracSfromStrange40008000->SetMarkerStyle(22);
1794 gfracSfromStrange40008000->Draw("p");
fa516524 1795
8fa56655 1796
fa516524 1797 TCanvas *c4 = new TCanvas("c4","c4");
1798 c4->Divide(3,1);
1799 c4_1->SetLogx();
1800 c4_2->SetLogx();
1801 c4_3->SetLogx();
1802
1803 c4->cd(1);
1804 Float_t intDatacut150200[15],intDataPcut150200[15],intDataPall150200[15];
1805 for(Int_t i=0;i<15;i++) {
1806 Int_t bin1=fHistd0rphiITSMIoneSPDInAcc150200->FindBin(-d0cut[i]);
1807 Int_t bin2=fHistd0rphiITSMIoneSPDInAcc150200->FindBin(+d0cut[i]);
1808 intDatacut150200[i]=fHistd0rphiITSMIoneSPDInAcc150200->Integral(bin1,bin2);
1809 intDataPcut150200[i]=intDatacut150200[i]*fracS150200[i];
1810 intDataPall150200[i]=intDataPcut150200[i]/fracP150200[i];
1811 }
1812 TGraph *gintDatacut150200=new TGraph(15,d0cut,intDatacut150200);
1813 gintDatacut150200->SetMarkerColor(1);
1814 gintDatacut150200->SetMarkerStyle(20);
1815 gintDatacut150200->Draw("ap");
1816 TGraph *gintDataPcut150200=new TGraph(15,d0cut,intDataPcut150200);
1817 gintDataPcut150200->SetMarkerColor(1);
1818 gintDataPcut150200->SetMarkerStyle(24);
1819 gintDataPcut150200->Draw("p");
1820 TGraph *gintDataPall150200=new TGraph(15,d0cut,intDataPall150200);
1821 gintDataPall150200->SetMarkerColor(2);
1822 gintDataPall150200->SetMarkerStyle(22);
1823 gintDataPall150200->Draw("p");
1824
1825 c4->cd(2);
1826 Float_t intDatacut500700[15],intDataPcut500700[15],intDataPall500700[15];
1827 for(Int_t i=0;i<15;i++) {
1828 Int_t bin1=fHistd0rphiITSMIoneSPDInAcc500700->FindBin(-d0cut[i]);
1829 Int_t bin2=fHistd0rphiITSMIoneSPDInAcc500700->FindBin(+d0cut[i]);
1830 intDatacut500700[i]=fHistd0rphiITSMIoneSPDInAcc500700->Integral(bin1,bin2);
1831 intDataPcut500700[i]=intDatacut500700[i]*fracS500700[i];
1832 intDataPall500700[i]=intDataPcut500700[i]/fracP500700[i];
1833 }
1834 TGraph *gintDatacut500700=new TGraph(15,d0cut,intDatacut500700);
1835 gintDatacut500700->SetMarkerColor(1);
1836 gintDatacut500700->SetMarkerStyle(20);
1837 gintDatacut500700->Draw("ap");
1838 TGraph *gintDataPcut500700=new TGraph(15,d0cut,intDataPcut500700);
1839 gintDataPcut500700->SetMarkerColor(1);
1840 gintDataPcut500700->SetMarkerStyle(24);
1841 gintDataPcut500700->Draw("p");
1842 TGraph *gintDataPall500700=new TGraph(15,d0cut,intDataPall500700);
1843 gintDataPall500700->SetMarkerColor(2);
1844 gintDataPall500700->SetMarkerStyle(22);
1845 gintDataPall500700->Draw("p");
1846
1847 c4->cd(3);
1848 Float_t intDatacut10001500[15],intDataPcut10001500[15],intDataPall10001500[15];
1849 for(Int_t i=0;i<15;i++) {
1850 Int_t bin1=fHistd0rphiITSMIoneSPDInAcc10001500->FindBin(-d0cut[i]);
1851 Int_t bin2=fHistd0rphiITSMIoneSPDInAcc10001500->FindBin(+d0cut[i]);
1852 intDatacut10001500[i]=fHistd0rphiITSMIoneSPDInAcc10001500->Integral(bin1,bin2);
1853 intDataPcut10001500[i]=intDatacut10001500[i]*fracS10001500[i];
1854 intDataPall10001500[i]=intDataPcut10001500[i]/fracP10001500[i];
1855 }
1856 TGraph *gintDatacut10001500=new TGraph(15,d0cut,intDatacut10001500);
1857 gintDatacut10001500->SetMarkerColor(1);
1858 gintDatacut10001500->SetMarkerStyle(20);
1859 gintDatacut10001500->Draw("ap");
1860 TGraph *gintDataPcut10001500=new TGraph(15,d0cut,intDataPcut10001500);
1861 gintDataPcut10001500->SetMarkerColor(1);
1862 gintDataPcut10001500->SetMarkerStyle(24);
1863 gintDataPcut10001500->Draw("p");
1864 TGraph *gintDataPall10001500=new TGraph(15,d0cut,intDataPall10001500);
1865 gintDataPall10001500->SetMarkerColor(2);
1866 gintDataPall10001500->SetMarkerStyle(22);
1867 gintDataPall10001500->Draw("p");
1868
1869
1870 TCanvas *c5 = new TCanvas("c5","c5");
1871 c5->Divide(3,1);
1872 c5_1->SetLogx();
1873 c5_2->SetLogx();
1874 c5_3->SetLogx();
1875
1876 c5->cd(1);
1877 Float_t intDataPall150200Strp30[15],intDataPall150200Strm30[15],intDataPall150200Matp10[15],intDataPall150200Matm10[15];
1878 for(Int_t i=0;i<15;i++) {
1879 intDataPall150200Strp30[i]=intDatacut150200[i]*fracS150200Strp30[i]/fracP150200[i];
1880 intDataPall150200Strm30[i]=intDatacut150200[i]*fracS150200Strm30[i]/fracP150200[i];
1881 intDataPall150200Matp10[i]=intDatacut150200[i]*fracS150200Matp10[i]/fracP150200[i];
1882 intDataPall150200Matm10[i]=intDatacut150200[i]*fracS150200Matm10[i]/fracP150200[i];
1883 }
1884 gintDataPall150200->Draw("ap");
1885
1886 TGraph *gintDataPall150200Strp30=new TGraph(15,d0cut,intDataPall150200Strp30);
1887 gintDataPall150200Strp30->SetMarkerColor(1);
1888 gintDataPall150200Strp30->SetMarkerStyle(22);
1889 gintDataPall150200Strp30->Draw("p");
1890 TGraph *gintDataPall150200Strm30=new TGraph(15,d0cut,intDataPall150200Strm30);
1891 gintDataPall150200Strm30->SetMarkerColor(1);
1892 gintDataPall150200Strm30->SetMarkerStyle(22);
1893 gintDataPall150200Strm30->Draw("p");
1894 TGraph *gintDataPall150200Matp10=new TGraph(15,d0cut,intDataPall150200Matp10);
1895 gintDataPall150200Matp10->SetMarkerColor(3);
1896 gintDataPall150200Matp10->SetMarkerStyle(22);
1897 gintDataPall150200Matp10->Draw("p");
1898 TGraph *gintDataPall150200Matm10=new TGraph(15,d0cut,intDataPall150200Matm10);
1899 gintDataPall150200Matm10->SetMarkerColor(3);
1900 gintDataPall150200Matm10->SetMarkerStyle(22);
1901 gintDataPall150200Matm10->Draw("p");
1902
1903 c5->cd(2);
1904 Float_t intDataPall500700Strp30[15],intDataPall500700Strm30[15],intDataPall500700Matp10[15],intDataPall500700Matm10[15];
1905 for(Int_t i=0;i<15;i++) {
1906 intDataPall500700Strp30[i]=intDatacut500700[i]*fracS500700Strp30[i]/fracP500700[i];
1907 intDataPall500700Strm30[i]=intDatacut500700[i]*fracS500700Strm30[i]/fracP500700[i];
1908 intDataPall500700Matp10[i]=intDatacut500700[i]*fracS500700Matp10[i]/fracP500700[i];
1909 intDataPall500700Matm10[i]=intDatacut500700[i]*fracS500700Matm10[i]/fracP500700[i];
1910 }
1911 gintDataPall500700->Draw("ap");
1912
1913 TGraph *gintDataPall500700Strp30=new TGraph(15,d0cut,intDataPall500700Strp30);
1914 gintDataPall500700Strp30->SetMarkerColor(1);
1915 gintDataPall500700Strp30->SetMarkerStyle(22);
1916 gintDataPall500700Strp30->Draw("p");
1917 TGraph *gintDataPall500700Strm30=new TGraph(15,d0cut,intDataPall500700Strm30);
1918 gintDataPall500700Strm30->SetMarkerColor(1);
1919 gintDataPall500700Strm30->SetMarkerStyle(22);
1920 gintDataPall500700Strm30->Draw("p");
1921 TGraph *gintDataPall500700Matp10=new TGraph(15,d0cut,intDataPall500700Matp10);
1922 gintDataPall500700Matp10->SetMarkerColor(3);
1923 gintDataPall500700Matp10->SetMarkerStyle(22);
1924 gintDataPall500700Matp10->Draw("p");
1925 TGraph *gintDataPall500700Matm10=new TGraph(15,d0cut,intDataPall500700Matm10);
1926 gintDataPall500700Matm10->SetMarkerColor(3);
1927 gintDataPall500700Matm10->SetMarkerStyle(22);
1928 gintDataPall500700Matm10->Draw("p");
1929
1930 c5->cd(3);
1931 Float_t intDataPall10001500Strp30[15],intDataPall10001500Strm30[15],intDataPall10001500Matp10[15],intDataPall10001500Matm10[15];
1932 for(Int_t i=0;i<15;i++) {
1933 intDataPall10001500Strp30[i]=intDatacut10001500[i]*fracS10001500Strp30[i]/fracP10001500[i];
1934 intDataPall10001500Strm30[i]=intDatacut10001500[i]*fracS10001500Strm30[i]/fracP10001500[i];
1935 intDataPall10001500Matp10[i]=intDatacut10001500[i]*fracS10001500Matp10[i]/fracP10001500[i];
1936 intDataPall10001500Matm10[i]=intDatacut10001500[i]*fracS10001500Matm10[i]/fracP10001500[i];
1937 }
1938 gintDataPall10001500->Draw("ap");
1939
1940 TGraph *gintDataPall10001500Strp30=new TGraph(15,d0cut,intDataPall10001500Strp30);
1941 gintDataPall10001500Strp30->SetMarkerColor(1);
1942 gintDataPall10001500Strp30->SetMarkerStyle(22);
1943 gintDataPall10001500Strp30->Draw("p");
1944 TGraph *gintDataPall10001500Strm30=new TGraph(15,d0cut,intDataPall10001500Strm30);
1945 gintDataPall10001500Strm30->SetMarkerColor(1);
1946 gintDataPall10001500Strm30->SetMarkerStyle(22);
1947 gintDataPall10001500Strm30->Draw("p");
1948 TGraph *gintDataPall10001500Matp10=new TGraph(15,d0cut,intDataPall10001500Matp10);
1949 gintDataPall10001500Matp10->SetMarkerColor(3);
1950 gintDataPall10001500Matp10->SetMarkerStyle(22);
1951 gintDataPall10001500Matp10->Draw("p");
1952 TGraph *gintDataPall10001500Matm10=new TGraph(15,d0cut,intDataPall10001500Matm10);
1953 gintDataPall10001500Matm10->SetMarkerColor(3);
1954 gintDataPall10001500Matm10->SetMarkerStyle(22);
1955 gintDataPall10001500Matm10->Draw("p");
8fa56655 1956
fa516524 1957 return;
1958}
1959//---------------------------------------------------------------------------
1960void PlotImpPar_z() {
1961
1962 TFile *fMC= new TFile("ITS.Performance_lhc10a8_100k.root");
1963
1964 TList *list=(TList*)fMC->Get("cOutputITS");
1965 TH1F *fHistd0zITSMIoneSPDInAccP150200MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP150200");
1966 TH1F *fHistd0zITSMIoneSPDInAccS150200MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS150200");
1967 TH1F *fHistd0zITSMIoneSPDInAccP500700MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP500700");
1968 TH1F *fHistd0zITSMIoneSPDInAccS500700MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS500700");
1969 TH1F *fHistd0zITSMIoneSPDInAccP10001500MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP10001500");
1970 TH1F *fHistd0zITSMIoneSPDInAccS10001500MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS10001500");
1971 TH1F *fHistd0zITSMIoneSPDInAcc150200MC=(TH1F*)fHistd0zITSMIoneSPDInAccP150200MC->Clone("fHistd0zITSMIoneSPDInAcc150200MC");
1972 fHistd0zITSMIoneSPDInAcc150200MC->Add(fHistd0zITSMIoneSPDInAccS150200MC);
1973 fHistd0zITSMIoneSPDInAcc150200MC->Scale(1./fHistd0zITSMIoneSPDInAcc150200MC->GetEntries());
1974 TH1F *fHistd0zITSMIoneSPDInAcc500700MC=(TH1F*)fHistd0zITSMIoneSPDInAccP500700MC->Clone("fHistd0zITSMIoneSPDInAcc500700MC");
1975 fHistd0zITSMIoneSPDInAcc500700MC->Add(fHistd0zITSMIoneSPDInAccS500700MC);
1976 fHistd0zITSMIoneSPDInAcc500700MC->Scale(1./fHistd0zITSMIoneSPDInAcc500700MC->GetEntries());
1977 TH1F *fHistd0zITSMIoneSPDInAcc10001500MC=(TH1F*)fHistd0zITSMIoneSPDInAccP10001500MC->Clone("fHistd0zITSMIoneSPDInAcc10001500MC");
1978 fHistd0zITSMIoneSPDInAcc10001500MC->Add(fHistd0zITSMIoneSPDInAccS10001500MC);
1979 fHistd0zITSMIoneSPDInAcc10001500MC->Scale(1./fHistd0zITSMIoneSPDInAcc10001500MC->GetEntries());
1980
1981
1982 TFile *f= new TFile("ITS.Performance_104892.root");
1983
1984 TList *list=(TList*)f->Get("cOutputITS");
1985 TH1F *fHistd0zITSMIoneSPDInAccP150200 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP150200");
1986 TH1F *fHistd0zITSMIoneSPDInAccS150200 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS150200");
1987 TH1F *fHistd0zITSMIoneSPDInAccP500700 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP500700");
1988 TH1F *fHistd0zITSMIoneSPDInAccS500700 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS500700");
1989 TH1F *fHistd0zITSMIoneSPDInAccP10001500 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP10001500");
1990 TH1F *fHistd0zITSMIoneSPDInAccS10001500 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS10001500");
1991 TH1F *fHistd0zITSMIoneSPDInAcc150200=(TH1F*)fHistd0zITSMIoneSPDInAccP150200->Clone("fHistd0zITSMIoneSPDInAcc150200");
1992 fHistd0zITSMIoneSPDInAcc150200->Add(fHistd0zITSMIoneSPDInAccS150200);
1993 fHistd0zITSMIoneSPDInAcc150200->Scale(1./fHistd0zITSMIoneSPDInAcc150200->GetEntries());
1994 TH1F *fHistd0zITSMIoneSPDInAcc500700=(TH1F*)fHistd0zITSMIoneSPDInAccP500700->Clone("fHistd0zITSMIoneSPDInAcc500700");
1995 fHistd0zITSMIoneSPDInAcc500700->Add(fHistd0zITSMIoneSPDInAccS500700);
1996 fHistd0zITSMIoneSPDInAcc500700->Scale(1./fHistd0zITSMIoneSPDInAcc500700->GetEntries());
1997 TH1F *fHistd0zITSMIoneSPDInAcc10001500=(TH1F*)fHistd0zITSMIoneSPDInAccP10001500->Clone("fHistd0zITSMIoneSPDInAcc10001500");
1998 fHistd0zITSMIoneSPDInAcc10001500->Add(fHistd0zITSMIoneSPDInAccS10001500);
1999 fHistd0zITSMIoneSPDInAcc10001500->Scale(1./fHistd0zITSMIoneSPDInAcc10001500->GetEntries());
2000
2001
2002 TCanvas *c1 = new TCanvas("c1","c1");
2003 c1->Divide(3,1);
2004 c1->cd(1);
2005 fHistd0zITSMIoneSPDInAcc150200MC->SetLineColor(2);
2006 fHistd0zITSMIoneSPDInAcc150200MC->Draw();
2007 fHistd0zITSMIoneSPDInAcc150200->SetLineColor(4);
2008 fHistd0zITSMIoneSPDInAcc150200->Draw("same");
2009 c1->cd(2);
2010 fHistd0zITSMIoneSPDInAcc500700MC->SetLineColor(2);
2011 fHistd0zITSMIoneSPDInAcc500700MC->Draw();
2012 fHistd0zITSMIoneSPDInAcc500700->SetLineColor(4);
2013 fHistd0zITSMIoneSPDInAcc500700->Draw("same");
2014 c1->cd(3);
2015 fHistd0zITSMIoneSPDInAcc10001500MC->SetLineColor(2);
2016 fHistd0zITSMIoneSPDInAcc10001500MC->Draw();
2017 fHistd0zITSMIoneSPDInAcc10001500->SetLineColor(4);
2018 fHistd0zITSMIoneSPDInAcc10001500->Draw("same");
2019
2020 TCanvas *c2 = new TCanvas("c2","c2");
2021 c2->Divide(3,1);
2022 c2_1->SetLogy();
2023 c2_2->SetLogy();
2024 c2_3->SetLogy();
2025 c2->cd(1);
2026 fHistd0zITSMIoneSPDInAccP150200MC->SetMinimum(1);
2027 fHistd0zITSMIoneSPDInAccS150200MC->SetMinimum(1);
2028 fHistd0zITSMIoneSPDInAccP150200MC->SetLineColor(1);
2029 fHistd0zITSMIoneSPDInAccP150200MC->Draw();
2030 fHistd0zITSMIoneSPDInAccS150200MC->SetLineColor(6);
2031 fHistd0zITSMIoneSPDInAccS150200MC->Draw("same");
2032 c2->cd(2);
2033 fHistd0zITSMIoneSPDInAccP500700MC->SetMinimum(1);
2034 fHistd0zITSMIoneSPDInAccS500700MC->SetMinimum(1);
2035 fHistd0zITSMIoneSPDInAccP500700MC->SetLineColor(1);
2036 fHistd0zITSMIoneSPDInAccP500700MC->Draw();
2037 fHistd0zITSMIoneSPDInAccS500700MC->SetLineColor(6);
2038 fHistd0zITSMIoneSPDInAccS500700MC->Draw("same");
2039 c2->cd(3);
2040 fHistd0zITSMIoneSPDInAccP10001500MC->SetMinimum(1);
2041 fHistd0zITSMIoneSPDInAccS10001500MC->SetMinimum(1);
2042 fHistd0zITSMIoneSPDInAccP10001500MC->SetLineColor(1);
2043 fHistd0zITSMIoneSPDInAccP10001500MC->Draw();
2044 fHistd0zITSMIoneSPDInAccS10001500MC->SetLineColor(6);
2045 fHistd0zITSMIoneSPDInAccS10001500MC->Draw("same");
2046
2047
2048 TCanvas *c3 = new TCanvas("c3","c3");
2049 c3->Divide(3,1);
2050 c3_1->SetLogx();
2051 c3_2->SetLogx();
2052 c3_3->SetLogx();
2053 Float_t d0cut[15]={0.05,0.06,0.08,0.10,0.15,0.2,0.25,0.30,0.35,0.40,0.50,0.60,0.80,1.,3.40};
2054
2055 c3->cd(1);
2056 Float_t fracP150200[15],fracS150200[15];
2057 Float_t intPtot150200=fHistd0zITSMIoneSPDInAccP150200MC->Integral(1,fHistd0zITSMIoneSPDInAccP150200MC->GetNbinsX());
2058 for(Int_t i=0;i<15;i++) {
2059 Int_t bin1=fHistd0zITSMIoneSPDInAccP150200MC->FindBin(-d0cut[i]);
2060 Int_t bin2=fHistd0zITSMIoneSPDInAccP150200MC->FindBin(+d0cut[i]);
2061 Float_t intPcut150200=fHistd0zITSMIoneSPDInAccP150200MC->Integral(bin1,bin2);
2062 Float_t intScut150200=fHistd0zITSMIoneSPDInAccS150200MC->Integral(bin1,bin2);
2063 fracP150200[i]=intPcut150200/intPtot150200;
2064 fracS150200[i]=1.-intScut150200/(intPcut150200+intScut150200);
2065 }
2066 TGraph *gfracP150200=new TGraph(15,d0cut,fracP150200);
2067 gfracP150200->SetMarkerColor(2);
2068 gfracP150200->SetMarkerStyle(20);
2069 gfracP150200->Draw("ap");
2070 TGraph *gfracS150200=new TGraph(15,d0cut,fracS150200);
2071 gfracS150200->SetMarkerColor(4);
2072 gfracS150200->SetMarkerStyle(21);
2073 gfracS150200->Draw("p");
2074
2075 c3->cd(2);
2076 Float_t fracP500700[15],fracS500700[15];
2077 Float_t intPtot500700=fHistd0zITSMIoneSPDInAccP500700MC->Integral(1,fHistd0zITSMIoneSPDInAccP500700MC->GetNbinsX());
2078 for(Int_t i=0;i<15;i++) {
2079 Int_t bin1=fHistd0zITSMIoneSPDInAccP500700MC->FindBin(-d0cut[i]);
2080 Int_t bin2=fHistd0zITSMIoneSPDInAccP500700MC->FindBin(+d0cut[i]);
2081 Float_t intPcut500700=fHistd0zITSMIoneSPDInAccP500700MC->Integral(bin1,bin2);
2082 Float_t intScut500700=fHistd0zITSMIoneSPDInAccS500700MC->Integral(bin1,bin2);
2083 fracP500700[i]=intPcut500700/intPtot500700;
2084 fracS500700[i]=1.-intScut500700/(intPcut500700+intScut500700);
2085 }
2086 TGraph *gfracP500700=new TGraph(15,d0cut,fracP500700);
2087 gfracP500700->SetMarkerColor(2);
2088 gfracP500700->SetMarkerStyle(20);
2089 gfracP500700->Draw("ap");
2090 TGraph *gfracS500700=new TGraph(15,d0cut,fracS500700);
2091 gfracS500700->SetMarkerColor(4);
2092 gfracS500700->SetMarkerStyle(21);
2093 gfracS500700->Draw("p");
2094
2095 c3->cd(3);
2096 Float_t fracP10001500[15],fracS10001500[15];
2097 Float_t intPtot10001500=fHistd0zITSMIoneSPDInAccP10001500MC->Integral(1,fHistd0zITSMIoneSPDInAccP10001500MC->GetNbinsX());
2098 for(Int_t i=0;i<15;i++) {
2099 Int_t bin1=fHistd0zITSMIoneSPDInAccP10001500MC->FindBin(-d0cut[i]);
2100 Int_t bin2=fHistd0zITSMIoneSPDInAccP10001500MC->FindBin(+d0cut[i]);
2101 Float_t intPcut10001500=fHistd0zITSMIoneSPDInAccP10001500MC->Integral(bin1,bin2);
2102 Float_t intScut10001500=fHistd0zITSMIoneSPDInAccS10001500MC->Integral(bin1,bin2);
2103 fracP10001500[i]=intPcut10001500/intPtot10001500;
2104 fracS10001500[i]=1.-intScut10001500/(intPcut10001500+intScut10001500);
2105 }
2106 TGraph *gfracP10001500=new TGraph(15,d0cut,fracP10001500);
2107 gfracP10001500->SetMarkerColor(2);
2108 gfracP10001500->SetMarkerStyle(20);
2109 gfracP10001500->Draw("ap");
2110 TGraph *gfracS10001500=new TGraph(15,d0cut,fracS10001500);
2111 gfracS10001500->SetMarkerColor(4);
2112 gfracS10001500->SetMarkerStyle(21);
2113 gfracS10001500->Draw("p");
2114
2115 return;
2116}
2117//----------------------------------------------------------------------------
2118void Corrections() {
2119
2120 TFile *f_cutits= new TFile("ITS.Performance_lhc09d10_500k_noz0cut_nodcaTPC.root");
2121
2122 TList *list=(TList*)f_cutits->Get("cOutputITS");
2123
2124 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
2125 Float_t nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
2126
2127 TH1F *fHistPtITSTPCsel = (TH1F*)list->FindObject("fHistPtITSTPCsel");
2128 NormalizePtHist(fHistPtITSTPCsel,nEventsSelWithVertex);
2129 TH1F *fHistPtITSTPCselP = (TH1F*)list->FindObject("fHistPtITSTPCselP");
2130 NormalizePtHist(fHistPtITSTPCselP,nEventsSelWithVertex);
2131 TH1F *fHistPtITSTPCselS = (TH1F*)list->FindObject("fHistPtITSTPCselS");
2132 NormalizePtHist(fHistPtITSTPCselS,nEventsSelWithVertex);
2133 TH1F *fHistPtITSTPCselPfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselPfromStrange");
2134 NormalizePtHist(fHistPtITSTPCselPfromStrange,nEventsSelWithVertex);
2135 TH1F *fHistPtITSTPCselSfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
2136 NormalizePtHist(fHistPtITSTPCselSfromStrange,nEventsSelWithVertex);
2137 TH1F *fHistPtITSTPCselSfromMat = (TH1F*)list->FindObject("fHistPtITSTPCselSfromMat");
2138 NormalizePtHist(fHistPtITSTPCselSfromMat,nEventsSelWithVertex);
2139
2140
2141 TFile *f_cuttpc= new TFile("ITS.Performance_lhc09d10_500k_noz0cut.root");
2142
2143 TList *list=(TList*)f_cuttpc->Get("cOutputITS");
2144
2145 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
2146 nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
2147
2148 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
2149 NormalizePtHist(fHistPtTPCInAcc,nEventsSelWithVertex);
2150 TH1F *fHistPtTPCInAccP = (TH1F*)list->FindObject("fHistPtTPCInAccP");
2151 NormalizePtHist(fHistPtTPCInAccP,nEventsSelWithVertex);
2152 TH1F *fHistPtTPCInAccS = (TH1F*)list->FindObject("fHistPtTPCInAccS");
2153 NormalizePtHist(fHistPtTPCInAccS,nEventsSelWithVertex);
2154 TH1F *fHistPtTPCInAccPfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccPfromStrange");
2155 NormalizePtHist(fHistPtTPCInAccPfromStrange,nEventsSelWithVertex);
2156 TH1F *fHistPtTPCInAccSfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccSfromStrange");
2157 NormalizePtHist(fHistPtTPCInAccSfromStrange,nEventsSelWithVertex);
2158 TH1F *fHistPtTPCInAccSfromMat = (TH1F*)list->FindObject("fHistPtTPCInAccSfromMat");
2159 NormalizePtHist(fHistPtTPCInAccSfromMat,nEventsSelWithVertex);
2160
2161
2162 TFile *f_nocut= new TFile("ITS.Performance_lhc09d10_500k_nod0z0cut.root");
2163
2164 TList *list=(TList*)f_nocut->Get("cOutputITS");
2165
2166 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
2167 nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
2168
2169 TH1F *fHistPtITSTPCselnocut = (TH1F*)list->FindObject("fHistPtITSTPCsel");
2170 NormalizePtHist(fHistPtITSTPCselnocut,nEventsSelWithVertex);
2171 TH1F *fHistPtITSTPCselnocutP = (TH1F*)list->FindObject("fHistPtITSTPCselP");
2172 NormalizePtHist(fHistPtITSTPCselnocutP,nEventsSelWithVertex);
2173 TH1F *fHistPtITSTPCselnocutS = (TH1F*)list->FindObject("fHistPtITSTPCselS");
2174 NormalizePtHist(fHistPtITSTPCselnocutS,nEventsSelWithVertex);
2175 TH1F *fHistPtITSTPCselnocutPfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselPfromStrange");
2176 NormalizePtHist(fHistPtITSTPCselnocutPfromStrange,nEventsSelWithVertex);
2177 TH1F *fHistPtITSTPCselnocutSfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
2178 NormalizePtHist(fHistPtITSTPCselnocutSfromStrange,nEventsSelWithVertex);
2179 TH1F *fHistPtITSTPCselnocutSfromMat = (TH1F*)list->FindObject("fHistPtITSTPCselSfromMat");
2180 NormalizePtHist(fHistPtITSTPCselnocutSfromMat,nEventsSelWithVertex);
2181
2182 TH1F *fHistPtTPCInAccnocut = (TH1F*)list->FindObject("fHistPtTPCInAcc");
2183 NormalizePtHist(fHistPtTPCInAccnocut,nEventsSelWithVertex);
2184 TH1F *fHistPtTPCInAccnocutP = (TH1F*)list->FindObject("fHistPtTPCInAccP");
2185 NormalizePtHist(fHistPtTPCInAccnocutP,nEventsSelWithVertex);
2186 TH1F *fHistPtTPCInAccnocutS = (TH1F*)list->FindObject("fHistPtTPCInAccS");
2187 NormalizePtHist(fHistPtTPCInAccnocutS,nEventsSelWithVertex);
2188 TH1F *fHistPtTPCInAccnocutPfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccPfromStrange");
2189 NormalizePtHist(fHistPtTPCInAccnocutPfromStrange,nEventsSelWithVertex);
2190 TH1F *fHistPtTPCInAccnocutSfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccSfromStrange");
2191 NormalizePtHist(fHistPtTPCInAccnocutSfromStrange,nEventsSelWithVertex);
2192 TH1F *fHistPtTPCInAccnocutSfromMat = (TH1F*)list->FindObject("fHistPtTPCInAccSfromMat");
2193 NormalizePtHist(fHistPtTPCInAccnocutSfromMat,nEventsSelWithVertex);
2194
2195 // cut efficiencies
2196
2197 TH1F *fHistPtITSTPCeffselP = (TH1F*)fHistPtITSTPCselP->Clone("fHistPtITSTPCeffselP");
2198 fHistPtITSTPCeffselP->Divide(fHistPtITSTPCselnocutP);
2199 TH1F *fHistPtITSTPCeffselS = (TH1F*)fHistPtITSTPCselS->Clone("fHistPtITSTPCeffselS");
2200 fHistPtITSTPCeffselS->Divide(fHistPtITSTPCselnocutS);
2201 TH1F *fHistPtITSTPCeffselSfromStrange = (TH1F*)fHistPtITSTPCselSfromStrange->Clone("fHistPtITSTPCeffselSfromStrange");
2202 fHistPtITSTPCeffselSfromStrange->Divide(fHistPtITSTPCselnocutSfromStrange);
2203 TH1F *fHistPtITSTPCeffselSfromMat = (TH1F*)fHistPtITSTPCselSfromMat->Clone("fHistPtITSTPCeffselSfromMat");
2204 fHistPtITSTPCeffselSfromMat->Divide(fHistPtITSTPCselnocutSfromMat);
2205
2206 TH1F *fHistPtTPCeffInAccP = (TH1F*)fHistPtTPCInAccP->Clone("fHistPtTPCeffInAccP");
2207 fHistPtTPCeffInAccP->Divide(fHistPtTPCInAccnocutP);
2208 TH1F *fHistPtTPCeffInAccS = (TH1F*)fHistPtTPCInAccS->Clone("fHistPtTPCeffInAccS");
2209 fHistPtTPCeffInAccS->Divide(fHistPtTPCInAccnocutS);
2210 TH1F *fHistPtTPCeffInAccSfromStrange = (TH1F*)fHistPtTPCInAccSfromStrange->Clone("fHistPtTPCeffInAccSfromStrange");
2211 fHistPtTPCeffInAccSfromStrange->Divide(fHistPtTPCInAccnocutSfromStrange);
2212 TH1F *fHistPtTPCeffInAccSfromMat = (TH1F*)fHistPtTPCInAccSfromMat->Clone("fHistPtTPCeffInAccSfromMat");
2213 fHistPtTPCeffInAccSfromMat->Divide(fHistPtTPCInAccnocutSfromMat);
2214
2215
2216 TLegend *l0=new TLegend(0.5,0.5,0.9,0.9);
2217 TCanvas *c0=new TCanvas("c0","c0");
2218 c0->Divide(2,1);
2219 c0->cd(1);
2220 TH1F *fHistPtITSTPCselFracS = (TH1F*)fHistPtITSTPCselS->Clone("fHistPtITSTPCselFracS");
2221 fHistPtITSTPCselFracS->Divide(fHistPtITSTPCsel);
2222 fHistPtITSTPCselFracS->SetLineColor(2);
2223 l0->AddEntry(fHistPtITSTPCselFracS,"secondaries","l");
2224 fHistPtITSTPCselFracS->Draw();
2225 TH1F *fHistPtITSTPCselFracSfromStrange = (TH1F*)fHistPtITSTPCselSfromStrange->Clone("fHistPtITSTPCselFracSfromStrange");
2226 fHistPtITSTPCselFracSfromStrange->Divide(fHistPtITSTPCsel);
2227 fHistPtITSTPCselFracSfromStrange->SetLineColor(6);
2228 l0->AddEntry(fHistPtITSTPCselFracSfromStrange,"sec. from strange","l");
2229 fHistPtITSTPCselFracSfromStrange->Draw("same");
2230 TH1F *fHistPtITSTPCselFracSfromMat = (TH1F*)fHistPtITSTPCselSfromMat->Clone("fHistPtITSTPCselFracSfromMat");
2231 fHistPtITSTPCselFracSfromMat->Divide(fHistPtITSTPCsel);
2232 fHistPtITSTPCselFracSfromMat->SetLineColor(kOrange+7);
2233 l0->AddEntry(fHistPtITSTPCselFracSfromMat,"sec. from material","l");
2234 fHistPtITSTPCselFracSfromMat->Draw("same");
2235 l0->Draw();
2236 c0->cd(2);
2237 TH1F *fHistPtTPCInAccFracS = (TH1F*)fHistPtTPCInAccS->Clone("fHistPtTPCInAccFracS");
2238 fHistPtTPCInAccFracS->Divide(fHistPtTPCInAcc);
2239 fHistPtTPCInAccFracS->SetLineColor(2);
2240 fHistPtTPCInAccFracS->Draw();
2241 TH1F *fHistPtTPCInAccFracSfromStrange = (TH1F*)fHistPtTPCInAccSfromStrange->Clone("fHistPtTPCInAccFracSfromStrange");
2242 fHistPtTPCInAccFracSfromStrange->Divide(fHistPtTPCInAcc);
2243 fHistPtTPCInAccFracSfromStrange->SetLineColor(6);
2244 fHistPtTPCInAccFracSfromStrange->Draw("same");
2245 TH1F *fHistPtTPCInAccFracSfromMat = (TH1F*)fHistPtTPCInAccSfromMat->Clone("fHistPtTPCInAccFracSfromMat");
2246 fHistPtTPCInAccFracSfromMat->Divide(fHistPtTPCInAcc);
2247 fHistPtTPCInAccFracSfromMat->SetLineColor(kOrange+7);
2248 fHistPtTPCInAccFracSfromMat->Draw("same");
2249
2250 TLegend *l1=new TLegend(0.5,0.5,0.9,0.9);
2251 TCanvas *c1=new TCanvas("c1","c1");
2252 c1->Divide(2,1);
2253 c1->cd(1);
2254 fHistPtITSTPCeffselP->SetYTitle("dca cut efficiency");
2255 fHistPtITSTPCeffselP->SetLineColor(1);
2256 l1->AddEntry(fHistPtITSTPCeffselP,"primaries","l");
2257 fHistPtITSTPCeffselP->Draw();
2258 fHistPtITSTPCeffselS->SetLineColor(2);
2259 l1->AddEntry(fHistPtITSTPCeffselS,"secondaries","l");
2260 fHistPtITSTPCeffselS->Draw("same");
2261 fHistPtITSTPCeffselSfromStrange->SetLineColor(6);
2262 l1->AddEntry(fHistPtITSTPCeffselSfromStrange,"sec. from strange","l");
2263 fHistPtITSTPCeffselSfromStrange->Draw("same");
2264 fHistPtITSTPCeffselSfromMat->SetLineColor(kOrange+7);
2265 l1->AddEntry(fHistPtITSTPCeffselSfromMat,"sec. from material","l");
2266 fHistPtITSTPCeffselSfromMat->Draw("same");
2267 l1->Draw();
2268 c1->cd(2);
2269 fHistPtTPCeffInAccP->SetYTitle("dca cut efficiency");
2270 fHistPtTPCeffInAccP->SetLineColor(1);
2271 fHistPtTPCeffInAccP->Draw();
2272 fHistPtTPCeffInAccS->SetLineColor(2);
2273 fHistPtTPCeffInAccS->Draw("same");
2274 fHistPtTPCeffInAccSfromStrange->SetLineColor(6);
2275 fHistPtTPCeffInAccSfromStrange->Draw("same");
2276 fHistPtTPCeffInAccSfromMat->SetLineColor(kOrange+7);
2277 fHistPtTPCeffInAccSfromMat->Draw("same");
2278
2279
2280
2281 Float_t weightStrange=1.;
2282 Float_t weightMat=1.;
2283
2284 TH1F *fHistPtITSTPCselCocktail = MakeCocktail(fHistPtITSTPCselP,
2285 fHistPtITSTPCselPfromStrange,
2286 fHistPtITSTPCselS,
2287 fHistPtITSTPCselSfromStrange,
2288 fHistPtITSTPCselSfromMat,
2289 weightStrange,
2290 weightMat,
2291 "fHistPtITSTPCselCocktail");
2292
2293 TH1F *fHistPtITSTPCselPCocktail = CorrectToP(fHistPtITSTPCselCocktail,
2294 fHistPtITSTPCselSfromStrange,
2295 fHistPtITSTPCselSfromMat,
2296 weightStrange,
2297 weightMat,
2298 "fHistPtITSTPCselPCocktail");
2299
2300 TH1F *fHistPtITSTPCallPCocktail = CorrectToPnocut(fHistPtITSTPCselPCocktail,
2301 fHistPtITSTPCeffselP,
2302 "fHistPtITSTPCallPCocktail");
2303
2304
2305 TH1F *fHistPtTPCInAccCocktail = MakeCocktail(fHistPtTPCInAccP,
2306 fHistPtTPCInAccPfromStrange,
2307 fHistPtTPCInAccS,
2308 fHistPtTPCInAccSfromStrange,
2309 fHistPtTPCInAccSfromMat,
2310 weightStrange,
2311 weightMat,
2312 "fHistPtTPCInAccCocktail");
2313
2314 TH1F *fHistPtTPCInAccPCocktail = CorrectToP(fHistPtTPCInAccCocktail,
2315 fHistPtTPCInAccSfromStrange,
2316 fHistPtTPCInAccSfromMat,
2317 weightStrange,
2318 weightMat,
2319 "fHistPtTPCInAccPCocktail");
2320
2321 TH1F *fHistPtTPCallPCocktail = CorrectToPnocut(fHistPtTPCInAccPCocktail,
2322 fHistPtTPCeffInAccP,
2323 "fHistPtTPCallPCocktail");
2324
2325
2326
2327 TLegend *l2=new TLegend(0.5,0.5,0.9,0.9);
2328 TCanvas *c2 = new TCanvas("c2","c2");
2329 c2->Divide(2,1);
2330 c2->cd(1);
2331 fHistPtITSTPCselCocktail->SetLineColor(4);
2332 l2->AddEntry(fHistPtITSTPCselCocktail,"all, with dca cut","l");
2333 fHistPtITSTPCselCocktail->Draw();
2334 fHistPtITSTPCselPCocktail->SetLineColor(7);
2335 l2->AddEntry(fHistPtITSTPCselPCocktail,"primaries, with dca cut","l");
2336 fHistPtITSTPCselPCocktail->Draw("same");
2337 fHistPtITSTPCallPCocktail->SetLineColor(6);
2338 l2->AddEntry(fHistPtITSTPCallPCocktail,"primaries, w/o dca cut","l");
2339 fHistPtITSTPCallPCocktail->Draw("same");
2340 l2->Draw();
2341 c2->cd(2);
2342 fHistPtTPCInAccCocktail->SetLineColor(4);
2343 fHistPtTPCInAccCocktail->Draw();
2344 fHistPtTPCInAccPCocktail->SetLineColor(7);
2345 fHistPtTPCInAccPCocktail->Draw("same");
2346 fHistPtTPCallPCocktail->SetLineColor(6);
2347 fHistPtTPCallPCocktail->Draw("same");
2348
2349
2350 // systematic from secondaries
2351 weightStrange=1.3;
2352 weightMat=1.0;
2353 TH1F *fHistPtITSTPCselPCocktail1 = CorrectToP(fHistPtITSTPCselCocktail,
2354 fHistPtITSTPCselSfromStrange,
2355 fHistPtITSTPCselSfromMat,
2356 weightStrange,
2357 weightMat,
2358 "fHistPtITSTPCselPCocktail1");
2359 TH1F *fHistPtITSTPCallPCocktail1 = CorrectToPnocut(fHistPtITSTPCselPCocktail1,
2360 fHistPtITSTPCeffselP,
2361 "fHistPtITSTPCallPCocktail1");
2362 weightStrange=0.7;
2363 weightMat=1.0;
2364 TH1F *fHistPtITSTPCselPCocktail2 = CorrectToP(fHistPtITSTPCselCocktail,
2365 fHistPtITSTPCselSfromStrange,
2366 fHistPtITSTPCselSfromMat,
2367 weightStrange,
2368 weightMat,
2369 "fHistPtITSTPCselPCocktail2");
2370 TH1F *fHistPtITSTPCallPCocktail2 = CorrectToPnocut(fHistPtITSTPCselPCocktail2,
2371 fHistPtITSTPCeffselP,
2372 "fHistPtITSTPCallPCocktail2");
2373 weightStrange=1.0;
2374 weightMat=1.1;
2375 TH1F *fHistPtITSTPCselPCocktail3 = CorrectToP(fHistPtITSTPCselCocktail,
2376 fHistPtITSTPCselSfromStrange,
2377 fHistPtITSTPCselSfromMat,
2378 weightStrange,
2379 weightMat,
2380 "fHistPtITSTPCselPCocktail3");
2381 TH1F *fHistPtITSTPCallPCocktail3 = CorrectToPnocut(fHistPtITSTPCselPCocktail3,
2382 fHistPtITSTPCeffselP,
2383 "fHistPtITSTPCallPCocktail3");
2384 weightStrange=1.0;
2385 weightMat=0.9;
2386 TH1F *fHistPtITSTPCselPCocktail4 = CorrectToP(fHistPtITSTPCselCocktail,
2387 fHistPtITSTPCselSfromStrange,
2388 fHistPtITSTPCselSfromMat,
2389 weightStrange,
2390 weightMat,
2391 "fHistPtITSTPCselPCocktail4");
2392 TH1F *fHistPtITSTPCallPCocktail4 = CorrectToPnocut(fHistPtITSTPCselPCocktail4,
2393 fHistPtITSTPCeffselP,
2394 "fHistPtITSTPCallPCocktail4");
2395 weightStrange=1.3;
2396 weightMat=1.1;
2397 TH1F *fHistPtITSTPCselPCocktail5 = CorrectToP(fHistPtITSTPCselCocktail,
2398 fHistPtITSTPCselSfromStrange,
2399 fHistPtITSTPCselSfromMat,
2400 weightStrange,
2401 weightMat,
2402 "fHistPtITSTPCselPCocktail5");
2403 TH1F *fHistPtITSTPCallPCocktail5 = CorrectToPnocut(fHistPtITSTPCselPCocktail5,
2404 fHistPtITSTPCeffselP,
2405 "fHistPtITSTPCallPCocktail5");
2406 weightStrange=0.7;
2407 weightMat=0.9;
2408 TH1F *fHistPtITSTPCselPCocktail6 = CorrectToP(fHistPtITSTPCselCocktail,
2409 fHistPtITSTPCselSfromStrange,
2410 fHistPtITSTPCselSfromMat,
2411 weightStrange,
2412 weightMat,
2413 "fHistPtITSTPCselPCocktail6");
2414 TH1F *fHistPtITSTPCallPCocktail6 = CorrectToPnocut(fHistPtITSTPCselPCocktail6,
2415 fHistPtITSTPCeffselP,
2416 "fHistPtITSTPCallPCocktail6");
2417
2418 weightStrange=1.3;
2419 weightMat=1.0;
2420 TH1F *fHistPtTPCInAccPCocktail1 = CorrectToP(fHistPtTPCInAccCocktail,
2421 fHistPtTPCInAccSfromStrange,
2422 fHistPtTPCInAccSfromMat,
2423 weightStrange,
2424 weightMat,
2425 "fHistPtTPCInAccPCocktail1");
2426 TH1F *fHistPtTPCallPCocktail1 = CorrectToPnocut(fHistPtTPCInAccPCocktail1,
2427 fHistPtTPCeffInAccP,
2428 "fHistPtTPCallPCocktail1");
2429 weightStrange=0.7;
2430 weightMat=1.0;
2431 TH1F *fHistPtTPCInAccPCocktail2 = CorrectToP(fHistPtTPCInAccCocktail,
2432 fHistPtTPCInAccSfromStrange,
2433 fHistPtTPCInAccSfromMat,
2434 weightStrange,
2435 weightMat,
2436 "fHistPtTPCInAccPCocktail2");
2437 TH1F *fHistPtTPCallPCocktail2 = CorrectToPnocut(fHistPtTPCInAccPCocktail2,
2438 fHistPtTPCeffInAccP,
2439 "fHistPtTPCallPCocktail2");
2440 weightStrange=1.0;
2441 weightMat=1.1;
2442 TH1F *fHistPtTPCInAccPCocktail3 = CorrectToP(fHistPtTPCInAccCocktail,
2443 fHistPtTPCInAccSfromStrange,
2444 fHistPtTPCInAccSfromMat,
2445 weightStrange,
2446 weightMat,
2447 "fHistPtTPCInAccPCocktail3");
2448 TH1F *fHistPtTPCallPCocktail3 = CorrectToPnocut(fHistPtTPCInAccPCocktail3,
2449 fHistPtTPCeffInAccP,
2450 "fHistPtTPCallPCocktail3");
2451 weightStrange=1.0;
2452 weightMat=0.9;
2453 TH1F *fHistPtTPCInAccPCocktail4 = CorrectToP(fHistPtTPCInAccCocktail,
2454 fHistPtTPCInAccSfromStrange,
2455 fHistPtTPCInAccSfromMat,
2456 weightStrange,
2457 weightMat,
2458 "fHistPtTPCInAccPCocktail4");
2459 TH1F *fHistPtTPCallPCocktail4 = CorrectToPnocut(fHistPtTPCInAccPCocktail4,
2460 fHistPtTPCeffInAccP,
2461 "fHistPtTPCallPCocktail4");
2462 weightStrange=1.3;
2463 weightMat=1.1;
2464 TH1F *fHistPtTPCInAccPCocktail5 = CorrectToP(fHistPtTPCInAccCocktail,
2465 fHistPtTPCInAccSfromStrange,
2466 fHistPtTPCInAccSfromMat,
2467 weightStrange,
2468 weightMat,
2469 "fHistPtTPCInAccPCocktail5");
2470 TH1F *fHistPtTPCallPCocktail5 = CorrectToPnocut(fHistPtTPCInAccPCocktail5,
2471 fHistPtTPCeffInAccP,
2472 "fHistPtTPCallPCocktail5");
2473 weightStrange=0.7;
2474 weightMat=0.9;
2475 TH1F *fHistPtTPCInAccPCocktail6 = CorrectToP(fHistPtTPCInAccCocktail,
2476 fHistPtTPCInAccSfromStrange,
2477 fHistPtTPCInAccSfromMat,
2478 weightStrange,
2479 weightMat,
2480 "fHistPtTPCInAccPCocktail6");
2481 TH1F *fHistPtTPCallPCocktail6 = CorrectToPnocut(fHistPtTPCInAccPCocktail6,
2482 fHistPtTPCeffInAccP,
2483 "fHistPtTPCallPCocktail6");
2484
2485 TLegend *l3=new TLegend(0.5,0.5,0.9,0.9);
2486 TCanvas *c3 = new TCanvas("c3","c3");
2487 c3->Divide(2,1);
2488 c3->cd(1);
2489 l3->AddEntry(fHistPtITSTPCallPCocktail,"central","l");
2490 fHistPtITSTPCallPCocktail->Draw();
2491 l3->AddEntry(fHistPtITSTPCallPCocktail1,"strange #pm 30%","l");
2492 fHistPtITSTPCallPCocktail1->SetLineColor(1);
2493 fHistPtITSTPCallPCocktail2->SetLineColor(1);
2494 fHistPtITSTPCallPCocktail1->Draw("same");
2495 fHistPtITSTPCallPCocktail2->Draw("same");
2496 l3->AddEntry(fHistPtITSTPCallPCocktail3,"material #pm 10%","l");
2497 fHistPtITSTPCallPCocktail3->SetLineColor(8);
2498 fHistPtITSTPCallPCocktail4->SetLineColor(8);
2499 fHistPtITSTPCallPCocktail3->Draw("same");
2500 fHistPtITSTPCallPCocktail4->Draw("same");
2501 l3->AddEntry(fHistPtITSTPCallPCocktail5,"strange #pm 30% & material #pm 10%","l");
2502 fHistPtITSTPCallPCocktail5->SetLineColor(4);
2503 fHistPtITSTPCallPCocktail6->SetLineColor(4);
2504 fHistPtITSTPCallPCocktail5->Draw("same");
2505 fHistPtITSTPCallPCocktail6->Draw("same");
2506 l3->Draw();
2507 c3->cd(2);
2508 fHistPtTPCallPCocktail->Draw();
2509 fHistPtTPCallPCocktail1->SetLineColor(1);
2510 fHistPtTPCallPCocktail2->SetLineColor(1);
2511 fHistPtTPCallPCocktail1->Draw("same");
2512 fHistPtTPCallPCocktail2->Draw("same");
2513 fHistPtTPCallPCocktail3->SetLineColor(8);
2514 fHistPtTPCallPCocktail4->SetLineColor(8);
2515 fHistPtTPCallPCocktail3->Draw("same");
2516 fHistPtTPCallPCocktail4->Draw("same");
2517 fHistPtTPCallPCocktail5->SetLineColor(4);
2518 fHistPtTPCallPCocktail6->SetLineColor(4);
2519 fHistPtTPCallPCocktail5->Draw("same");
2520 fHistPtTPCallPCocktail6->Draw("same");
2521
2522 /*
2523 TCanvas *c4 = new TCanvas("c4","c4");
2524 c4->Divide(2,1);
2525 c4->cd(1);
2526 fHistPtITSTPCallPCocktail2->Add(fHistPtITSTPCallPCocktail,-1);
2527 fHistPtITSTPCallPCocktail2->Divide(fHistPtITSTPCallPCocktail);
2528 fHistPtITSTPCallPCocktail2->Draw();
2529 fHistPtITSTPCallPCocktail4->Add(fHistPtITSTPCallPCocktail,-1);
2530 fHistPtITSTPCallPCocktail4->Divide(fHistPtITSTPCallPCocktail);
2531 fHistPtITSTPCallPCocktail4->Draw("same");
2532 fHistPtITSTPCallPCocktail6->Add(fHistPtITSTPCallPCocktail,-1);
2533 fHistPtITSTPCallPCocktail6->Divide(fHistPtITSTPCallPCocktail);
2534 fHistPtITSTPCallPCocktail6->Draw("same");
2535 l3->Draw();
2536 c4->cd(2);
2537 fHistPtTPCallPCocktail2->Add(fHistPtTPCallPCocktail,-1);
2538 fHistPtTPCallPCocktail2->Divide(fHistPtTPCallPCocktail);
2539 fHistPtTPCallPCocktail2->Draw();
2540 fHistPtTPCallPCocktail4->Add(fHistPtTPCallPCocktail,-1);
2541 fHistPtTPCallPCocktail4->Divide(fHistPtTPCallPCocktail);
2542 fHistPtTPCallPCocktail4->Draw("same");
2543 fHistPtTPCallPCocktail6->Add(fHistPtTPCallPCocktail,-1);
2544 fHistPtTPCallPCocktail6->Divide(fHistPtTPCallPCocktail);
2545 fHistPtTPCallPCocktail6->Draw("same");
2546 */
2547
2548 return;
2549}
2550//-------------------------------------------------------------
2551void NormalizePtHist(TH1F *h,Float_t n) {
2552
2553 for(Int_t i=1; i<=h->GetNbinsX(); i++) {
2554 if(h->GetBinContent(i)) {
2555 h->SetBinContent(i,h->GetBinContent(i)/h->GetBinWidth(i));
2556 h->SetBinError(i,h->GetBinError(i)/h->GetBinWidth(i));
2557 }
2558 }
2559 h->Scale(1./n);
2560
2561 return;
2562}
2563//----------------------------------------------------------------------
2564TH1F *MakeCocktail(TH1F *fHistPtselP,
2565 TH1F *fHistPtselPfromStrange,
2566 TH1F *fHistPtselS,
2567 TH1F *fHistPtselSfromStrange,
2568 TH1F *fHistPtselSfromMat,
2569 Float_t weightStrange,
2570 Float_t weightMat,TString name) {
2571
2572 // build histogram with yield passing cut
2573 TH1F *fHistPtselCocktail = (TH1F*)fHistPtselP->Clone(name.Data());
2574 fHistPtselCocktail->Add(fHistPtselS);
2575 fHistPtselCocktail->Add(fHistPtselPfromStrange,weightStrange-1.);
2576 fHistPtselCocktail->Add(fHistPtselSfromStrange,weightStrange-1.);
2577 fHistPtselCocktail->Add(fHistPtselSfromMat,weightMat-1.);
2578
2579
2580 return fHistPtselCocktail;
2581}
2582//----------------------------------------------------------------------
2583TH1F *CorrectToP(TH1F *fHistPtselCocktail,
2584 TH1F *fHistPtselSfromStrange,
2585 TH1F *fHistPtselSfromMat,
2586 Float_t weightStrange,
2587 Float_t weightMat,TString name) {
2588
2589
2590 // build histogram with yield of primaries passing cut
2591 TH1F* fHistPtselPCocktail = (TH1F*)fHistPtselCocktail->Clone(name.Data());
2592 fHistPtselPCocktail->Add(fHistPtselSfromStrange,-weightStrange);
2593 fHistPtselPCocktail->Add(fHistPtselSfromMat,-weightMat);
2594
2595 return fHistPtselPCocktail;
2596}
2597//----------------------------------------------------------------------
2598TH1F *CorrectToPnocut(TH1F *fHistPtselPCocktail,
2599 TH1F *fHistPteffselP,TString name) {
2600
2601 // build histogram with yield of primaries without cut
2602 TH1F *fHistPtallPCocktail = (TH1F*)fHistPtselPCocktail->Clone(name.Data());
2603 //printf("%f\n",fHistPtallPCocktail->GetBinContent(6));
2604 fHistPtallPCocktail->Divide(fHistPteffselP);
2605 //printf("%f\n",fHistPtallPCocktail->GetBinContent(6));
2606
2607 return fHistPtallPCocktail;
2608}
2609//----------------------------------------------------------------------
2610void CompareTPCPt() {
2611
2612 TFile *f_tpcrefit= new TFile("ITS.Performance_104892pass2_etale05.root");
2613
2614 TList *list=(TList*)f_tpcrefit->Get("cOutputITS");
2615
2616 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
2617 Float_t nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
2618
2619 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
2620 NormalizePtHist(fHistPtTPCInAcc,nEventsSelWithVertex);
2621 TH1F *fHistPtTPCInAccNoTRDout = (TH1F*)list->FindObject("fHistPtTPCInAccNoTRDout");
2622 NormalizePtHist(fHistPtTPCInAccNoTRDout,nEventsSelWithVertex);
2623 TH1F *fHistPtTPCInAccNoTOFout = (TH1F*)list->FindObject("fHistPtTPCInAccNoTOFout");
2624 NormalizePtHist(fHistPtTPCInAccNoTOFout,nEventsSelWithVertex);
2625 TH1F *fHistPtTPCInAccWithPtTPCAtVtx = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtVtx");
2626 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtVtx,nEventsSelWithVertex);
2627 TH1F *fHistPtTPCInAccWithPtTPCAtInnerWall = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtInnerWall");
2628 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtInnerWall,nEventsSelWithVertex);
2629
2630
2631 TFile *f_notpcrefit= new TFile("ITS.Performance_104892pass3_etale05.root");
2632
2633 TList *list=(TList*)f_notpcrefit->Get("cOutputITS");
2634
2635 TH1F *fHistNEventsNoRefit = (TH1F*)list->FindObject("fHistNEvents");
2636 nEventsSelWithVertex = fHistNEventsNoRefit->GetBinContent(5);
2637
2638 TH1F *fHistPtTPCInAccNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAcc");
2639 NormalizePtHist(fHistPtTPCInAccNoRefit,nEventsSelWithVertex);
2640 TH1F *fHistPtTPCInAccNoTRDoutNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccNoTRDout");
2641 NormalizePtHist(fHistPtTPCInAccNoTRDoutNoRefit,nEventsSelWithVertex);
2642 TH1F *fHistPtTPCInAccNoTOFoutNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccNoTOFout");
2643 NormalizePtHist(fHistPtTPCInAccNoTOFoutNoRefit,nEventsSelWithVertex);
2644 TH1F *fHistPtTPCInAccWithPtTPCAtVtxNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtVtx");
2645 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtVtxNoRefit,nEventsSelWithVertex);
2646 TH1F *fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtInnerWall");
2647 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit,nEventsSelWithVertex);
2648
2649
2650
2651 TLegend *l1 = new TLegend(0.5,0.5,0.9,0.9);
2652 TCanvas *c1=new TCanvas("c1","c1",0,0,1000,600);
2653 c1->Divide(2,1);
2654 c1->cd(1);
2655 fHistPtTPCInAcc->SetMarkerStyle(20);
2656 l1->AddEntry(fHistPtTPCInAcc,"p_{t} measurement:","");
2657 l1->AddEntry(fHistPtTPCInAcc,"at vtx, from main params","p");
2658 fHistPtTPCInAcc->Draw();
2659 fHistPtTPCInAccWithPtTPCAtVtx->SetMarkerStyle(21);
2660 fHistPtTPCInAccWithPtTPCAtVtx->SetMarkerColor(2);
2661 l1->AddEntry(fHistPtTPCInAccWithPtTPCAtVtx,"at vtx, from TPC params","p");
2662 fHistPtTPCInAccWithPtTPCAtVtx->Draw("same");
2663 fHistPtTPCInAccWithPtTPCAtInnerWall->SetMarkerStyle(25);
2664 fHistPtTPCInAccWithPtTPCAtInnerWall->SetMarkerColor(2);
2665 l1->AddEntry(fHistPtTPCInAccWithPtTPCAtInnerWall,"at r=80cm, from TPC params","p");
2666 fHistPtTPCInAccWithPtTPCAtInnerWall->Draw("same");
2667 fHistPtTPCInAccNoTOFout->SetMarkerStyle(23);
2668 fHistPtTPCInAccNoTOFout->SetMarkerColor(3);
2669 fHistPtTPCInAccNoTOFout->Draw("same");
2670 fHistPtTPCInAccNoTRDout->SetMarkerStyle(24);
2671 fHistPtTPCInAccNoTRDout->SetMarkerColor(4);
2672 fHistPtTPCInAccNoTRDout->Draw("same");
2673 l1->Draw();
2674 c1->cd(2);
2675 fHistPtTPCInAccNoRefit->SetMarkerStyle(20);
2676 fHistPtTPCInAccNoRefit->Draw();
2677 fHistPtTPCInAccWithPtTPCAtVtxNoRefit->SetMarkerStyle(21);
2678 fHistPtTPCInAccWithPtTPCAtVtxNoRefit->SetMarkerColor(2);
2679 fHistPtTPCInAccWithPtTPCAtVtxNoRefit->Draw("same");
2680 fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit->SetMarkerStyle(25);
2681 fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit->SetMarkerColor(2);
2682 fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit->Draw("same");
2683 fHistPtTPCInAccNoTOFoutNoRefit->SetMarkerStyle(23);
2684 fHistPtTPCInAccNoTOFoutNoRefit->SetMarkerColor(3);
2685 fHistPtTPCInAccNoTOFoutNoRefit->Draw("same");
2686 fHistPtTPCInAccNoTRDoutNoRefit->SetMarkerStyle(24);
2687 fHistPtTPCInAccNoTRDoutNoRefit->SetMarkerColor(4);
2688 fHistPtTPCInAccNoTRDoutNoRefit->Draw("same");
2689 //l1->Draw();
2690
2691 return;
2692}
2693
2694TH1F* gHistTemplateStrange;
2695TH1F* gHistTemplatePandMat;
2696Float_t gd0min=0.15;
2697Float_t gd0max=1.5;
2698
2699
2700Double_t FitMCTemplate(Double_t *x, Double_t *par) {
2701 Double_t xx = x[0];
2702 if(TMath::Abs(xx)<gd0min) return 0;
2703 Int_t bin = gHistTemplateStrange->GetXaxis()->FindBin(xx);
2704 Double_t pandmat = par[0]*gHistTemplatePandMat->GetBinContent(bin);
2705 Double_t strange = par[1]*gHistTemplateStrange->GetBinContent(bin);
2706 return strange + pandmat;
2707}
2708
2709void FitImpPar_rphi() {
2710
2711 Int_t rebin=1;
2712
2713 TFile *fMC= new TFile("ITS.Performance_lhc10a8a12_trkvtx_noTPCdca_120310.root");
2714
2715 TList *list=(TList*)fMC->Get("cOutputITS");
2716 TH1F *fHistd0rphiITSMIoneSPDInAccPMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
2717 TH1F *fHistd0rphiITSMIoneSPDInAccSMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
2718 TH1F *fHistd0rphiITSMIoneSPDInAccSfromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200fromMat");
2719 fHistd0rphiITSMIoneSPDInAccPMC->Rebin(rebin);
2720 fHistd0rphiITSMIoneSPDInAccSMC->Rebin(rebin);
2721 fHistd0rphiITSMIoneSPDInAccSfromMatMC->Rebin(rebin);
2722
2723 TH1F *fHistd0rphiITSMIoneSPDInAccSfromStrangeMC=(TH1F*)fHistd0rphiITSMIoneSPDInAccSMC->Clone("fHistd0rphiITSMIoneSPDInAccSfromStrangeMC");
2724 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Add(fHistd0rphiITSMIoneSPDInAccSfromMatMC,-1.);
2725 TH1F *fHistd0rphiITSMIoneSPDInAccMC=(TH1F*)fHistd0rphiITSMIoneSPDInAccPMC->Clone("fHistd0rphiITSMIoneSPDInAccMC");
2726 fHistd0rphiITSMIoneSPDInAccMC->Add(fHistd0rphiITSMIoneSPDInAccSMC);
2727 //fHistd0rphiITSMIoneSPDInAccMC->Scale(1./fHistd0rphiITSMIoneSPDInAccMC->GetEntries());
2728 TH1F *fHistd0rphiITSMIoneSPDInAccPandMatMC=(TH1F*)fHistd0rphiITSMIoneSPDInAccPMC->Clone("fHistd0rphiITSMIoneSPDInAccPandMatMC");
2729 fHistd0rphiITSMIoneSPDInAccPandMatMC->Add(fHistd0rphiITSMIoneSPDInAccSfromMatMC);
2730
2731 // templates
2732 gHistTemplateStrange = (TH1F*)fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Clone("gHistTemplateStrange");
2733 gHistTemplatePandMat = (TH1F*)fHistd0rphiITSMIoneSPDInAccPandMatMC->Clone("gHistTemplatePandMat");
2734
2735
2736
2737
2738 TFile *f= new TFile("ITS.Performance_104065_104892_pass4_trkvtx_noTPCdca.root");
2739 //TFile *f= new TFile("ITS.Performance_lhc10a8_300k_trkvtx_noTPCdca.root");
2740
2741 TList *list=(TList*)f->Get("cOutputITS");
2742 TH1F *fHistd0rphiITSMIoneSPDInAccP = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
2743 TH1F *fHistd0rphiITSMIoneSPDInAccS = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
2744 fHistd0rphiITSMIoneSPDInAccP->Rebin(rebin);
2745 fHistd0rphiITSMIoneSPDInAccS->Rebin(rebin);
2746
2747 TH1F *fHistd0rphiITSMIoneSPDInAcc=(TH1F*)fHistd0rphiITSMIoneSPDInAccP->Clone("fHistd0rphiITSMIoneSPDInAcc");
2748 fHistd0rphiITSMIoneSPDInAcc->Add(fHistd0rphiITSMIoneSPDInAccS);
2749
2750
2751 TF1 *fitMCtemplate = new TF1("fitMCtemplate",FitMCTemplate,-gd0max,+gd0max,2);
2752
2753 TH1F* fHistd0rphiITSMIoneSPDInAccCopy = (TH1F*)fHistd0rphiITSMIoneSPDInAcc->Clone("fHistd0rphiITSMIoneSPDInAccCopy");
2754 for(Int_t ib=1;ib<=fHistd0rphiITSMIoneSPDInAccCopy->GetNbinsX();ib++) {
2755 if(TMath::Abs(fHistd0rphiITSMIoneSPDInAccCopy->GetBinCenter(ib))<gd0min) {
2756 fHistd0rphiITSMIoneSPDInAccCopy->SetBinContent(ib,0);
2757 fHistd0rphiITSMIoneSPDInAccCopy->SetBinError(ib,0);
2758 }
2759 }
2760 fHistd0rphiITSMIoneSPDInAccCopy->Fit("fitMCtemplate","RL");
2761
2762
2763 fHistd0rphiITSMIoneSPDInAccPandMatMC->Scale(fitMCtemplate->GetParameter(0));
2764 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Scale(fitMCtemplate->GetParameter(1));
2765 TH1F *fHistd0rphiITSMIoneSPDInAccMCFit=(TH1F*)fHistd0rphiITSMIoneSPDInAccPandMatMC->Clone("fHistd0rphiITSMIoneSPDInAccMCFit");
2766 fHistd0rphiITSMIoneSPDInAccMCFit->Add(fHistd0rphiITSMIoneSPDInAccSfromStrangeMC);
2767
2768 TCanvas *c1 = new TCanvas("c1","c1",1000,1000);
2769 c1->SetLogy();
2770 fHistd0rphiITSMIoneSPDInAccCopy->SetMinimum(1);
2771 fHistd0rphiITSMIoneSPDInAccCopy->SetMaximum(2*fHistd0rphiITSMIoneSPDInAcc->GetMaximum());
2772 fHistd0rphiITSMIoneSPDInAcc->SetMinimum(1);
2773 fHistd0rphiITSMIoneSPDInAccMC->SetMinimum(1);
2774 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->SetMinimum(1);
2775 fHistd0rphiITSMIoneSPDInAccPandMatMC->SetMinimum(1);
2776 fHistd0rphiITSMIoneSPDInAccCopy->SetLineColor(4);
2777 fHistd0rphiITSMIoneSPDInAccCopy->Draw();
2778 fHistd0rphiITSMIoneSPDInAcc->SetLineColor(4);
2779 fHistd0rphiITSMIoneSPDInAcc->Draw("same");
2780 fHistd0rphiITSMIoneSPDInAccMCFit->SetLineColor(2);
2781 fHistd0rphiITSMIoneSPDInAccMCFit->Draw("same");
2782 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->SetLineColor(3);
2783 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Draw("same");
2784 fHistd0rphiITSMIoneSPDInAccPandMatMC->SetLineColor(15);
2785 fHistd0rphiITSMIoneSPDInAccPandMatMC->Draw("same");
2786
2787
2788
2789 TCanvas *c3 = new TCanvas("c3","c3");
2790 c3->SetLogx();
2791 c3->SetGridx();
2792 c3->SetGridy();
2793 Float_t d0cut[15]={0.0301,0.0401,0.0501,0.0601,0.0801,0.101,0.151,0.201,0.251,0.301,0.401,0.601,0.801,1.001,1.401};
2794
2795 Float_t fracSfromStrange[15],fracSfromStrangeMC[15],fracSfromMatMC[15],fracS[15],fracSMC[15],fracSRatio[15];
2796 for(Int_t i=0;i<15;i++) {
2797 Int_t bin1=fHistd0rphiITSMIoneSPDInAccPandMatMC->FindBin(-d0cut[i]);
2798 Int_t bin2=fHistd0rphiITSMIoneSPDInAccPandMatMC->FindBin(+d0cut[i]);
2799 Float_t intDatacut=fHistd0rphiITSMIoneSPDInAcc->Integral(bin1,bin2);
2800 Float_t intMCcut=fHistd0rphiITSMIoneSPDInAccMC->Integral(bin1,bin2);
2801 Float_t intSfromStrangecut=fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Integral(bin1,bin2);
2802 Float_t intSfromMatcut=fHistd0rphiITSMIoneSPDInAccSfromMatMC->Integral(bin1,bin2);
2803 fracSfromStrange[i]=intSfromStrangecut/intDatacut;
2804 fracSfromStrangeMC[i]=intSfromStrangecut/fitMCtemplate->GetParameter(1)/intMCcut;
2805 fracSfromMatMC[i]=intSfromMatcut/intMCcut;
2806
2807 fracS[i]=fracSfromStrange[i]+fracSfromMatMC[i];
2808 fracSMC[i]=fracSfromStrangeMC[i]+fracSfromMatMC[i];
2809
2810 fracSRatio[i]=fracS[i]/fracSMC[i];
2811 }
2812 TGraph *gfracSfromStrange=new TGraph(15,d0cut,fracSfromStrange);
2813 gfracSfromStrange->SetMarkerColor(2);
2814 gfracSfromStrange->SetMarkerStyle(20);
2815 gfracSfromStrange->Draw("ap");
2816
2817 TGraph *gfracSfromStrangeMC=new TGraph(15,d0cut,fracSfromStrangeMC);
2818 gfracSfromStrangeMC->SetMarkerColor(1);
2819 gfracSfromStrangeMC->SetMarkerStyle(22);
2820 gfracSfromStrangeMC->Draw("p");
2821
2822 TGraph *gfracSfromMatMC=new TGraph(15,d0cut,fracSfromMatMC);
2823 gfracSfromMatMC->SetMarkerColor(4);
2824 gfracSfromMatMC->SetMarkerStyle(21);
2825 gfracSfromMatMC->Draw("p");
2826
2827 TGraph *gfracS=new TGraph(15,d0cut,fracS);
2828 gfracS->SetMarkerColor(2);
2829 gfracS->SetMarkerStyle(24);
2830 gfracS->Draw("p");
2831
2832 TGraph *gfracSMC=new TGraph(15,d0cut,fracSMC);
2833 gfracSMC->SetMarkerColor(1);
2834 gfracSMC->SetMarkerStyle(25);
2835 gfracSMC->Draw("p");
2836
2837
2838 TCanvas *c4 = new TCanvas("c4","c4");
2839 c4->SetLogx();
2840 c4->SetGridx();
2841 c4->SetGridy();
2842
2843 TGraph *gfracSRatio=new TGraph(15,d0cut,fracSRatio);
2844 gfracSRatio->SetMarkerColor(9);
2845 gfracSRatio->SetMarkerStyle(21);
2846 gfracSRatio->Draw("ap");
2847
2848
2849 return;
2850}
2851//-----------------------------------------------------------------------------
2852void PlotEffVSErrors() {
2853
2854 TFile *f1= new TFile("eff_104892_pd_testnew1.root");
2855
2856 TH1F *fHistPtITSMI6InAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMI6InAcc");
2857 TH1F *fHistPtITSMIge2InAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMIge2InAcc");
2858 TH1F *fHistPtITSMISPDInAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMISPDInAcc");
2859 TH1F *fHistPtITSMIoneSPDInAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMIoneSPDInAcc");
2860
2861 fHistPtITSMI6InAcctestnew1->SetMarkerStyle(20);
2862 fHistPtITSMIge2InAcctestnew1->SetMarkerStyle(20);
2863 fHistPtITSMISPDInAcctestnew1->SetMarkerStyle(20);
2864 fHistPtITSMIoneSPDInAcctestnew1->SetMarkerStyle(20);
2865 fHistPtITSMI6InAcctestnew1->SetMarkerColor(2);
2866 fHistPtITSMIge2InAcctestnew1->SetMarkerColor(1);
2867 fHistPtITSMISPDInAcctestnew1->SetMarkerColor(9);
2868 fHistPtITSMIoneSPDInAcctestnew1->SetMarkerColor(15);
2869
2870
2871 TFile *f2= new TFile("eff_104892_pd_testnew2.root");
2872
2873 TH1F *fHistPtITSMI6InAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMI6InAcc");
2874 TH1F *fHistPtITSMIge2InAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMIge2InAcc");
2875 TH1F *fHistPtITSMISPDInAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMISPDInAcc");
2876 TH1F *fHistPtITSMIoneSPDInAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMIoneSPDInAcc");
2877
2878 fHistPtITSMI6InAcctestnew2->SetMarkerStyle(21);
2879 fHistPtITSMIge2InAcctestnew2->SetMarkerStyle(21);
2880 fHistPtITSMISPDInAcctestnew2->SetMarkerStyle(21);
2881 fHistPtITSMIoneSPDInAcctestnew2->SetMarkerStyle(21);
2882 fHistPtITSMI6InAcctestnew2->SetMarkerColor(2);
2883 fHistPtITSMIge2InAcctestnew2->SetMarkerColor(1);
2884 fHistPtITSMISPDInAcctestnew2->SetMarkerColor(9);
2885 fHistPtITSMIoneSPDInAcctestnew2->SetMarkerColor(15);
2886
2887
2888 TFile *f3= new TFile("eff_104892_pd_testnew3.root");
2889
2890 TH1F *fHistPtITSMI6InAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMI6InAcc");
2891 TH1F *fHistPtITSMIge2InAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMIge2InAcc");
2892 TH1F *fHistPtITSMISPDInAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMISPDInAcc");
2893 TH1F *fHistPtITSMIoneSPDInAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMIoneSPDInAcc");
2894
2895 fHistPtITSMI6InAcctestnew3->SetMarkerStyle(22);
2896 fHistPtITSMIge2InAcctestnew3->SetMarkerStyle(22);
2897 fHistPtITSMISPDInAcctestnew3->SetMarkerStyle(22);
2898 fHistPtITSMIoneSPDInAcctestnew3->SetMarkerStyle(22);
2899 fHistPtITSMI6InAcctestnew3->SetMarkerColor(2);
2900 fHistPtITSMIge2InAcctestnew3->SetMarkerColor(1);
2901 fHistPtITSMISPDInAcctestnew3->SetMarkerColor(9);
2902 fHistPtITSMIoneSPDInAcctestnew3->SetMarkerColor(15);
2903
2904
2905 TFile *f5= new TFile("eff_104892_pd_testnew5.root");
2906
2907 TH1F *fHistPtITSMI6InAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMI6InAcc");
2908 TH1F *fHistPtITSMIge2InAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMIge2InAcc");
2909 TH1F *fHistPtITSMISPDInAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMISPDInAcc");
2910 TH1F *fHistPtITSMIoneSPDInAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMIoneSPDInAcc");
2911
2912 fHistPtITSMI6InAcctestnew5->SetMarkerStyle(24);
2913 fHistPtITSMIge2InAcctestnew5->SetMarkerStyle(24);
2914 fHistPtITSMISPDInAcctestnew5->SetMarkerStyle(24);
2915 fHistPtITSMIoneSPDInAcctestnew5->SetMarkerStyle(24);
2916 fHistPtITSMI6InAcctestnew5->SetMarkerColor(2);
2917 fHistPtITSMIge2InAcctestnew5->SetMarkerColor(1);
2918 fHistPtITSMISPDInAcctestnew5->SetMarkerColor(9);
2919 fHistPtITSMIoneSPDInAcctestnew5->SetMarkerColor(15);
2920
2921
2922 TFile *f6= new TFile("eff_104892_pd_testnew6.root");
2923
2924 TH1F *fHistPtITSMI6InAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMI6InAcc");
2925 TH1F *fHistPtITSMIge2InAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMIge2InAcc");
2926 TH1F *fHistPtITSMISPDInAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMISPDInAcc");
2927 TH1F *fHistPtITSMIoneSPDInAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMIoneSPDInAcc");
2928
2929 fHistPtITSMI6InAcctestnew6->SetMarkerStyle(25);
2930 fHistPtITSMIge2InAcctestnew6->SetMarkerStyle(25);
2931 fHistPtITSMISPDInAcctestnew6->SetMarkerStyle(25);
2932 fHistPtITSMIoneSPDInAcctestnew6->SetMarkerStyle(25);
2933 fHistPtITSMI6InAcctestnew6->SetMarkerColor(2);
2934 fHistPtITSMIge2InAcctestnew6->SetMarkerColor(1);
2935 fHistPtITSMISPDInAcctestnew6->SetMarkerColor(9);
2936 fHistPtITSMIoneSPDInAcctestnew6->SetMarkerColor(15);
2937
2938
2939 TFile *f7= new TFile("eff_104892_pd_testnew7.root");
2940
2941 TH1F *fHistPtITSMI6InAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMI6InAcc");
2942 TH1F *fHistPtITSMIge2InAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMIge2InAcc");
2943 TH1F *fHistPtITSMISPDInAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMISPDInAcc");
2944 TH1F *fHistPtITSMIoneSPDInAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMIoneSPDInAcc");
2945
2946 fHistPtITSMI6InAcctestnew7->SetMarkerStyle(23);
2947 fHistPtITSMIge2InAcctestnew7->SetMarkerStyle(23);
2948 fHistPtITSMISPDInAcctestnew7->SetMarkerStyle(23);
2949 fHistPtITSMIoneSPDInAcctestnew7->SetMarkerStyle(23);
2950 fHistPtITSMI6InAcctestnew7->SetMarkerColor(2);
2951 fHistPtITSMIge2InAcctestnew7->SetMarkerColor(1);
2952 fHistPtITSMISPDInAcctestnew7->SetMarkerColor(9);
2953 fHistPtITSMIoneSPDInAcctestnew7->SetMarkerColor(15);
2954
2955
2956 TFile *f8= new TFile("eff_104892_pd_testnew8.root");
2957
2958 TH1F *fHistPtITSMI6InAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMI6InAcc");
2959 TH1F *fHistPtITSMIge2InAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMIge2InAcc");
2960 TH1F *fHistPtITSMISPDInAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMISPDInAcc");
2961 TH1F *fHistPtITSMIoneSPDInAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMIoneSPDInAcc");
2962
2963 fHistPtITSMI6InAcctestnew8->SetMarkerStyle(28);
2964 fHistPtITSMIge2InAcctestnew8->SetMarkerStyle(28);
2965 fHistPtITSMISPDInAcctestnew8->SetMarkerStyle(28);
2966 fHistPtITSMIoneSPDInAcctestnew8->SetMarkerStyle(28);
2967 fHistPtITSMI6InAcctestnew8->SetMarkerColor(2);
2968 fHistPtITSMIge2InAcctestnew8->SetMarkerColor(1);
2969 fHistPtITSMISPDInAcctestnew8->SetMarkerColor(9);
2970 fHistPtITSMIoneSPDInAcctestnew8->SetMarkerColor(15);
2971
2972
2973 TCanvas *c= new TCanvas("c","c");
2974 fHistPtITSMI6InAcctestnew1->Draw();
2975 fHistPtITSMIge2InAcctestnew1->Draw("same");
2976 fHistPtITSMISPDInAcctestnew1->Draw("same");
2977 fHistPtITSMIoneSPDInAcctestnew1->Draw("same");
2978 fHistPtITSMI6InAcctestnew2->Draw("same");
2979 fHistPtITSMIge2InAcctestnew2->Draw("same");
2980 fHistPtITSMISPDInAcctestnew2->Draw("same");
2981 fHistPtITSMIoneSPDInAcctestnew2->Draw("same");
2982 fHistPtITSMI6InAcctestnew3->Draw("same");
2983 fHistPtITSMIge2InAcctestnew3->Draw("same");
2984 fHistPtITSMISPDInAcctestnew3->Draw("same");
2985 fHistPtITSMIoneSPDInAcctestnew3->Draw("same");
2986 fHistPtITSMI6InAcctestnew5->Draw("same");
2987 fHistPtITSMIge2InAcctestnew5->Draw("same");
2988 fHistPtITSMISPDInAcctestnew5->Draw("same");
2989 fHistPtITSMIoneSPDInAcctestnew5->Draw("same");
2990 fHistPtITSMI6InAcctestnew6->Draw("same");
2991 fHistPtITSMIge2InAcctestnew6->Draw("same");
2992 fHistPtITSMISPDInAcctestnew6->Draw("same");
2993 fHistPtITSMIoneSPDInAcctestnew6->Draw("same");
2994 fHistPtITSMI6InAcctestnew7->Draw("same");
2995 fHistPtITSMIge2InAcctestnew7->Draw("same");
2996 fHistPtITSMISPDInAcctestnew7->Draw("same");
2997 fHistPtITSMIoneSPDInAcctestnew7->Draw("same");
2998 fHistPtITSMI6InAcctestnew8->Draw("same");
2999 fHistPtITSMIge2InAcctestnew8->Draw("same");
3000 fHistPtITSMISPDInAcctestnew8->Draw("same");
3001 fHistPtITSMIoneSPDInAcctestnew8->Draw("same");
3002
3003
8fa56655 3004 return;
3005}
3006//--------------------------------------------------------------------------
3007void ReweightStrange(TH1F *hPt,TH1F* hPtPfromStrange,TH1F* hPtSfromStrange) {
3008
3009 for(Int_t i=1;i<=hPt->GetNbinsX();i++) {
3010 Double_t pt=hPt->GetBinCenter(i);
3011 Double_t weight=1.;
3012 Double_t weightP=1.;
3013 if(pt<.25) {
3014 weight=1.3;//1.;
3015 } else if(pt>=0.25 && pt<0.5) {
3016 weight=1.5;//1.3;
3017 } else if(pt>=0.5 && pt<1.1) {
3018 weight=1.7;//1.44;
3019 } else if(pt>=1.1) {
3020 weight=2;//1.7;
3021 }
3022
3023 weightP*=1.3;
3024 hPt->SetBinContent(i,hPt->GetBinContent(i)+(weight-1.)*hPtSfromStrange->GetBinContent(i));
3025 if(hPtPfromStrange) hPt->SetBinContent(i,hPt->GetBinContent(i)+(weightP-1.)*hPtPfromStrange->GetBinContent(i));
3026 }
3027
2308c814 3028 return;
3029}
b51e971f 3030//--------------------------------------------------------------------------
e75c9f34 3031void ITSTrackingTrending(Int_t firstrun=158508,Int_t lastrun=158511,
3032 //TString pathBeforeRun="/alice/sim/LHC11a10a_bis/",
3033 //TString pathAfterRun="/QA67/Stage_4/004/QAresults.root",
3034 TString pathBeforeRun="/alice/data/2011/LHC11d/000",
3035 TString pathAfterRun="/ESDs/pass1/QAresults.root",
71f9d353 3036 //TString pathBeforeRun="/alice/data/2010/LHC10b/000",
3037 //TString pathAfterRun="/ESDs/pass2/QA9/QAresults.root",
af874c53 3038 TString pathAfterRun2="")
3039{
b51e971f 3040 //
3041 // Make ITS efficiency trending plots from QAresults.root files
3042 //
3043 gStyle->SetOptStat(0);
3044
d388aa11 3045 TGrid::Connect("alien://");
3046
b51e971f 3047 Float_t ioValues[100],ioErrors[100];
3048 Int_t index;
3049 Int_t nruns=lastrun-firstrun+1;
3050
3051
3052 TH1F *hFracSPD1 = new TH1F("hFracSPD1","SPD inner; run number; Fraction of HSs",nruns,firstrun-0.5,lastrun+0.5);
3053 hFracSPD1->SetLineColor(3);
3054 hFracSPD1->SetMarkerColor(3);
e75c9f34 3055 hFracSPD1->SetMarkerStyle(20);
b51e971f 3056 TH1F *hFracSPD2 = new TH1F("hFracSPD2","SPD outer; run number; Fraction of HSs",nruns,firstrun-0.5,lastrun+0.5);
3057 hFracSPD2->SetLineColor(8);
3058 hFracSPD2->SetMarkerColor(8);
3059
3060
3061 TH1F *hEffge2Pt02 = new TH1F("hEffge2Pt02","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3062 hEffge2Pt02->SetLineWidth(2);
3063 hEffge2Pt02->SetLineColor(1);
3064 hEffge2Pt02->SetMarkerColor(1);
d388aa11 3065 hEffge2Pt02->SetMarkerStyle(20);
b51e971f 3066 TH1F *hEffge2Pt1 = new TH1F("hEffge2Pt1","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3067 hEffge2Pt1->SetLineWidth(2);
3068 hEffge2Pt1->SetLineColor(1);
3069 hEffge2Pt1->SetMarkerColor(1);
d388aa11 3070 hEffge2Pt1->SetMarkerStyle(20);
b51e971f 3071 TH1F *hEffge2Pt10 = new TH1F("hEffge2Pt10","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3072 hEffge2Pt10->SetLineWidth(2);
3073 hEffge2Pt10->SetLineColor(1);
3074 hEffge2Pt10->SetMarkerColor(1);
d388aa11 3075 hEffge2Pt10->SetMarkerStyle(20);
b51e971f 3076
3077 TH1F *hEff6Pt02 = new TH1F("hEff6Pt02","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3078 hEff6Pt02->SetLineWidth(2);
3079 hEff6Pt02->SetLineColor(2);
3080 hEff6Pt02->SetMarkerColor(2);
d388aa11 3081 hEff6Pt02->SetMarkerStyle(20);
b51e971f 3082 TH1F *hEff6Pt1 = new TH1F("hEff6Pt1","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3083 hEff6Pt1->SetLineWidth(2);
3084 hEff6Pt1->SetLineColor(2);
3085 hEff6Pt1->SetMarkerColor(2);
d388aa11 3086 hEff6Pt1->SetMarkerStyle(20);
b51e971f 3087 TH1F *hEff6Pt10 = new TH1F("hEff6Pt10","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3088 hEff6Pt10->SetLineWidth(2);
3089 hEff6Pt10->SetLineColor(2);
3090 hEff6Pt10->SetMarkerColor(2);
d388aa11 3091 hEff6Pt10->SetMarkerStyle(20);
b51e971f 3092
3093 TH1F *hEffSPDPt02 = new TH1F("hEffSPDPt02","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3094 hEffSPDPt02->SetLineWidth(2);
3095 hEffSPDPt02->SetLineColor(9);
3096 hEffSPDPt02->SetMarkerColor(9);
d388aa11 3097 hEffSPDPt02->SetMarkerStyle(20);
b51e971f 3098 TH1F *hEffSPDPt1 = new TH1F("hEffSPDPt1","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3099 hEffSPDPt1->SetLineWidth(2);
3100 hEffSPDPt1->SetLineColor(9);
3101 hEffSPDPt1->SetMarkerColor(9);
d388aa11 3102 hEffSPDPt1->SetMarkerStyle(20);
b51e971f 3103 TH1F *hEffSPDPt10 = new TH1F("hEffSPDPt10","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3104 hEffSPDPt10->SetLineWidth(2);
3105 hEffSPDPt10->SetLineColor(9);
3106 hEffSPDPt10->SetMarkerColor(9);
d388aa11 3107 hEffSPDPt10->SetMarkerStyle(20);
b51e971f 3108
3109 TH1F *hEffoneSPDPt02 = new TH1F("hEffoneSPDPt02","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3110 hEffoneSPDPt02->SetLineWidth(2);
3111 hEffoneSPDPt02->SetLineColor(15);
3112 hEffoneSPDPt02->SetMarkerColor(15);
d388aa11 3113 hEffoneSPDPt02->SetMarkerStyle(20);
b51e971f 3114 TH1F *hEffoneSPDPt1 = new TH1F("hEffoneSPDPt1","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3115 hEffoneSPDPt1->SetLineWidth(2);
3116 hEffoneSPDPt1->SetLineColor(15);
3117 hEffoneSPDPt1->SetMarkerColor(15);
d388aa11 3118 hEffoneSPDPt1->SetMarkerStyle(20);
b51e971f 3119 TH1F *hEffoneSPDPt10 = new TH1F("hEffoneSPDPt10","Efficiency; run number; TPC+ITS / TPC",nruns,firstrun-0.5,lastrun+0.5);
3120 hEffoneSPDPt10->SetLineWidth(2);
3121 hEffoneSPDPt10->SetLineColor(15);
3122 hEffoneSPDPt10->SetMarkerColor(15);
d388aa11 3123 hEffoneSPDPt10->SetMarkerStyle(20);
3124
b51e971f 3125 // loop on runs
3126 for(Int_t irun=firstrun; irun<=lastrun; irun++) {
3127 if(!KeepRun(irun)) continue;
3128
3129 TString path=pathBeforeRun;
3130 path+=irun;
3131 path.Append(pathAfterRun.Data());
3132 path.Prepend("alien://");
3133
e75c9f34 3134 printf("%s\n",path.Data());
3135
af874c53 3136 if(!PlotITSTrackingHists(path.Data(),ioValues,ioErrors)) {
3137 if(pathAfterRun2.Data()=="") continue;
3138 TString path2=pathBeforeRun;
3139 path2+=irun;
3140 path2.Append(pathAfterRun2.Data());
3141 path2.Prepend("alien://");
3142 if(!PlotITSTrackingHists(path2.Data(),ioValues,ioErrors)) continue;
3143 }
b51e971f 3144
3145 Int_t bin=hEffge2Pt1->FindBin(irun);
3146
3147
3148 // fill histos
3149 index=0;
3150 hFracSPD1->SetBinContent(bin,ioValues[index]);
3151 hFracSPD1->SetBinError(bin,.01);
3152 index=1;
3153 hFracSPD2->SetBinContent(bin,ioValues[index]);
3154 hFracSPD2->SetBinError(bin,.01);
3155
3156 index=2;
3157 hEffge2Pt02->SetBinContent(bin,ioValues[index]);
3158 hEffge2Pt02->SetBinError(bin,ioErrors[index]);
3159 index=3;
3160 hEffge2Pt1->SetBinContent(bin,ioValues[index]);
3161 hEffge2Pt1->SetBinError(bin,ioErrors[index]);
3162 index=4;
3163 hEffge2Pt10->SetBinContent(bin,ioValues[index]);
3164 hEffge2Pt10->SetBinError(bin,ioErrors[index]);
3165
3166 index=5;
3167 hEff6Pt02->SetBinContent(bin,ioValues[index]);
3168 hEff6Pt02->SetBinError(bin,ioErrors[index]);
3169 index=6;
3170 hEff6Pt1->SetBinContent(bin,ioValues[index]);
3171 hEff6Pt1->SetBinError(bin,ioErrors[index]);
3172 index=7;
3173 hEff6Pt10->SetBinContent(bin,ioValues[index]);
3174 hEff6Pt10->SetBinError(bin,ioErrors[index]);
3175
3176 index=8;
3177 hEffSPDPt02->SetBinContent(bin,ioValues[index]);
3178 hEffSPDPt02->SetBinError(bin,ioErrors[index]);
3179 index=9;
3180 hEffSPDPt1->SetBinContent(bin,ioValues[index]);
3181 hEffSPDPt1->SetBinError(bin,ioErrors[index]);
3182 index=10;
3183 hEffSPDPt10->SetBinContent(bin,ioValues[index]);
3184 hEffSPDPt10->SetBinError(bin,ioErrors[index]);
3185
3186 index=11;
3187 hEffoneSPDPt02->SetBinContent(bin,ioValues[index]);
3188 hEffoneSPDPt02->SetBinError(bin,ioErrors[index]);
3189 index=12;
3190 hEffoneSPDPt1->SetBinContent(bin,ioValues[index]);
3191 hEffoneSPDPt1->SetBinError(bin,ioErrors[index]);
3192 index=13;
3193 hEffoneSPDPt10->SetBinContent(bin,ioValues[index]);
3194 hEffoneSPDPt10->SetBinError(bin,ioErrors[index]);
3195
3196 }
3197
3198 TCanvas *cSPD = new TCanvas("cSPD","cSPD",0,0,1000,300);
3199 cSPD->SetGridy();
3200 hFracSPD1->SetMaximum(1.2);
3201 hFracSPD1->SetMaximum(0);
3202 hFracSPD1->Draw("p");
3203 hFracSPD2->Draw("same,p");
3204 TLegend* lSPD=new TLegend(0.9,0.8,1,1);
3205 lSPD->AddEntry(hFracSPD1,"Frac. SPD1 ON","l");
3206 lSPD->AddEntry(hFracSPD2,"Frac. SPD2 ON","l");
3207 lSPD->Draw();
3208
3209 TCanvas *cPt02 = new TCanvas("cPt02","cPt02",0,0,1000,300);
3210 cPt02->SetGridy();
3211 hEffge2Pt02->SetMaximum(1.2);
3212 hEffge2Pt02->SetMaximum(0);
3213 hEffge2Pt02->Draw();
3214 hEff6Pt02->Draw("same");
3215 hEffSPDPt02->Draw("same");
3216 hEffoneSPDPt02->Draw("same");
3217 TLegend* lPt02=new TLegend(0.9,0.8,1,1);
3218 lPt02->AddEntry(hEffge2Pt02,">=2 cls","l");
3219 lPt02->AddEntry(hEffoneSPDPt02,">=1 SPD + any","l");
3220 lPt02->AddEntry(hEffSPDPt02,"2 SPD + any","l");
3221 lPt02->AddEntry(hEff6Pt02,"6 cls","l");
3222 lPt02->Draw();
3223
3224 TCanvas *cPt1 = new TCanvas("cPt1","cPt1",0,0,1000,300);
3225 cPt1->SetGridy();
3226 hEffge2Pt1->SetMaximum(1.2);
3227 hEffge2Pt1->SetMaximum(0);
3228 hEffge2Pt1->Draw();
3229 hEff6Pt1->Draw("same");
3230 hEffSPDPt1->Draw("same");
3231 hEffoneSPDPt1->Draw("same");
3232 lPt02->Draw();
3233
3234 TCanvas *cPt10 = new TCanvas("cPt10","cPt10",0,0,1000,300);
3235 cPt10->SetGridy();
3236 hEffge2Pt10->SetMaximum(1.2);
3237 hEffge2Pt10->SetMaximum(0);
3238 hEffge2Pt10->Draw("p");
3239 hEff6Pt10->Draw("same,p");
3240 hEffSPDPt10->Draw("same,p");
3241 hEffoneSPDPt10->Draw("same,p");
3242 lPt02->Draw();
3243
3244 return;
3245}
3246//------------------------------------------------------------------------
3247Bool_t KeepRun(Int_t irun) {
3248
3249
71f9d353 3250
3251 // LHC10c good runs
3252 Int_t nruns10b=33;
3253 Int_t goodruns10b[33]={117222, 117220, 117116, 117112, 117109, 117099, 117092, 117086, 117077, 117065, 117063, 117060, 117059, 117054, 117053, 117052, 117050, 117048, 116645, 116643, 116574, 116571, 116562, 116403, 116402, 116288, 116102, 115414, 115401, 115393, 115193, 115186, 114931};
3254
b51e971f 3255 // LHC10c good runs
af874c53 3256 Int_t nruns10c=46;
3257 Int_t goodruns10c[46]={121040, 121039, 120829, 120825, 120824, 120823, 120822, 120821, 120820, 120758, 120750, 120741, 120671, 120617, 120616, 120505, 120504, 120503, 120244, 120079, 120076, 120073, 120072, 120069, 120067, 119862, 119859, 119856, 119853, 119849, 119846, 119845, 119844, 119842, 119841, 119163, 119161, 119159, 118561, 118560, 118558, 118556, 118518, 118512, 118507, 118506};
3258 //
3259
3260 // LHC10d good runs
3261 Int_t nruns10d=80;
3262 Int_t goodruns10d[80]={126437, 126432, 126425, 126424, 126422, 126409, 126408, 126407, 126406, 126405, 126404, 126403, 126359, 126352, 126351, 126350, 126285, 126284, 126283, 126168, 126167, 126160, 126158, 126097, 126090, 126088, 126082, 126081, 126078, 126073, 126008, 126007, 126004, 125855, 125851, 125850, 125849, 125848, 125847, 125844, 125843, 125842, 125633, 125632, 125630, 125628, 125296, 125295, 125186, 125156, 125140, 125139, 125134, 125133, 125101, 125100, 125097, 125085, 125083, 125023, 124751, 124750, 124746, 124702, 124608, 124607, 124606, 124605, 124604, 124381, 124380, 124378, 124367, 124362, 124358, 124355, 124191, 124187, 122375, 122374};
b51e971f 3263 //
3264
af874c53 3265 // LHC10e good runs
e75c9f34 3266 Int_t nruns10e=158;
3267 Int_t goodruns10e[158]={130369, 130365, 130360, 130358, 130356, 130354, 130353, 130348, 130343, 130342, 130179, 130178, 130172, 130170, 130168, 130158, 130157, 130156, 130151, 130149, 130148, 129983, 129966, 129962, 129961, 129960, 129959, 129763, 129760, 129750, 129748, 129747, 129745, 129744, 129742, 129738, 129736, 129735, 129734, 129731, 129729, 129726, 129725, 129723, 129667, 129666, 129659, 129655, 129654, 129653, 129652, 129651, 129650, 129649, 129648, 129647, 129641, 129639, 129599, 129598, 129597, 129587, 129586, 129541, 129540, 129536, 129529, 129528, 129527, 129526, 129525, 129524, 129523, 129522, 129521, 129520, 129519, 129516, 129515, 129514, 129513, 129512, 129510, 129508, 129042, 129041, 128913, 128912, 128911, 128910, 128855, 128853, 128850, 128849, 128843, 128836, 128835, 128833, 128824, 128823, 128820, 128819, 128817, 128813, 128777, 128776, 128678, 128677, 128621, 128615, 128611, 128609, 128605, 128596, 128594, 128592, 128590, 128582, 128581, 128507, 128506, 128505, 128504, 128503, 128498, 128495, 128494, 128486, 128483, 128452, 128366, 128260, 128257, 128192, 128191, 128189, 128186, 128185, 128182, 128180, 128175, 128053, 127942, 127941, 127937, 127936, 127935, 127933, 127932, 127931, 127822, 127817, 127815, 127814, 127730, 127729, 127724, 127719};
3268
d388aa11 3269
3270
3271 // LHC10h good runs
e75c9f34 3272 Int_t nruns10h=49;
3273 Int_t goodruns10h[49]={138624, 138637, 138666, 138730, 138731, 138732, 138737, 138740, 138742, 138830, 138836, 138871, 138872, 138924, 138965, 138972, 138977, 138978, 138980, 138982, 138983, 139029, 139036, 139037, 139038, 139042, 139105, 139107, 139309, 139310, 139314, 139360, 139438, 139439, 139440, 139441, 139465, 139466, 139467, 139470, 139471, 139503, 139505, 139507, 139510, 139511, 139513, 139514, 139517};
d388aa11 3274
e75c9f34 3275 /*
3276 Int_t nruns10h=91;
3277 Int_t goodruns10h[91]={139510, 139507, 139505, 139503, 139465, 139438, 139437, 139360, 139329, 139328, 139314, 139310, 139309, 139173, 139107, 139105, 139038, 139037, 139036, 139029, 139028, 138872, 138871, 138870, 138837, 138732, 138730, 138666, 138662, 138653, 138652, 138638, 138624, 138621, 138583, 138582, 138579, 138578, 138534, 138469, 138442, 138439, 138438, 138396, 138364, 138275, 138225, 138201, 138197, 138192, 138190, 137848, 137844, 137752, 137751, 137724, 137722, 137718, 137704, 137693, 137692, 137691, 137686, 137685, 137639, 137638, 137608, 137595, 137549, 137546, 137544, 137541, 137539, 137531, 137530, 137443, 137441, 137440, 137439, 137434, 137432, 137431, 137430, 137366, 137243, 137236, 137235, 137232, 137231, 137162, 137161};
3278 */
d388aa11 3279
e75c9f34 3280 // LHC11b good runs
3281 Int_t nruns11b=83;
3282 Int_t goodruns11b[83]={148531, 148541, 148544, 148648, 148800, 148847, 148857, 149102, 149115, 149122, 149123, 149367, 149382, 149389, 149395, 149398, 149432, 149434, 149442, 149449, 149452, 149453, 149456, 149461, 149462, 149472, 149484, 149487, 149499, 149529, 149530, 149531, 149533, 149534, 149540, 149543, 149544, 149545, 149547, 149549, 149582, 149584, 149585, 149588, 149590, 149592, 149594, 149596, 149598, 149600, 149604, 149657, 149669, 149680, 149682, 149712, 149723, 149726, 149727, 149734, 149735, 149740, 149746, 149749, 149757, 149761, 149769, 149771, 149772, 149773, 149774, 149775, 149776, 149781, 149785, 149786, 149789, 149790, 150375, 150421, 150586, 150589, 150705};
d388aa11 3283
e75c9f34 3284 // LHC11c good runs
3285 Int_t nruns11c=244;
3286 Int_t goodruns11c[244]={154570, 154495, 154485, 154483, 154482, 154480, 154478, 154448, 154385, 154383, 154382, 154296, 154293, 154289, 154286, 154283, 154281, 154273, 154270, 154269, 154266, 154264, 154261, 154257, 154252, 154222, 154221, 154220, 154219, 154211, 154207, 154163, 154158, 154151, 154145, 154143, 154141, 154138, 154136, 154132, 154131, 154130, 154129, 154126, 154125, 154091, 154083, 154081, 154039, 154031, 154030, 154026, 154024, 154018, 153954, 153946, 153944, 153939, 153938, 153935, 153929, 153924, 153916, 153911, 153909, 153906, 153876, 153875, 153873, 153812, 153808, 153807, 153805, 153798, 153796, 153794, 153784, 153781, 153779, 153778, 153777, 153776, 153738, 153733, 153728, 153727, 153726, 153725, 153718, 153709, 153702, 153594, 153591, 153589, 153588, 153587, 153583, 153578, 153571, 153570, 153566, 153560, 153558, 153552, 153548, 153544, 153542, 153541, 153539, 153536, 153533, 153513, 153465, 153415, 153413, 153373, 153371, 153369, 153363, 153362, 153360, 153296, 153234, 153232, 153227, 153223, 153127, 153121, 153120, 153117, 153116, 153115, 153059, 153056, 152935, 152934, 152907, 152823, 152822, 152821, 152820, 152819, 152817, 152780, 152773, 152772, 152751, 152750, 152718, 152717, 152716, 152715, 152708, 152702, 152701, 152698, 152697, 152696, 152695, 152658, 152581, 152513, 152512, 152488, 152455, 152377, 152371, 152369, 152368, 152367, 152334, 152332, 152323, 152322, 152321, 152320, 152319, 152318, 152314, 152313, 152312, 152311, 152310, 152309, 152307, 152306, 152304, 152285, 152284, 152257, 152256, 152214, 152209, 152208, 152207, 152206, 152146, 152138, 152137, 152136, 152091, 152090, 152089, 152087, 152086, 152083, 152082, 152081, 152079, 152078, 152046, 152015, 152011, 152008, 152007, 152005, 152003, 152002, 151852, 151851, 151850, 151849, 151810, 151809, 151752, 151751, 151732, 151724, 151689, 151681, 151680, 151678, 151674, 151672, 151671, 151669, 151666, 151665, 151664, 151661, 151660, 151655, 151638, 151636};
d388aa11 3287
e75c9f34 3288 // LHC11d good runs
3289 Int_t nruns11d=81;
3290 Int_t goodruns11d[81]={158508, 158509, 158511, 158518, 158520, 158521, 158526, 158528, 158531, 158592, 158602, 158604, 158611, 158615, 158622, 158626, 158706, 158717, 158722, 158745, 158777, 158780, 158784, 158788, 158791, 158793, 158794, 158848, 158868, 158876, 158877, 158879, 159028, 159042, 159076, 159090, 159120, 159121, 159128, 159147, 159154, 159162, 159168, 159173, 159185, 159191, 159193, 159199, 159201, 159205, 159207, 159214, 159216, 159218, 159221, 159223, 159255, 159258, 159260, 159285, 159286, 159318, 159379, 159451, 159503, 159517, 159521, 159535, 159538, 159571, 159577, 159580, 159581, 159582, 159586, 159593, 159595, 159599, 159602, 159606, 159635};
af874c53 3291
b51e971f 3292
3293 Bool_t found=kFALSE;
af874c53 3294 Int_t i=0;
3295
71f9d353 3296 for(i=0; i<nruns10b; i++) {
3297 if(irun==goodruns10b[i]) {found=kTRUE; break;}
3298 }
3299 if(found) return kTRUE;
3300
af874c53 3301 for(i=0; i<nruns10c; i++) {
3302 if(irun==goodruns10c[i]) {found=kTRUE; break;}
3303 }
3304 if(found) return kTRUE;
3305
3306 for(i=0; i<nruns10d; i++) {
3307 if(irun==goodruns10d[i]) {found=kTRUE; break;}
3308 }
b51e971f 3309 if(found) return kTRUE;
af874c53 3310
3311 for(i=0; i<nruns10e; i++) {
3312 if(irun==goodruns10e[i]) {found=kTRUE; break;}
3313 }
3314 if(found) return kTRUE;
3315
d388aa11 3316 for(i=0; i<nruns10h; i++) {
3317 if(irun==goodruns10h[i]) {found=kTRUE; break;}
3318 }
3319 if(found) return kTRUE;
3320
e75c9f34 3321 for(i=0; i<nruns11b; i++) {
3322 if(irun==goodruns11b[i]) {found=kTRUE; break;}
3323 }
3324 if(found) return kTRUE;
3325
3326 for(i=0; i<nruns11c; i++) {
3327 if(irun==goodruns11c[i]) {found=kTRUE; break;}
3328 }
3329 if(found) return kTRUE;
3330
3331 for(i=0; i<nruns11d; i++) {
3332 if(irun==goodruns11d[i]) {found=kTRUE; break;}
d388aa11 3333 }
3334 if(found) return kTRUE;
3335
3336
b51e971f 3337 return kFALSE;
3338}
af874c53 3339//-------------------------------------------------------------------------
3340void FakesWithChi2Cut(TString fname) {
3341
3342 gStyle->SetOptStat(0);
3343
3344 TFile *f= TFile::Open(fname.Data());
3345 if(!f) return kFALSE;
3346
3347 TList *list=(TList*)f->Get("cOutputITS");
3348 TDirectoryFile *dir=0;
3349 if(!list) {
3350 dir=(TDirectoryFile*)f->GetDirectory("ITS_Performance");
3351 if(dir) list = (TList*)dir->Get("cOutputITS");
3352 }
3353
3354 if(!list) return kFALSE;
3355
3356
3357
3358
3359 TH1F *fHistRedChi2AllPt02 = (TH1F*)list->FindObject("fHistITSRedChi2NonFakePt02");
3360 TH1F *fHistRedChi2FakesPt02 = (TH1F*)list->FindObject("fHistITSRedChi2FakePt02");
3361 TH1F *fHistRedChi2NonFakesPt02 = (TH1F*)fHistRedChi2AllPt02->Clone("fHistITSRedChi2NonFakePt02");
3362 fHistRedChi2NonFakesPt02->Add(fHistRedChi2FakesPt02,-1.);
3363
3364 TH1F *fHistFakeFracVSChi2CutPt02 = (TH1F*)fHistRedChi2AllPt02->Clone("fHistFakeFracVSChi2CutPt02");
3365 fHistFakeFracVSChi2CutPt02->SetLineColor(2);
3366 fHistFakeFracVSChi2CutPt02->SetLineWidth(2);
3367 TH1F *fHistNonFakeEffVSChi2CutPt02 = (TH1F*)fHistRedChi2AllPt02->Clone("fHistFakeEffVSChi2CutPt02");
3368 fHistNonFakeEffVSChi2CutPt02->SetLineColor(4);
3369 fHistNonFakeEffVSChi2CutPt02->SetLineWidth(2);
3370
3371 for(Int_t bin=1;bin<=fHistFakeFracVSChi2CutPt02->GetNbinsX();bin++) {
3372 Float_t fakesBelowCut=fHistRedChi2FakesPt02->Integral(1,bin);
3373 Float_t nonfakesBelowCut=fHistRedChi2NonFakesPt02->Integral(1,bin);
3374 Float_t allBelowCut=fHistRedChi2AllPt02->Integral(1,bin);
3375 Float_t nonfakesNoCut=fHistRedChi2NonFakesPt02->Integral();
3376 fHistFakeFracVSChi2CutPt02->SetBinContent(bin,fakesBelowCut/allBelowCut);
3377 fHistNonFakeEffVSChi2CutPt02->SetBinContent(bin,nonfakesBelowCut/nonfakesNoCut);
3378 }
3379
3380 TH1F *fHistRedChi2AllPt05 = (TH1F*)list->FindObject("fHistITSRedChi2NonFakePt05");
3381 TH1F *fHistRedChi2FakesPt05 = (TH1F*)list->FindObject("fHistITSRedChi2FakePt05");
3382 TH1F *fHistRedChi2NonFakesPt05 = (TH1F*)fHistRedChi2AllPt05->Clone("fHistITSRedChi2NonFakePt05");
3383 fHistRedChi2NonFakesPt05->Add(fHistRedChi2FakesPt05,-1.);
3384
3385 TH1F *fHistFakeFracVSChi2CutPt05 = (TH1F*)fHistRedChi2AllPt05->Clone("fHistFakeFracVSChi2CutPt05");
3386 fHistFakeFracVSChi2CutPt05->SetLineColor(2);
3387 fHistFakeFracVSChi2CutPt05->SetLineWidth(2);
3388 TH1F *fHistNonFakeEffVSChi2CutPt05 = (TH1F*)fHistRedChi2AllPt05->Clone("fHistFakeEffVSChi2CutPt05");
3389 fHistNonFakeEffVSChi2CutPt05->SetLineColor(4);
3390 fHistNonFakeEffVSChi2CutPt05->SetLineWidth(2);
3391
3392 for(Int_t bin=1;bin<=fHistFakeFracVSChi2CutPt05->GetNbinsX();bin++) {
3393 Float_t fakesBelowCut=fHistRedChi2FakesPt05->Integral(1,bin);
3394 Float_t nonfakesBelowCut=fHistRedChi2NonFakesPt05->Integral(1,bin);
3395 Float_t allBelowCut=fHistRedChi2AllPt05->Integral(1,bin);
3396 Float_t nonfakesNoCut=fHistRedChi2NonFakesPt05->Integral();
3397 fHistFakeFracVSChi2CutPt05->SetBinContent(bin,fakesBelowCut/allBelowCut);
3398 fHistNonFakeEffVSChi2CutPt05->SetBinContent(bin,nonfakesBelowCut/nonfakesNoCut);
3399 }
3400
3401 TH1F *fHistRedChi2AllPt1 = (TH1F*)list->FindObject("fHistITSRedChi2NonFakePt1");
3402 TH1F *fHistRedChi2FakesPt1 = (TH1F*)list->FindObject("fHistITSRedChi2FakePt1");
3403 TH1F *fHistRedChi2NonFakesPt1 = (TH1F*)fHistRedChi2AllPt1->Clone("fHistITSRedChi2NonFakePt1");
3404 fHistRedChi2NonFakesPt1->Add(fHistRedChi2FakesPt1,-1.);
3405
3406 TH1F *fHistFakeFracVSChi2CutPt1 = (TH1F*)fHistRedChi2AllPt1->Clone("fHistFakeFracVSChi2CutPt1");
3407 fHistFakeFracVSChi2CutPt1->SetLineColor(2);
3408 fHistFakeFracVSChi2CutPt1->SetLineWidth(2);
3409 TH1F *fHistNonFakeEffVSChi2CutPt1 = (TH1F*)fHistRedChi2AllPt1->Clone("fHistFakeEffVSChi2CutPt1");
3410 fHistNonFakeEffVSChi2CutPt1->SetLineColor(4);
3411 fHistNonFakeEffVSChi2CutPt1->SetLineWidth(2);
3412
3413 for(Int_t bin=1;bin<=fHistFakeFracVSChi2CutPt1->GetNbinsX();bin++) {
3414 Float_t fakesBelowCut=fHistRedChi2FakesPt1->Integral(1,bin);
3415 Float_t nonfakesBelowCut=fHistRedChi2NonFakesPt1->Integral(1,bin);
3416 Float_t allBelowCut=fHistRedChi2AllPt1->Integral(1,bin);
3417 Float_t nonfakesNoCut=fHistRedChi2NonFakesPt1->Integral();
3418 fHistFakeFracVSChi2CutPt1->SetBinContent(bin,fakesBelowCut/allBelowCut);
3419 fHistNonFakeEffVSChi2CutPt1->SetBinContent(bin,nonfakesBelowCut/nonfakesNoCut);
3420 }
3421
3422
3423 TCanvas *c=new TCanvas("c","c",0,0,1000,500);
3424 c->Divide(3,1);
3425 c->cd(1);
3426 fHistRedChi2NonFakesPt02->Scale(1./fHistRedChi2NonFakesPt02->Integral());
3427 fHistRedChi2NonFakesPt02->SetLineColor(4);
3428 fHistRedChi2NonFakesPt02->Draw();
3429 fHistRedChi2FakesPt02->Scale(1./fHistRedChi2FakesPt02->Integral());
3430 fHistRedChi2FakesPt02->SetLineColor(2);
3431 fHistRedChi2FakesPt02->Draw("same");
3432 c->cd(2);
3433 fHistRedChi2NonFakesPt05->Scale(1./fHistRedChi2NonFakesPt05->Integral());
3434 fHistRedChi2NonFakesPt05->SetLineColor(4);
3435 fHistRedChi2NonFakesPt05->Draw();
3436 fHistRedChi2FakesPt05->Scale(1./fHistRedChi2FakesPt05->Integral());
3437 fHistRedChi2FakesPt05->SetLineColor(2);
3438 fHistRedChi2FakesPt05->Draw("same");
3439 c->cd(3);
3440 fHistRedChi2NonFakesPt1->Scale(1./fHistRedChi2NonFakesPt1->Integral());
3441 fHistRedChi2NonFakesPt1->SetLineColor(4);
3442 fHistRedChi2NonFakesPt1->Draw();
3443 fHistRedChi2FakesPt1->Scale(1./fHistRedChi2FakesPt1->Integral());
3444 fHistRedChi2FakesPt1->SetLineColor(2);
3445 fHistRedChi2FakesPt1->Draw("same");
3446
3447 TCanvas *cc=new TCanvas("cc","cc",0,0,1000,500);
3448 cc->Divide(3,1);
3449 cc_1->SetLogx();
3450 cc_2->SetLogx();
3451 cc_3->SetLogx();
3452 cc->cd(1);
3453 fHistFakeFracVSChi2CutPt02->SetXTitle("maximum ITS #chi^{2}/nclusters");
3454 fHistFakeFracVSChi2CutPt02->SetMinimum(0);
3455 fHistFakeFracVSChi2CutPt02->SetMaximum(1);
3456 fHistFakeFracVSChi2CutPt02->Draw();
3457 fHistNonFakeEffVSChi2CutPt02->Draw("same");
3458 cc->cd(2);
3459 fHistFakeFracVSChi2CutPt05->SetXTitle("maximum ITS #chi^{2}/nclusters");
3460 fHistFakeFracVSChi2CutPt05->SetMinimum(0);
3461 fHistFakeFracVSChi2CutPt05->SetMaximum(1);
3462 fHistFakeFracVSChi2CutPt05->Draw();
3463 fHistNonFakeEffVSChi2CutPt05->Draw("same");
3464 cc->cd(3);
3465 fHistFakeFracVSChi2CutPt1->SetXTitle("maximum ITS #chi^{2}/nclusters");
3466 fHistFakeFracVSChi2CutPt1->SetMinimum(0);
3467 fHistFakeFracVSChi2CutPt1->SetMaximum(1);
3468 fHistFakeFracVSChi2CutPt1->Draw();
3469 fHistNonFakeEffVSChi2CutPt1->Draw("same");
3470
3471
3472
3473 return;
3474}