]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG1/macros/PlotITSTrackingHists.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Added possibility to recalibrate CaloClusters during the analysis
[u/mrichter/AliRoot.git] / PWG1 / macros / PlotITSTrackingHists.C
CommitLineData
c2893e82 1void PlotITSTrackingHists(TString fname="ITS.Performance.root") {
2308c814 2 //
3 // Macro to plot the histos from the task AliAnalysisTaskITSTrackingCheck
4 // A. Dainese 28.11.09
5 //
6
fa516524 7 //gStyle->SetOptStat(0);
2308c814 8
c2893e82 9 TFile *f= new TFile(fname.Data());
2308c814 10
fa516524 11 TList *list=(TList*)f->Get("cOutputITS");
2308c814 12 TH1F *fHistNclsITSSA = (TH1F*)list->FindObject("fHistNclsITSSA");
13 TH1F *fHistClusterMapITSSA = (TH1F*)list->FindObject("fHistClusterMapITSSA");
14 TH1F *fHistClusterMapITSSAok = (TH1F*)list->FindObject("fHistClusterMapITSSAok");
15 TH1F *fHistClusterMapITSSAbad = (TH1F*)list->FindObject("fHistClusterMapITSSAbad");
16 TH1F *fHistClusterMapITSSAskipped = (TH1F*)list->FindObject("fHistClusterMapITSSAskipped");
17 TH1F *fHistClusterMapITSSAoutinz = (TH1F*)list->FindObject("fHistClusterMapITSSAoutinz");
18 TH1F *fHistClusterMapITSSAokoutinzbad = (TH1F*)list->FindObject("fHistClusterMapITSSAokoutinzbad");
19 TH1F *fHistClusterMapITSSAnorefit = (TH1F*)list->FindObject("fHistClusterMapITSSAnorefit");
20 TH1F *fHistClusterMapITSSAnocls = (TH1F*)list->FindObject("fHistClusterMapITSSAnocls");
21 TH1F *fHistNclsITSSAInAcc = (TH1F*)list->FindObject("fHistNclsITSSAInAcc");
22 TH1F *fHistClusterMapITSSAInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAInAcc");
23 TH1F *fHistClusterMapITSSAokInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAokInAcc");
24 TH1F *fHistClusterMapITSSAbadInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAbadInAcc");
fa516524 25 TH1F *fHistClusterMapModuleITSMIokInAcc = (TH1F*)list->FindObject("fHistClusterMapModuleITSMIokInAcc");
26 TH1F *fHistClusterMapModuleITSMIbadInAcc = (TH1F*)list->FindObject("fHistClusterMapModuleITSMIbadInAcc");
2308c814 27 TH1F *fHistClusterMapITSSAskippedInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAskippedInAcc");
28 TH1F *fHistClusterMapITSSAoutinzInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAoutinzInAcc");
29 TH1F *fHistClusterMapITSSAokoutinzbadInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAokoutinzbadInAcc");
30 TH1F *fHistClusterMapITSSAnorefitInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAnorefitInAcc");
31 TH1F *fHistClusterMapITSSAnoclsInAcc = (TH1F*)list->FindObject("fHistClusterMapITSSAnoclsInAcc");
fa516524 32 TH1F *fHistClusterMapModuleITSMInoclsInAcc = (TH1F*)list->FindObject("fHistClusterMapModuleITSMInoclsInAcc");
2308c814 33 TH1F *fHistNclsITSMI = (TH1F*)list->FindObject("fHistNclsITSMI");
34 TH1F *fHistClusterMapITSMI = (TH1F*)list->FindObject("fHistClusterMapITSMI");
35 TH1F *fHistClusterMapITSMIok = (TH1F*)list->FindObject("fHistClusterMapITSMIok");
36 TH1F *fHistClusterMapITSMIbad = (TH1F*)list->FindObject("fHistClusterMapITSMIbad");
37 TH1F *fHistClusterMapITSMIskipped = (TH1F*)list->FindObject("fHistClusterMapITSMIskipped");
38 TH1F *fHistClusterMapITSMIoutinz = (TH1F*)list->FindObject("fHistClusterMapITSMIoutinz");
39 TH1F *fHistClusterMapITSMIokoutinzbad = (TH1F*)list->FindObject("fHistClusterMapITSMIokoutinzbad");
40 TH1F *fHistClusterMapITSMInorefit = (TH1F*)list->FindObject("fHistClusterMapITSMInorefit");
41 TH1F *fHistClusterMapITSMInocls = (TH1F*)list->FindObject("fHistClusterMapITSMInocls");
42
c2893e82 43 TH1F *fHistPhiTPCInAcc = (TH1F*)list->FindObject("fHistPhiTPCInAcc");
44 TH1F *fHistPhiITSMIokbadoutinz6InAcc = (TH1F*)list->FindObject("fHistPhiITSMIokbadoutinz6InAcc");
45
46 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
fa516524 47 TH1F *fHistPtTPCInAccMCtwoSPD = (TH1F*)list->FindObject("fHistPtTPCInAccMCtwoSPD");
48 TH1F *fHistPtTPCInAccMConeSPD = (TH1F*)list->FindObject("fHistPtTPCInAccMConeSPD");
c2893e82 49 TH1F *fHistPtITSMI6InAcc = (TH1F*)list->FindObject("fHistPtITSMI6InAcc");
50 TH1F *fHistPtITSMI5InAcc = (TH1F*)list->FindObject("fHistPtITSMI5InAcc");
51 TH1F *fHistPtITSMI4InAcc = (TH1F*)list->FindObject("fHistPtITSMI4InAcc");
52 TH1F *fHistPtITSMI3InAcc = (TH1F*)list->FindObject("fHistPtITSMI3InAcc");
53 TH1F *fHistPtITSMI2InAcc = (TH1F*)list->FindObject("fHistPtITSMI2InAcc");
54 TH1F *fHistPtITSMISPDInAcc = (TH1F*)list->FindObject("fHistPtITSMISPDInAcc");
fa516524 55 TH1F *fHistPtITSMIoneSPDInAcc = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAcc");
c2893e82 56 TH1F *fHistPtITSMIokbadoutinz6InAcc = (TH1F*)list->FindObject("fHistPtITSMIokbadoutinz6InAcc");
57 TH1F *fHistPtITSMIokbadoutinz5InAcc = (TH1F*)list->FindObject("fHistPtITSMIokbadoutinz5InAcc");
58 TH1F *fHistPtITSMIokbadoutinz4InAcc = (TH1F*)list->FindObject("fHistPtITSMIokbadoutinz4InAcc");
fa516524 59 TH1F *fHistPtTPCInAccP = (TH1F*)list->FindObject("fHistPtTPCInAccP");
60 TH1F *fHistPtITSMI6InAccP = (TH1F*)list->FindObject("fHistPtITSMI6InAccP");
61 TH1F *fHistPtITSMI5InAccP = (TH1F*)list->FindObject("fHistPtITSMI5InAccP");
62 TH1F *fHistPtITSMI4InAccP = (TH1F*)list->FindObject("fHistPtITSMI4InAccP");
63 TH1F *fHistPtITSMI3InAccP = (TH1F*)list->FindObject("fHistPtITSMI3InAccP");
64 TH1F *fHistPtITSMI2InAccP = (TH1F*)list->FindObject("fHistPtITSMI2InAccP");
65 TH1F *fHistPtITSMISPDInAccP = (TH1F*)list->FindObject("fHistPtITSMISPDInAccP");
66 TH1F *fHistPtITSMIoneSPDInAccP = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAccP");
67 TH1F *fHistPtTPCInAccS = (TH1F*)list->FindObject("fHistPtTPCInAccS");
68 TH1F *fHistPtITSMI6InAccS = (TH1F*)list->FindObject("fHistPtITSMI6InAccS");
69 TH1F *fHistPtITSMI5InAccS = (TH1F*)list->FindObject("fHistPtITSMI5InAccS");
70 TH1F *fHistPtITSMI4InAccS = (TH1F*)list->FindObject("fHistPtITSMI4InAccS");
71 TH1F *fHistPtITSMI3InAccS = (TH1F*)list->FindObject("fHistPtITSMI3InAccS");
72 TH1F *fHistPtITSMI2InAccS = (TH1F*)list->FindObject("fHistPtITSMI2InAccS");
73 TH1F *fHistPtITSMISPDInAccS = (TH1F*)list->FindObject("fHistPtITSMISPDInAccS");
74 TH1F *fHistPtITSMIoneSPDInAccS = (TH1F*)list->FindObject("fHistPtITSMIoneSPDInAccS");
75
76 TH1F *fHistRProdVtxInAccP = (TH1F*)list->FindObject("fHistRProdVtxInAccP");
77 TH1F *fHistRProdVtxInAccS = (TH1F*)list->FindObject("fHistRProdVtxInAccS");
78
79 TH1F *fHistxlocSDDok = (TH1F*)list->FindObject("fHistxlocSDDok");
80 TH1F *fHistxlocSDDall = (TH1F*)list->FindObject("fHistxlocSDDall");
81 TH1F *fHistzlocSDDok = (TH1F*)list->FindObject("fHistzlocSDDok");
82 TH1F *fHistzlocSDDall = (TH1F*)list->FindObject("fHistzlocSDDall");
83
84 TH1F *fHistPtITSTPCsel = (TH1F*)list->FindObject("fHistPtITSTPCsel");
85 TH1F *fHistPtITSTPCselP = (TH1F*)list->FindObject("fHistPtITSTPCselP");
86 TH1F *fHistPtITSTPCselS = (TH1F*)list->FindObject("fHistPtITSTPCselS");
87 TH1F *fHistPtITSTPCselPfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
88 TH1F *fHistPtITSTPCselSfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
89 TH1F *fHistPtITSTPCselSfromMat = (TH1F*)list->FindObject("fHistPtITSTPCselSfromMat");
90
91 TH1F *fHistPtTPCInAccPfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccSfromStrange");
92 TH1F *fHistPtTPCInAccSfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccSfromStrange");
93 TH1F *fHistPtTPCInAccSfromMat = (TH1F*)list->FindObject("fHistPtTPCInAccSfromMat");
94
95
c2893e82 96
97 //---------------------------------------------------------------
fa516524 98 TH1F *fHistPtITSMIge2InAcc = (TH1F*)fHistPtITSMI6InAcc->Clone("fHistPtITSMIge2InAcc");
99 fHistPtITSMIge2InAcc->Add(fHistPtITSMI5InAcc);
100 fHistPtITSMIge2InAcc->Add(fHistPtITSMI4InAcc);
101 fHistPtITSMIge2InAcc->Add(fHistPtITSMI3InAcc);
102 fHistPtITSMIge2InAcc->Add(fHistPtITSMI2InAcc);
103 TH1F *fHistPtITSMIge2InAccP = (TH1F*)fHistPtITSMI6InAccP->Clone("fHistPtITSMIge2InAccP");
104 fHistPtITSMIge2InAccP->Add(fHistPtITSMI5InAccP);
105 fHistPtITSMIge2InAccP->Add(fHistPtITSMI4InAccP);
106 fHistPtITSMIge2InAccP->Add(fHistPtITSMI3InAccP);
107 fHistPtITSMIge2InAccP->Add(fHistPtITSMI2InAccP);
108 TH1F *fHistPtITSMIge2InAccS = (TH1F*)fHistPtITSMI6InAccS->Clone("fHistPtITSMIge2InAccS");
109 fHistPtITSMIge2InAccS->Add(fHistPtITSMI5InAccS);
110 fHistPtITSMIge2InAccS->Add(fHistPtITSMI4InAccS);
111 fHistPtITSMIge2InAccS->Add(fHistPtITSMI3InAccS);
112 fHistPtITSMIge2InAccS->Add(fHistPtITSMI2InAccS);
113
114 TH1F* fHistPtITSMISPDInAccMC=(TH1F*)fHistPtITSMISPDInAcc->Clone("fHistPtITSMISPDInAccMC");
115 TH1F* fHistPtITSMIoneSPDInAccMC=(TH1F*)fHistPtITSMIoneSPDInAcc->Clone("fHistPtITSMIoneSPDInAccMC");
116
c2893e82 117
2308c814 118 TLegend *l1=new TLegend(0.5,0.5,0.9,0.9);
119 TLegend *l2=new TLegend(0.5,0.5,0.9,0.9);
120
121 TCanvas *c1= new TCanvas("c1","c1",10,10,600,500);
122 fHistNclsITSSA->SetMinimum(0);
123 fHistNclsITSSA->SetLineColor(1);
124 l1->AddEntry(fHistNclsITSSA,"ITS-SA","l");
125 fHistNclsITSSA->Draw();
126 fHistNclsITSSAInAcc->SetLineColor(4);
127 l1->AddEntry(fHistNclsITSSAInAcc,"ITS-SA in acc.","l");
128 fHistNclsITSSAInAcc->Draw("same");
129 fHistNclsITSMI->SetLineColor(2);
130 l1->AddEntry(fHistNclsITSMI,"ITS from TPC","l");
131 fHistNclsITSMI->Draw("same");
132 l1->Draw();
133
134
135 TCanvas *c2 =new TCanvas("c2","c2",10,10,1200,800);
136 c2->Divide(3,2);
137 c2->cd(1);
138 //
139 fHistClusterMapITSSAokoutinzbad->SetLineColor(1);
140 fHistClusterMapITSSAokoutinzbad->SetMarkerColor(1);
141 fHistClusterMapITSSAokoutinzbad->SetMarkerStyle(20);
142 fHistClusterMapITSSAokoutinzbad->Draw();
c2893e82 143 fHistClusterMapITSMIokoutinzbad->SetLineColor(2);
144 fHistClusterMapITSMIokoutinzbad->SetMarkerColor(2);
2308c814 145 fHistClusterMapITSMIokoutinzbad->SetMarkerStyle(20);
146 fHistClusterMapITSMIokoutinzbad->Draw("same");
147 l1->Draw();
148 //
149 c2->cd(2);
150 fHistClusterMapITSSAok->SetLineColor(1);
151 fHistClusterMapITSSAok->SetMarkerColor(1);
152 fHistClusterMapITSSAok->SetMarkerStyle(21);
153 fHistClusterMapITSSAok->Draw();
154 fHistClusterMapITSMIok->SetLineColor(2);
155 fHistClusterMapITSMIok->SetMarkerColor(2);
156 fHistClusterMapITSMIok->SetMarkerStyle(21);
157 fHistClusterMapITSMIok->Draw("same");
158 //
159 c2->cd(3);
160 fHistClusterMapITSSAoutinz->SetLineColor(1);
161 fHistClusterMapITSSAoutinz->SetMarkerColor(1);
162 fHistClusterMapITSSAoutinz->SetMarkerStyle(22);
163 fHistClusterMapITSSAoutinz->Draw();
164 fHistClusterMapITSMIoutinz->SetLineColor(2);
165 fHistClusterMapITSMIoutinz->SetMarkerColor(2);
166 fHistClusterMapITSMIoutinz->SetMarkerStyle(22);
167 fHistClusterMapITSMIoutinz->Draw("same");
168 //
169 c2->cd(4);
170 fHistClusterMapITSSAbad->SetLineColor(1);
171 fHistClusterMapITSSAbad->SetMarkerColor(1);
172 fHistClusterMapITSSAbad->SetMarkerStyle(23);
173 fHistClusterMapITSSAbad->Draw();
174 fHistClusterMapITSMIbad->SetLineColor(2);
175 fHistClusterMapITSMIbad->SetMarkerColor(2);
176 fHistClusterMapITSMIbad->SetMarkerStyle(23);
177 fHistClusterMapITSMIbad->Draw("same");
178 //
179 c2->cd(5);
180 fHistClusterMapITSSAnocls->SetLineColor(1);
181 fHistClusterMapITSSAnocls->SetMarkerColor(1);
182 fHistClusterMapITSSAnocls->SetMarkerStyle(24);
183 fHistClusterMapITSSAnocls->Draw();
184 fHistClusterMapITSMInocls->SetLineColor(2);
185 fHistClusterMapITSMInocls->SetMarkerColor(2);
186 fHistClusterMapITSMInocls->SetMarkerStyle(24);
187 fHistClusterMapITSMInocls->Draw("same");
188
fa516524 189 TCanvas *c3 =new TCanvas("c3","c3",10,10,1200,800);
190 c3->Divide(1,2);
191 c3->cd(1);
192 fHistClusterMapModuleITSMIokInAcc->SetLineColor(1);
193 fHistClusterMapModuleITSMIokInAcc->Draw();
194 l2->AddEntry(fHistClusterMapModuleITSMIokInAcc,"ok","l");
195 fHistClusterMapModuleITSMIbadInAcc->SetLineColor(3);
196 fHistClusterMapModuleITSMIbadInAcc->Draw("same");
197 l2->AddEntry(fHistClusterMapModuleITSMIbadInAcc,"bad","l");
198 fHistClusterMapModuleITSMInoclsInAcc->SetLineColor(2);
199 fHistClusterMapModuleITSMInoclsInAcc->Draw("same");
200 l2->AddEntry(fHistClusterMapModuleITSMInoclsInAcc,"no cls","l");
2308c814 201 l2->Draw();
fa516524 202 c3->cd(2);
203 TH1F *fHistClusterMapModuleITSMIallInAcc=(TH1F*)fHistClusterMapModuleITSMIokInAcc->Clone("fHistClusterMapModuleITSMIallInAcc");
204 fHistClusterMapModuleITSMIallInAcc->Add(fHistClusterMapModuleITSMIbadInAcc);
205 fHistClusterMapModuleITSMIallInAcc->Add(fHistClusterMapModuleITSMInoclsInAcc);
206 TH1F *fHistClusterMapModuleITSMIokRatioInAcc=(TH1F*)fHistClusterMapModuleITSMIokInAcc->Clone("fHistClusterMapModuleITSMIokRatioInAcc");
207 fHistClusterMapModuleITSMIokRatioInAcc->Divide(fHistClusterMapModuleITSMIallInAcc);
208 for(Int_t ib=1;ib<=fHistClusterMapModuleITSMIokRatioInAcc->GetNbinsX();ib++) {
209 fHistClusterMapModuleITSMIokRatioInAcc->SetBinError(ib,0);
210 if(fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)) fHistClusterMapModuleITSMIokRatioInAcc->SetBinError(ib,TMath::Sqrt(fHistClusterMapModuleITSMIokRatioInAcc->GetBinContent(ib)*(1.-fHistClusterMapModuleITSMIokRatioInAcc->GetBinContent(ib))/fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)));
211 }
212 fHistClusterMapModuleITSMIokRatioInAcc->Draw();
213 TH1F *fHistClusterMapModuleITSMIbadRatioInAcc=(TH1F*)fHistClusterMapModuleITSMIbadInAcc->Clone("fHistClusterMapModuleITSMIbadRatioInAcc");
214 fHistClusterMapModuleITSMIbadRatioInAcc->Divide(fHistClusterMapModuleITSMIallInAcc);
215 for(Int_t ib=1;ib<=fHistClusterMapModuleITSMIokRatioInAcc->GetNbinsX();ib++) {
216 fHistClusterMapModuleITSMIbadRatioInAcc->SetBinError(ib,0);
217 if(fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)) fHistClusterMapModuleITSMIbadRatioInAcc->SetBinError(ib,TMath::Sqrt(fHistClusterMapModuleITSMIbadRatioInAcc->GetBinContent(ib)*(1.-fHistClusterMapModuleITSMIbadRatioInAcc->GetBinContent(ib))/fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)));
218 }
219 fHistClusterMapModuleITSMIbadRatioInAcc->Draw("same");
220 TH1F *fHistClusterMapModuleITSMInoclsRatioInAcc=(TH1F*)fHistClusterMapModuleITSMInoclsInAcc->Clone("fHistClusterMapModuleITSMInoclsRatioInAcc");
221 fHistClusterMapModuleITSMInoclsRatioInAcc->Divide(fHistClusterMapModuleITSMIallInAcc);
222 for(Int_t ib=1;ib<=fHistClusterMapModuleITSMIokRatioInAcc->GetNbinsX();ib++) {
223 fHistClusterMapModuleITSMInoclsRatioInAcc->SetBinError(ib,0);
224 if(fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)) fHistClusterMapModuleITSMInoclsRatioInAcc->SetBinError(ib,TMath::Sqrt(fHistClusterMapModuleITSMInoclsRatioInAcc->GetBinContent(ib)*(1.-fHistClusterMapModuleITSMInoclsRatioInAcc->GetBinContent(ib))/fHistClusterMapModuleITSMIallInAcc->GetBinContent(ib)));
225 }
226 fHistClusterMapModuleITSMInoclsRatioInAcc->Draw("same");
2308c814 227
fa516524 228 TCanvas *c3b = new TCanvas("c3b","c3b");
229 c3b->Divide(2,1);
230 c3b->cd(1);
231 fHistxlocSDDok->Divide(fHistxlocSDDall);
232 fHistxlocSDDok->Draw();
233 c3b->cd(2);
234 fHistzlocSDDok->Divide(fHistzlocSDDall);
235 fHistzlocSDDok->Draw();
2308c814 236
c2893e82 237 TCanvas *c4 =new TCanvas("c4","c4",10,10,500,500);
238 c4->SetGridy();
239 fHistPhiITSMIokbadoutinz6InAcc->Divide(fHistPhiTPCInAcc);
240 fHistPhiITSMIokbadoutinz6InAcc->SetMinimum(0);
241 fHistPhiITSMIokbadoutinz6InAcc->SetMaximum(1.5);
242 fHistPhiITSMIokbadoutinz6InAcc->SetYTitle("ITS+TPC / TPC");
243 fHistPhiITSMIokbadoutinz6InAcc->SetTitle("Fraction of tracks with 6 layers ok");
244 fHistPhiITSMIokbadoutinz6InAcc->Draw();
245
fa516524 246 TCanvas *c5c =new TCanvas("c5c","c5c",10,10,600,600);
247 c5c->SetGridy();
248 TH1F *fHistPtTPCInAccSecFrac = (TH1F*)fHistPtTPCInAccS->Clone("fHistPtTPCInAccSecFrac");
249 fHistPtTPCInAccSecFrac->Divide(fHistPtTPCInAcc);
250 fHistPtTPCInAccSecFrac->SetLineColor(kOrange+7);
251 fHistPtTPCInAccSecFrac->SetTitle("Fraction of secondaries");
252 fHistPtTPCInAccSecFrac->SetYTitle("sec/all");
253 fHistPtTPCInAccSecFrac->Draw();
254 TH1F *fHistPtITSMI2InAccSecFrac = (TH1F*)fHistPtITSMI2InAccS->Clone("fHistPtITSMI2InAccSecFrac");
255 fHistPtITSMI2InAccSecFrac->Divide(fHistPtITSMI2InAcc);
256 fHistPtITSMI2InAccSecFrac->SetLineColor(7);
257 fHistPtITSMI2InAccSecFrac->Draw("same");
258 TH1F *fHistPtITSMIge2InAccSecFrac = (TH1F*)fHistPtITSMIge2InAccS->Clone("fHistPtITSMIge2InAccSecFrac");
259 fHistPtITSMIge2InAccSecFrac->Divide(fHistPtITSMIge2InAcc);
260 fHistPtITSMIge2InAccSecFrac->Draw("same");
261 TH1F *fHistPtITSMI3InAccSecFrac = (TH1F*)fHistPtITSMI3InAccS->Clone("fHistPtITSMI3InAccSecFrac");
262 fHistPtITSMI3InAccSecFrac->Divide(fHistPtITSMI3InAcc);
263 fHistPtITSMI3InAccSecFrac->SetLineColor(6);
264 fHistPtITSMI3InAccSecFrac->Draw("same");
265 TH1F *fHistPtITSMI4InAccSecFrac = (TH1F*)fHistPtITSMI4InAccS->Clone("fHistPtITSMI4InAccSecFrac");
266 fHistPtITSMI4InAccSecFrac->Divide(fHistPtITSMI4InAcc);
267 fHistPtITSMI4InAccSecFrac->SetLineColor(4);
268 fHistPtITSMI4InAccSecFrac->Draw("same");
269 TH1F *fHistPtITSMI5InAccSecFrac = (TH1F*)fHistPtITSMI5InAccS->Clone("fHistPtITSMI5InAccSecFrac");
270 fHistPtITSMI5InAccSecFrac->Divide(fHistPtITSMI5InAcc);
271 fHistPtITSMI5InAccSecFrac->SetLineColor(3);
272 fHistPtITSMI5InAccSecFrac->Draw("same");
273 TH1F *fHistPtITSMI6InAccSecFrac = (TH1F*)fHistPtITSMI6InAccS->Clone("fHistPtITSMI6InAccSecFrac");
274 fHistPtITSMI6InAccSecFrac->Divide(fHistPtITSMI6InAcc);
275 fHistPtITSMI6InAccSecFrac->SetLineColor(2);
276 fHistPtITSMI6InAccSecFrac->SetLineColor(2);
277 fHistPtITSMI6InAccSecFrac->Draw("same");
278 TH1F *fHistPtITSMISPDInAccSecFrac = (TH1F*)fHistPtITSMISPDInAccS->Clone("fHistPtITSMISPDInAccSecFrac");
279 fHistPtITSMISPDInAccSecFrac->Divide(fHistPtITSMISPDInAcc);
280 fHistPtITSMISPDInAccSecFrac->SetLineColor(9);
281 fHistPtITSMISPDInAccSecFrac->Draw("same");
282 TH1F *fHistPtITSMIoneSPDInAccSecFrac = (TH1F*)fHistPtITSMIoneSPDInAccS->Clone("fHistPtITSMIoneSPDInAccSecFrac");
283 fHistPtITSMIoneSPDInAccSecFrac->Divide(fHistPtITSMIoneSPDInAcc);
284 fHistPtITSMIoneSPDInAccSecFrac->SetLineColor(15);
285 fHistPtITSMIoneSPDInAccSecFrac->Draw("same");
286 TH1F *fHistPtITSTPCselSecFrac = (TH1F*)fHistPtITSTPCselS->Clone("fHistPtITSTPCselSecFrac");
287 fHistPtITSTPCselSecFrac->Divide(fHistPtITSTPCsel);
288 fHistPtITSTPCselSecFrac->SetLineColor(15);
289 fHistPtITSTPCselSecFrac->Draw("same");
290
291
c2893e82 292 TLegend *l3=new TLegend(0.5,0.5,0.9,0.9);
293 TLegend *l4=new TLegend(0.5,0.5,0.9,0.9);
294 TCanvas *c5 =new TCanvas("c5","c5",10,10,1200,600);
295 c5->Divide(2,1);
296 c5_1->SetGridy();
297 c5_2->SetGridy();
298 c5->cd(1);
c2893e82 299 fHistPtITSMIge2InAcc->SetMaximum(1.5);
300 fHistPtITSMIge2InAcc->SetMinimum(0);
301 fHistPtITSMIge2InAcc->SetTitle("Fraction of prolonged tracks with N ITS points");
302 fHistPtITSMIge2InAcc->SetYTitle("ITS+TPC / TPC");
fa516524 303 fHistPtITSMIge2InAcc->Divide(fHistPtTPCInAcc);
c2893e82 304 fHistPtITSMIge2InAcc->Draw();
305 l3->AddEntry(fHistPtITSMIge2InAcc,">=2 cls","l");
306 fHistPtITSMI6InAcc->Divide(fHistPtTPCInAcc);
307 fHistPtITSMI6InAcc->SetLineColor(2);
308 l3->AddEntry(fHistPtITSMI6InAcc,"6 cls","l");
309 fHistPtITSMI6InAcc->Draw("same");
310 fHistPtITSMI5InAcc->Divide(fHistPtTPCInAcc);
311 fHistPtITSMI5InAcc->SetLineColor(3);
312 l3->AddEntry(fHistPtITSMI5InAcc,"5 cls","l");
313 fHistPtITSMI5InAcc->Draw("same");
314 fHistPtITSMI4InAcc->Divide(fHistPtTPCInAcc);
315 fHistPtITSMI4InAcc->SetLineColor(4);
316 l3->AddEntry(fHistPtITSMI4InAcc,"4 cls","l");
317 fHistPtITSMI4InAcc->Draw("same");
318 fHistPtITSMI3InAcc->Divide(fHistPtTPCInAcc);
319 fHistPtITSMI3InAcc->SetLineColor(6);
320 l3->AddEntry(fHistPtITSMI3InAcc,"3 cls","l");
321 fHistPtITSMI3InAcc->Draw("same");
322 fHistPtITSMI2InAcc->Divide(fHistPtTPCInAcc);
323 fHistPtITSMI2InAcc->SetLineColor(7);
324 l3->AddEntry(fHistPtITSMI2InAcc,"2 cls","l");
325 fHistPtITSMI2InAcc->Draw("same");
326 fHistPtITSMISPDInAcc->Divide(fHistPtTPCInAcc);
327 fHistPtITSMISPDInAcc->SetLineColor(9);
328 l3->AddEntry(fHistPtITSMISPDInAcc,"2SPD + any","l");
329 fHistPtITSMISPDInAcc->Draw("same");
fa516524 330 fHistPtITSMIoneSPDInAcc->Divide(fHistPtTPCInAcc);
331 fHistPtITSMIoneSPDInAcc->SetLineColor(15);
332 l3->AddEntry(fHistPtITSMIoneSPDInAcc,">=1SPD + any","l");
333 fHistPtITSMIoneSPDInAcc->Draw("same");
334 fHistPtITSTPCsel->Divide(fHistPtTPCInAcc);
335 fHistPtITSTPCsel->SetLineColor(kAzure+1);
336 l3->AddEntry(fHistPtITSTPCsel,">=1SPD + any + d_{0} cut","l");
337 fHistPtITSTPCsel->Draw("same");
c2893e82 338 fHistPtITSMIge2InAcc->Draw("same");
339 l3->Draw();
340 c5->cd(2);
341 TH1F *fHistPtITSMIokbadoutinzge4InAcc = (TH1F*)fHistPtITSMIokbadoutinz6InAcc->Clone("fHistPtITSMIokbadoutinzge4InAcc");
342 fHistPtITSMIokbadoutinzge4InAcc->Add(fHistPtITSMIokbadoutinz5InAcc);
343 fHistPtITSMIokbadoutinzge4InAcc->Add(fHistPtITSMIokbadoutinz4InAcc);
344 fHistPtITSMIokbadoutinzge4InAcc->SetMaximum(1.5);
345 fHistPtITSMIokbadoutinzge4InAcc->SetMinimum(0);
346 fHistPtITSMIokbadoutinzge4InAcc->SetTitle("Fraction of prolonged tracks with N ITS layers \"ok\"");
347 fHistPtITSMIokbadoutinzge4InAcc->SetYTitle("ITS+TPC / TPC");
348 fHistPtITSMIokbadoutinzge4InAcc->Divide(fHistPtTPCInAcc);
349 fHistPtITSMIokbadoutinzge4InAcc->SetLineColor(1);
350 fHistPtITSMIokbadoutinzge4InAcc->Draw();
351 fHistPtITSMIokbadoutinz6InAcc->Divide(fHistPtTPCInAcc);
352 fHistPtITSMIokbadoutinz6InAcc->SetLineColor(2);
353 fHistPtITSMIokbadoutinz6InAcc->Draw("same");
354 fHistPtITSMIokbadoutinz5InAcc->Divide(fHistPtTPCInAcc);
355 fHistPtITSMIokbadoutinz5InAcc->SetLineColor(3);
356 fHistPtITSMIokbadoutinz5InAcc->Draw("same");
357 fHistPtITSMIokbadoutinz4InAcc->Divide(fHistPtTPCInAcc);
358 fHistPtITSMIokbadoutinz4InAcc->SetLineColor(4);
359 fHistPtITSMIokbadoutinz4InAcc->Draw("same");
360 fHistPtITSMIokbadoutinzge4InAcc->Draw("same");
361 l4->AddEntry(fHistPtITSMIokbadoutinzge4InAcc,">=4 layers","l");
362 l4->AddEntry(fHistPtITSMIokbadoutinz6InAcc,"6 layers","l");
363 l4->AddEntry(fHistPtITSMIokbadoutinz5InAcc,"5 layers","l");
364 l4->AddEntry(fHistPtITSMIokbadoutinz4InAcc,"4 layers","l");
365 l4->Draw();
366
367
fa516524 368 TLegend *l4d=new TLegend(0.5,0.5,0.9,0.9);
369 TCanvas *c5d =new TCanvas("c5d","c5d",10,10,600,600);
370 c5d->SetGridy();
371 c5d->SetGridy();
372 fHistPtITSMISPDInAccMC->SetMaximum(1.5);
373 fHistPtITSMISPDInAccMC->SetMinimum(0);
374 fHistPtITSMISPDInAccMC->SetTitle("Fraction of prolonged tracks with N ITS points");
375 fHistPtITSMISPDInAccMC->SetYTitle("ITS+TPC / TPC");
376 fHistPtITSMISPDInAccMC->Divide(fHistPtTPCInAccMCtwoSPD);
377 fHistPtITSMISPDInAccMC->SetLineColor(9);
378 l4d->AddEntry(fHistPtITSMISPDInAccMC,"2SPD + any","l");
379 fHistPtITSMISPDInAccMC->Draw();
380 fHistPtITSMIoneSPDInAccMC->Divide(fHistPtTPCInAccMConeSPD);
381 fHistPtITSMIoneSPDInAccMC->SetLineColor(15);
382 l4d->AddEntry(fHistPtITSMIoneSPDInAccMC,">=1SPD + any","l");
383 fHistPtITSMIoneSPDInAccMC->Draw("same");
384 l4d->Draw();
385
386 TFile *eff=new TFile("eff.root","recreate");
387 fHistPtITSMIge2InAcc->Write();
388 fHistPtITSMI6InAcc->Write();
389 fHistPtITSMI5InAcc->Write();
390 fHistPtITSMI4InAcc->Write();
391 fHistPtITSMI3InAcc->Write();
392 fHistPtITSMI2InAcc->Write();
393 fHistPtITSMISPDInAcc->Write();
394 fHistPtITSMIoneSPDInAcc->Write();
395 fHistPtITSTPCsel->Write();
396 fHistClusterMapModuleITSMIokRatioInAcc->Write();
397 fHistClusterMapModuleITSMIbadRatioInAcc->Write();
398 fHistClusterMapModuleITSMInoclsRatioInAcc->Write();
399 eff->Close();
400
401 c5c->cd();
402 l3->Draw();
403
404
405 TCanvas *c5b =new TCanvas("c5b","c5b",10,10,600,600);
406 c5b->SetLogy();
407 fHistRProdVtxInAccP->SetLineColor(2);
408 fHistRProdVtxInAccP->SetMinimum(1);
409 fHistRProdVtxInAccS->SetMinimum(1);
410 fHistRProdVtxInAccP->Draw();
411 fHistRProdVtxInAccS->Draw("sames");
412
413 TCanvas *c5a =new TCanvas("c5a","c5a",10,10,1200,600);
414 c5a->Divide(2,1);
415 c5a_1->SetGridy();
416 c5a_2->SetGridy();
417 c5a->cd(1);
418 fHistPtITSMIge2InAccP->SetMaximum(1.5);
419 fHistPtITSMIge2InAccP->SetMinimum(0);
420 fHistPtITSMIge2InAccP->SetTitle("Fraction of prolonged tracks with N ITS points");
421 fHistPtITSMIge2InAccP->SetYTitle("ITS+TPC / TPC");
422 fHistPtITSMIge2InAccP->Draw();
423 fHistPtITSMIge2InAccP->Divide(fHistPtTPCInAccP);
424 fHistPtITSMI6InAccP->Divide(fHistPtTPCInAccP);
425 fHistPtITSMI6InAccP->SetLineColor(2);
426 fHistPtITSMI6InAccP->Draw("same");
427 fHistPtITSMI5InAccP->Divide(fHistPtTPCInAccP);
428 fHistPtITSMI5InAccP->SetLineColor(3);
429 fHistPtITSMI5InAccP->Draw("same");
430 fHistPtITSMI4InAccP->Divide(fHistPtTPCInAccP);
431 fHistPtITSMI4InAccP->SetLineColor(4);
432 fHistPtITSMI4InAccP->Draw("same");
433 fHistPtITSMI3InAccP->Divide(fHistPtTPCInAccP);
434 fHistPtITSMI3InAccP->SetLineColor(6);
435 fHistPtITSMI3InAccP->Draw("same");
436 fHistPtITSMI2InAccP->Divide(fHistPtTPCInAccP);
437 fHistPtITSMI2InAccP->SetLineColor(7);
438 fHistPtITSMI2InAccP->Draw("same");
439 fHistPtITSMISPDInAccP->Divide(fHistPtTPCInAccP);
440 fHistPtITSMISPDInAccP->SetLineColor(9);
441 fHistPtITSMISPDInAccP->Draw("same");
442 fHistPtITSMIoneSPDInAccP->Divide(fHistPtTPCInAccP);
443 fHistPtITSMIoneSPDInAccP->SetLineColor(15);
444 fHistPtITSMIoneSPDInAccP->Draw("same");
445 fHistPtITSTPCselP->Divide(fHistPtTPCInAccP);
446 fHistPtITSTPCselP->SetLineColor(kAzure+1);
447 fHistPtITSTPCselP->Draw("same");
448 fHistPtITSMIge2InAccP->Draw("same");
449 l3->Draw();
450 c5a->cd(2);
451 fHistPtITSMIge2InAccS->SetMaximum(1.5);
452 fHistPtITSMIge2InAccS->SetMinimum(0);
453 fHistPtITSMIge2InAccS->SetTitle("Fraction of prolonged tracks with N ITS points");
454 fHistPtITSMIge2InAccS->SetYTitle("ITS+TPC / TPC");
455 fHistPtITSMIge2InAccS->Divide(fHistPtTPCInAccS);
456 fHistPtITSMIge2InAccS->Draw();
457 fHistPtITSMI6InAccS->Divide(fHistPtTPCInAccS);
458 fHistPtITSMI6InAccS->SetLineColor(2);
459 fHistPtITSMI6InAccS->Draw("same");
460 fHistPtITSMI5InAccS->Divide(fHistPtTPCInAccS);
461 fHistPtITSMI5InAccS->SetLineColor(3);
462 fHistPtITSMI5InAccS->Draw("same");
463 fHistPtITSMI4InAccS->Divide(fHistPtTPCInAccS);
464 fHistPtITSMI4InAccS->SetLineColor(4);
465 fHistPtITSMI4InAccS->Draw("same");
466 fHistPtITSMI3InAccS->Divide(fHistPtTPCInAccS);
467 fHistPtITSMI3InAccS->SetLineColor(6);
468 fHistPtITSMI3InAccS->Draw("same");
469 fHistPtITSMI2InAccS->Divide(fHistPtTPCInAccS);
470 fHistPtITSMI2InAccS->SetLineColor(7);
471 fHistPtITSMI2InAccS->Draw("same");
472 fHistPtITSMISPDInAccS->Divide(fHistPtTPCInAccS);
473 fHistPtITSMISPDInAccS->SetLineColor(9);
474 fHistPtITSMISPDInAccS->Draw("same");
475 fHistPtITSMIoneSPDInAccS->Divide(fHistPtTPCInAccS);
476 fHistPtITSMIoneSPDInAccS->SetLineColor(15);
477 fHistPtITSMIoneSPDInAccS->Draw("same");
478 fHistPtITSTPCselS->Divide(fHistPtTPCInAccS);
479 fHistPtITSTPCselS->SetLineColor(kAzure+1);
480 fHistPtITSTPCselS->Draw("same");
481 fHistPtITSMIge2InAccS->Draw("same");
482 l3->Draw();
483
484
485
486
c2893e82 487 // PLOT ALIGNMENT CHECKS
488 //
489 TH1F *hSPDTrackletsTBdxy = new TH1F("hSPDTrackletsTBdxy","SPD tracklets; SPD tracklet to SPD vertex distance in (x,y) [cm]; entries",100,-0.1,0.1);
490 TH1F *hSPDTrackletsLRdxy = new TH1F("hSPDTrackletsLRdxy","SPD tracklets; SPD tracklet to SPD vertex distance in (x,y) [cm]; entries",100,-0.1,0.1);
491 TH1F *hSPDTrackletsTBdz = new TH1F("hSPDTrackletsTBdz","SPD tracklets; SPD tracklet to SPD vertex distance in z [cm]; entries",100,-0.1,0.1);
492 TH1F *hSPDTrackletsLRdz = new TH1F("hSPDTrackletsLRdz","SPD tracklets; SPD tracklet to SPD vertex distance in z [cm]; entries",100,-0.1,0.1);
493
494 TNtuple *fNtupleITSAlignSPDTracklets = (TNtuple*)list->FindObject("fNtupleITSAlignSPDTracklets");
495 Float_t dxy,dz,phi,pt;
496 fNtupleITSAlignSPDTracklets->SetBranchAddress("pt",&pt);
497 fNtupleITSAlignSPDTracklets->SetBranchAddress("phi",&phi);
498 fNtupleITSAlignSPDTracklets->SetBranchAddress("dxy",&dxy);
499 fNtupleITSAlignSPDTracklets->SetBranchAddress("dz",&dz);
500
501 for(Int_t i=0;i<fNtupleITSAlignSPDTracklets->GetEntries();i++) {
502 fNtupleITSAlignSPDTracklets->GetEvent(i);
fa516524 503 if(pt<.01) continue;
c2893e82 504 if(TMath::Abs(TMath::Abs(phi)-0.5*TMath::Pi())<0.25*TMath::Pi()) {
505 // top-bottom
506 hSPDTrackletsTBdxy->Fill(dxy);
507 hSPDTrackletsTBdz->Fill(dz);
508 } else {
509 // left-right
510 hSPDTrackletsLRdxy->Fill(dxy);
511 hSPDTrackletsLRdz->Fill(dz);
512 }
513 }
514
515 TLegend *l6=new TLegend(0.5,0.5,0.9,0.9);
516
517 TCanvas *c6 = new TCanvas("c6","c6",0,0,1000,500);
518 c6->Divide(2,1);
519 c6->cd(1);
520 hSPDTrackletsTBdxy->SetLineColor(4);
521 hSPDTrackletsTBdxy->Draw();
522 hSPDTrackletsLRdxy->SetLineColor(2);
fa516524 523 hSPDTrackletsLRdxy->Draw("sames");
c2893e82 524 l6->AddEntry(hSPDTrackletsTBdxy,"top-bottom","l");
525 l6->AddEntry(hSPDTrackletsLRdxy,"left-right","l");
526 l6->Draw();
527 c6->cd(2);
528 hSPDTrackletsTBdz->SetLineColor(4);
529 hSPDTrackletsTBdz->Draw();
530 hSPDTrackletsLRdz->SetLineColor(2);
fa516524 531 hSPDTrackletsLRdz->Draw("sames");
c2893e82 532 l6->Draw();
533
534
535 TH1F *hSPDExtraClsTBdxy = new TH1F("hSPDExtraClsTBdxy","SPD extra clusters; track-to-point distance in (x,y) [cm]; entries",100,-0.1,0.1);
536 TH1F *hSPDExtraClsLRdxy = new TH1F("hSPDExtraClsLRdxy","SPD extra clusters; track-to-point distance in (x,y) [cm]; entries",100,-0.1,0.1);
537 TH1F *hSPDExtraClsTBdz = new TH1F("hSPDExtraClsTBdz","SPD extra clusters; track-to-point distance in z [cm]; entries",100,-0.1,0.1);
538 TH1F *hSPDExtraClsLRdz = new TH1F("hSPDExtraClsLRdz","SPD extra clusters; track-to-point distance in z [cm]; entries",100,-0.1,0.1);
539
540 TNtuple *fNtupleITSAlignExtra = (TNtuple*)list->FindObject("fNtupleITSAlignExtra");
541 Float_t layer,npoints,x,y,z;
542 fNtupleITSAlignExtra->SetBranchAddress("layer",&layer);
543 fNtupleITSAlignExtra->SetBranchAddress("npoints",&npoints);
544 fNtupleITSAlignExtra->SetBranchAddress("x",&x);
545 fNtupleITSAlignExtra->SetBranchAddress("y",&y);
546 fNtupleITSAlignExtra->SetBranchAddress("z",&z);
547 fNtupleITSAlignExtra->SetBranchAddress("dxy",&dxy);
548 fNtupleITSAlignExtra->SetBranchAddress("dz",&dz);
549 fNtupleITSAlignExtra->SetBranchAddress("pt",&pt);
550
551 for(Int_t i=0;i<fNtupleITSAlignExtra->GetEntries();i++) {
552 fNtupleITSAlignExtra->GetEvent(i);
fa516524 553 if(pt<0.1) continue;
c2893e82 554 if(layer!=0 && layer!=1) continue;
fa516524 555 if(npoints<3) continue;
c2893e82 556 phi = TMath::ATan2(y,x);
557 if(TMath::Abs(TMath::Abs(phi)-0.5*TMath::Pi())<0.25*TMath::Pi()) {
558 // top-bottom
559 hSPDExtraClsTBdxy->Fill(dxy);
560 hSPDExtraClsTBdz->Fill(dz);
561 } else {
562 // left-right
563 hSPDExtraClsLRdxy->Fill(dxy);
564 hSPDExtraClsLRdz->Fill(dz);
565 }
566 }
567
568 TCanvas *c7 = new TCanvas("c7","c7",0,0,1000,500);
569 c7->Divide(2,1);
570 c7->cd(1);
571 hSPDExtraClsTBdxy->SetLineColor(4);
572 hSPDExtraClsTBdxy->Draw();
573 hSPDExtraClsLRdxy->SetLineColor(2);
574 hSPDExtraClsLRdxy->Draw("same");
575 l6->Draw();
576 c7->cd(2);
577 hSPDExtraClsTBdz->SetLineColor(4);
578 hSPDExtraClsTBdz->Draw();
579 hSPDExtraClsLRdz->SetLineColor(2);
580 hSPDExtraClsLRdz->Draw("same");
581 l6->Draw();
582
583
fa516524 584 return;
585}
586//-----------------------------------------------------------------------------
587void PlotEffRatio() {
588
589 TFile *f= new TFile("eff_lhc10a8_100k_trkvtx.root");
590
591 TH1F *fHistPtITSMI6InAcc = (TH1F*)f->Get("fHistPtITSMI6InAcc");
592 TH1F *fHistPtITSMI5InAcc = (TH1F*)f->Get("fHistPtITSMI5InAcc");
593 TH1F *fHistPtITSMI4InAcc = (TH1F*)f->Get("fHistPtITSMI4InAcc");
594 TH1F *fHistPtITSMI3InAcc = (TH1F*)f->Get("fHistPtITSMI3InAcc");
595 TH1F *fHistPtITSMI2InAcc = (TH1F*)f->Get("fHistPtITSMI2InAcc");
596 TH1F *fHistPtITSMIge2InAcc = (TH1F*)f->Get("fHistPtITSMIge2InAcc");
597 TH1F *fHistPtITSMISPDInAcc = (TH1F*)f->Get("fHistPtITSMISPDInAcc");
598 TH1F *fHistPtITSMIoneSPDInAcc = (TH1F*)f->Get("fHistPtITSMIoneSPDInAcc");
599 TH1F *fHistPtITSTPCsel = (TH1F*)f->Get("fHistPtITSTPCsel");
600 TH1F *fHistClusterMapModuleITSMIokRatioInAcc = (TH1F*)f->Get("fHistClusterMapModuleITSMIokRatioInAcc");
601
602 TFile *fMC= new TFile("eff_lhc10a8_100k_spdvtx.root");
603
604 TH1F *fHistPtITSMI6InAccMC = (TH1F*)fMC->Get("fHistPtITSMI6InAcc");
605 TH1F *fHistPtITSMI5InAccMC = (TH1F*)fMC->Get("fHistPtITSMI5InAcc");
606 TH1F *fHistPtITSMI4InAccMC = (TH1F*)fMC->Get("fHistPtITSMI4InAcc");
607 TH1F *fHistPtITSMI3InAccMC = (TH1F*)fMC->Get("fHistPtITSMI3InAcc");
608 TH1F *fHistPtITSMI2InAccMC = (TH1F*)fMC->Get("fHistPtITSMI2InAcc");
609 TH1F *fHistPtITSMIge2InAccMC = (TH1F*)fMC->Get("fHistPtITSMIge2InAcc");
610 TH1F *fHistPtITSMISPDInAccMC = (TH1F*)fMC->Get("fHistPtITSMISPDInAcc");
611 TH1F *fHistPtITSMIoneSPDInAccMC = (TH1F*)fMC->Get("fHistPtITSMIoneSPDInAcc");
612 TH1F *fHistPtITSTPCselMC = (TH1F*)fMC->Get("fHistPtITSTPCsel");
613 TH1F *fHistClusterMapModuleITSMIokRatioInAccMC = (TH1F*)fMC->Get("fHistClusterMapModuleITSMIokRatioInAcc");
614
615 TCanvas *c = new TCanvas("c","c",0,0,600,500);
616 c->SetGridy();
617 fHistPtITSMI6InAcc->Divide(fHistPtITSMI6InAccMC);
618 fHistPtITSMI6InAcc->Draw();
619 fHistPtITSMI5InAcc->Divide(fHistPtITSMI5InAccMC);
620 fHistPtITSMI5InAcc->Draw("same");
621 fHistPtITSMI4InAcc->Divide(fHistPtITSMI4InAccMC);
622 fHistPtITSMI4InAcc->Draw("same");
623 fHistPtITSMI3InAcc->Divide(fHistPtITSMI3InAccMC);
624 fHistPtITSMI3InAcc->Draw("same");
625 fHistPtITSMI2InAcc->Divide(fHistPtITSMI2InAccMC);
626 fHistPtITSMI2InAcc->Draw("same");
627 fHistPtITSMIge2InAcc->Divide(fHistPtITSMIge2InAccMC);
628 fHistPtITSMIge2InAcc->Draw("same");
629 fHistPtITSMISPDInAcc->Divide(fHistPtITSMISPDInAccMC);
630 fHistPtITSMISPDInAcc->Draw("same");
631 fHistPtITSMIoneSPDInAcc->Divide(fHistPtITSMIoneSPDInAccMC);
632 fHistPtITSMIoneSPDInAcc->Draw("same");
633 fHistPtITSTPCsel->Divide(fHistPtITSTPCselMC);
634 fHistPtITSTPCsel->Draw("same");
635 TLegend *l3=new TLegend(0.5,0.5,0.9,0.9);
636 l3->AddEntry(fHistPtITSMIge2InAcc,">=2 cls","l");
637 l3->AddEntry(fHistPtITSMI6InAcc,"6 cls","l");
638 l3->AddEntry(fHistPtITSMI5InAcc,"5 cls","l");
639 l3->AddEntry(fHistPtITSMI4InAcc,"4 cls","l");
640 l3->AddEntry(fHistPtITSMI3InAcc,"3 cls","l");
641 l3->AddEntry(fHistPtITSMI2InAcc,"2 cls","l");
642 l3->AddEntry(fHistPtITSMISPDInAcc,"2SPD + any","l");
643 l3->AddEntry(fHistPtITSMIoneSPDInAcc,">=1SPD + any","l");
644 l3->AddEntry(fHistPtITSTPCsel,">=1SPD + any + d_{0} cut","l");
645 l3->Draw();
646
647 TCanvas *cc = new TCanvas("cc","cc",0,0,600,500);
648 cc->Divide(1,2);
649 cc->cd(1);
650 fHistClusterMapModuleITSMIokRatioInAccMC->SetLineColor(4);
651 fHistClusterMapModuleITSMIokRatioInAccMC->Draw();
652 fHistClusterMapModuleITSMIokRatioInAcc->Draw("same");
653 cc->cd(2);
654 for(Int_t i=1;i<=fHistClusterMapModuleITSMIokRatioInAccMC->GetNbinsX();i++) {
655 if(fHistClusterMapModuleITSMIokRatioInAccMC->GetBinContent(i)<0.02) {
656 fHistClusterMapModuleITSMIokRatioInAccMC->SetBinContent(i,1);
657 if(fHistClusterMapModuleITSMIokRatioInAcc->GetBinContent(i)<0.02) {
658 fHistClusterMapModuleITSMIokRatioInAcc->SetBinContent(i,1);
659 }
660 }
661 }
662 TH1F* fHistClusterMapDataOverMC=(TH1F*)fHistClusterMapModuleITSMIokRatioInAcc->Clone("fHistClusterMapDataOverMC");
663 fHistClusterMapDataOverMC->Divide(fHistClusterMapModuleITSMIokRatioInAccMC);
664 fHistClusterMapDataOverMC->Draw();
665
666 return;
667}
668
669void PlotImpPar_rphi(Int_t rebin=1) {
670
671
672 TFile *fMC= new TFile("ITS.Performance.root");
673
674 TList *list=(TList*)fMC->Get("cOutputITS");
675 TH1F *fHistd0rphiITSMIoneSPDInAccP150200MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
676 TH1F *fHistd0rphiITSMIoneSPDInAccS150200MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
677 TH1F *fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200fromStrange");
678 TH1F *fHistd0rphiITSMIoneSPDInAccS150200fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200fromMat");
679 /* TH1F *fHistd0rphiITSMIoneSPDInAccP350450MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP350450");
680 TH1F *fHistd0rphiITSMIoneSPDInAccS350450MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450");
681 TH1F *fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450fromStrange");
682 TH1F *fHistd0rphiITSMIoneSPDInAccS350450fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450fromMat");*/
683 TH1F *fHistd0rphiITSMIoneSPDInAccP500700MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP500700");
684 TH1F *fHistd0rphiITSMIoneSPDInAccS500700MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700");
685 TH1F *fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700fromStrange");
686 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from211MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from211");
687 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from22MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from22");
688 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from310MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from310");
689 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from321MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from321");
690 TH1F *fHistd0rphiITSMIoneSPDInAccS500700from3122MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700from3122");
691 TH1F *fHistd0rphiITSMIoneSPDInAccS500700fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700fromMat");
692 TH1F *fHistd0rphiITSMIoneSPDInAccP10001500MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP10001500");
693 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500");
694 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500fromStrange");
695 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500fromMat");
696 /* TH1F *fHistd0rphiITSMIoneSPDInAccP25004000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP25004000");
697 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000");
698 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000fromStrange");
699 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000fromMat");
700 TH1F *fHistd0rphiITSMIoneSPDInAccP40008000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP40008000");
701 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000MC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000");
702 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000fromStrange");
703 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000fromMat"); */
704
705 TH1F *fHistd0rphiITSMIoneSPDInAcc150200MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP150200MC->Clone("fHistd0rphiITSMIoneSPDInAcc150200MC");
706 fHistd0rphiITSMIoneSPDInAcc150200MC->Add(fHistd0rphiITSMIoneSPDInAccS150200MC);
707 fHistd0rphiITSMIoneSPDInAcc150200MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc150200MC->GetEntries());
708 /* TH1F *fHistd0rphiITSMIoneSPDInAcc350450MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP350450MC->Clone("fHistd0rphiITSMIoneSPDInAcc350450MC");
709 fHistd0rphiITSMIoneSPDInAcc350450MC->Add(fHistd0rphiITSMIoneSPDInAccS350450MC);
710 fHistd0rphiITSMIoneSPDInAcc350450MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc350450MC->GetEntries());*/
711 TH1F *fHistd0rphiITSMIoneSPDInAcc500700MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP500700MC->Clone("fHistd0rphiITSMIoneSPDInAcc500700MC");
712 fHistd0rphiITSMIoneSPDInAcc500700MC->Add(fHistd0rphiITSMIoneSPDInAccS500700MC);
713 fHistd0rphiITSMIoneSPDInAcc500700MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc500700MC->GetEntries());
714 TH1F *fHistd0rphiITSMIoneSPDInAcc10001500MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP10001500MC->Clone("fHistd0rphiITSMIoneSPDInAcc10001500MC");
715 fHistd0rphiITSMIoneSPDInAcc10001500MC->Add(fHistd0rphiITSMIoneSPDInAccS10001500MC);
716 fHistd0rphiITSMIoneSPDInAcc10001500MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc10001500MC->GetEntries());
717 /* TH1F *fHistd0rphiITSMIoneSPDInAcc25004000MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP25004000MC->Clone("fHistd0rphiITSMIoneSPDInAcc25004000MC");
718 fHistd0rphiITSMIoneSPDInAcc25004000MC->Add(fHistd0rphiITSMIoneSPDInAccS25004000MC);
719 fHistd0rphiITSMIoneSPDInAcc25004000MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc25004000MC->GetEntries());
720 TH1F *fHistd0rphiITSMIoneSPDInAcc40008000MC=(TH1F*)fHistd0rphiITSMIoneSPDInAccP40008000MC->Clone("fHistd0rphiITSMIoneSPDInAcc40008000MC");
721 fHistd0rphiITSMIoneSPDInAcc40008000MC->Add(fHistd0rphiITSMIoneSPDInAccS40008000MC);
722 fHistd0rphiITSMIoneSPDInAcc40008000MC->Scale(1./fHistd0rphiITSMIoneSPDInAcc40008000MC->GetEntries());
723 */
724
725 TFile *f= new TFile("ITS.Performance_104065_104892_pass4_trkvtx_noTPCdca.root");
726
727 TList *list=(TList*)f->Get("cOutputITS");
728 TH1F *fHistd0rphiITSMIoneSPDInAccP150200 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
729 TH1F *fHistd0rphiITSMIoneSPDInAccS150200 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
730 /*TH1F *fHistd0rphiITSMIoneSPDInAccP350450 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP350450");
731 TH1F *fHistd0rphiITSMIoneSPDInAccS350450 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS350450");*/
732 TH1F *fHistd0rphiITSMIoneSPDInAccP500700 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP500700");
733 TH1F *fHistd0rphiITSMIoneSPDInAccS500700 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS500700");
734 TH1F *fHistd0rphiITSMIoneSPDInAccP10001500 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP10001500");
735 TH1F *fHistd0rphiITSMIoneSPDInAccS10001500 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS10001500");
736 /*TH1F *fHistd0rphiITSMIoneSPDInAccP25004000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP25004000");
737 TH1F *fHistd0rphiITSMIoneSPDInAccS25004000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS25004000");
738 TH1F *fHistd0rphiITSMIoneSPDInAccP40008000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP40008000");
739 TH1F *fHistd0rphiITSMIoneSPDInAccS40008000 = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS40008000");*/
740 TH1F *fHistd0rphiITSMIoneSPDInAcc150200=(TH1F*)fHistd0rphiITSMIoneSPDInAccP150200->Clone("fHistd0rphiITSMIoneSPDInAcc150200");
741 fHistd0rphiITSMIoneSPDInAcc150200->Add(fHistd0rphiITSMIoneSPDInAccS150200);
742 fHistd0rphiITSMIoneSPDInAcc150200->Scale(1./fHistd0rphiITSMIoneSPDInAcc150200->GetEntries());
743 /*TH1F *fHistd0rphiITSMIoneSPDInAcc350450=(TH1F*)fHistd0rphiITSMIoneSPDInAccP350450->Clone("fHistd0rphiITSMIoneSPDInAcc350450");
744 fHistd0rphiITSMIoneSPDInAcc350450->Add(fHistd0rphiITSMIoneSPDInAccS350450);
745 fHistd0rphiITSMIoneSPDInAcc350450->Scale(1./fHistd0rphiITSMIoneSPDInAcc350450->GetEntries());*/
746 TH1F *fHistd0rphiITSMIoneSPDInAcc500700=(TH1F*)fHistd0rphiITSMIoneSPDInAccP500700->Clone("fHistd0rphiITSMIoneSPDInAcc500700");
747 fHistd0rphiITSMIoneSPDInAcc500700->Add(fHistd0rphiITSMIoneSPDInAccS500700);
748 fHistd0rphiITSMIoneSPDInAcc500700->Scale(1./fHistd0rphiITSMIoneSPDInAcc500700->GetEntries());
749 TH1F *fHistd0rphiITSMIoneSPDInAcc10001500=(TH1F*)fHistd0rphiITSMIoneSPDInAccP10001500->Clone("fHistd0rphiITSMIoneSPDInAcc10001500");
750 fHistd0rphiITSMIoneSPDInAcc10001500->Add(fHistd0rphiITSMIoneSPDInAccS10001500);
751 fHistd0rphiITSMIoneSPDInAcc10001500->Scale(1./fHistd0rphiITSMIoneSPDInAcc10001500->GetEntries());
752 /*TH1F *fHistd0rphiITSMIoneSPDInAcc25004000=(TH1F*)fHistd0rphiITSMIoneSPDInAccP25004000->Clone("fHistd0rphiITSMIoneSPDInAcc25004000");
753 fHistd0rphiITSMIoneSPDInAcc25004000->Add(fHistd0rphiITSMIoneSPDInAccS25004000);
754 fHistd0rphiITSMIoneSPDInAcc25004000->Scale(1./fHistd0rphiITSMIoneSPDInAcc25004000->GetEntries());
755 TH1F *fHistd0rphiITSMIoneSPDInAcc40008000=(TH1F*)fHistd0rphiITSMIoneSPDInAccP40008000->Clone("fHistd0rphiITSMIoneSPDInAcc40008000");
756 fHistd0rphiITSMIoneSPDInAcc40008000->Add(fHistd0rphiITSMIoneSPDInAccS40008000);
757 fHistd0rphiITSMIoneSPDInAcc40008000->Scale(1./fHistd0rphiITSMIoneSPDInAcc40008000->GetEntries()); */
758
759
760 TCanvas *c1 = new TCanvas("c1","c1");
761 c1->Divide(3,2);
762 c1->cd(1);
763 fHistd0rphiITSMIoneSPDInAcc150200MC->Rebin(rebin);
764 fHistd0rphiITSMIoneSPDInAcc150200->Rebin(rebin);
765 fHistd0rphiITSMIoneSPDInAcc150200MC->SetLineColor(2);
766 fHistd0rphiITSMIoneSPDInAcc150200MC->Draw();
767 fHistd0rphiITSMIoneSPDInAcc150200->SetLineColor(4);
768 fHistd0rphiITSMIoneSPDInAcc150200->Draw("same");
769 /* c1->cd(2);
770 fHistd0rphiITSMIoneSPDInAcc350450MC->Rebin(rebin);
771 fHistd0rphiITSMIoneSPDInAcc350450->Rebin(rebin);
772 fHistd0rphiITSMIoneSPDInAcc350450MC->SetLineColor(2);
773 fHistd0rphiITSMIoneSPDInAcc350450MC->Draw();
774 fHistd0rphiITSMIoneSPDInAcc350450->SetLineColor(4);
775 fHistd0rphiITSMIoneSPDInAcc350450->Draw("same"); */
776 c1->cd(3);
777 fHistd0rphiITSMIoneSPDInAcc500700MC->Rebin(rebin);
778 fHistd0rphiITSMIoneSPDInAcc500700->Rebin(rebin);
779 fHistd0rphiITSMIoneSPDInAcc500700MC->SetLineColor(2);
780 fHistd0rphiITSMIoneSPDInAcc500700MC->Draw();
781 fHistd0rphiITSMIoneSPDInAcc500700->SetLineColor(4);
782 fHistd0rphiITSMIoneSPDInAcc500700->Draw("same");
783 c1->cd(4);
784 fHistd0rphiITSMIoneSPDInAcc10001500MC->Rebin(rebin);
785 fHistd0rphiITSMIoneSPDInAcc10001500->Rebin(rebin);
786 fHistd0rphiITSMIoneSPDInAcc10001500MC->SetLineColor(2);
787 fHistd0rphiITSMIoneSPDInAcc10001500MC->Draw();
788 fHistd0rphiITSMIoneSPDInAcc10001500->SetLineColor(4);
789 fHistd0rphiITSMIoneSPDInAcc10001500->Draw("same");
790 /* c1->cd(5);
791 fHistd0rphiITSMIoneSPDInAcc25004000MC->Rebin(rebin);
792 fHistd0rphiITSMIoneSPDInAcc25004000->Rebin(rebin);
793 fHistd0rphiITSMIoneSPDInAcc25004000MC->SetLineColor(2);
794 fHistd0rphiITSMIoneSPDInAcc25004000MC->Draw();
795 fHistd0rphiITSMIoneSPDInAcc25004000->SetLineColor(4);
796 fHistd0rphiITSMIoneSPDInAcc25004000->Draw("same");
797 c1->cd(6);
798 fHistd0rphiITSMIoneSPDInAcc40008000MC->Rebin(rebin);
799 fHistd0rphiITSMIoneSPDInAcc40008000->Rebin(rebin);
800 fHistd0rphiITSMIoneSPDInAcc40008000MC->SetLineColor(2);
801 fHistd0rphiITSMIoneSPDInAcc40008000MC->Draw();
802 fHistd0rphiITSMIoneSPDInAcc40008000->SetLineColor(4);
803 fHistd0rphiITSMIoneSPDInAcc40008000->Draw("same");
804 */
805
806 TCanvas *c1a = new TCanvas("c1a","c1a");
807 c1a->Divide(3,1);
808 c1a->cd(1);
809 TH1F* fHistd0rphiITSMIoneSPDInAcc150200Ratio=(TH1F*)fHistd0rphiITSMIoneSPDInAcc150200->Clone("fHistd0rphiITSMIoneSPDInAcc150200Ratio");
810 fHistd0rphiITSMIoneSPDInAcc150200Ratio->Divide(fHistd0rphiITSMIoneSPDInAcc150200MC);
811 fHistd0rphiITSMIoneSPDInAcc150200Ratio->Draw();
812 c1a->cd(2);
813 TH1F* fHistd0rphiITSMIoneSPDInAcc500700Ratio=(TH1F*)fHistd0rphiITSMIoneSPDInAcc500700->Clone("fHistd0rphiITSMIoneSPDInAcc500700Ratio");
814 fHistd0rphiITSMIoneSPDInAcc500700Ratio->Divide(fHistd0rphiITSMIoneSPDInAcc500700MC);
815 fHistd0rphiITSMIoneSPDInAcc500700Ratio->Draw();
816 c1a->cd(3);
817 TH1F* fHistd0rphiITSMIoneSPDInAcc10001500Ratio=(TH1F*)fHistd0rphiITSMIoneSPDInAcc10001500->Clone("fHistd0rphiITSMIoneSPDInAcc10001500Ratio");
818 fHistd0rphiITSMIoneSPDInAcc10001500Ratio->Divide(fHistd0rphiITSMIoneSPDInAcc10001500MC);
819 fHistd0rphiITSMIoneSPDInAcc10001500Ratio->Draw();
820
821
822
823 TCanvas *c2 = new TCanvas("c2","c2");
824 c2->Divide(3,2);
825 c2_1->SetLogy();
826 c2_2->SetLogy();
827 c2_3->SetLogy();
828 c2_4->SetLogy();
829 c2_5->SetLogy();
830 c2_6->SetLogy();
831 c2->cd(1);
832 fHistd0rphiITSMIoneSPDInAccP150200MC->SetMinimum(1);
833 fHistd0rphiITSMIoneSPDInAccS150200MC->SetMinimum(1);
834 fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->SetMinimum(1);
835 fHistd0rphiITSMIoneSPDInAccS150200fromMatMC->SetMinimum(1);
836 fHistd0rphiITSMIoneSPDInAccP150200MC->SetLineColor(1);
837 fHistd0rphiITSMIoneSPDInAccP150200MC->Draw();
838 fHistd0rphiITSMIoneSPDInAccS150200MC->SetLineColor(6);
839 fHistd0rphiITSMIoneSPDInAccS150200MC->Draw("same");
840 fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->SetLineColor(8);
841 fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->Draw("same");
842 fHistd0rphiITSMIoneSPDInAccS150200fromMatMC->SetLineColor(9);
843 fHistd0rphiITSMIoneSPDInAccS150200fromMatMC->Draw("same");
844 /* c2->cd(2);
845 fHistd0rphiITSMIoneSPDInAccP350450MC->SetMinimum(1);
846 fHistd0rphiITSMIoneSPDInAccS350450MC->SetMinimum(1);
847 fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->SetMinimum(1);
848 fHistd0rphiITSMIoneSPDInAccS350450fromMatMC->SetMinimum(1);
849 fHistd0rphiITSMIoneSPDInAccP350450MC->SetLineColor(1);
850 fHistd0rphiITSMIoneSPDInAccP350450MC->Draw();
851 fHistd0rphiITSMIoneSPDInAccS350450MC->SetLineColor(6);
852 fHistd0rphiITSMIoneSPDInAccS350450MC->Draw("same");
853 fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->SetLineColor(8);
854 fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->Draw("same");
855 fHistd0rphiITSMIoneSPDInAccS350450fromMatMC->SetLineColor(9);
856 fHistd0rphiITSMIoneSPDInAccS350450fromMatMC->Draw("same");
857*/
858 c2->cd(3);
859 fHistd0rphiITSMIoneSPDInAccP500700MC->SetMinimum(1);
860 fHistd0rphiITSMIoneSPDInAccS500700MC->SetMinimum(1);
861 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetMinimum(1);
862 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->SetMinimum(1);
863 fHistd0rphiITSMIoneSPDInAccP500700MC->SetLineColor(1);
864 fHistd0rphiITSMIoneSPDInAccP500700MC->Draw();
865 fHistd0rphiITSMIoneSPDInAccS500700MC->SetLineColor(6);
866 fHistd0rphiITSMIoneSPDInAccS500700MC->Draw("same");
867 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetLineColor(8);
868 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->Draw("same");
869 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->SetLineColor(9);
870 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->Draw("same");
871 c2->cd(4);
872 fHistd0rphiITSMIoneSPDInAccP10001500MC->SetMinimum(1);
873 fHistd0rphiITSMIoneSPDInAccS10001500MC->SetMinimum(1);
874 fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->SetMinimum(1);
875 fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC->SetMinimum(1);
876 fHistd0rphiITSMIoneSPDInAccP10001500MC->SetLineColor(1);
877 fHistd0rphiITSMIoneSPDInAccP10001500MC->Draw();
878 fHistd0rphiITSMIoneSPDInAccS10001500MC->SetLineColor(6);
879 fHistd0rphiITSMIoneSPDInAccS10001500MC->Draw("same");
880 fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->SetLineColor(8);
881 fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->Draw("same");
882 fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC->SetLineColor(9);
883 fHistd0rphiITSMIoneSPDInAccS10001500fromMatMC->Draw("same");
884 /*
885 c2->cd(5);
886 fHistd0rphiITSMIoneSPDInAccP25004000MC->SetMinimum(1);
887 fHistd0rphiITSMIoneSPDInAccS25004000MC->SetMinimum(1);
888 fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->SetMinimum(1);
889 fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC->SetMinimum(1);
890 fHistd0rphiITSMIoneSPDInAccP25004000MC->SetLineColor(1);
891 fHistd0rphiITSMIoneSPDInAccP25004000MC->Draw();
892 fHistd0rphiITSMIoneSPDInAccS25004000MC->SetLineColor(6);
893 fHistd0rphiITSMIoneSPDInAccS25004000MC->Draw("same");
894 fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->SetLineColor(8);
895 fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->Draw("same");
896 fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC->SetLineColor(9);
897 fHistd0rphiITSMIoneSPDInAccS25004000fromMatMC->Draw("same");
898 c2->cd(6);
899 fHistd0rphiITSMIoneSPDInAccP40008000MC->SetMinimum(1);
900 fHistd0rphiITSMIoneSPDInAccS40008000MC->SetMinimum(1);
901 fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->SetMinimum(1);
902 fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC->SetMinimum(1);
903 fHistd0rphiITSMIoneSPDInAccP40008000MC->SetLineColor(1);
904 fHistd0rphiITSMIoneSPDInAccP40008000MC->Draw();
905 fHistd0rphiITSMIoneSPDInAccS40008000MC->SetLineColor(6);
906 fHistd0rphiITSMIoneSPDInAccS40008000MC->Draw("same");
907 fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->SetLineColor(8);
908 fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->Draw("same");
909 fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC->SetLineColor(9);
910 fHistd0rphiITSMIoneSPDInAccS40008000fromMatMC->Draw("same");
911*/
912 TCanvas *c2b = new TCanvas("c2b","c2b");
913 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetMinimum(1);
914 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->SetLineColor(8);
915 fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->Draw();
916 fHistd0rphiITSMIoneSPDInAccS500700from310MC->SetMinimum(1);
917 fHistd0rphiITSMIoneSPDInAccS500700from310MC->SetLineColor(2);
918 fHistd0rphiITSMIoneSPDInAccS500700from310MC->Draw("same");
919 fHistd0rphiITSMIoneSPDInAccS500700from321MC->SetMinimum(1);
920 fHistd0rphiITSMIoneSPDInAccS500700from321MC->SetLineColor(5);
921 fHistd0rphiITSMIoneSPDInAccS500700from321MC->Draw("same");
922 fHistd0rphiITSMIoneSPDInAccS500700from3122MC->SetMinimum(1);
923 fHistd0rphiITSMIoneSPDInAccS500700from3122MC->SetLineColor(1);
924 fHistd0rphiITSMIoneSPDInAccS500700from3122MC->Draw("same");
925
926 TCanvas *c2c = new TCanvas("c2c","c2c");
927 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->SetMinimum(1);
928 fHistd0rphiITSMIoneSPDInAccS500700fromMatMC->Draw();
929 fHistd0rphiITSMIoneSPDInAccS500700from211MC->SetMinimum(1);
930 fHistd0rphiITSMIoneSPDInAccS500700from211MC->SetLineColor(2);
931 fHistd0rphiITSMIoneSPDInAccS500700from211MC->Draw("same");
932 fHistd0rphiITSMIoneSPDInAccS500700from22MC->SetMinimum(1);
933 fHistd0rphiITSMIoneSPDInAccS500700from22MC->SetLineColor(1);
934 fHistd0rphiITSMIoneSPDInAccS500700from22MC->Draw("same");
935
936
937
938 TCanvas *c3 = new TCanvas("c3","c3");
939 c3->Divide(3,2);
940 c3_1->SetLogx();
941 c3_2->SetLogx();
942 c3_3->SetLogx();
943 c3_4->SetLogx();
944 c3_5->SetLogx();
945 c3_6->SetLogx();
946 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};
947 /*
948 c3->cd(1);
949 Float_t fracP150200[15],fracS150200[15],fracSfromStrange150200[15],fracS150200Strp30[15],fracS150200Strm30[15],fracS150200Matm10[15],fracS150200Matp10[15];
950 Float_t intPtot150200=fHistd0rphiITSMIoneSPDInAccP150200MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP150200MC->GetNbinsX());
951 for(Int_t i=0;i<15;i++) {
952 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP150200MC->FindBin(-d0cut[i]);
953 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP150200MC->FindBin(+d0cut[i]);
954 Float_t intPcut150200=fHistd0rphiITSMIoneSPDInAccP150200MC->Integral(bin1,bin2);
955 Float_t intScut150200=fHistd0rphiITSMIoneSPDInAccS150200MC->Integral(bin1,bin2);
956 Float_t intSfromStrangecut150200=fHistd0rphiITSMIoneSPDInAccS150200fromStrangeMC->Integral(bin1,bin2);
957 Float_t intScut150200Strp30 = intScut150200 + 0.3*intSfromStrangecut150200;
958 Float_t intScut150200Strm30 = intScut150200 - 0.3*intSfromStrangecut150200;
959 Float_t intScut150200Matp10 = intScut150200 + 0.1*(intScut150200-intSfromStrangecut150200);
960 Float_t intScut150200Matm10 = intScut150200 - 0.1*(intScut150200-intSfromStrangecut150200);
961 fracP150200[i]=intPcut150200/intPtot150200;
962 fracS150200[i]=1.-intScut150200/(intPcut150200+intScut150200);
963 fracS150200Strp30[i]=1.-intScut150200Strp30/(intPcut150200+intScut150200Strp30);
964 fracS150200Strm30[i]=1.-intScut150200Strm30/(intPcut150200+intScut150200Strm30);
965 fracS150200Matp10[i]=1.-intScut150200Matp10/(intPcut150200+intScut150200Matp10);
966 fracS150200Matm10[i]=1.-intScut150200Matm10/(intPcut150200+intScut150200Matm10);
967 fracSfromStrange150200[i]=1.-intSfromStrangecut150200/(intPcut150200+intScut150200);
968 }
969 TGraph *gfracP150200=new TGraph(15,d0cut,fracP150200);
970 gfracP150200->SetMarkerColor(2);
971 gfracP150200->SetMarkerStyle(20);
972 gfracP150200->Draw("ap");
973 TGraph *gfracS150200=new TGraph(15,d0cut,fracS150200);
974 gfracS150200->SetMarkerColor(4);
975 gfracS150200->SetMarkerStyle(21);
976 gfracS150200->Draw("p");
977 TGraph *gfracSfromStrange150200=new TGraph(15,d0cut,fracSfromStrange150200);
978 gfracSfromStrange150200->SetMarkerColor(8);
979 gfracSfromStrange150200->SetMarkerStyle(22);
980 gfracSfromStrange150200->Draw("p"); */
981 /*
982 c3->cd(2);
983 Float_t fracP350450[15],fracS350450[15],fracSfromStrange350450[15],fracS350450Strp30[15],fracS350450Strm30[15],fracS350450Matm10[15],fracS350450Matp10[15];
984 Float_t intPtot350450=fHistd0rphiITSMIoneSPDInAccP350450MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP350450MC->GetNbinsX());
985 for(Int_t i=0;i<15;i++) {
986 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP350450MC->FindBin(-d0cut[i]);
987 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP350450MC->FindBin(+d0cut[i]);
988 Float_t intPcut350450=fHistd0rphiITSMIoneSPDInAccP350450MC->Integral(bin1,bin2);
989 Float_t intScut350450=fHistd0rphiITSMIoneSPDInAccS350450MC->Integral(bin1,bin2);
990 Float_t intSfromStrangecut350450=fHistd0rphiITSMIoneSPDInAccS350450fromStrangeMC->Integral(bin1,bin2);
991 Float_t intScut350450Strp30 = intScut350450 + 0.3*intSfromStrangecut350450;
992 Float_t intScut350450Strm30 = intScut350450 - 0.3*intSfromStrangecut350450;
993 Float_t intScut350450Matp10 = intScut350450 + 0.1*(intScut350450-intSfromStrangecut350450);
994 Float_t intScut350450Matm10 = intScut350450 - 0.1*(intScut350450-intSfromStrangecut350450);
995 fracP350450[i]=intPcut350450/intPtot350450;
996 fracS350450[i]=1.-intScut350450/(intPcut350450+intScut350450);
997 fracS350450Strp30[i]=1.-intScut350450Strp30/(intPcut350450+intScut350450Strp30);
998 fracS350450Strm30[i]=1.-intScut350450Strm30/(intPcut350450+intScut350450Strm30);
999 fracS350450Matp10[i]=1.-intScut350450Matp10/(intPcut350450+intScut350450Matp10);
1000 fracS350450Matm10[i]=1.-intScut350450Matm10/(intPcut350450+intScut350450Matm10);
1001 fracSfromStrange350450[i]=1.-intSfromStrangecut350450/(intPcut350450+intScut350450);
1002 }
1003 TGraph *gfracP350450=new TGraph(15,d0cut,fracP350450);
1004 gfracP350450->SetMarkerColor(2);
1005 gfracP350450->SetMarkerStyle(20);
1006 gfracP350450->Draw("ap");
1007 TGraph *gfracS350450=new TGraph(15,d0cut,fracS350450);
1008 gfracS350450->SetMarkerColor(4);
1009 gfracS350450->SetMarkerStyle(21);
1010 gfracS350450->Draw("p");
1011 TGraph *gfracSfromStrange350450=new TGraph(15,d0cut,fracSfromStrange350450);
1012 gfracSfromStrange350450->SetMarkerColor(8);
1013 gfracSfromStrange350450->SetMarkerStyle(22);
1014 gfracSfromStrange350450->Draw("p");
1015*/
1016 c3->cd(3);
1017 Float_t fracP500700[15],fracS500700[15],fracSfromStrange500700[15],fracS500700Strp30[15],fracS500700Strm30[15],fracS500700Matm10[15],fracS500700Matp10[15];
1018 Float_t intPtot500700=fHistd0rphiITSMIoneSPDInAccP500700MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP500700MC->GetNbinsX());
1019 for(Int_t i=0;i<15;i++) {
1020 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP500700MC->FindBin(-d0cut[i]);
1021 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP500700MC->FindBin(+d0cut[i]);
1022 Float_t intPcut500700=fHistd0rphiITSMIoneSPDInAccP500700MC->Integral(bin1,bin2);
1023 Float_t intScut500700=fHistd0rphiITSMIoneSPDInAccS500700MC->Integral(bin1,bin2);
1024 Float_t intSfromStrangecut500700=fHistd0rphiITSMIoneSPDInAccS500700fromStrangeMC->Integral(bin1,bin2);
1025 Float_t intScut500700Strp30 = intScut500700 + 0.3*intSfromStrangecut500700;
1026 Float_t intScut500700Strm30 = intScut500700 - 0.3*intSfromStrangecut500700;
1027 Float_t intScut500700Matp10 = intScut500700 + 0.1*(intScut500700-intSfromStrangecut500700);
1028 Float_t intScut500700Matm10 = intScut500700 - 0.1*(intScut500700-intSfromStrangecut500700);
1029 fracP500700[i]=intPcut500700/intPtot500700;
1030 fracS500700[i]=1.-intScut500700/(intPcut500700+intScut500700);
1031 fracS500700Strp30[i]=1.-intScut500700Strp30/(intPcut500700+intScut500700Strp30);
1032 fracS500700Strm30[i]=1.-intScut500700Strm30/(intPcut500700+intScut500700Strm30);
1033 fracS500700Matp10[i]=1.-intScut500700Matp10/(intPcut500700+intScut500700Matp10);
1034 fracS500700Matm10[i]=1.-intScut500700Matm10/(intPcut500700+intScut500700Matm10);
1035 fracSfromStrange500700[i]=1.-intSfromStrangecut500700/(intPcut500700+intScut500700);
1036 }
1037 TGraph *gfracP500700=new TGraph(15,d0cut,fracP500700);
1038 gfracP500700->SetMarkerColor(2);
1039 gfracP500700->SetMarkerStyle(20);
1040 gfracP500700->Draw("ap");
1041 TGraph *gfracS500700=new TGraph(15,d0cut,fracS500700);
1042 gfracS500700->SetMarkerColor(4);
1043 gfracS500700->SetMarkerStyle(21);
1044 gfracS500700->Draw("p");
1045 TGraph *gfracSfromStrange500700=new TGraph(15,d0cut,fracSfromStrange500700);
1046 gfracSfromStrange500700->SetMarkerColor(8);
1047 gfracSfromStrange500700->SetMarkerStyle(22);
1048 gfracSfromStrange500700->Draw("p");
1049
1050 c3->cd(4);
1051 Float_t fracP10001500[15],fracS10001500[15],fracSfromStrange10001500[15],fracS10001500Strp30[15],fracS10001500Strm30[15],fracS10001500Matm10[15],fracS10001500Matp10[15];
1052 Float_t intPtot10001500=fHistd0rphiITSMIoneSPDInAccP10001500MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP10001500MC->GetNbinsX());
1053 for(Int_t i=0;i<15;i++) {
1054 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP10001500MC->FindBin(-d0cut[i]);
1055 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP10001500MC->FindBin(+d0cut[i]);
1056 Float_t intPcut10001500=fHistd0rphiITSMIoneSPDInAccP10001500MC->Integral(bin1,bin2);
1057 Float_t intScut10001500=fHistd0rphiITSMIoneSPDInAccS10001500MC->Integral(bin1,bin2);
1058 Float_t intSfromStrangecut10001500=fHistd0rphiITSMIoneSPDInAccS10001500fromStrangeMC->Integral(bin1,bin2);
1059 Float_t intScut10001500Strp30 = intScut10001500 + 0.3*intSfromStrangecut10001500;
1060 Float_t intScut10001500Strm30 = intScut10001500 - 0.3*intSfromStrangecut10001500;
1061 Float_t intScut10001500Matp10 = intScut10001500 + 0.1*(intScut10001500-intSfromStrangecut10001500);
1062 Float_t intScut10001500Matm10 = intScut10001500 - 0.1*(intScut10001500-intSfromStrangecut10001500);
1063 fracP10001500[i]=intPcut10001500/intPtot10001500;
1064 fracS10001500[i]=1.-intScut10001500/(intPcut10001500+intScut10001500);
1065 fracS10001500Strp30[i]=1.-intScut10001500Strp30/(intPcut10001500+intScut10001500Strp30);
1066 fracS10001500Strm30[i]=1.-intScut10001500Strm30/(intPcut10001500+intScut10001500Strm30);
1067 fracS10001500Matp10[i]=1.-intScut10001500Matp10/(intPcut10001500+intScut10001500Matp10);
1068 fracS10001500Matm10[i]=1.-intScut10001500Matm10/(intPcut10001500+intScut10001500Matm10);
1069 fracSfromStrange10001500[i]=1.-intSfromStrangecut10001500/(intPcut10001500+intScut10001500);
1070 }
1071 TGraph *gfracP10001500=new TGraph(15,d0cut,fracP10001500);
1072 gfracP10001500->SetMarkerColor(2);
1073 gfracP10001500->SetMarkerStyle(20);
1074 gfracP10001500->Draw("ap");
1075 TGraph *gfracS10001500=new TGraph(15,d0cut,fracS10001500);
1076 gfracS10001500->SetMarkerColor(4);
1077 gfracS10001500->SetMarkerStyle(21);
1078 gfracS10001500->Draw("p");
1079 TGraph *gfracSfromStrange10001500=new TGraph(15,d0cut,fracSfromStrange10001500);
1080 gfracSfromStrange10001500->SetMarkerColor(8);
1081 gfracSfromStrange10001500->SetMarkerStyle(22);
1082 gfracSfromStrange10001500->Draw("p");
1083 /*
1084 c3->cd(5);
1085 Float_t fracP25004000[15],fracS25004000[15],fracSfromStrange25004000[15],fracS25004000Strp30[15],fracS25004000Strm30[15],fracS25004000Matm10[15],fracS25004000Matp10[15];
1086 Float_t intPtot25004000=fHistd0rphiITSMIoneSPDInAccP25004000MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP25004000MC->GetNbinsX());
1087 for(Int_t i=0;i<15;i++) {
1088 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP25004000MC->FindBin(-d0cut[i]);
1089 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP25004000MC->FindBin(+d0cut[i]);
1090 Float_t intPcut25004000=fHistd0rphiITSMIoneSPDInAccP25004000MC->Integral(bin1,bin2);
1091 Float_t intScut25004000=fHistd0rphiITSMIoneSPDInAccS25004000MC->Integral(bin1,bin2);
1092 Float_t intSfromStrangecut25004000=fHistd0rphiITSMIoneSPDInAccS25004000fromStrangeMC->Integral(bin1,bin2);
1093 Float_t intScut25004000Strp30 = intScut25004000 + 0.3*intSfromStrangecut25004000;
1094 Float_t intScut25004000Strm30 = intScut25004000 - 0.3*intSfromStrangecut25004000;
1095 Float_t intScut25004000Matp10 = intScut25004000 + 0.1*(intScut25004000-intSfromStrangecut25004000);
1096 Float_t intScut25004000Matm10 = intScut25004000 - 0.1*(intScut25004000-intSfromStrangecut25004000);
1097 fracP25004000[i]=intPcut25004000/intPtot25004000;
1098 fracS25004000[i]=1.-intScut25004000/(intPcut25004000+intScut25004000);
1099 fracS25004000Strp30[i]=1.-intScut25004000Strp30/(intPcut25004000+intScut25004000Strp30);
1100 fracS25004000Strm30[i]=1.-intScut25004000Strm30/(intPcut25004000+intScut25004000Strm30);
1101 fracS25004000Matp10[i]=1.-intScut25004000Matp10/(intPcut25004000+intScut25004000Matp10);
1102 fracS25004000Matm10[i]=1.-intScut25004000Matm10/(intPcut25004000+intScut25004000Matm10);
1103 fracSfromStrange25004000[i]=1.-intSfromStrangecut25004000/(intPcut25004000+intScut25004000);
1104 }
1105 TGraph *gfracP25004000=new TGraph(15,d0cut,fracP25004000);
1106 gfracP25004000->SetMarkerColor(2);
1107 gfracP25004000->SetMarkerStyle(20);
1108 gfracP25004000->Draw("ap");
1109 TGraph *gfracS25004000=new TGraph(15,d0cut,fracS25004000);
1110 gfracS25004000->SetMarkerColor(4);
1111 gfracS25004000->SetMarkerStyle(21);
1112 gfracS25004000->Draw("p");
1113 TGraph *gfracSfromStrange25004000=new TGraph(15,d0cut,fracSfromStrange25004000);
1114 gfracSfromStrange25004000->SetMarkerColor(8);
1115 gfracSfromStrange25004000->SetMarkerStyle(22);
1116 gfracSfromStrange25004000->Draw("p");
1117
1118 c3->cd(6);
1119 Float_t fracP40008000[15],fracS40008000[15],fracSfromStrange40008000[15],fracS40008000Strp30[15],fracS40008000Strm30[15],fracS40008000Matm10[15],fracS40008000Matp10[15];
1120 Float_t intPtot40008000=fHistd0rphiITSMIoneSPDInAccP40008000MC->Integral(1,fHistd0rphiITSMIoneSPDInAccP40008000MC->GetNbinsX());
1121 for(Int_t i=0;i<15;i++) {
1122 Int_t bin1=fHistd0rphiITSMIoneSPDInAccP40008000MC->FindBin(-d0cut[i]);
1123 Int_t bin2=fHistd0rphiITSMIoneSPDInAccP40008000MC->FindBin(+d0cut[i]);
1124 Float_t intPcut40008000=fHistd0rphiITSMIoneSPDInAccP40008000MC->Integral(bin1,bin2);
1125 Float_t intScut40008000=fHistd0rphiITSMIoneSPDInAccS40008000MC->Integral(bin1,bin2);
1126 Float_t intSfromStrangecut40008000=fHistd0rphiITSMIoneSPDInAccS40008000fromStrangeMC->Integral(bin1,bin2);
1127 Float_t intScut40008000Strp30 = intScut40008000 + 0.3*intSfromStrangecut40008000;
1128 Float_t intScut40008000Strm30 = intScut40008000 - 0.3*intSfromStrangecut40008000;
1129 Float_t intScut40008000Matp10 = intScut40008000 + 0.1*(intScut40008000-intSfromStrangecut40008000);
1130 Float_t intScut40008000Matm10 = intScut40008000 - 0.1*(intScut40008000-intSfromStrangecut40008000);
1131 fracP40008000[i]=intPcut40008000/intPtot40008000;
1132 fracS40008000[i]=1.-intScut40008000/(intPcut40008000+intScut40008000);
1133 fracS40008000Strp30[i]=1.-intScut40008000Strp30/(intPcut40008000+intScut40008000Strp30);
1134 fracS40008000Strm30[i]=1.-intScut40008000Strm30/(intPcut40008000+intScut40008000Strm30);
1135 fracS40008000Matp10[i]=1.-intScut40008000Matp10/(intPcut40008000+intScut40008000Matp10);
1136 fracS40008000Matm10[i]=1.-intScut40008000Matm10/(intPcut40008000+intScut40008000Matm10);
1137 fracSfromStrange40008000[i]=1.-intSfromStrangecut40008000/(intPcut40008000+intScut40008000);
1138 }
1139 TGraph *gfracP40008000=new TGraph(15,d0cut,fracP40008000);
1140 gfracP40008000->SetMarkerColor(2);
1141 gfracP40008000->SetMarkerStyle(20);
1142 gfracP40008000->Draw("ap");
1143 TGraph *gfracS40008000=new TGraph(15,d0cut,fracS40008000);
1144 gfracS40008000->SetMarkerColor(4);
1145 gfracS40008000->SetMarkerStyle(21);
1146 gfracS40008000->Draw("p");
1147 TGraph *gfracSfromStrange40008000=new TGraph(15,d0cut,fracSfromStrange40008000);
1148 gfracSfromStrange40008000->SetMarkerColor(8);
1149 gfracSfromStrange40008000->SetMarkerStyle(22);
1150 gfracSfromStrange40008000->Draw("p");
1151*/
1152
1153 /*
1154 TCanvas *c4 = new TCanvas("c4","c4");
1155 c4->Divide(3,1);
1156 c4_1->SetLogx();
1157 c4_2->SetLogx();
1158 c4_3->SetLogx();
1159
1160 c4->cd(1);
1161 Float_t intDatacut150200[15],intDataPcut150200[15],intDataPall150200[15];
1162 for(Int_t i=0;i<15;i++) {
1163 Int_t bin1=fHistd0rphiITSMIoneSPDInAcc150200->FindBin(-d0cut[i]);
1164 Int_t bin2=fHistd0rphiITSMIoneSPDInAcc150200->FindBin(+d0cut[i]);
1165 intDatacut150200[i]=fHistd0rphiITSMIoneSPDInAcc150200->Integral(bin1,bin2);
1166 intDataPcut150200[i]=intDatacut150200[i]*fracS150200[i];
1167 intDataPall150200[i]=intDataPcut150200[i]/fracP150200[i];
1168 }
1169 TGraph *gintDatacut150200=new TGraph(15,d0cut,intDatacut150200);
1170 gintDatacut150200->SetMarkerColor(1);
1171 gintDatacut150200->SetMarkerStyle(20);
1172 gintDatacut150200->Draw("ap");
1173 TGraph *gintDataPcut150200=new TGraph(15,d0cut,intDataPcut150200);
1174 gintDataPcut150200->SetMarkerColor(1);
1175 gintDataPcut150200->SetMarkerStyle(24);
1176 gintDataPcut150200->Draw("p");
1177 TGraph *gintDataPall150200=new TGraph(15,d0cut,intDataPall150200);
1178 gintDataPall150200->SetMarkerColor(2);
1179 gintDataPall150200->SetMarkerStyle(22);
1180 gintDataPall150200->Draw("p");
1181
1182 c4->cd(2);
1183 Float_t intDatacut500700[15],intDataPcut500700[15],intDataPall500700[15];
1184 for(Int_t i=0;i<15;i++) {
1185 Int_t bin1=fHistd0rphiITSMIoneSPDInAcc500700->FindBin(-d0cut[i]);
1186 Int_t bin2=fHistd0rphiITSMIoneSPDInAcc500700->FindBin(+d0cut[i]);
1187 intDatacut500700[i]=fHistd0rphiITSMIoneSPDInAcc500700->Integral(bin1,bin2);
1188 intDataPcut500700[i]=intDatacut500700[i]*fracS500700[i];
1189 intDataPall500700[i]=intDataPcut500700[i]/fracP500700[i];
1190 }
1191 TGraph *gintDatacut500700=new TGraph(15,d0cut,intDatacut500700);
1192 gintDatacut500700->SetMarkerColor(1);
1193 gintDatacut500700->SetMarkerStyle(20);
1194 gintDatacut500700->Draw("ap");
1195 TGraph *gintDataPcut500700=new TGraph(15,d0cut,intDataPcut500700);
1196 gintDataPcut500700->SetMarkerColor(1);
1197 gintDataPcut500700->SetMarkerStyle(24);
1198 gintDataPcut500700->Draw("p");
1199 TGraph *gintDataPall500700=new TGraph(15,d0cut,intDataPall500700);
1200 gintDataPall500700->SetMarkerColor(2);
1201 gintDataPall500700->SetMarkerStyle(22);
1202 gintDataPall500700->Draw("p");
1203
1204 c4->cd(3);
1205 Float_t intDatacut10001500[15],intDataPcut10001500[15],intDataPall10001500[15];
1206 for(Int_t i=0;i<15;i++) {
1207 Int_t bin1=fHistd0rphiITSMIoneSPDInAcc10001500->FindBin(-d0cut[i]);
1208 Int_t bin2=fHistd0rphiITSMIoneSPDInAcc10001500->FindBin(+d0cut[i]);
1209 intDatacut10001500[i]=fHistd0rphiITSMIoneSPDInAcc10001500->Integral(bin1,bin2);
1210 intDataPcut10001500[i]=intDatacut10001500[i]*fracS10001500[i];
1211 intDataPall10001500[i]=intDataPcut10001500[i]/fracP10001500[i];
1212 }
1213 TGraph *gintDatacut10001500=new TGraph(15,d0cut,intDatacut10001500);
1214 gintDatacut10001500->SetMarkerColor(1);
1215 gintDatacut10001500->SetMarkerStyle(20);
1216 gintDatacut10001500->Draw("ap");
1217 TGraph *gintDataPcut10001500=new TGraph(15,d0cut,intDataPcut10001500);
1218 gintDataPcut10001500->SetMarkerColor(1);
1219 gintDataPcut10001500->SetMarkerStyle(24);
1220 gintDataPcut10001500->Draw("p");
1221 TGraph *gintDataPall10001500=new TGraph(15,d0cut,intDataPall10001500);
1222 gintDataPall10001500->SetMarkerColor(2);
1223 gintDataPall10001500->SetMarkerStyle(22);
1224 gintDataPall10001500->Draw("p");
1225
1226
1227 TCanvas *c5 = new TCanvas("c5","c5");
1228 c5->Divide(3,1);
1229 c5_1->SetLogx();
1230 c5_2->SetLogx();
1231 c5_3->SetLogx();
1232
1233 c5->cd(1);
1234 Float_t intDataPall150200Strp30[15],intDataPall150200Strm30[15],intDataPall150200Matp10[15],intDataPall150200Matm10[15];
1235 for(Int_t i=0;i<15;i++) {
1236 intDataPall150200Strp30[i]=intDatacut150200[i]*fracS150200Strp30[i]/fracP150200[i];
1237 intDataPall150200Strm30[i]=intDatacut150200[i]*fracS150200Strm30[i]/fracP150200[i];
1238 intDataPall150200Matp10[i]=intDatacut150200[i]*fracS150200Matp10[i]/fracP150200[i];
1239 intDataPall150200Matm10[i]=intDatacut150200[i]*fracS150200Matm10[i]/fracP150200[i];
1240 }
1241 gintDataPall150200->Draw("ap");
1242
1243 TGraph *gintDataPall150200Strp30=new TGraph(15,d0cut,intDataPall150200Strp30);
1244 gintDataPall150200Strp30->SetMarkerColor(1);
1245 gintDataPall150200Strp30->SetMarkerStyle(22);
1246 gintDataPall150200Strp30->Draw("p");
1247 TGraph *gintDataPall150200Strm30=new TGraph(15,d0cut,intDataPall150200Strm30);
1248 gintDataPall150200Strm30->SetMarkerColor(1);
1249 gintDataPall150200Strm30->SetMarkerStyle(22);
1250 gintDataPall150200Strm30->Draw("p");
1251 TGraph *gintDataPall150200Matp10=new TGraph(15,d0cut,intDataPall150200Matp10);
1252 gintDataPall150200Matp10->SetMarkerColor(3);
1253 gintDataPall150200Matp10->SetMarkerStyle(22);
1254 gintDataPall150200Matp10->Draw("p");
1255 TGraph *gintDataPall150200Matm10=new TGraph(15,d0cut,intDataPall150200Matm10);
1256 gintDataPall150200Matm10->SetMarkerColor(3);
1257 gintDataPall150200Matm10->SetMarkerStyle(22);
1258 gintDataPall150200Matm10->Draw("p");
1259
1260 c5->cd(2);
1261 Float_t intDataPall500700Strp30[15],intDataPall500700Strm30[15],intDataPall500700Matp10[15],intDataPall500700Matm10[15];
1262 for(Int_t i=0;i<15;i++) {
1263 intDataPall500700Strp30[i]=intDatacut500700[i]*fracS500700Strp30[i]/fracP500700[i];
1264 intDataPall500700Strm30[i]=intDatacut500700[i]*fracS500700Strm30[i]/fracP500700[i];
1265 intDataPall500700Matp10[i]=intDatacut500700[i]*fracS500700Matp10[i]/fracP500700[i];
1266 intDataPall500700Matm10[i]=intDatacut500700[i]*fracS500700Matm10[i]/fracP500700[i];
1267 }
1268 gintDataPall500700->Draw("ap");
1269
1270 TGraph *gintDataPall500700Strp30=new TGraph(15,d0cut,intDataPall500700Strp30);
1271 gintDataPall500700Strp30->SetMarkerColor(1);
1272 gintDataPall500700Strp30->SetMarkerStyle(22);
1273 gintDataPall500700Strp30->Draw("p");
1274 TGraph *gintDataPall500700Strm30=new TGraph(15,d0cut,intDataPall500700Strm30);
1275 gintDataPall500700Strm30->SetMarkerColor(1);
1276 gintDataPall500700Strm30->SetMarkerStyle(22);
1277 gintDataPall500700Strm30->Draw("p");
1278 TGraph *gintDataPall500700Matp10=new TGraph(15,d0cut,intDataPall500700Matp10);
1279 gintDataPall500700Matp10->SetMarkerColor(3);
1280 gintDataPall500700Matp10->SetMarkerStyle(22);
1281 gintDataPall500700Matp10->Draw("p");
1282 TGraph *gintDataPall500700Matm10=new TGraph(15,d0cut,intDataPall500700Matm10);
1283 gintDataPall500700Matm10->SetMarkerColor(3);
1284 gintDataPall500700Matm10->SetMarkerStyle(22);
1285 gintDataPall500700Matm10->Draw("p");
1286
1287 c5->cd(3);
1288 Float_t intDataPall10001500Strp30[15],intDataPall10001500Strm30[15],intDataPall10001500Matp10[15],intDataPall10001500Matm10[15];
1289 for(Int_t i=0;i<15;i++) {
1290 intDataPall10001500Strp30[i]=intDatacut10001500[i]*fracS10001500Strp30[i]/fracP10001500[i];
1291 intDataPall10001500Strm30[i]=intDatacut10001500[i]*fracS10001500Strm30[i]/fracP10001500[i];
1292 intDataPall10001500Matp10[i]=intDatacut10001500[i]*fracS10001500Matp10[i]/fracP10001500[i];
1293 intDataPall10001500Matm10[i]=intDatacut10001500[i]*fracS10001500Matm10[i]/fracP10001500[i];
1294 }
1295 gintDataPall10001500->Draw("ap");
1296
1297 TGraph *gintDataPall10001500Strp30=new TGraph(15,d0cut,intDataPall10001500Strp30);
1298 gintDataPall10001500Strp30->SetMarkerColor(1);
1299 gintDataPall10001500Strp30->SetMarkerStyle(22);
1300 gintDataPall10001500Strp30->Draw("p");
1301 TGraph *gintDataPall10001500Strm30=new TGraph(15,d0cut,intDataPall10001500Strm30);
1302 gintDataPall10001500Strm30->SetMarkerColor(1);
1303 gintDataPall10001500Strm30->SetMarkerStyle(22);
1304 gintDataPall10001500Strm30->Draw("p");
1305 TGraph *gintDataPall10001500Matp10=new TGraph(15,d0cut,intDataPall10001500Matp10);
1306 gintDataPall10001500Matp10->SetMarkerColor(3);
1307 gintDataPall10001500Matp10->SetMarkerStyle(22);
1308 gintDataPall10001500Matp10->Draw("p");
1309 TGraph *gintDataPall10001500Matm10=new TGraph(15,d0cut,intDataPall10001500Matm10);
1310 gintDataPall10001500Matm10->SetMarkerColor(3);
1311 gintDataPall10001500Matm10->SetMarkerStyle(22);
1312 gintDataPall10001500Matm10->Draw("p");
1313 */
1314 return;
1315}
1316//---------------------------------------------------------------------------
1317void PlotImpPar_z() {
1318
1319 TFile *fMC= new TFile("ITS.Performance_lhc10a8_100k.root");
1320
1321 TList *list=(TList*)fMC->Get("cOutputITS");
1322 TH1F *fHistd0zITSMIoneSPDInAccP150200MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP150200");
1323 TH1F *fHistd0zITSMIoneSPDInAccS150200MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS150200");
1324 TH1F *fHistd0zITSMIoneSPDInAccP500700MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP500700");
1325 TH1F *fHistd0zITSMIoneSPDInAccS500700MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS500700");
1326 TH1F *fHistd0zITSMIoneSPDInAccP10001500MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP10001500");
1327 TH1F *fHistd0zITSMIoneSPDInAccS10001500MC = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS10001500");
1328 TH1F *fHistd0zITSMIoneSPDInAcc150200MC=(TH1F*)fHistd0zITSMIoneSPDInAccP150200MC->Clone("fHistd0zITSMIoneSPDInAcc150200MC");
1329 fHistd0zITSMIoneSPDInAcc150200MC->Add(fHistd0zITSMIoneSPDInAccS150200MC);
1330 fHistd0zITSMIoneSPDInAcc150200MC->Scale(1./fHistd0zITSMIoneSPDInAcc150200MC->GetEntries());
1331 TH1F *fHistd0zITSMIoneSPDInAcc500700MC=(TH1F*)fHistd0zITSMIoneSPDInAccP500700MC->Clone("fHistd0zITSMIoneSPDInAcc500700MC");
1332 fHistd0zITSMIoneSPDInAcc500700MC->Add(fHistd0zITSMIoneSPDInAccS500700MC);
1333 fHistd0zITSMIoneSPDInAcc500700MC->Scale(1./fHistd0zITSMIoneSPDInAcc500700MC->GetEntries());
1334 TH1F *fHistd0zITSMIoneSPDInAcc10001500MC=(TH1F*)fHistd0zITSMIoneSPDInAccP10001500MC->Clone("fHistd0zITSMIoneSPDInAcc10001500MC");
1335 fHistd0zITSMIoneSPDInAcc10001500MC->Add(fHistd0zITSMIoneSPDInAccS10001500MC);
1336 fHistd0zITSMIoneSPDInAcc10001500MC->Scale(1./fHistd0zITSMIoneSPDInAcc10001500MC->GetEntries());
1337
1338
1339 TFile *f= new TFile("ITS.Performance_104892.root");
1340
1341 TList *list=(TList*)f->Get("cOutputITS");
1342 TH1F *fHistd0zITSMIoneSPDInAccP150200 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP150200");
1343 TH1F *fHistd0zITSMIoneSPDInAccS150200 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS150200");
1344 TH1F *fHistd0zITSMIoneSPDInAccP500700 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP500700");
1345 TH1F *fHistd0zITSMIoneSPDInAccS500700 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS500700");
1346 TH1F *fHistd0zITSMIoneSPDInAccP10001500 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccP10001500");
1347 TH1F *fHistd0zITSMIoneSPDInAccS10001500 = (TH1F*)list->FindObject("fHistd0zITSMIoneSPDInAccS10001500");
1348 TH1F *fHistd0zITSMIoneSPDInAcc150200=(TH1F*)fHistd0zITSMIoneSPDInAccP150200->Clone("fHistd0zITSMIoneSPDInAcc150200");
1349 fHistd0zITSMIoneSPDInAcc150200->Add(fHistd0zITSMIoneSPDInAccS150200);
1350 fHistd0zITSMIoneSPDInAcc150200->Scale(1./fHistd0zITSMIoneSPDInAcc150200->GetEntries());
1351 TH1F *fHistd0zITSMIoneSPDInAcc500700=(TH1F*)fHistd0zITSMIoneSPDInAccP500700->Clone("fHistd0zITSMIoneSPDInAcc500700");
1352 fHistd0zITSMIoneSPDInAcc500700->Add(fHistd0zITSMIoneSPDInAccS500700);
1353 fHistd0zITSMIoneSPDInAcc500700->Scale(1./fHistd0zITSMIoneSPDInAcc500700->GetEntries());
1354 TH1F *fHistd0zITSMIoneSPDInAcc10001500=(TH1F*)fHistd0zITSMIoneSPDInAccP10001500->Clone("fHistd0zITSMIoneSPDInAcc10001500");
1355 fHistd0zITSMIoneSPDInAcc10001500->Add(fHistd0zITSMIoneSPDInAccS10001500);
1356 fHistd0zITSMIoneSPDInAcc10001500->Scale(1./fHistd0zITSMIoneSPDInAcc10001500->GetEntries());
1357
1358
1359 TCanvas *c1 = new TCanvas("c1","c1");
1360 c1->Divide(3,1);
1361 c1->cd(1);
1362 fHistd0zITSMIoneSPDInAcc150200MC->SetLineColor(2);
1363 fHistd0zITSMIoneSPDInAcc150200MC->Draw();
1364 fHistd0zITSMIoneSPDInAcc150200->SetLineColor(4);
1365 fHistd0zITSMIoneSPDInAcc150200->Draw("same");
1366 c1->cd(2);
1367 fHistd0zITSMIoneSPDInAcc500700MC->SetLineColor(2);
1368 fHistd0zITSMIoneSPDInAcc500700MC->Draw();
1369 fHistd0zITSMIoneSPDInAcc500700->SetLineColor(4);
1370 fHistd0zITSMIoneSPDInAcc500700->Draw("same");
1371 c1->cd(3);
1372 fHistd0zITSMIoneSPDInAcc10001500MC->SetLineColor(2);
1373 fHistd0zITSMIoneSPDInAcc10001500MC->Draw();
1374 fHistd0zITSMIoneSPDInAcc10001500->SetLineColor(4);
1375 fHistd0zITSMIoneSPDInAcc10001500->Draw("same");
1376
1377 TCanvas *c2 = new TCanvas("c2","c2");
1378 c2->Divide(3,1);
1379 c2_1->SetLogy();
1380 c2_2->SetLogy();
1381 c2_3->SetLogy();
1382 c2->cd(1);
1383 fHistd0zITSMIoneSPDInAccP150200MC->SetMinimum(1);
1384 fHistd0zITSMIoneSPDInAccS150200MC->SetMinimum(1);
1385 fHistd0zITSMIoneSPDInAccP150200MC->SetLineColor(1);
1386 fHistd0zITSMIoneSPDInAccP150200MC->Draw();
1387 fHistd0zITSMIoneSPDInAccS150200MC->SetLineColor(6);
1388 fHistd0zITSMIoneSPDInAccS150200MC->Draw("same");
1389 c2->cd(2);
1390 fHistd0zITSMIoneSPDInAccP500700MC->SetMinimum(1);
1391 fHistd0zITSMIoneSPDInAccS500700MC->SetMinimum(1);
1392 fHistd0zITSMIoneSPDInAccP500700MC->SetLineColor(1);
1393 fHistd0zITSMIoneSPDInAccP500700MC->Draw();
1394 fHistd0zITSMIoneSPDInAccS500700MC->SetLineColor(6);
1395 fHistd0zITSMIoneSPDInAccS500700MC->Draw("same");
1396 c2->cd(3);
1397 fHistd0zITSMIoneSPDInAccP10001500MC->SetMinimum(1);
1398 fHistd0zITSMIoneSPDInAccS10001500MC->SetMinimum(1);
1399 fHistd0zITSMIoneSPDInAccP10001500MC->SetLineColor(1);
1400 fHistd0zITSMIoneSPDInAccP10001500MC->Draw();
1401 fHistd0zITSMIoneSPDInAccS10001500MC->SetLineColor(6);
1402 fHistd0zITSMIoneSPDInAccS10001500MC->Draw("same");
1403
1404
1405 TCanvas *c3 = new TCanvas("c3","c3");
1406 c3->Divide(3,1);
1407 c3_1->SetLogx();
1408 c3_2->SetLogx();
1409 c3_3->SetLogx();
1410 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};
1411
1412 c3->cd(1);
1413 Float_t fracP150200[15],fracS150200[15];
1414 Float_t intPtot150200=fHistd0zITSMIoneSPDInAccP150200MC->Integral(1,fHistd0zITSMIoneSPDInAccP150200MC->GetNbinsX());
1415 for(Int_t i=0;i<15;i++) {
1416 Int_t bin1=fHistd0zITSMIoneSPDInAccP150200MC->FindBin(-d0cut[i]);
1417 Int_t bin2=fHistd0zITSMIoneSPDInAccP150200MC->FindBin(+d0cut[i]);
1418 Float_t intPcut150200=fHistd0zITSMIoneSPDInAccP150200MC->Integral(bin1,bin2);
1419 Float_t intScut150200=fHistd0zITSMIoneSPDInAccS150200MC->Integral(bin1,bin2);
1420 fracP150200[i]=intPcut150200/intPtot150200;
1421 fracS150200[i]=1.-intScut150200/(intPcut150200+intScut150200);
1422 }
1423 TGraph *gfracP150200=new TGraph(15,d0cut,fracP150200);
1424 gfracP150200->SetMarkerColor(2);
1425 gfracP150200->SetMarkerStyle(20);
1426 gfracP150200->Draw("ap");
1427 TGraph *gfracS150200=new TGraph(15,d0cut,fracS150200);
1428 gfracS150200->SetMarkerColor(4);
1429 gfracS150200->SetMarkerStyle(21);
1430 gfracS150200->Draw("p");
1431
1432 c3->cd(2);
1433 Float_t fracP500700[15],fracS500700[15];
1434 Float_t intPtot500700=fHistd0zITSMIoneSPDInAccP500700MC->Integral(1,fHistd0zITSMIoneSPDInAccP500700MC->GetNbinsX());
1435 for(Int_t i=0;i<15;i++) {
1436 Int_t bin1=fHistd0zITSMIoneSPDInAccP500700MC->FindBin(-d0cut[i]);
1437 Int_t bin2=fHistd0zITSMIoneSPDInAccP500700MC->FindBin(+d0cut[i]);
1438 Float_t intPcut500700=fHistd0zITSMIoneSPDInAccP500700MC->Integral(bin1,bin2);
1439 Float_t intScut500700=fHistd0zITSMIoneSPDInAccS500700MC->Integral(bin1,bin2);
1440 fracP500700[i]=intPcut500700/intPtot500700;
1441 fracS500700[i]=1.-intScut500700/(intPcut500700+intScut500700);
1442 }
1443 TGraph *gfracP500700=new TGraph(15,d0cut,fracP500700);
1444 gfracP500700->SetMarkerColor(2);
1445 gfracP500700->SetMarkerStyle(20);
1446 gfracP500700->Draw("ap");
1447 TGraph *gfracS500700=new TGraph(15,d0cut,fracS500700);
1448 gfracS500700->SetMarkerColor(4);
1449 gfracS500700->SetMarkerStyle(21);
1450 gfracS500700->Draw("p");
1451
1452 c3->cd(3);
1453 Float_t fracP10001500[15],fracS10001500[15];
1454 Float_t intPtot10001500=fHistd0zITSMIoneSPDInAccP10001500MC->Integral(1,fHistd0zITSMIoneSPDInAccP10001500MC->GetNbinsX());
1455 for(Int_t i=0;i<15;i++) {
1456 Int_t bin1=fHistd0zITSMIoneSPDInAccP10001500MC->FindBin(-d0cut[i]);
1457 Int_t bin2=fHistd0zITSMIoneSPDInAccP10001500MC->FindBin(+d0cut[i]);
1458 Float_t intPcut10001500=fHistd0zITSMIoneSPDInAccP10001500MC->Integral(bin1,bin2);
1459 Float_t intScut10001500=fHistd0zITSMIoneSPDInAccS10001500MC->Integral(bin1,bin2);
1460 fracP10001500[i]=intPcut10001500/intPtot10001500;
1461 fracS10001500[i]=1.-intScut10001500/(intPcut10001500+intScut10001500);
1462 }
1463 TGraph *gfracP10001500=new TGraph(15,d0cut,fracP10001500);
1464 gfracP10001500->SetMarkerColor(2);
1465 gfracP10001500->SetMarkerStyle(20);
1466 gfracP10001500->Draw("ap");
1467 TGraph *gfracS10001500=new TGraph(15,d0cut,fracS10001500);
1468 gfracS10001500->SetMarkerColor(4);
1469 gfracS10001500->SetMarkerStyle(21);
1470 gfracS10001500->Draw("p");
1471
1472 return;
1473}
1474//----------------------------------------------------------------------------
1475void Corrections() {
1476
1477 TFile *f_cutits= new TFile("ITS.Performance_lhc09d10_500k_noz0cut_nodcaTPC.root");
1478
1479 TList *list=(TList*)f_cutits->Get("cOutputITS");
1480
1481 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
1482 Float_t nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
1483
1484 TH1F *fHistPtITSTPCsel = (TH1F*)list->FindObject("fHistPtITSTPCsel");
1485 NormalizePtHist(fHistPtITSTPCsel,nEventsSelWithVertex);
1486 TH1F *fHistPtITSTPCselP = (TH1F*)list->FindObject("fHistPtITSTPCselP");
1487 NormalizePtHist(fHistPtITSTPCselP,nEventsSelWithVertex);
1488 TH1F *fHistPtITSTPCselS = (TH1F*)list->FindObject("fHistPtITSTPCselS");
1489 NormalizePtHist(fHistPtITSTPCselS,nEventsSelWithVertex);
1490 TH1F *fHistPtITSTPCselPfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselPfromStrange");
1491 NormalizePtHist(fHistPtITSTPCselPfromStrange,nEventsSelWithVertex);
1492 TH1F *fHistPtITSTPCselSfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
1493 NormalizePtHist(fHistPtITSTPCselSfromStrange,nEventsSelWithVertex);
1494 TH1F *fHistPtITSTPCselSfromMat = (TH1F*)list->FindObject("fHistPtITSTPCselSfromMat");
1495 NormalizePtHist(fHistPtITSTPCselSfromMat,nEventsSelWithVertex);
1496
1497
1498 TFile *f_cuttpc= new TFile("ITS.Performance_lhc09d10_500k_noz0cut.root");
1499
1500 TList *list=(TList*)f_cuttpc->Get("cOutputITS");
1501
1502 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
1503 nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
1504
1505 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
1506 NormalizePtHist(fHistPtTPCInAcc,nEventsSelWithVertex);
1507 TH1F *fHistPtTPCInAccP = (TH1F*)list->FindObject("fHistPtTPCInAccP");
1508 NormalizePtHist(fHistPtTPCInAccP,nEventsSelWithVertex);
1509 TH1F *fHistPtTPCInAccS = (TH1F*)list->FindObject("fHistPtTPCInAccS");
1510 NormalizePtHist(fHistPtTPCInAccS,nEventsSelWithVertex);
1511 TH1F *fHistPtTPCInAccPfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccPfromStrange");
1512 NormalizePtHist(fHistPtTPCInAccPfromStrange,nEventsSelWithVertex);
1513 TH1F *fHistPtTPCInAccSfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccSfromStrange");
1514 NormalizePtHist(fHistPtTPCInAccSfromStrange,nEventsSelWithVertex);
1515 TH1F *fHistPtTPCInAccSfromMat = (TH1F*)list->FindObject("fHistPtTPCInAccSfromMat");
1516 NormalizePtHist(fHistPtTPCInAccSfromMat,nEventsSelWithVertex);
1517
1518
1519 TFile *f_nocut= new TFile("ITS.Performance_lhc09d10_500k_nod0z0cut.root");
1520
1521 TList *list=(TList*)f_nocut->Get("cOutputITS");
1522
1523 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
1524 nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
1525
1526 TH1F *fHistPtITSTPCselnocut = (TH1F*)list->FindObject("fHistPtITSTPCsel");
1527 NormalizePtHist(fHistPtITSTPCselnocut,nEventsSelWithVertex);
1528 TH1F *fHistPtITSTPCselnocutP = (TH1F*)list->FindObject("fHistPtITSTPCselP");
1529 NormalizePtHist(fHistPtITSTPCselnocutP,nEventsSelWithVertex);
1530 TH1F *fHistPtITSTPCselnocutS = (TH1F*)list->FindObject("fHistPtITSTPCselS");
1531 NormalizePtHist(fHistPtITSTPCselnocutS,nEventsSelWithVertex);
1532 TH1F *fHistPtITSTPCselnocutPfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselPfromStrange");
1533 NormalizePtHist(fHistPtITSTPCselnocutPfromStrange,nEventsSelWithVertex);
1534 TH1F *fHistPtITSTPCselnocutSfromStrange = (TH1F*)list->FindObject("fHistPtITSTPCselSfromStrange");
1535 NormalizePtHist(fHistPtITSTPCselnocutSfromStrange,nEventsSelWithVertex);
1536 TH1F *fHistPtITSTPCselnocutSfromMat = (TH1F*)list->FindObject("fHistPtITSTPCselSfromMat");
1537 NormalizePtHist(fHistPtITSTPCselnocutSfromMat,nEventsSelWithVertex);
1538
1539 TH1F *fHistPtTPCInAccnocut = (TH1F*)list->FindObject("fHistPtTPCInAcc");
1540 NormalizePtHist(fHistPtTPCInAccnocut,nEventsSelWithVertex);
1541 TH1F *fHistPtTPCInAccnocutP = (TH1F*)list->FindObject("fHistPtTPCInAccP");
1542 NormalizePtHist(fHistPtTPCInAccnocutP,nEventsSelWithVertex);
1543 TH1F *fHistPtTPCInAccnocutS = (TH1F*)list->FindObject("fHistPtTPCInAccS");
1544 NormalizePtHist(fHistPtTPCInAccnocutS,nEventsSelWithVertex);
1545 TH1F *fHistPtTPCInAccnocutPfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccPfromStrange");
1546 NormalizePtHist(fHistPtTPCInAccnocutPfromStrange,nEventsSelWithVertex);
1547 TH1F *fHistPtTPCInAccnocutSfromStrange = (TH1F*)list->FindObject("fHistPtTPCInAccSfromStrange");
1548 NormalizePtHist(fHistPtTPCInAccnocutSfromStrange,nEventsSelWithVertex);
1549 TH1F *fHistPtTPCInAccnocutSfromMat = (TH1F*)list->FindObject("fHistPtTPCInAccSfromMat");
1550 NormalizePtHist(fHistPtTPCInAccnocutSfromMat,nEventsSelWithVertex);
1551
1552 // cut efficiencies
1553
1554 TH1F *fHistPtITSTPCeffselP = (TH1F*)fHistPtITSTPCselP->Clone("fHistPtITSTPCeffselP");
1555 fHistPtITSTPCeffselP->Divide(fHistPtITSTPCselnocutP);
1556 TH1F *fHistPtITSTPCeffselS = (TH1F*)fHistPtITSTPCselS->Clone("fHistPtITSTPCeffselS");
1557 fHistPtITSTPCeffselS->Divide(fHistPtITSTPCselnocutS);
1558 TH1F *fHistPtITSTPCeffselSfromStrange = (TH1F*)fHistPtITSTPCselSfromStrange->Clone("fHistPtITSTPCeffselSfromStrange");
1559 fHistPtITSTPCeffselSfromStrange->Divide(fHistPtITSTPCselnocutSfromStrange);
1560 TH1F *fHistPtITSTPCeffselSfromMat = (TH1F*)fHistPtITSTPCselSfromMat->Clone("fHistPtITSTPCeffselSfromMat");
1561 fHistPtITSTPCeffselSfromMat->Divide(fHistPtITSTPCselnocutSfromMat);
1562
1563 TH1F *fHistPtTPCeffInAccP = (TH1F*)fHistPtTPCInAccP->Clone("fHistPtTPCeffInAccP");
1564 fHistPtTPCeffInAccP->Divide(fHistPtTPCInAccnocutP);
1565 TH1F *fHistPtTPCeffInAccS = (TH1F*)fHistPtTPCInAccS->Clone("fHistPtTPCeffInAccS");
1566 fHistPtTPCeffInAccS->Divide(fHistPtTPCInAccnocutS);
1567 TH1F *fHistPtTPCeffInAccSfromStrange = (TH1F*)fHistPtTPCInAccSfromStrange->Clone("fHistPtTPCeffInAccSfromStrange");
1568 fHistPtTPCeffInAccSfromStrange->Divide(fHistPtTPCInAccnocutSfromStrange);
1569 TH1F *fHistPtTPCeffInAccSfromMat = (TH1F*)fHistPtTPCInAccSfromMat->Clone("fHistPtTPCeffInAccSfromMat");
1570 fHistPtTPCeffInAccSfromMat->Divide(fHistPtTPCInAccnocutSfromMat);
1571
1572
1573 TLegend *l0=new TLegend(0.5,0.5,0.9,0.9);
1574 TCanvas *c0=new TCanvas("c0","c0");
1575 c0->Divide(2,1);
1576 c0->cd(1);
1577 TH1F *fHistPtITSTPCselFracS = (TH1F*)fHistPtITSTPCselS->Clone("fHistPtITSTPCselFracS");
1578 fHistPtITSTPCselFracS->Divide(fHistPtITSTPCsel);
1579 fHistPtITSTPCselFracS->SetLineColor(2);
1580 l0->AddEntry(fHistPtITSTPCselFracS,"secondaries","l");
1581 fHistPtITSTPCselFracS->Draw();
1582 TH1F *fHistPtITSTPCselFracSfromStrange = (TH1F*)fHistPtITSTPCselSfromStrange->Clone("fHistPtITSTPCselFracSfromStrange");
1583 fHistPtITSTPCselFracSfromStrange->Divide(fHistPtITSTPCsel);
1584 fHistPtITSTPCselFracSfromStrange->SetLineColor(6);
1585 l0->AddEntry(fHistPtITSTPCselFracSfromStrange,"sec. from strange","l");
1586 fHistPtITSTPCselFracSfromStrange->Draw("same");
1587 TH1F *fHistPtITSTPCselFracSfromMat = (TH1F*)fHistPtITSTPCselSfromMat->Clone("fHistPtITSTPCselFracSfromMat");
1588 fHistPtITSTPCselFracSfromMat->Divide(fHistPtITSTPCsel);
1589 fHistPtITSTPCselFracSfromMat->SetLineColor(kOrange+7);
1590 l0->AddEntry(fHistPtITSTPCselFracSfromMat,"sec. from material","l");
1591 fHistPtITSTPCselFracSfromMat->Draw("same");
1592 l0->Draw();
1593 c0->cd(2);
1594 TH1F *fHistPtTPCInAccFracS = (TH1F*)fHistPtTPCInAccS->Clone("fHistPtTPCInAccFracS");
1595 fHistPtTPCInAccFracS->Divide(fHistPtTPCInAcc);
1596 fHistPtTPCInAccFracS->SetLineColor(2);
1597 fHistPtTPCInAccFracS->Draw();
1598 TH1F *fHistPtTPCInAccFracSfromStrange = (TH1F*)fHistPtTPCInAccSfromStrange->Clone("fHistPtTPCInAccFracSfromStrange");
1599 fHistPtTPCInAccFracSfromStrange->Divide(fHistPtTPCInAcc);
1600 fHistPtTPCInAccFracSfromStrange->SetLineColor(6);
1601 fHistPtTPCInAccFracSfromStrange->Draw("same");
1602 TH1F *fHistPtTPCInAccFracSfromMat = (TH1F*)fHistPtTPCInAccSfromMat->Clone("fHistPtTPCInAccFracSfromMat");
1603 fHistPtTPCInAccFracSfromMat->Divide(fHistPtTPCInAcc);
1604 fHistPtTPCInAccFracSfromMat->SetLineColor(kOrange+7);
1605 fHistPtTPCInAccFracSfromMat->Draw("same");
1606
1607 TLegend *l1=new TLegend(0.5,0.5,0.9,0.9);
1608 TCanvas *c1=new TCanvas("c1","c1");
1609 c1->Divide(2,1);
1610 c1->cd(1);
1611 fHistPtITSTPCeffselP->SetYTitle("dca cut efficiency");
1612 fHistPtITSTPCeffselP->SetLineColor(1);
1613 l1->AddEntry(fHistPtITSTPCeffselP,"primaries","l");
1614 fHistPtITSTPCeffselP->Draw();
1615 fHistPtITSTPCeffselS->SetLineColor(2);
1616 l1->AddEntry(fHistPtITSTPCeffselS,"secondaries","l");
1617 fHistPtITSTPCeffselS->Draw("same");
1618 fHistPtITSTPCeffselSfromStrange->SetLineColor(6);
1619 l1->AddEntry(fHistPtITSTPCeffselSfromStrange,"sec. from strange","l");
1620 fHistPtITSTPCeffselSfromStrange->Draw("same");
1621 fHistPtITSTPCeffselSfromMat->SetLineColor(kOrange+7);
1622 l1->AddEntry(fHistPtITSTPCeffselSfromMat,"sec. from material","l");
1623 fHistPtITSTPCeffselSfromMat->Draw("same");
1624 l1->Draw();
1625 c1->cd(2);
1626 fHistPtTPCeffInAccP->SetYTitle("dca cut efficiency");
1627 fHistPtTPCeffInAccP->SetLineColor(1);
1628 fHistPtTPCeffInAccP->Draw();
1629 fHistPtTPCeffInAccS->SetLineColor(2);
1630 fHistPtTPCeffInAccS->Draw("same");
1631 fHistPtTPCeffInAccSfromStrange->SetLineColor(6);
1632 fHistPtTPCeffInAccSfromStrange->Draw("same");
1633 fHistPtTPCeffInAccSfromMat->SetLineColor(kOrange+7);
1634 fHistPtTPCeffInAccSfromMat->Draw("same");
1635
1636
1637
1638 Float_t weightStrange=1.;
1639 Float_t weightMat=1.;
1640
1641 TH1F *fHistPtITSTPCselCocktail = MakeCocktail(fHistPtITSTPCselP,
1642 fHistPtITSTPCselPfromStrange,
1643 fHistPtITSTPCselS,
1644 fHistPtITSTPCselSfromStrange,
1645 fHistPtITSTPCselSfromMat,
1646 weightStrange,
1647 weightMat,
1648 "fHistPtITSTPCselCocktail");
1649
1650 TH1F *fHistPtITSTPCselPCocktail = CorrectToP(fHistPtITSTPCselCocktail,
1651 fHistPtITSTPCselSfromStrange,
1652 fHistPtITSTPCselSfromMat,
1653 weightStrange,
1654 weightMat,
1655 "fHistPtITSTPCselPCocktail");
1656
1657 TH1F *fHistPtITSTPCallPCocktail = CorrectToPnocut(fHistPtITSTPCselPCocktail,
1658 fHistPtITSTPCeffselP,
1659 "fHistPtITSTPCallPCocktail");
1660
1661
1662 TH1F *fHistPtTPCInAccCocktail = MakeCocktail(fHistPtTPCInAccP,
1663 fHistPtTPCInAccPfromStrange,
1664 fHistPtTPCInAccS,
1665 fHistPtTPCInAccSfromStrange,
1666 fHistPtTPCInAccSfromMat,
1667 weightStrange,
1668 weightMat,
1669 "fHistPtTPCInAccCocktail");
1670
1671 TH1F *fHistPtTPCInAccPCocktail = CorrectToP(fHistPtTPCInAccCocktail,
1672 fHistPtTPCInAccSfromStrange,
1673 fHistPtTPCInAccSfromMat,
1674 weightStrange,
1675 weightMat,
1676 "fHistPtTPCInAccPCocktail");
1677
1678 TH1F *fHistPtTPCallPCocktail = CorrectToPnocut(fHistPtTPCInAccPCocktail,
1679 fHistPtTPCeffInAccP,
1680 "fHistPtTPCallPCocktail");
1681
1682
1683
1684 TLegend *l2=new TLegend(0.5,0.5,0.9,0.9);
1685 TCanvas *c2 = new TCanvas("c2","c2");
1686 c2->Divide(2,1);
1687 c2->cd(1);
1688 fHistPtITSTPCselCocktail->SetLineColor(4);
1689 l2->AddEntry(fHistPtITSTPCselCocktail,"all, with dca cut","l");
1690 fHistPtITSTPCselCocktail->Draw();
1691 fHistPtITSTPCselPCocktail->SetLineColor(7);
1692 l2->AddEntry(fHistPtITSTPCselPCocktail,"primaries, with dca cut","l");
1693 fHistPtITSTPCselPCocktail->Draw("same");
1694 fHistPtITSTPCallPCocktail->SetLineColor(6);
1695 l2->AddEntry(fHistPtITSTPCallPCocktail,"primaries, w/o dca cut","l");
1696 fHistPtITSTPCallPCocktail->Draw("same");
1697 l2->Draw();
1698 c2->cd(2);
1699 fHistPtTPCInAccCocktail->SetLineColor(4);
1700 fHistPtTPCInAccCocktail->Draw();
1701 fHistPtTPCInAccPCocktail->SetLineColor(7);
1702 fHistPtTPCInAccPCocktail->Draw("same");
1703 fHistPtTPCallPCocktail->SetLineColor(6);
1704 fHistPtTPCallPCocktail->Draw("same");
1705
1706
1707 // systematic from secondaries
1708 weightStrange=1.3;
1709 weightMat=1.0;
1710 TH1F *fHistPtITSTPCselPCocktail1 = CorrectToP(fHistPtITSTPCselCocktail,
1711 fHistPtITSTPCselSfromStrange,
1712 fHistPtITSTPCselSfromMat,
1713 weightStrange,
1714 weightMat,
1715 "fHistPtITSTPCselPCocktail1");
1716 TH1F *fHistPtITSTPCallPCocktail1 = CorrectToPnocut(fHistPtITSTPCselPCocktail1,
1717 fHistPtITSTPCeffselP,
1718 "fHistPtITSTPCallPCocktail1");
1719 weightStrange=0.7;
1720 weightMat=1.0;
1721 TH1F *fHistPtITSTPCselPCocktail2 = CorrectToP(fHistPtITSTPCselCocktail,
1722 fHistPtITSTPCselSfromStrange,
1723 fHistPtITSTPCselSfromMat,
1724 weightStrange,
1725 weightMat,
1726 "fHistPtITSTPCselPCocktail2");
1727 TH1F *fHistPtITSTPCallPCocktail2 = CorrectToPnocut(fHistPtITSTPCselPCocktail2,
1728 fHistPtITSTPCeffselP,
1729 "fHistPtITSTPCallPCocktail2");
1730 weightStrange=1.0;
1731 weightMat=1.1;
1732 TH1F *fHistPtITSTPCselPCocktail3 = CorrectToP(fHistPtITSTPCselCocktail,
1733 fHistPtITSTPCselSfromStrange,
1734 fHistPtITSTPCselSfromMat,
1735 weightStrange,
1736 weightMat,
1737 "fHistPtITSTPCselPCocktail3");
1738 TH1F *fHistPtITSTPCallPCocktail3 = CorrectToPnocut(fHistPtITSTPCselPCocktail3,
1739 fHistPtITSTPCeffselP,
1740 "fHistPtITSTPCallPCocktail3");
1741 weightStrange=1.0;
1742 weightMat=0.9;
1743 TH1F *fHistPtITSTPCselPCocktail4 = CorrectToP(fHistPtITSTPCselCocktail,
1744 fHistPtITSTPCselSfromStrange,
1745 fHistPtITSTPCselSfromMat,
1746 weightStrange,
1747 weightMat,
1748 "fHistPtITSTPCselPCocktail4");
1749 TH1F *fHistPtITSTPCallPCocktail4 = CorrectToPnocut(fHistPtITSTPCselPCocktail4,
1750 fHistPtITSTPCeffselP,
1751 "fHistPtITSTPCallPCocktail4");
1752 weightStrange=1.3;
1753 weightMat=1.1;
1754 TH1F *fHistPtITSTPCselPCocktail5 = CorrectToP(fHistPtITSTPCselCocktail,
1755 fHistPtITSTPCselSfromStrange,
1756 fHistPtITSTPCselSfromMat,
1757 weightStrange,
1758 weightMat,
1759 "fHistPtITSTPCselPCocktail5");
1760 TH1F *fHistPtITSTPCallPCocktail5 = CorrectToPnocut(fHistPtITSTPCselPCocktail5,
1761 fHistPtITSTPCeffselP,
1762 "fHistPtITSTPCallPCocktail5");
1763 weightStrange=0.7;
1764 weightMat=0.9;
1765 TH1F *fHistPtITSTPCselPCocktail6 = CorrectToP(fHistPtITSTPCselCocktail,
1766 fHistPtITSTPCselSfromStrange,
1767 fHistPtITSTPCselSfromMat,
1768 weightStrange,
1769 weightMat,
1770 "fHistPtITSTPCselPCocktail6");
1771 TH1F *fHistPtITSTPCallPCocktail6 = CorrectToPnocut(fHistPtITSTPCselPCocktail6,
1772 fHistPtITSTPCeffselP,
1773 "fHistPtITSTPCallPCocktail6");
1774
1775 weightStrange=1.3;
1776 weightMat=1.0;
1777 TH1F *fHistPtTPCInAccPCocktail1 = CorrectToP(fHistPtTPCInAccCocktail,
1778 fHistPtTPCInAccSfromStrange,
1779 fHistPtTPCInAccSfromMat,
1780 weightStrange,
1781 weightMat,
1782 "fHistPtTPCInAccPCocktail1");
1783 TH1F *fHistPtTPCallPCocktail1 = CorrectToPnocut(fHistPtTPCInAccPCocktail1,
1784 fHistPtTPCeffInAccP,
1785 "fHistPtTPCallPCocktail1");
1786 weightStrange=0.7;
1787 weightMat=1.0;
1788 TH1F *fHistPtTPCInAccPCocktail2 = CorrectToP(fHistPtTPCInAccCocktail,
1789 fHistPtTPCInAccSfromStrange,
1790 fHistPtTPCInAccSfromMat,
1791 weightStrange,
1792 weightMat,
1793 "fHistPtTPCInAccPCocktail2");
1794 TH1F *fHistPtTPCallPCocktail2 = CorrectToPnocut(fHistPtTPCInAccPCocktail2,
1795 fHistPtTPCeffInAccP,
1796 "fHistPtTPCallPCocktail2");
1797 weightStrange=1.0;
1798 weightMat=1.1;
1799 TH1F *fHistPtTPCInAccPCocktail3 = CorrectToP(fHistPtTPCInAccCocktail,
1800 fHistPtTPCInAccSfromStrange,
1801 fHistPtTPCInAccSfromMat,
1802 weightStrange,
1803 weightMat,
1804 "fHistPtTPCInAccPCocktail3");
1805 TH1F *fHistPtTPCallPCocktail3 = CorrectToPnocut(fHistPtTPCInAccPCocktail3,
1806 fHistPtTPCeffInAccP,
1807 "fHistPtTPCallPCocktail3");
1808 weightStrange=1.0;
1809 weightMat=0.9;
1810 TH1F *fHistPtTPCInAccPCocktail4 = CorrectToP(fHistPtTPCInAccCocktail,
1811 fHistPtTPCInAccSfromStrange,
1812 fHistPtTPCInAccSfromMat,
1813 weightStrange,
1814 weightMat,
1815 "fHistPtTPCInAccPCocktail4");
1816 TH1F *fHistPtTPCallPCocktail4 = CorrectToPnocut(fHistPtTPCInAccPCocktail4,
1817 fHistPtTPCeffInAccP,
1818 "fHistPtTPCallPCocktail4");
1819 weightStrange=1.3;
1820 weightMat=1.1;
1821 TH1F *fHistPtTPCInAccPCocktail5 = CorrectToP(fHistPtTPCInAccCocktail,
1822 fHistPtTPCInAccSfromStrange,
1823 fHistPtTPCInAccSfromMat,
1824 weightStrange,
1825 weightMat,
1826 "fHistPtTPCInAccPCocktail5");
1827 TH1F *fHistPtTPCallPCocktail5 = CorrectToPnocut(fHistPtTPCInAccPCocktail5,
1828 fHistPtTPCeffInAccP,
1829 "fHistPtTPCallPCocktail5");
1830 weightStrange=0.7;
1831 weightMat=0.9;
1832 TH1F *fHistPtTPCInAccPCocktail6 = CorrectToP(fHistPtTPCInAccCocktail,
1833 fHistPtTPCInAccSfromStrange,
1834 fHistPtTPCInAccSfromMat,
1835 weightStrange,
1836 weightMat,
1837 "fHistPtTPCInAccPCocktail6");
1838 TH1F *fHistPtTPCallPCocktail6 = CorrectToPnocut(fHistPtTPCInAccPCocktail6,
1839 fHistPtTPCeffInAccP,
1840 "fHistPtTPCallPCocktail6");
1841
1842 TLegend *l3=new TLegend(0.5,0.5,0.9,0.9);
1843 TCanvas *c3 = new TCanvas("c3","c3");
1844 c3->Divide(2,1);
1845 c3->cd(1);
1846 l3->AddEntry(fHistPtITSTPCallPCocktail,"central","l");
1847 fHistPtITSTPCallPCocktail->Draw();
1848 l3->AddEntry(fHistPtITSTPCallPCocktail1,"strange #pm 30%","l");
1849 fHistPtITSTPCallPCocktail1->SetLineColor(1);
1850 fHistPtITSTPCallPCocktail2->SetLineColor(1);
1851 fHistPtITSTPCallPCocktail1->Draw("same");
1852 fHistPtITSTPCallPCocktail2->Draw("same");
1853 l3->AddEntry(fHistPtITSTPCallPCocktail3,"material #pm 10%","l");
1854 fHistPtITSTPCallPCocktail3->SetLineColor(8);
1855 fHistPtITSTPCallPCocktail4->SetLineColor(8);
1856 fHistPtITSTPCallPCocktail3->Draw("same");
1857 fHistPtITSTPCallPCocktail4->Draw("same");
1858 l3->AddEntry(fHistPtITSTPCallPCocktail5,"strange #pm 30% & material #pm 10%","l");
1859 fHistPtITSTPCallPCocktail5->SetLineColor(4);
1860 fHistPtITSTPCallPCocktail6->SetLineColor(4);
1861 fHistPtITSTPCallPCocktail5->Draw("same");
1862 fHistPtITSTPCallPCocktail6->Draw("same");
1863 l3->Draw();
1864 c3->cd(2);
1865 fHistPtTPCallPCocktail->Draw();
1866 fHistPtTPCallPCocktail1->SetLineColor(1);
1867 fHistPtTPCallPCocktail2->SetLineColor(1);
1868 fHistPtTPCallPCocktail1->Draw("same");
1869 fHistPtTPCallPCocktail2->Draw("same");
1870 fHistPtTPCallPCocktail3->SetLineColor(8);
1871 fHistPtTPCallPCocktail4->SetLineColor(8);
1872 fHistPtTPCallPCocktail3->Draw("same");
1873 fHistPtTPCallPCocktail4->Draw("same");
1874 fHistPtTPCallPCocktail5->SetLineColor(4);
1875 fHistPtTPCallPCocktail6->SetLineColor(4);
1876 fHistPtTPCallPCocktail5->Draw("same");
1877 fHistPtTPCallPCocktail6->Draw("same");
1878
1879 /*
1880 TCanvas *c4 = new TCanvas("c4","c4");
1881 c4->Divide(2,1);
1882 c4->cd(1);
1883 fHistPtITSTPCallPCocktail2->Add(fHistPtITSTPCallPCocktail,-1);
1884 fHistPtITSTPCallPCocktail2->Divide(fHistPtITSTPCallPCocktail);
1885 fHistPtITSTPCallPCocktail2->Draw();
1886 fHistPtITSTPCallPCocktail4->Add(fHistPtITSTPCallPCocktail,-1);
1887 fHistPtITSTPCallPCocktail4->Divide(fHistPtITSTPCallPCocktail);
1888 fHistPtITSTPCallPCocktail4->Draw("same");
1889 fHistPtITSTPCallPCocktail6->Add(fHistPtITSTPCallPCocktail,-1);
1890 fHistPtITSTPCallPCocktail6->Divide(fHistPtITSTPCallPCocktail);
1891 fHistPtITSTPCallPCocktail6->Draw("same");
1892 l3->Draw();
1893 c4->cd(2);
1894 fHistPtTPCallPCocktail2->Add(fHistPtTPCallPCocktail,-1);
1895 fHistPtTPCallPCocktail2->Divide(fHistPtTPCallPCocktail);
1896 fHistPtTPCallPCocktail2->Draw();
1897 fHistPtTPCallPCocktail4->Add(fHistPtTPCallPCocktail,-1);
1898 fHistPtTPCallPCocktail4->Divide(fHistPtTPCallPCocktail);
1899 fHistPtTPCallPCocktail4->Draw("same");
1900 fHistPtTPCallPCocktail6->Add(fHistPtTPCallPCocktail,-1);
1901 fHistPtTPCallPCocktail6->Divide(fHistPtTPCallPCocktail);
1902 fHistPtTPCallPCocktail6->Draw("same");
1903 */
1904
1905 return;
1906}
1907//-------------------------------------------------------------
1908void NormalizePtHist(TH1F *h,Float_t n) {
1909
1910 for(Int_t i=1; i<=h->GetNbinsX(); i++) {
1911 if(h->GetBinContent(i)) {
1912 h->SetBinContent(i,h->GetBinContent(i)/h->GetBinWidth(i));
1913 h->SetBinError(i,h->GetBinError(i)/h->GetBinWidth(i));
1914 }
1915 }
1916 h->Scale(1./n);
1917
1918 return;
1919}
1920//----------------------------------------------------------------------
1921TH1F *MakeCocktail(TH1F *fHistPtselP,
1922 TH1F *fHistPtselPfromStrange,
1923 TH1F *fHistPtselS,
1924 TH1F *fHistPtselSfromStrange,
1925 TH1F *fHistPtselSfromMat,
1926 Float_t weightStrange,
1927 Float_t weightMat,TString name) {
1928
1929 // build histogram with yield passing cut
1930 TH1F *fHistPtselCocktail = (TH1F*)fHistPtselP->Clone(name.Data());
1931 fHistPtselCocktail->Add(fHistPtselS);
1932 fHistPtselCocktail->Add(fHistPtselPfromStrange,weightStrange-1.);
1933 fHistPtselCocktail->Add(fHistPtselSfromStrange,weightStrange-1.);
1934 fHistPtselCocktail->Add(fHistPtselSfromMat,weightMat-1.);
1935
1936
1937 return fHistPtselCocktail;
1938}
1939//----------------------------------------------------------------------
1940TH1F *CorrectToP(TH1F *fHistPtselCocktail,
1941 TH1F *fHistPtselSfromStrange,
1942 TH1F *fHistPtselSfromMat,
1943 Float_t weightStrange,
1944 Float_t weightMat,TString name) {
1945
1946
1947 // build histogram with yield of primaries passing cut
1948 TH1F* fHistPtselPCocktail = (TH1F*)fHistPtselCocktail->Clone(name.Data());
1949 fHistPtselPCocktail->Add(fHistPtselSfromStrange,-weightStrange);
1950 fHistPtselPCocktail->Add(fHistPtselSfromMat,-weightMat);
1951
1952 return fHistPtselPCocktail;
1953}
1954//----------------------------------------------------------------------
1955TH1F *CorrectToPnocut(TH1F *fHistPtselPCocktail,
1956 TH1F *fHistPteffselP,TString name) {
1957
1958 // build histogram with yield of primaries without cut
1959 TH1F *fHistPtallPCocktail = (TH1F*)fHistPtselPCocktail->Clone(name.Data());
1960 //printf("%f\n",fHistPtallPCocktail->GetBinContent(6));
1961 fHistPtallPCocktail->Divide(fHistPteffselP);
1962 //printf("%f\n",fHistPtallPCocktail->GetBinContent(6));
1963
1964 return fHistPtallPCocktail;
1965}
1966//----------------------------------------------------------------------
1967void CompareTPCPt() {
1968
1969 TFile *f_tpcrefit= new TFile("ITS.Performance_104892pass2_etale05.root");
1970
1971 TList *list=(TList*)f_tpcrefit->Get("cOutputITS");
1972
1973 TH1F *fHistNEvents = (TH1F*)list->FindObject("fHistNEvents");
1974 Float_t nEventsSelWithVertex = fHistNEvents->GetBinContent(5);
1975
1976 TH1F *fHistPtTPCInAcc = (TH1F*)list->FindObject("fHistPtTPCInAcc");
1977 NormalizePtHist(fHistPtTPCInAcc,nEventsSelWithVertex);
1978 TH1F *fHistPtTPCInAccNoTRDout = (TH1F*)list->FindObject("fHistPtTPCInAccNoTRDout");
1979 NormalizePtHist(fHistPtTPCInAccNoTRDout,nEventsSelWithVertex);
1980 TH1F *fHistPtTPCInAccNoTOFout = (TH1F*)list->FindObject("fHistPtTPCInAccNoTOFout");
1981 NormalizePtHist(fHistPtTPCInAccNoTOFout,nEventsSelWithVertex);
1982 TH1F *fHistPtTPCInAccWithPtTPCAtVtx = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtVtx");
1983 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtVtx,nEventsSelWithVertex);
1984 TH1F *fHistPtTPCInAccWithPtTPCAtInnerWall = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtInnerWall");
1985 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtInnerWall,nEventsSelWithVertex);
1986
1987
1988 TFile *f_notpcrefit= new TFile("ITS.Performance_104892pass3_etale05.root");
1989
1990 TList *list=(TList*)f_notpcrefit->Get("cOutputITS");
1991
1992 TH1F *fHistNEventsNoRefit = (TH1F*)list->FindObject("fHistNEvents");
1993 nEventsSelWithVertex = fHistNEventsNoRefit->GetBinContent(5);
1994
1995 TH1F *fHistPtTPCInAccNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAcc");
1996 NormalizePtHist(fHistPtTPCInAccNoRefit,nEventsSelWithVertex);
1997 TH1F *fHistPtTPCInAccNoTRDoutNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccNoTRDout");
1998 NormalizePtHist(fHistPtTPCInAccNoTRDoutNoRefit,nEventsSelWithVertex);
1999 TH1F *fHistPtTPCInAccNoTOFoutNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccNoTOFout");
2000 NormalizePtHist(fHistPtTPCInAccNoTOFoutNoRefit,nEventsSelWithVertex);
2001 TH1F *fHistPtTPCInAccWithPtTPCAtVtxNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtVtx");
2002 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtVtxNoRefit,nEventsSelWithVertex);
2003 TH1F *fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit = (TH1F*)list->FindObject("fHistPtTPCInAccWithPtTPCAtInnerWall");
2004 NormalizePtHist(fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit,nEventsSelWithVertex);
2005
2006
2007
2008 TLegend *l1 = new TLegend(0.5,0.5,0.9,0.9);
2009 TCanvas *c1=new TCanvas("c1","c1",0,0,1000,600);
2010 c1->Divide(2,1);
2011 c1->cd(1);
2012 fHistPtTPCInAcc->SetMarkerStyle(20);
2013 l1->AddEntry(fHistPtTPCInAcc,"p_{t} measurement:","");
2014 l1->AddEntry(fHistPtTPCInAcc,"at vtx, from main params","p");
2015 fHistPtTPCInAcc->Draw();
2016 fHistPtTPCInAccWithPtTPCAtVtx->SetMarkerStyle(21);
2017 fHistPtTPCInAccWithPtTPCAtVtx->SetMarkerColor(2);
2018 l1->AddEntry(fHistPtTPCInAccWithPtTPCAtVtx,"at vtx, from TPC params","p");
2019 fHistPtTPCInAccWithPtTPCAtVtx->Draw("same");
2020 fHistPtTPCInAccWithPtTPCAtInnerWall->SetMarkerStyle(25);
2021 fHistPtTPCInAccWithPtTPCAtInnerWall->SetMarkerColor(2);
2022 l1->AddEntry(fHistPtTPCInAccWithPtTPCAtInnerWall,"at r=80cm, from TPC params","p");
2023 fHistPtTPCInAccWithPtTPCAtInnerWall->Draw("same");
2024 fHistPtTPCInAccNoTOFout->SetMarkerStyle(23);
2025 fHistPtTPCInAccNoTOFout->SetMarkerColor(3);
2026 fHistPtTPCInAccNoTOFout->Draw("same");
2027 fHistPtTPCInAccNoTRDout->SetMarkerStyle(24);
2028 fHistPtTPCInAccNoTRDout->SetMarkerColor(4);
2029 fHistPtTPCInAccNoTRDout->Draw("same");
2030 l1->Draw();
2031 c1->cd(2);
2032 fHistPtTPCInAccNoRefit->SetMarkerStyle(20);
2033 fHistPtTPCInAccNoRefit->Draw();
2034 fHistPtTPCInAccWithPtTPCAtVtxNoRefit->SetMarkerStyle(21);
2035 fHistPtTPCInAccWithPtTPCAtVtxNoRefit->SetMarkerColor(2);
2036 fHistPtTPCInAccWithPtTPCAtVtxNoRefit->Draw("same");
2037 fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit->SetMarkerStyle(25);
2038 fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit->SetMarkerColor(2);
2039 fHistPtTPCInAccWithPtTPCAtInnerWallNoRefit->Draw("same");
2040 fHistPtTPCInAccNoTOFoutNoRefit->SetMarkerStyle(23);
2041 fHistPtTPCInAccNoTOFoutNoRefit->SetMarkerColor(3);
2042 fHistPtTPCInAccNoTOFoutNoRefit->Draw("same");
2043 fHistPtTPCInAccNoTRDoutNoRefit->SetMarkerStyle(24);
2044 fHistPtTPCInAccNoTRDoutNoRefit->SetMarkerColor(4);
2045 fHistPtTPCInAccNoTRDoutNoRefit->Draw("same");
2046 //l1->Draw();
2047
2048 return;
2049}
2050
2051TH1F* gHistTemplateStrange;
2052TH1F* gHistTemplatePandMat;
2053Float_t gd0min=0.15;
2054Float_t gd0max=1.5;
2055
2056
2057Double_t FitMCTemplate(Double_t *x, Double_t *par) {
2058 Double_t xx = x[0];
2059 if(TMath::Abs(xx)<gd0min) return 0;
2060 Int_t bin = gHistTemplateStrange->GetXaxis()->FindBin(xx);
2061 Double_t pandmat = par[0]*gHistTemplatePandMat->GetBinContent(bin);
2062 Double_t strange = par[1]*gHistTemplateStrange->GetBinContent(bin);
2063 return strange + pandmat;
2064}
2065
2066void FitImpPar_rphi() {
2067
2068 Int_t rebin=1;
2069
2070 TFile *fMC= new TFile("ITS.Performance_lhc10a8a12_trkvtx_noTPCdca_120310.root");
2071
2072 TList *list=(TList*)fMC->Get("cOutputITS");
2073 TH1F *fHistd0rphiITSMIoneSPDInAccPMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
2074 TH1F *fHistd0rphiITSMIoneSPDInAccSMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
2075 TH1F *fHistd0rphiITSMIoneSPDInAccSfromMatMC = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200fromMat");
2076 fHistd0rphiITSMIoneSPDInAccPMC->Rebin(rebin);
2077 fHistd0rphiITSMIoneSPDInAccSMC->Rebin(rebin);
2078 fHistd0rphiITSMIoneSPDInAccSfromMatMC->Rebin(rebin);
2079
2080 TH1F *fHistd0rphiITSMIoneSPDInAccSfromStrangeMC=(TH1F*)fHistd0rphiITSMIoneSPDInAccSMC->Clone("fHistd0rphiITSMIoneSPDInAccSfromStrangeMC");
2081 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Add(fHistd0rphiITSMIoneSPDInAccSfromMatMC,-1.);
2082 TH1F *fHistd0rphiITSMIoneSPDInAccMC=(TH1F*)fHistd0rphiITSMIoneSPDInAccPMC->Clone("fHistd0rphiITSMIoneSPDInAccMC");
2083 fHistd0rphiITSMIoneSPDInAccMC->Add(fHistd0rphiITSMIoneSPDInAccSMC);
2084 //fHistd0rphiITSMIoneSPDInAccMC->Scale(1./fHistd0rphiITSMIoneSPDInAccMC->GetEntries());
2085 TH1F *fHistd0rphiITSMIoneSPDInAccPandMatMC=(TH1F*)fHistd0rphiITSMIoneSPDInAccPMC->Clone("fHistd0rphiITSMIoneSPDInAccPandMatMC");
2086 fHistd0rphiITSMIoneSPDInAccPandMatMC->Add(fHistd0rphiITSMIoneSPDInAccSfromMatMC);
2087
2088 // templates
2089 gHistTemplateStrange = (TH1F*)fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Clone("gHistTemplateStrange");
2090 gHistTemplatePandMat = (TH1F*)fHistd0rphiITSMIoneSPDInAccPandMatMC->Clone("gHistTemplatePandMat");
2091
2092
2093
2094
2095 TFile *f= new TFile("ITS.Performance_104065_104892_pass4_trkvtx_noTPCdca.root");
2096 //TFile *f= new TFile("ITS.Performance_lhc10a8_300k_trkvtx_noTPCdca.root");
2097
2098 TList *list=(TList*)f->Get("cOutputITS");
2099 TH1F *fHistd0rphiITSMIoneSPDInAccP = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccP150200");
2100 TH1F *fHistd0rphiITSMIoneSPDInAccS = (TH1F*)list->FindObject("fHistd0rphiITSMIoneSPDInAccS150200");
2101 fHistd0rphiITSMIoneSPDInAccP->Rebin(rebin);
2102 fHistd0rphiITSMIoneSPDInAccS->Rebin(rebin);
2103
2104 TH1F *fHistd0rphiITSMIoneSPDInAcc=(TH1F*)fHistd0rphiITSMIoneSPDInAccP->Clone("fHistd0rphiITSMIoneSPDInAcc");
2105 fHistd0rphiITSMIoneSPDInAcc->Add(fHistd0rphiITSMIoneSPDInAccS);
2106
2107
2108 TF1 *fitMCtemplate = new TF1("fitMCtemplate",FitMCTemplate,-gd0max,+gd0max,2);
2109
2110 TH1F* fHistd0rphiITSMIoneSPDInAccCopy = (TH1F*)fHistd0rphiITSMIoneSPDInAcc->Clone("fHistd0rphiITSMIoneSPDInAccCopy");
2111 for(Int_t ib=1;ib<=fHistd0rphiITSMIoneSPDInAccCopy->GetNbinsX();ib++) {
2112 if(TMath::Abs(fHistd0rphiITSMIoneSPDInAccCopy->GetBinCenter(ib))<gd0min) {
2113 fHistd0rphiITSMIoneSPDInAccCopy->SetBinContent(ib,0);
2114 fHistd0rphiITSMIoneSPDInAccCopy->SetBinError(ib,0);
2115 }
2116 }
2117 fHistd0rphiITSMIoneSPDInAccCopy->Fit("fitMCtemplate","RL");
2118
2119
2120 fHistd0rphiITSMIoneSPDInAccPandMatMC->Scale(fitMCtemplate->GetParameter(0));
2121 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Scale(fitMCtemplate->GetParameter(1));
2122 TH1F *fHistd0rphiITSMIoneSPDInAccMCFit=(TH1F*)fHistd0rphiITSMIoneSPDInAccPandMatMC->Clone("fHistd0rphiITSMIoneSPDInAccMCFit");
2123 fHistd0rphiITSMIoneSPDInAccMCFit->Add(fHistd0rphiITSMIoneSPDInAccSfromStrangeMC);
2124
2125 TCanvas *c1 = new TCanvas("c1","c1",1000,1000);
2126 c1->SetLogy();
2127 fHistd0rphiITSMIoneSPDInAccCopy->SetMinimum(1);
2128 fHistd0rphiITSMIoneSPDInAccCopy->SetMaximum(2*fHistd0rphiITSMIoneSPDInAcc->GetMaximum());
2129 fHistd0rphiITSMIoneSPDInAcc->SetMinimum(1);
2130 fHistd0rphiITSMIoneSPDInAccMC->SetMinimum(1);
2131 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->SetMinimum(1);
2132 fHistd0rphiITSMIoneSPDInAccPandMatMC->SetMinimum(1);
2133 fHistd0rphiITSMIoneSPDInAccCopy->SetLineColor(4);
2134 fHistd0rphiITSMIoneSPDInAccCopy->Draw();
2135 fHistd0rphiITSMIoneSPDInAcc->SetLineColor(4);
2136 fHistd0rphiITSMIoneSPDInAcc->Draw("same");
2137 fHistd0rphiITSMIoneSPDInAccMCFit->SetLineColor(2);
2138 fHistd0rphiITSMIoneSPDInAccMCFit->Draw("same");
2139 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->SetLineColor(3);
2140 fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Draw("same");
2141 fHistd0rphiITSMIoneSPDInAccPandMatMC->SetLineColor(15);
2142 fHistd0rphiITSMIoneSPDInAccPandMatMC->Draw("same");
2143
2144
2145
2146 TCanvas *c3 = new TCanvas("c3","c3");
2147 c3->SetLogx();
2148 c3->SetGridx();
2149 c3->SetGridy();
2150 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};
2151
2152 Float_t fracSfromStrange[15],fracSfromStrangeMC[15],fracSfromMatMC[15],fracS[15],fracSMC[15],fracSRatio[15];
2153 for(Int_t i=0;i<15;i++) {
2154 Int_t bin1=fHistd0rphiITSMIoneSPDInAccPandMatMC->FindBin(-d0cut[i]);
2155 Int_t bin2=fHistd0rphiITSMIoneSPDInAccPandMatMC->FindBin(+d0cut[i]);
2156 Float_t intDatacut=fHistd0rphiITSMIoneSPDInAcc->Integral(bin1,bin2);
2157 Float_t intMCcut=fHistd0rphiITSMIoneSPDInAccMC->Integral(bin1,bin2);
2158 Float_t intSfromStrangecut=fHistd0rphiITSMIoneSPDInAccSfromStrangeMC->Integral(bin1,bin2);
2159 Float_t intSfromMatcut=fHistd0rphiITSMIoneSPDInAccSfromMatMC->Integral(bin1,bin2);
2160 fracSfromStrange[i]=intSfromStrangecut/intDatacut;
2161 fracSfromStrangeMC[i]=intSfromStrangecut/fitMCtemplate->GetParameter(1)/intMCcut;
2162 fracSfromMatMC[i]=intSfromMatcut/intMCcut;
2163
2164 fracS[i]=fracSfromStrange[i]+fracSfromMatMC[i];
2165 fracSMC[i]=fracSfromStrangeMC[i]+fracSfromMatMC[i];
2166
2167 fracSRatio[i]=fracS[i]/fracSMC[i];
2168 }
2169 TGraph *gfracSfromStrange=new TGraph(15,d0cut,fracSfromStrange);
2170 gfracSfromStrange->SetMarkerColor(2);
2171 gfracSfromStrange->SetMarkerStyle(20);
2172 gfracSfromStrange->Draw("ap");
2173
2174 TGraph *gfracSfromStrangeMC=new TGraph(15,d0cut,fracSfromStrangeMC);
2175 gfracSfromStrangeMC->SetMarkerColor(1);
2176 gfracSfromStrangeMC->SetMarkerStyle(22);
2177 gfracSfromStrangeMC->Draw("p");
2178
2179 TGraph *gfracSfromMatMC=new TGraph(15,d0cut,fracSfromMatMC);
2180 gfracSfromMatMC->SetMarkerColor(4);
2181 gfracSfromMatMC->SetMarkerStyle(21);
2182 gfracSfromMatMC->Draw("p");
2183
2184 TGraph *gfracS=new TGraph(15,d0cut,fracS);
2185 gfracS->SetMarkerColor(2);
2186 gfracS->SetMarkerStyle(24);
2187 gfracS->Draw("p");
2188
2189 TGraph *gfracSMC=new TGraph(15,d0cut,fracSMC);
2190 gfracSMC->SetMarkerColor(1);
2191 gfracSMC->SetMarkerStyle(25);
2192 gfracSMC->Draw("p");
2193
2194
2195 TCanvas *c4 = new TCanvas("c4","c4");
2196 c4->SetLogx();
2197 c4->SetGridx();
2198 c4->SetGridy();
2199
2200 TGraph *gfracSRatio=new TGraph(15,d0cut,fracSRatio);
2201 gfracSRatio->SetMarkerColor(9);
2202 gfracSRatio->SetMarkerStyle(21);
2203 gfracSRatio->Draw("ap");
2204
2205
2206 return;
2207}
2208//-----------------------------------------------------------------------------
2209void PlotEffVSErrors() {
2210
2211 TFile *f1= new TFile("eff_104892_pd_testnew1.root");
2212
2213 TH1F *fHistPtITSMI6InAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMI6InAcc");
2214 TH1F *fHistPtITSMIge2InAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMIge2InAcc");
2215 TH1F *fHistPtITSMISPDInAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMISPDInAcc");
2216 TH1F *fHistPtITSMIoneSPDInAcctestnew1 = (TH1F*)f1->Get("fHistPtITSMIoneSPDInAcc");
2217
2218 fHistPtITSMI6InAcctestnew1->SetMarkerStyle(20);
2219 fHistPtITSMIge2InAcctestnew1->SetMarkerStyle(20);
2220 fHistPtITSMISPDInAcctestnew1->SetMarkerStyle(20);
2221 fHistPtITSMIoneSPDInAcctestnew1->SetMarkerStyle(20);
2222 fHistPtITSMI6InAcctestnew1->SetMarkerColor(2);
2223 fHistPtITSMIge2InAcctestnew1->SetMarkerColor(1);
2224 fHistPtITSMISPDInAcctestnew1->SetMarkerColor(9);
2225 fHistPtITSMIoneSPDInAcctestnew1->SetMarkerColor(15);
2226
2227
2228 TFile *f2= new TFile("eff_104892_pd_testnew2.root");
2229
2230 TH1F *fHistPtITSMI6InAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMI6InAcc");
2231 TH1F *fHistPtITSMIge2InAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMIge2InAcc");
2232 TH1F *fHistPtITSMISPDInAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMISPDInAcc");
2233 TH1F *fHistPtITSMIoneSPDInAcctestnew2 = (TH1F*)f2->Get("fHistPtITSMIoneSPDInAcc");
2234
2235 fHistPtITSMI6InAcctestnew2->SetMarkerStyle(21);
2236 fHistPtITSMIge2InAcctestnew2->SetMarkerStyle(21);
2237 fHistPtITSMISPDInAcctestnew2->SetMarkerStyle(21);
2238 fHistPtITSMIoneSPDInAcctestnew2->SetMarkerStyle(21);
2239 fHistPtITSMI6InAcctestnew2->SetMarkerColor(2);
2240 fHistPtITSMIge2InAcctestnew2->SetMarkerColor(1);
2241 fHistPtITSMISPDInAcctestnew2->SetMarkerColor(9);
2242 fHistPtITSMIoneSPDInAcctestnew2->SetMarkerColor(15);
2243
2244
2245 TFile *f3= new TFile("eff_104892_pd_testnew3.root");
2246
2247 TH1F *fHistPtITSMI6InAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMI6InAcc");
2248 TH1F *fHistPtITSMIge2InAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMIge2InAcc");
2249 TH1F *fHistPtITSMISPDInAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMISPDInAcc");
2250 TH1F *fHistPtITSMIoneSPDInAcctestnew3 = (TH1F*)f3->Get("fHistPtITSMIoneSPDInAcc");
2251
2252 fHistPtITSMI6InAcctestnew3->SetMarkerStyle(22);
2253 fHistPtITSMIge2InAcctestnew3->SetMarkerStyle(22);
2254 fHistPtITSMISPDInAcctestnew3->SetMarkerStyle(22);
2255 fHistPtITSMIoneSPDInAcctestnew3->SetMarkerStyle(22);
2256 fHistPtITSMI6InAcctestnew3->SetMarkerColor(2);
2257 fHistPtITSMIge2InAcctestnew3->SetMarkerColor(1);
2258 fHistPtITSMISPDInAcctestnew3->SetMarkerColor(9);
2259 fHistPtITSMIoneSPDInAcctestnew3->SetMarkerColor(15);
2260
2261
2262 TFile *f5= new TFile("eff_104892_pd_testnew5.root");
2263
2264 TH1F *fHistPtITSMI6InAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMI6InAcc");
2265 TH1F *fHistPtITSMIge2InAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMIge2InAcc");
2266 TH1F *fHistPtITSMISPDInAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMISPDInAcc");
2267 TH1F *fHistPtITSMIoneSPDInAcctestnew5 = (TH1F*)f5->Get("fHistPtITSMIoneSPDInAcc");
2268
2269 fHistPtITSMI6InAcctestnew5->SetMarkerStyle(24);
2270 fHistPtITSMIge2InAcctestnew5->SetMarkerStyle(24);
2271 fHistPtITSMISPDInAcctestnew5->SetMarkerStyle(24);
2272 fHistPtITSMIoneSPDInAcctestnew5->SetMarkerStyle(24);
2273 fHistPtITSMI6InAcctestnew5->SetMarkerColor(2);
2274 fHistPtITSMIge2InAcctestnew5->SetMarkerColor(1);
2275 fHistPtITSMISPDInAcctestnew5->SetMarkerColor(9);
2276 fHistPtITSMIoneSPDInAcctestnew5->SetMarkerColor(15);
2277
2278
2279 TFile *f6= new TFile("eff_104892_pd_testnew6.root");
2280
2281 TH1F *fHistPtITSMI6InAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMI6InAcc");
2282 TH1F *fHistPtITSMIge2InAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMIge2InAcc");
2283 TH1F *fHistPtITSMISPDInAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMISPDInAcc");
2284 TH1F *fHistPtITSMIoneSPDInAcctestnew6 = (TH1F*)f6->Get("fHistPtITSMIoneSPDInAcc");
2285
2286 fHistPtITSMI6InAcctestnew6->SetMarkerStyle(25);
2287 fHistPtITSMIge2InAcctestnew6->SetMarkerStyle(25);
2288 fHistPtITSMISPDInAcctestnew6->SetMarkerStyle(25);
2289 fHistPtITSMIoneSPDInAcctestnew6->SetMarkerStyle(25);
2290 fHistPtITSMI6InAcctestnew6->SetMarkerColor(2);
2291 fHistPtITSMIge2InAcctestnew6->SetMarkerColor(1);
2292 fHistPtITSMISPDInAcctestnew6->SetMarkerColor(9);
2293 fHistPtITSMIoneSPDInAcctestnew6->SetMarkerColor(15);
2294
2295
2296 TFile *f7= new TFile("eff_104892_pd_testnew7.root");
2297
2298 TH1F *fHistPtITSMI6InAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMI6InAcc");
2299 TH1F *fHistPtITSMIge2InAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMIge2InAcc");
2300 TH1F *fHistPtITSMISPDInAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMISPDInAcc");
2301 TH1F *fHistPtITSMIoneSPDInAcctestnew7 = (TH1F*)f7->Get("fHistPtITSMIoneSPDInAcc");
2302
2303 fHistPtITSMI6InAcctestnew7->SetMarkerStyle(23);
2304 fHistPtITSMIge2InAcctestnew7->SetMarkerStyle(23);
2305 fHistPtITSMISPDInAcctestnew7->SetMarkerStyle(23);
2306 fHistPtITSMIoneSPDInAcctestnew7->SetMarkerStyle(23);
2307 fHistPtITSMI6InAcctestnew7->SetMarkerColor(2);
2308 fHistPtITSMIge2InAcctestnew7->SetMarkerColor(1);
2309 fHistPtITSMISPDInAcctestnew7->SetMarkerColor(9);
2310 fHistPtITSMIoneSPDInAcctestnew7->SetMarkerColor(15);
2311
2312
2313 TFile *f8= new TFile("eff_104892_pd_testnew8.root");
2314
2315 TH1F *fHistPtITSMI6InAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMI6InAcc");
2316 TH1F *fHistPtITSMIge2InAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMIge2InAcc");
2317 TH1F *fHistPtITSMISPDInAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMISPDInAcc");
2318 TH1F *fHistPtITSMIoneSPDInAcctestnew8 = (TH1F*)f8->Get("fHistPtITSMIoneSPDInAcc");
2319
2320 fHistPtITSMI6InAcctestnew8->SetMarkerStyle(28);
2321 fHistPtITSMIge2InAcctestnew8->SetMarkerStyle(28);
2322 fHistPtITSMISPDInAcctestnew8->SetMarkerStyle(28);
2323 fHistPtITSMIoneSPDInAcctestnew8->SetMarkerStyle(28);
2324 fHistPtITSMI6InAcctestnew8->SetMarkerColor(2);
2325 fHistPtITSMIge2InAcctestnew8->SetMarkerColor(1);
2326 fHistPtITSMISPDInAcctestnew8->SetMarkerColor(9);
2327 fHistPtITSMIoneSPDInAcctestnew8->SetMarkerColor(15);
2328
2329
2330 TCanvas *c= new TCanvas("c","c");
2331 fHistPtITSMI6InAcctestnew1->Draw();
2332 fHistPtITSMIge2InAcctestnew1->Draw("same");
2333 fHistPtITSMISPDInAcctestnew1->Draw("same");
2334 fHistPtITSMIoneSPDInAcctestnew1->Draw("same");
2335 fHistPtITSMI6InAcctestnew2->Draw("same");
2336 fHistPtITSMIge2InAcctestnew2->Draw("same");
2337 fHistPtITSMISPDInAcctestnew2->Draw("same");
2338 fHistPtITSMIoneSPDInAcctestnew2->Draw("same");
2339 fHistPtITSMI6InAcctestnew3->Draw("same");
2340 fHistPtITSMIge2InAcctestnew3->Draw("same");
2341 fHistPtITSMISPDInAcctestnew3->Draw("same");
2342 fHistPtITSMIoneSPDInAcctestnew3->Draw("same");
2343 fHistPtITSMI6InAcctestnew5->Draw("same");
2344 fHistPtITSMIge2InAcctestnew5->Draw("same");
2345 fHistPtITSMISPDInAcctestnew5->Draw("same");
2346 fHistPtITSMIoneSPDInAcctestnew5->Draw("same");
2347 fHistPtITSMI6InAcctestnew6->Draw("same");
2348 fHistPtITSMIge2InAcctestnew6->Draw("same");
2349 fHistPtITSMISPDInAcctestnew6->Draw("same");
2350 fHistPtITSMIoneSPDInAcctestnew6->Draw("same");
2351 fHistPtITSMI6InAcctestnew7->Draw("same");
2352 fHistPtITSMIge2InAcctestnew7->Draw("same");
2353 fHistPtITSMISPDInAcctestnew7->Draw("same");
2354 fHistPtITSMIoneSPDInAcctestnew7->Draw("same");
2355 fHistPtITSMI6InAcctestnew8->Draw("same");
2356 fHistPtITSMIge2InAcctestnew8->Draw("same");
2357 fHistPtITSMISPDInAcctestnew8->Draw("same");
2358 fHistPtITSMIoneSPDInAcctestnew8->Draw("same");
2359
2360
2308c814 2361 return;
2362}