]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG1/macros/plotSparseAL.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
moving AddTaskMeanVertexCalib.C
[u/mrichter/AliRoot.git] / PWG1 / macros / plotSparseAL.C
CommitLineData
852b4956 1{
2
3//TFile f("output/pp-trunk010709-pdc1-00XXX.root")
4TFile f("results_V0QA_v416rev08_090709.root");
5
6
7THnSparseF* sparseAL;
8
9sparseAL=(THnSparseF*)(outputTRD->FindObject("sparseAL"));
10
11sparseAL->GetAxis(5)->SetRangeUser(0.,3.9);
12sparseAL->GetAxis(9)->SetRangeUser(0.,3.9);
13
14
15TH1* ptgeantPPhys=sparseAL->Projection(5);
16TH1* etageantPPhys=sparseAL->Projection(6);
17TH1* phigeantPPhys=sparseAL->Projection(7);
18
19TH1* ptResSinglePPhys=sparseAL->Projection(13);
20TH1* etasinglePPhys=sparseAL->Projection(14);
21TH1* phisinglePPhys=sparseAL->Projection(15);
22
23
24TH1* ptgeantNPhys=sparseAL->Projection(9);
25TH1* etageantNPhys=sparseAL->Projection(10);
26TH1* phigeantNPhys=sparseAL->Projection(11);
27
28TH1* ptResSingleNPhys=sparseAL->Projection(21);
29TH1* etasingleNPhys=sparseAL->Projection(22);
30TH1* phisingleNPhys=sparseAL->Projection(23);
31
32
33TH1* efiEtaSingPPhys=(TH1*) etageantPPhys->Clone();
34efiEtaSingPPhys->Sumw2();
35efiEtaSingPPhys->Divide(etasinglePPhys,etageantPPhys);
36efiEtaSingPPhys->SetMinimum(0.);
37efiEtaSingPPhys->SetMaximum(1.2);
38efiEtaSingPPhys->SetAxisRange(-1.5,1.5);
39efiEtaSingPPhys->SetMarkerStyle(20);
40efiEtaSingPPhys->SetMarkerColor(2);
41
42
43
44TH1* efiEtaSingNPhys=(TH1*) etageantNPhys->Clone();
45efiEtaSingNPhys->Sumw2();
46efiEtaSingNPhys->Divide(etasingleNPhys,etageantNPhys);
47efiEtaSingNPhys->SetMinimum(0.);
48efiEtaSingNPhys->SetMaximum(1.2);
49efiEtaSingNPhys->SetAxisRange(-1.5,1.5);
50efiEtaSingNPhys->SetMarkerStyle(20);
51efiEtaSingNPhys->SetMarkerColor(2);
52
53
54TH1* efiPhiSingPPhys=(TH1*) phigeantPPhys->Clone();
55efiPhiSingPPhys->Sumw2();
56efiPhiSingPPhys->Divide(phisinglePPhys,phigeantPPhys);
57efiPhiSingPPhys->SetMinimum(0.);
58efiPhiSingPPhys->SetMaximum(1.2);
59efiPhiSingPPhys->SetAxisRange(-TMath::Pi(),TMath::Pi());
60efiPhiSingPPhys->SetMarkerStyle(20);
61efiPhiSingPPhys->SetMarkerColor(2);
62
63
64TH1* efiPhiSingNPhys=(TH1*) phigeantNPhys->Clone();
65efiPhiSingNPhys->Sumw2();
66efiPhiSingNPhys->Divide(phisingleNPhys,phigeantNPhys);
67efiPhiSingNPhys->SetMinimum(0.);
68efiPhiSingNPhys->SetMaximum(1.2);
69efiPhiSingNPhys->SetAxisRange(-TMath::Pi(),TMath::Pi());
70efiPhiSingNPhys->SetMarkerStyle(20);
71efiPhiSingNPhys->SetMarkerColor(2);
72
73
74
75// reconstructable tracks
76
77// track Lenght
78sparseAL->GetAxis(8)->SetRangeUser(50,200);
79sparseAL->GetAxis(12)->SetRangeUser(50,200);
80// momentum
81sparseAL->GetAxis(5)->SetRangeUser(0.,3.9);
82sparseAL->GetAxis(9)->SetRangeUser(0.,3.9);
83
84
85TH1* rgeantTrk=sparseAL->Projection(3);
86TH1* ptgeantPTrk=sparseAL->Projection(5);
87TH1* etageantPTrk=sparseAL->Projection(6);
88TH1* phigeantPTrk=sparseAL->Projection(7);
89
90
91TH1* ptgeantNTrk=sparseAL->Projection(9);
92TH1* etageantNTrk=sparseAL->Projection(10);
93TH1* phigeantNTrk=sparseAL->Projection(11);
94
95 sparseAL->GetAxis(20)->SetRangeUser(0,2.);
96 sparseAL->GetAxis(28)->SetRangeUser(0,2.);
97
98TH1* rsingleTrk=sparseAL->Projection(3);
99TH1* ptsinglePTrk=sparseAL->Projection(5);
100TH1* etasinglePTrk=sparseAL->Projection(6);
101TH1* phisinglePTrk=sparseAL->Projection(7);
102
103
104
105
106TH1* ptsingleNTrk=sparseAL->Projection(9);
107TH1* etasingleNTrk=sparseAL->Projection(10);
108TH1* phisingleNTrk=sparseAL->Projection(11);
109
110
111 sparseAL->GetAxis(32)->SetRangeUser(0,2.);
112 sparseAL->GetAxis(36)->SetRangeUser(0,2.);
113
114TH1* rV0Trk=sparseAL->Projection(3);
115TH1* ptV0PTrk=sparseAL->Projection(5);
116TH1* etaV0PTrk=sparseAL->Projection(6);
117TH1* phiV0PTrk=sparseAL->Projection(7);
118
119
120
121TH1* ptV0NTrk=sparseAL->Projection(9);
122TH1* etaV0NTrk=sparseAL->Projection(10);
123TH1* phiV0NTrk=sparseAL->Projection(11);
124
125
126TH1* efiSingPPhys=(TH1*) ptgeantPPhys->Clone();
127efiSingPPhys->Sumw2();
128efiSingPPhys->Divide(ptgeantPTrk,ptgeantPPhys);
129efiSingPPhys->SetAxisRange(0.,5);
130efiSingPPhys->SetMinimum(0.);
131efiSingPPhys->SetMaximum(1.2);
132efiSingPPhys->SetMarkerStyle(20);
133efiSingPPhys->SetMarkerColor(2);
134
135
136TH1* efiSingNPhys=(TH1*) ptgeantNPhys->Clone();
137efiSingNPhys->Sumw2();
138efiSingNPhys->Divide(ptgeantNTrk,ptgeantNPhys);
139efiSingNPhys->SetAxisRange(0.,5);
140efiSingNPhys->SetMinimum(0.);
141efiSingNPhys->SetMaximum(1.2);
142efiSingNPhys->SetMarkerStyle(20);
143efiSingNPhys->SetMarkerColor(2);
144
145
146
147
148TH1* efiSingPTrk=(TH1*) ptgeantPTrk->Clone();
149efiSingPTrk->Sumw2();
150efiSingPTrk->Divide(ptsinglePTrk,ptgeantPTrk);
151efiSingPTrk->SetAxisRange(0.,5);
152efiSingPTrk->SetMinimum(0.);
153efiSingPTrk->SetMaximum(1.2);
154efiSingPTrk->SetMarkerStyle(21);
155efiSingPTrk->SetMarkerColor(4);
156
157
158TH1* efiSingNTrk=(TH1*) ptgeantNTrk->Clone();
159efiSingNTrk->Sumw2();
160efiSingNTrk->Divide(ptsingleNTrk,ptgeantNTrk);
161efiSingNTrk->SetAxisRange(0.,5);
162efiSingNTrk->SetMinimum(0.);
163efiSingNTrk->SetMaximum(1.2);
164efiSingNTrk->SetMarkerStyle(21);
165efiSingNTrk->SetMarkerColor(4);
166
167
168TH1* efiV0PTrk=(TH1*) ptsinglePTrk->Clone();
169efiV0PTrk->Sumw2();
170efiV0PTrk->Divide(ptV0PTrk,ptsinglePTrk);
171efiV0PTrk->SetAxisRange(0.,5);
172efiV0PTrk->SetMinimum(0.);
173efiV0PTrk->SetMaximum(1.2);
174efiV0PTrk->SetMarkerStyle(21);
175efiV0PTrk->SetMarkerColor(3);
176
177
178TH1* efiV0NTrk=(TH1*) ptsingleNTrk->Clone();
179efiV0NTrk->Sumw2();
180efiV0NTrk->Divide(ptV0NTrk,ptsingleNTrk);
181efiV0NTrk->SetAxisRange(0.,5);
182efiV0NTrk->SetMinimum(0.);
183efiV0NTrk->SetMaximum(1.2);
184efiV0NTrk->SetMarkerStyle(21);
185efiV0NTrk->SetMarkerColor(3);
186
187
188
189TH1* efiEtaSingPTrk=(TH1*) etageantPTrk->Clone();
190efiEtaSingPTrk->Sumw2();
191efiEtaSingPTrk->Divide(etasinglePTrk,etageantPTrk);
192efiEtaSingPTrk->SetMinimum(0.);
193efiEtaSingPTrk->SetMaximum(1.2);
194efiEtaSingPTrk->SetAxisRange(-1.5,1.5);
195efiEtaSingPTrk->SetMarkerStyle(21);
196efiEtaSingPTrk->SetMarkerColor(4);
197
198
199TH1* efiEtaSingNTrk=(TH1*) etageantNTrk->Clone();
200efiEtaSingNTrk->Sumw2();
201efiEtaSingNTrk->Divide(etasingleNTrk,etageantNTrk);
202efiEtaSingNTrk->SetMinimum(0.);
203efiEtaSingNTrk->SetMaximum(1.2);
204efiEtaSingNTrk->SetAxisRange(-1.5,1.5);
205efiEtaSingNTrk->SetMarkerStyle(21);
206efiEtaSingNTrk->SetMarkerColor(4);
207
208
209TH1* efiEtaV0PTrk=(TH1*) etasinglePTrk->Clone();
210efiEtaV0PTrk->Sumw2();
211efiEtaV0PTrk->Divide(etaV0PTrk,etasinglePTrk);
212efiEtaV0PTrk->SetMinimum(0.);
213efiEtaV0PTrk->SetMaximum(1.2);
214efiEtaV0PTrk->SetMarkerStyle(21);
215efiEtaV0PTrk->SetMarkerColor(3);
216
217
218TH1* efiEtaV0NTrk=(TH1*) etasingleNTrk->Clone();
219efiEtaV0NTrk->Sumw2();
220efiEtaV0NTrk->Divide(etaV0NTrk,etasingleNTrk);
221efiEtaV0NTrk->SetMinimum(0.);
222efiEtaV0NTrk->SetMaximum(1.2);
223efiEtaV0NTrk->SetMarkerStyle(21);
224efiEtaV0NTrk->SetMarkerColor(3);
225
226
227
228TH1* efiPhiSingPTrk=(TH1*) phigeantPTrk->Clone();
229efiPhiSingPTrk->Sumw2();
230efiPhiSingPTrk->Divide(phisinglePTrk,phigeantPTrk);
231efiPhiSingPTrk->SetMinimum(0.);
232efiPhiSingPTrk->SetMaximum(1.2);
233efiPhiSingPTrk->SetAxisRange(-TMath::Pi(),TMath::Pi());
234efiPhiSingPTrk->SetMarkerStyle(20);
235efiPhiSingPTrk->SetMarkerColor(4);
236
237
238TH1* efiPhiSingNTrk=(TH1*) phigeantNTrk->Clone();
239efiPhiSingNTrk->Sumw2();
240efiPhiSingNTrk->Divide(phisingleNTrk,phigeantNTrk);
241efiPhiSingNTrk->SetMinimum(0.);
242efiPhiSingNTrk->SetMaximum(1.2);
243efiPhiSingNTrk->SetAxisRange(-TMath::Pi(),TMath::Pi());
244efiPhiSingNTrk->SetMarkerStyle(20);
245efiPhiSingNTrk->SetMarkerColor(4);
246
247
248
249
250TH1* efiPhiV0PTrk=(TH1*) phisinglePTrk->Clone();
251efiPhiV0PTrk->Sumw2();
252efiPhiV0PTrk->Divide(phiV0PTrk,phisinglePTrk);
253efiPhiV0PTrk->SetMinimum(0.);
254efiPhiV0PTrk->SetMaximum(1.2);
255efiPhiV0PTrk->SetMarkerStyle(21);
256efiPhiV0PTrk->SetMarkerColor(3);
257
258
259TH1* efiPhiV0NTrk=(TH1*) phisingleNTrk->Clone();
260efiPhiV0NTrk->Sumw2();
261efiPhiV0NTrk->Divide(phiV0NTrk,phisingleNTrk);
262efiPhiV0NTrk->SetMinimum(0.);
263efiPhiV0NTrk->SetMaximum(1.2);
264efiPhiV0NTrk->SetMarkerStyle(21);
265efiPhiV0NTrk->SetMarkerColor(3);
266
267
268
269sparseAL->GetAxis(32)->SetRangeUser(-2,0.);
270sparseAL->GetAxis(36)->SetRangeUser(-2,0.);
271
272
273
274TH1* sigbXYZEPlusRecSingle=sparseAL->Projection(17);
275TH1* sigbXYZEMinusRecSingle=sparseAL->Projection(25);
276
277TH1* ptsinglePTrkNo=sparseAL->Projection(5);
278TH1* etasinglePTrkNo=sparseAL->Projection(6);
279TH1* phisinglePTrkNo=sparseAL->Projection(7);
280
281TH1* ptsingleNTrkNo=sparseAL->Projection(9);
282TH1* etasingleNTrkNo=sparseAL->Projection(10);
283TH1* phisingleNTrkNo=sparseAL->Projection(11);
284
285
286
287
288
289TH1* efiV0PTrkNo=(TH1*) ptsinglePTrkNo->Clone();
290efiV0PTrkNo->Sumw2();
291efiV0PTrkNo->Divide(ptsinglePTrkNo,ptsinglePTrk);
292efiV0PTrkNo->SetAxisRange(0.,5);
293efiV0PTrkNo->SetMinimum(0.);
294efiV0PTrkNo->SetMaximum(1.2);
295efiV0PTrkNo->SetMarkerStyle(21);
296efiV0PTrkNo->SetMarkerColor(6);
297
298
299
300TH1* efiV0NTrkNo=(TH1*) ptsingleNTrkNo->Clone();
301efiV0NTrkNo->Sumw2();
302efiV0NTrkNo->Divide(ptsingleNTrkNo,ptsingleNTrk);
303efiV0NTrkNo->SetAxisRange(0.,5);
304efiV0NTrkNo->SetMinimum(0.);
305efiV0NTrkNo->SetMaximum(1.2);
306efiV0NTrkNo->SetMarkerStyle(21);
307efiV0NTrkNo->SetMarkerColor(6);
308
309
310TH1* efiEtaV0PTrkNo=(TH1*) etasinglePTrkNo->Clone();
311efiEtaV0PTrkNo->Sumw2();
312efiEtaV0PTrkNo->Divide(etasinglePTrkNo,etasinglePTrk);
313efiEtaV0PTrkNo->SetAxisRange(-1.5,1.5);
314efiEtaV0PTrkNo->SetMinimum(0.);
315efiEtaV0PTrkNo->SetMaximum(1.2);
316efiEtaV0PTrkNo->SetMarkerStyle(21);
317efiEtaV0PTrkNo->SetMarkerColor(6);
318
319
320
321TH1* efiEtaV0NTrkNo=(TH1*) etasingleNTrkNo->Clone();
322efiEtaV0NTrkNo->Sumw2();
323efiEtaV0NTrkNo->Divide(etasingleNTrkNo,etasingleNTrk);
324efiEtaV0NTrkNo->SetAxisRange(-1.5,1.5);
325efiEtaV0NTrkNo->SetMinimum(0.);
326efiEtaV0NTrkNo->SetMaximum(1.2);
327efiEtaV0NTrkNo->SetMarkerStyle(21);
328efiEtaV0NTrkNo->SetMarkerColor(6);
329
330
331TH1* efiPhiV0PTrkNo=(TH1*) phisinglePTrkNo->Clone();
332efiPhiV0PTrkNo->Sumw2();
333efiPhiV0PTrkNo->Divide(phisinglePTrkNo,phisinglePTrk);
334efiPhiV0PTrkNo->SetAxisRange(-3.5,3.5);
335efiPhiV0PTrkNo->SetMinimum(0.);
336efiPhiV0PTrkNo->SetMaximum(1.2);
337efiPhiV0PTrkNo->SetMarkerStyle(21);
338efiPhiV0PTrkNo->SetMarkerColor(6);
339
340
341
342TH1* efiPhiV0NTrkNo=(TH1*) phisingleNTrkNo->Clone();
343efiPhiV0NTrkNo->Sumw2();
344efiPhiV0NTrkNo->Divide(phisingleNTrkNo,phisingleNTrk);
345efiPhiV0NTrkNo->SetAxisRange(-3.5,3.5);
346efiPhiV0NTrkNo->SetMinimum(0.);
347efiPhiV0NTrkNo->SetMaximum(1.2);
348efiPhiV0NTrkNo->SetMarkerStyle(21);
349efiPhiV0NTrkNo->SetMarkerColor(6);
350
351
352
353
354
355 TCanvas* c1 = new TCanvas("c1","MC Info",200,10,700,900);
356 c1->Divide(2,4);
357 c1->SetFillColor(0);
358 c1->GetFrame()->SetFillColor(0);
359 c1->SetBorderMode(0);
360 TPad *c1_1 = (TPad*)(c1->GetPrimitive("c1_1"));
361 TPad *c1_2 = (TPad*)(c1->GetPrimitive("c1_2"));
362 TPad *c1_3 = (TPad*)(c1->GetPrimitive("c1_3"));
363 TPad *c1_4 = (TPad*)(c1->GetPrimitive("c1_4"));
364 TPad *c1_5 = (TPad*)(c1->GetPrimitive("c1_5"));
365 TPad *c1_6 = (TPad*)(c1->GetPrimitive("c1_6"));
366
367
368 c1_1->SetGridx(1);
369 c1_3->SetGridx(1);
370 c1_2->SetGridx(1);
371 c1_4->SetGridx(1);
372 c1_5->SetGridx(1);
373 c1_6->SetGridx(1);
374
375 c1_1->SetGridy(1);
376 c1_3->SetGridy(1);
377 c1_2->SetGridy(1);
378 c1_4->SetGridy(1);
379 c1_5->SetGridy(1);
380 c1_6->SetGridy(1);
381
382
383c1->cd(1);
384efiSingPTrk->Draw();
385efiSingPTrk->SetXTitle("1/sqrt(pt) pos (sqrt(GeV))");
386efiSingPPhys->Draw("same");
387efiV0PTrk->Draw("same");
388efiV0PTrkNo->Draw("same");
389
390TLegend *leg1 = new TLegend(0.6,0.75,0.9,0.9);
391leg1->AddEntry(efiSingPTrk,"Single tracking eff");
392leg1->AddEntry(efiV0PTrk,"V0 tracking eff");
393leg1->AddEntry(efiV0PTrkNo,"V0 tracking ineff");
394leg1->Draw();
395c1->cd(2);
396efiSingNTrk->Draw();
397efiSingNTrk->SetXTitle("1/sqrt(pt) neg (sqrt(GeV))");
398efiSingNPhys->Draw("same");
399efiV0NTrk->Draw("same");
400efiV0NTrkNo->Draw("same");
401
402c1->cd(3);
403efiEtaSingPTrk->Draw();
404//efiEtaSingPPhys->Draw("same");
405efiEtaV0PTrk->Draw("same");
406efiEtaV0PTrkNo->Draw("same");
407
408c1->cd(4);
409efiEtaSingNTrk->Draw();
410//efiEtaSingNPhys->Draw("same");
411efiEtaV0NTrk->Draw("same");
412efiEtaV0NTrkNo->Draw("same");
413
414c1->cd(5);
415efiPhiSingPTrk->Draw();
416//efiPhiSingPPhys->Draw("same");
417efiPhiV0PTrk->Draw("same");
418efiPhiV0PTrkNo->Draw("same");
419
420c1->cd(6);
421efiPhiSingNTrk->Draw();
422//efiPhiSingNPhys->Draw("same");
423efiPhiV0NTrk->Draw("same");
424efiPhiV0NTrkNo->Draw("same");
425
426c1->cd(7);
427TH1* efiRSingTrk=(TH1*) rgeantTrk->Clone();
428efiRSingTrk->Sumw2();
429efiRSingTrk->Divide(rsingleTrk,rgeantTrk);
430efiRSingTrk->SetMinimum(0.);
431efiRSingTrk->SetMaximum(1.2);
432efiRSingTrk->SetMarkerStyle(21);
433efiRSingTrk->SetMarkerColor(4);
434efiRSingTrk->Draw();
435
436TH1* efiRV0Trk=(TH1*) rgeantTrk->Clone();
437efiRV0Trk->Sumw2();
438efiRV0Trk->Divide(rV0Trk,rsingleTrk);
439efiRV0Trk->SetMinimum(0.);
440efiRV0Trk->SetMaximum(1.2);
441efiRV0Trk->SetMarkerStyle(21);
442efiRV0Trk->SetMarkerColor(3);
443efiRV0Trk->Draw("same");
444
445
446c1->Print("Effi-AL-norefit-v416rev08_090709.ps");
447//c1->Print("Effi-AL-norefit-trunk010709.ps");
448
449}