AliPhysics  vAN-20150822 (d56cf94)
 All Classes Namespaces Files Functions Variables Enumerations Enumerator Macros
ReadCFHeavyFlavourOutput.C
Go to the documentation of this file.
1 #include <Riostream.h>
2 
3 extern TRandom *gRandom;
4 extern TBenchmark *gBenchmark;
5 extern TSystem *gSystem;
6 
8 
9  // example macro to read the CF produced for D->Kpi
10  // 8 steps considered
11  // 12 variables
12 
13  gROOT->SetStyle("Plain");
14  gStyle->SetPalette(1);
15  gStyle->SetOptStat(0);
16  gStyle->SetPalette(1);
17  gStyle->SetCanvasColor(0);
18  gStyle->SetFrameFillColor(0);
19  gStyle->SetOptTitle(0);
20 
21  gSystem->SetIncludePath("-I. -I$ALICE_ROOT/include -I$ROOTSYS/include");
22  gSystem->Load("libANALYSIS");
23  gSystem->Load("libANALYSISalice");
24  gSystem->Load("$ALICE_ROOT/CORRFW/libCORRFW") ;
25 
26  // Setting up the container grid...
27 
28  // Flag the sel steps
29  Int_t stepGen=0;
30  Int_t stepGenAcc=1;
31  Int_t stepGenVertex=2;
32  Int_t stepGenRefit=3;
33  Int_t stepRec=4;
34  Int_t stepRecAcc=5;
35  Int_t stepRecITSClus=6;
36  Int_t stepRecPPR=7;
37 
38  // The sensitive variables, their indeces
39  Int_t ipt =0;
40  Int_t iy=1;
41  UInt_t icTS = 2; //cTS stands for cosThetaStar
42  UInt_t iptPi = 3;
43  UInt_t iptK = 4;
44  UInt_t icT = 5;
45  UInt_t idca = 6;
46  UInt_t id0pi = 7;
47  UInt_t id0K = 8;
48  UInt_t id0xd0 = 9;
49  UInt_t ipointing = 10;
50  UInt_t iphi = 11;
51 
52  // Get the container
53  TFile* f = new TFile("AnalysisResults.root");
54  TDirectoryFile* d = (TDirectoryFile*)f->Get("PWG3_D2H_CFtaskD0toKpi");
55  AliCFContainer *data = (AliCFContainer*) (d->Get("container"));
56 
57  // Make some 1 & 2-D projections.
58  // MC level
59  TCanvas *cmc1 =new TCanvas("cmc1","The MC distributions",0,0,1600,1200); // canvas for all steps, 1st three variables
60  TCanvas *cmc2 =new TCanvas("cmc2","The MC distributions",0,0,1600,1200); // canvas for all steps, 2nd three variables
61 
62  TCanvas *cmc3 =new TCanvas("cmc3","The MC distributions",0,0,1600,1200); // canvas for all steps, 2nd three variables
63  TCanvas *cmc4 =new TCanvas("cmc4","The MC distributions",0,0,1600,765); // canvas for all steps, 2nd three variables
64  TCanvas *cmcpt = new TCanvas("cmcpt","pt distribution from MC",50,50,550,550);
65  TCanvas *cmcy = new TCanvas("cmcy","y distribution from MC",50,50,550,550);
66  TCanvas *cmccTS = new TCanvas("cmcTS","cosThetaStar distribution from MC",50,50,550,550); // cTS stands for cosThetaStar
67  TCanvas *cmcptPi = new TCanvas("cmcptPi","pt_pi distribution from MC",50,50,550,550);
68  TCanvas *cmcptK = new TCanvas("cmcptK","pt_K distribution from MC",50,50,550,550);
69  TCanvas *cmccT = new TCanvas("cmccT","cT distribution from MC",50,50,550,550);
70  TCanvas *cmcdca = new TCanvas("cmcdca","dca distribution from MC",50,50,550,550);
71  TCanvas *cmcd0pi = new TCanvas("cmcd0pi","d0pi distribution from MC",50,50,550,550);
72  TCanvas *cmcd0K = new TCanvas("cmcd0K","d0K distribution from MC",50,50,550,550);
73  TCanvas *cmcd0xd0 = new TCanvas("cmcd0xd0","d0xd0 distribution from MC",50,50,550,550);
74  TCanvas *cmcpointing = new TCanvas("cmcpointing","cosPointingAngle distribution from MC",50,50,550,550);
75  TCanvas *cmcphi = new TCanvas("cmcphi","Phi distribution from MC",50,50,550,550);
76 
77  // MC Acceptance level
78  TCanvas *cmcAccpt = new TCanvas("cmcAccpt","pt distribution from MC - Acceptance cuts",50,50,550,550);
79  TCanvas *cmcAccy = new TCanvas("cmcAccy","y distribution from MC - Acceptance cuts",50,50,550,550);
80  TCanvas *cmcAcccTS = new TCanvas("cmcAccTS","cosThetaStar distribution from MC - Acceptance cuts",50,50,550,550); // cTS stands for cosThetaStar
81  TCanvas *cmcAccptPi = new TCanvas("cmcAccptPi","pt_pi distribution from MC - Acceptance cuts",50,50,550,550);
82  TCanvas *cmcAccptK = new TCanvas("cmcAccptK","pt_K distribution from MC - Acceptance cuts",50,50,550,550);
83  TCanvas *cmcAcccT = new TCanvas("cmcAcccT","cT distribution from MC - Acceptance cuts",50,50,550,550);
84  TCanvas *cmcAccdca = new TCanvas("cmcAccdca","dca distribution from MC - Acceptance cuts",50,50,550,550);
85  TCanvas *cmcAccd0pi = new TCanvas("cmcAccd0pi","d0pi distribution from MC - Acceptance cuts",50,50,550,550);
86  TCanvas *cmcAccd0K = new TCanvas("cmcAccd0K","d0K distribution from MC - Acceptance cuts",50,50,550,550);
87  TCanvas *cmcAccd0xd0 = new TCanvas("cmcAccd0xd0","d0xd0 distribution from MC - Acceptance cuts",50,50,550,550);
88  TCanvas *cmcAccpointing = new TCanvas("cmcAccpointing","cosPointingAngle distribution from MC - Acceptance cuts",50,50,550,550);
89  TCanvas *cmcAccphi = new TCanvas("cmcAccphi","Phi distribution from MC - Acceptance cuts",50,50,550,550);
90 
91  // MC Vertex level
92  TCanvas *cmcVertexpt = new TCanvas("cmcVertexpt","pt distribution from MC - Vertex cuts",50,50,550,550);
93  TCanvas *cmcVertexy = new TCanvas("cmcVertexy","y distribution from MC - Vertex cuts",50,50,550,550);
94  TCanvas *cmcVertexcTS = new TCanvas("cmcVertexTS","cosThetaStar distribution from MC - Vertex cuts",50,50,550,550); // cTS stands for cosThetaStar
95  TCanvas *cmcVertexptPi = new TCanvas("cmcVertexptPi","pt_pi distribution from MC - Vertex cuts",50,50,550,550);
96  TCanvas *cmcVertexptK = new TCanvas("cmcVertexptK","pt_K distribution from MC - Vertex cuts",50,50,550,550);
97  TCanvas *cmcVertexcT = new TCanvas("cmcVertexcT","cT distribution from MC - Vertex cuts",50,50,550,550);
98  TCanvas *cmcVertexdca = new TCanvas("cmcVertexdca","dca distribution from MC - Vertex cuts",50,50,550,550);
99  TCanvas *cmcVertexd0pi = new TCanvas("cmcVertexd0pi","d0pi distribution from MC - Vertex cuts",50,50,550,550);
100  TCanvas *cmcVertexd0K = new TCanvas("cmcVertexd0K","d0K distribution from MC - Vertex cuts",50,50,550,550);
101  TCanvas *cmcVertexd0xd0 = new TCanvas("cmcVertexd0xd0","d0xd0 distribution from MC - Vertex cuts",50,50,550,550);
102  TCanvas *cmcVertexpointing = new TCanvas("cmcVertexpointing","cosPointingAngle distribution from MC - Vertex cuts",50,50,550,550);
103  TCanvas *cmcVertexphi = new TCanvas("cmcVertexphi","Phi distribution from MC - Vertex cuts",50,50,550,550);
104 
105  // MC Refit level
106  TCanvas *cmcRefitpt = new TCanvas("cmcRefitpt","pt distribution from MC - Refit cuts",50,50,550,550);
107  TCanvas *cmcRefity = new TCanvas("cmcRefity","y distribution from MC - Refit cuts",50,50,550,550);
108  TCanvas *cmcRefitcTS = new TCanvas("cmcRefitTS","cosThetaStar distribution from MC - Refit cuts",50,50,550,550); // cTS stands for cosThetaStar
109  TCanvas *cmcRefitptPi = new TCanvas("cmcRefitptPi","pt_pi distribution from MC - Refit cuts",50,50,550,550);
110  TCanvas *cmcRefitptK = new TCanvas("cmcRefitptK","pt_K distribution from MC - Refit cuts",50,50,550,550);
111  TCanvas *cmcRefitcT = new TCanvas("cmcRefitcT","cT distribution from MC - Refit cuts",50,50,550,550);
112  TCanvas *cmcRefitdca = new TCanvas("cmcRefitdca","dca distribution from MC - Refit cuts",50,50,550,550);
113  TCanvas *cmcRefitd0pi = new TCanvas("cmcRefitd0pi","d0pi distribution from MC - Refit cuts",50,50,550,550);
114  TCanvas *cmcRefitd0K = new TCanvas("cmcRefitd0K","d0K distribution from MC - Refit cuts",50,50,550,550);
115  TCanvas *cmcRefitd0xd0 = new TCanvas("cmcRefitd0xd0","d0xd0 distribution from MC - Refit cuts",50,50,550,550);
116  TCanvas *cmcRefitpointing = new TCanvas("cmcRefitpointing","cosPointingAngle distribution from MC - Refit cuts",50,50,550,550);
117  TCanvas *cmcRefitphi = new TCanvas("cmcRefitphi","Phi distribution from MC - Refit cuts",50,50,550,550);
118 
119  // Reco-aod level
120  TCanvas *cpt = new TCanvas("cpt","pt distribution from reco aod",50,50,550,550);
121  TCanvas *cy = new TCanvas("cy","y distribution from reco aod",50,50,550,550);
122  TCanvas *ccTS = new TCanvas("cTS","cosThetaStar distribution from reco aod",50,50,550,550);
123  TCanvas *cptPi = new TCanvas("cptPi","pt_pi distribution from reco aod",50,50,550,550);
124  TCanvas *cptK = new TCanvas("cptK","pt_K distribution from reco aod",50,50,550,550);
125  TCanvas *ccT = new TCanvas("ccT","cT distribution from reco aod",50,50,550,550);
126  TCanvas *cdca = new TCanvas("cdca","dca distribution from reco aod",50,50,550,550);
127  TCanvas *cd0pi = new TCanvas("cd0pi","d0pi distribution from reco aod",50,50,550,550);
128  TCanvas *cd0K = new TCanvas("cd0K","d0K distribution from reco aod",50,50,550,550);
129  TCanvas *cd0xd0 = new TCanvas("cd0xd0","d0xd0 distribution from reco aod",50,50,550,550);
130  TCanvas *cpointing = new TCanvas("cpointing","cosPointingAngle distribution from reco aod",50,50,550,550);
131  TCanvas *cphi = new TCanvas("cphi","Phi distribution from reco aod",50,50,550,550);
132 
133  // Reco-aod Acceptance level
134  TCanvas *cAccpt = new TCanvas("cAccpt","pt distribution from reco aod - Acceptance cuts",50,50,550,550);
135  TCanvas *cAccy = new TCanvas("cAccy","y distribution from reco aod - Acceptance cuts",50,50,550,550);
136  TCanvas *cAcccTS = new TCanvas("cAccTS","cosThetaStar distribution from reco aod - Acceptance cuts",50,50,550,550);
137  TCanvas *cAccptPi = new TCanvas("cAccptPi","pt_pi distribution from reco aod - Acceptance cuts",50,50,550,550);
138  TCanvas *cAccptK = new TCanvas("cAccptK","pt_K distribution from reco aod - Acceptance cuts",50,50,550,550);
139  TCanvas *cAccdca = new TCanvas("cAccdca","dca distribution from reco aod - Acceptance cuts",50,50,550,550);
140  TCanvas *cAccd0pi = new TCanvas("cAccd0pi","d0pi distribution from reco aod - Acceptance cuts",50,50,550,550);
141  TCanvas *cAccd0K = new TCanvas("cAccd0K","d0K distribution from reco aod - Acceptance cuts",50,50,550,550);
142  TCanvas *cAccd0xd0 = new TCanvas("cAccd0xd0","d0xd0 distribution from reco aod - Acceptance cuts",50,50,550,550);
143  TCanvas *cAccpointing = new TCanvas("cAccpointing","cosPointingAngle distribution from reco aod - Acceptance cuts",50,50,550,550);
144  TCanvas *cAcccT = new TCanvas("cAcccT","cT distribution from reco aod - Acceptance cuts",50,50,550,550);
145  TCanvas *cAccphi = new TCanvas("cAccphi","Phi distribution from reco aod - Acceptance cuts",50,50,550,550);
146 
147  // Reco-aod ITS Clusters level
148  TCanvas *cITSClpt = new TCanvas("cITSClpt","pt distribution from reco aod - ITS Clusters cuts",50,50,550,550);
149  TCanvas *cITSCly = new TCanvas("cITSCly","y distribution from reco aod - ITS Clusters cuts",50,50,550,550);
150  TCanvas *cITSClcTS = new TCanvas("cITSClTS","cosThetaStar distribution from reco aod - ITS Clusters cuts",50,50,550,550);
151  TCanvas *cITSClptPi = new TCanvas("cITSClptPi","pt_pi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
152  TCanvas *cITSClptK = new TCanvas("cITSClptK","pt_K distribution from reco aod - ITS Clusters cuts",50,50,550,550);
153  TCanvas *cITSClcT = new TCanvas("cITSClcT","cT distribution from reco aod - ITS Clusters cuts",50,50,550,550);
154  TCanvas *cITSCldca = new TCanvas("cITSCldca","dca distribution from reco aod - ITS Clusters cuts",50,50,550,550);
155  TCanvas *cITSCld0pi = new TCanvas("cITSCld0pi","d0pi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
156  TCanvas *cITSCld0K = new TCanvas("cITSCld0K","d0K distribution from reco aod - ITS Clusters cuts",50,50,550,550);
157  TCanvas *cITSCld0xd0 = new TCanvas("cITSCld0xd0","d0xd0 distribution from reco aod - ITS Clusters cuts",50,50,550,550);
158  TCanvas *cITSClpointing = new TCanvas("cITSClpointing","cosPointingAngle distribution from reco aod - ITS Clusters cuts",50,50,550,550);
159  TCanvas *cITSClphi = new TCanvas("cITSClphi","Phi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
160 
161  // Reco-aod PPR level
162  TCanvas *cPPRpt = new TCanvas("cPPRpt","pt distribution from reco aod - PPR cuts",50,50,550,550);
163  TCanvas *cPPRy = new TCanvas("cPPRy","y distribution from reco aod - PPR cuts",50,50,550,550);
164  TCanvas *cPPRcTS = new TCanvas("cPPRTS","cosThetaStar distribution from reco aod - PPR cuts",50,50,550,550);
165  TCanvas *cPPRptPi = new TCanvas("cPPRptPi","pt_pi distribution from reco aod - PPR cuts",50,50,550,550);
166  TCanvas *cPPRptK = new TCanvas("cPPRptK","pt_K distribution from reco aod - PPR cuts",50,50,550,550);
167  TCanvas *cPPRcT = new TCanvas("cPPRcT","cT distribution from reco aod - PPR cuts",50,50,550,550);
168  TCanvas *cPPRdca = new TCanvas("cPPRdca","dca distribution from reco aod - PPR cuts",50,50,550,550);
169  TCanvas *cPPRd0pi = new TCanvas("cPPRd0pi","d0pi distribution from reco aod - PPR cuts",50,50,550,550);
170  TCanvas *cPPRd0K = new TCanvas("cPPRd0K","d0K distribution from reco aod - PPR cuts",50,50,550,550);
171  TCanvas *cPPRd0xd0 = new TCanvas("cPPRd0xd0","d0xd0 distribution from reco aod - PPR cuts",50,50,550,550);
172  TCanvas *cPPRpointing = new TCanvas("cPPRpointing","cosPointingAngle distribution from reco aod - PPR cuts",50,50,550,550);
173  TCanvas *cPPRphi = new TCanvas("cPPRphi","Phi distribution from reco aod - PPR cuts",50,50,550,550);
174 
175  // pt, y, cosThetaStar for all steps
176  cmc1->Divide(8,3);
177 
178  cmc1->cd(1);
179  TH1D *hMCpt1D = data->ShowProjection(ipt, stepGen);
180  Double_t maxpt = hMCpt1D->GetMaximum();
181  hMCpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
182  hMCpt1D->Sumw2();
183  hMCpt1D->SetMinimum(0.01);
184  SetHistoProperties(hMCpt1D,2,3,2,20,2,3005);
185  hMCpt1D->GetXaxis()->SetTitleOffset(1.2);
186  hMCpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC data");
187  hMCpt1D->Draw("hist");
188  cmcpt->cd();
189  cmcpt->SetLeftMargin(0.15);
190  cmcpt->SetRightMargin(0.05);
191  hMCpt1D->Draw("hist");
192  hMCpt1D->Draw("err same");
193  cmcpt->Update();
194 
195  cmc1->cd(2);
196  TH1D *hMCAccpt1D = data->ShowProjection(ipt, stepGenAcc);
197  maxpt = hMCAccpt1D->GetMaximum();
198  hMCAccpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
199  SetHistoProperties(hMCAccpt1D,8,3,8,20,8,3006);
200  hMCAccpt1D->GetXaxis()->SetTitleOffset(1.2);;
201  hMCAccpt1D->Sumw2();
202  hMCAccpt1D->SetMinimum(0.01);
203  hMCAccpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Acc cut");
204  hMCAccpt1D->Draw("hist");
205  cmcAccpt->cd();
206  cmcAccpt->SetLeftMargin(0.15);
207  cmcAccpt->SetRightMargin(0.05);
208  hMCAccpt1D->Draw("hist");
209  hMCAccpt1D->Draw("err same");
210  cmcAccpt->Update();
211 
212  cmc1->cd(3);
213  TH1D *hMCVertexpt1D = data->ShowProjection(ipt, stepGenVertex);
214  maxpt = hMCVertexpt1D->GetMaximum();
215  hMCVertexpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
216  SetHistoProperties(hMCVertexpt1D,49,3,49,20,49,3006);
217  hMCVertexpt1D->GetXaxis()->SetTitleOffset(1.2);;
218  hMCVertexpt1D->Sumw2();
219  hMCVertexpt1D->SetMinimum(0.01);
220  hMCVertexpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Vertex cut");
221  hMCVertexpt1D->Draw("hist");
222  cmcVertexpt->cd();
223  cmcVertexpt->SetLeftMargin(0.15);
224  cmcVertexpt->SetRightMargin(0.05);
225  hMCVertexpt1D->Draw("hist");
226  hMCVertexpt1D->Draw("err same");
227  cmcVertexpt->Update();
228 
229  cmc1->cd(4);
230  TH1D *hMCRefitpt1D = data->ShowProjection(ipt, stepGenRefit);
231  maxpt = hMCRefitpt1D->GetMaximum();
232  hMCRefitpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
233  SetHistoProperties(hMCRefitpt1D,42,3,42,20,42,3006);
234  hMCRefitpt1D->GetXaxis()->SetTitleOffset(1.2);;
235  hMCRefitpt1D->Sumw2();
236  hMCRefitpt1D->SetMinimum(0.01);
237  hMCRefitpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Refit cut");
238  hMCRefitpt1D->Draw("hist");
239  cmcRefitpt->cd();
240  cmcRefitpt->SetLeftMargin(0.15);
241  cmcRefitpt->SetRightMargin(0.05);
242  hMCRefitpt1D->Draw("hist");
243  hMCRefitpt1D->Draw("err same");
244  cmcRefitpt->Update();
245 
246  cmc1->cd(5);
247  TH1D *hRECpt1D = data->ShowProjection(ipt, stepRec);
248  maxpt = hRECpt1D->GetMaximum();
249  hRECpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
250  SetHistoProperties(hRECpt1D,4,3,4,20,4,3004);
251  hRECpt1D->GetXaxis()->SetTitleOffset(1.2);;
252  hRECpt1D->Sumw2();
253  hRECpt1D->SetMinimum(0.01);
254  hRECpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD");
255  hRECpt1D->Draw("hist");
256  cpt->cd();
257  cpt->SetLeftMargin(0.15);
258  cpt->SetRightMargin(0.05);
259  hRECpt1D->Draw("hist");
260  hRECpt1D->Draw("err same");
261  cpt->Update();
262 
263  cmc1->cd(6);
264  TH1D *hRECAccpt1D = data->ShowProjection(ipt, stepRecAcc);
265  maxpt = hRECAccpt1D->GetMaximum();
266  hRECAccpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
267  SetHistoProperties(hRECAccpt1D,38,3,38,20,38,3007);
268  hRECAccpt1D->GetXaxis()->SetTitleOffset(1.2);;
269  hRECAccpt1D->Sumw2();
270  hRECAccpt1D->SetMinimum(0.01);
271  hRECAccpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - Acc cuts");
272  hRECAccpt1D->Draw("hist");
273  cAccpt->cd();
274  cAccpt->SetLeftMargin(0.15);
275  cAccpt->SetRightMargin(0.05);
276  hRECAccpt1D->Draw("hist");
277  hRECAccpt1D->Draw("err same");
278  cAccpt->Update();
279 
280  cmc1->cd(7);
281  TH1D *hRECITSClpt1D = data->ShowProjection(ipt, stepRecITSClus);
282  maxpt = hRECITSClpt1D->GetMaximum();
283  hRECITSClpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
284  SetHistoProperties(hRECITSClpt1D,13,3,13,20,13,3019);
285  hRECITSClpt1D->GetXaxis()->SetTitleOffset(1.2);;
286  hRECITSClpt1D->Sumw2();
287  hRECITSClpt1D->SetMinimum(0.01);
288  hRECITSClpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - ITS Cl cuts");
289  hRECITSClpt1D->Draw("hist");
290  cITSClpt->cd();
291  cITSClpt->SetLeftMargin(0.15);
292  cITSClpt->SetRightMargin(0.05);
293  hRECITSClpt1D->Draw("hist");
294  hRECITSClpt1D->Draw("err same");
295  cITSClpt->Update();
296 
297  cmc1->cd(8);
298  TH1D *hRECPPRpt1D = data->ShowProjection(ipt, stepRecPPR);
299  maxpt = hRECPPRpt1D->GetMaximum();
300  hRECPPRpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
301  SetHistoProperties(hRECPPRpt1D,6,3,6,20,6,3019);
302  hRECPPRpt1D->GetXaxis()->SetTitleOffset(1.2);;
303  hRECPPRpt1D->Sumw2();
304  hRECPPRpt1D->SetMinimum(0.01);
305  hRECPPRpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - PPR cuts");
306  hRECPPRpt1D->Draw("hist");
307  cPPRpt->cd();
308  cPPRpt->SetLeftMargin(0.15);
309  cPPRpt->SetRightMargin(0.05);
310  hRECPPRpt1D->Draw("hist");
311  hRECPPRpt1D->Draw("err same");
312  cPPRpt->Update();
313 
314  cmc1->cd(9);
315  TH1D *hMCy1D = data->ShowProjection(iy, stepGen);
316  Double_t maxy = hMCy1D->GetMaximum();
317  hMCy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
318  SetHistoProperties(hMCy1D,2,3,2,20,2,3005);
319  hMCy1D->SetLineColor(2);
320  hMCy1D->SetLineWidth(3);
321  hMCy1D->SetMarkerColor(2);
322  hMCy1D->SetFillColor(2);
323  hMCy1D->SetFillStyle(3005);
324  hMCy1D->SetMarkerStyle(20);
325  hMCy1D->GetXaxis()->SetTitleOffset(1.2);;
326  hMCy1D->Sumw2();
327  hMCy1D->SetMinimum(0.01);
328  hMCy1D->GetXaxis()->SetTitle("y, MC data");
329  hMCy1D->Draw("hist");
330  cmcy->cd();
331  cmcy->SetLeftMargin(0.15);
332  cmcy->SetRightMargin(0.05);
333  hMCy1D->Draw("hist");
334  hMCy1D->Draw("err same");
335  cmcy->Update();
336 
337  cmc1->cd(10);
338  TH1D *hMCAccy1D = data->ShowProjection(iy, stepGenAcc);
339  maxy = hMCAccy1D->GetMaximum();
340  hMCAccy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
341  SetHistoProperties(hMCAccy1D,8,3,8,20,8,3006);
342  hMCAccy1D->GetXaxis()->SetTitleOffset(1.2);
343  hMCAccy1D->Sumw2();
344  hMCAccy1D->SetMinimum(0.01);
345  hMCAccy1D->GetXaxis()->SetTitle("y, MC data - Acc cuts");
346  hMCAccy1D->Draw("hist");
347  cmcAccy->cd();
348  cmcAccy->SetLeftMargin(0.15);
349  cmcAccy->SetRightMargin(0.05);
350  hMCAccy1D->Draw("hist");
351  hMCAccy1D->Draw("err same");
352  cmcy->Update();
353 
354  cmc1->cd(11);
355  TH1D *hMCVertexy1D = data->ShowProjection(iy, stepGenVertex);
356  maxy = hMCVertexy1D->GetMaximum();
357  hMCVertexy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
358  SetHistoProperties(hMCVertexy1D,49,3,49,20,49,3006);
359  hMCVertexy1D->GetXaxis()->SetTitleOffset(1.2);
360  hMCVertexy1D->Sumw2();
361  hMCVertexy1D->SetMinimum(0.01);
362  hMCVertexy1D->GetXaxis()->SetTitle("y, MC data - Vertex cuts");
363  hMCVertexy1D->Draw("hist");
364  cmcVertexy->cd();
365  cmcVertexy->SetLeftMargin(0.15);
366  cmcVertexy->SetRightMargin(0.05);
367  hMCVertexy1D->Draw("hist");
368  hMCVertexy1D->Draw("err same");
369  cmcy->Update();
370 
371  cmc1->cd(12);
372  TH1D *hMCRefity1D = data->ShowProjection(iy, stepGenRefit);
373  maxy = hMCRefity1D->GetMaximum();
374  hMCRefity1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
375  SetHistoProperties(hMCRefity1D,42,3,42,20,42,3006);
376  hMCRefity1D->GetXaxis()->SetTitleOffset(1.2);
377  hMCRefity1D->Sumw2();
378  hMCRefity1D->SetMinimum(0.01);
379  hMCRefity1D->GetXaxis()->SetTitle("y, MC data - Refit cuts");
380  hMCRefity1D->Draw("hist");
381  cmcRefity->cd();
382  cmcRefity->SetLeftMargin(0.15);
383  cmcRefity->SetRightMargin(0.05);
384  hMCRefity1D->Draw("hist");
385  hMCRefity1D->Draw("err same");
386  cmcy->Update();
387 
388  cmc1->cd(13);
389  TH1D *hRECy1D = data->ShowProjection(iy, stepRec);
390  maxy = hRECy1D->GetMaximum();
391  hRECy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
392  SetHistoProperties(hRECy1D,4,3,4,20,4,3004);
393  hRECy1D->GetXaxis()->SetTitleOffset(1.2);;
394  hRECy1D->Sumw2();
395  hRECy1D->SetMinimum(0.01);
396  hRECy1D->Draw("hist");
397  cy->cd();
398  cy->SetLeftMargin(0.15);
399  cy->SetRightMargin(0.05);
400  hRECy1D->GetXaxis()->SetTitle("y, AOD");
401  hRECy1D->Draw("hist");
402  hRECy1D->Draw("err same");
403  cy->Update();
404 
405  cmc1->cd(14);
406  TH1D *hRECAccy1D = data->ShowProjection(iy, stepRecAcc);
407  maxy = hRECAccy1D->GetMaximum();
408  hRECAccy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
409  SetHistoProperties(hRECAccy1D,38,3,38,20,38,3007);
410  hRECAccy1D->GetXaxis()->SetTitleOffset(1.2);;
411  hRECAccy1D->Sumw2();
412  hRECAccy1D->SetMinimum(0.01);
413  hRECAccy1D->Draw("hist");
414  cAccy->cd();
415  cAccy->SetLeftMargin(0.15);
416  cAccy->SetRightMargin(0.05);
417  hRECAccy1D->GetXaxis()->SetTitle("y, AOD - Acc cuts");
418  hRECAccy1D->Draw("hist");
419  hRECAccy1D->Draw("err same");
420  cAccy->Update();
421 
422  cmc1->cd(15);
423  TH1D *hRECITSCly1D = data->ShowProjection(iy, stepRecITSClus);
424  maxy = hRECITSCly1D->GetMaximum();
425  hRECITSCly1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
426  SetHistoProperties(hRECITSCly1D,13,3,13,20,13,3019);
427  hRECITSCly1D->GetXaxis()->SetTitleOffset(1.2);;
428  hRECITSCly1D->Sumw2();
429  hRECITSCly1D->SetMinimum(0.01);
430  hRECITSCly1D->Draw("hist");
431  cITSCly->cd();
432  cITSCly->SetLeftMargin(0.15);
433  cITSCly->SetRightMargin(0.05);
434  hRECITSCly1D->GetXaxis()->SetTitle("y, AOD - ITS Cl cuts");
435  hRECITSCly1D->Draw("hist");
436  hRECITSCly1D->Draw("err same");
437  cITSCly->Update();
438 
439  cmc1->cd(16);
440  TH1D *hRECPPRy1D = data->ShowProjection(iy, stepRecPPR);
441  maxy = hRECPPRy1D->GetMaximum();
442  hRECPPRy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
443  SetHistoProperties(hRECPPRy1D,6,3,6,20,6,3019);
444  hRECPPRy1D->SetMarkerStyle(20);
445  hRECPPRy1D->GetXaxis()->SetTitleOffset(1.2);;
446  hRECPPRy1D->Sumw2();
447  hRECPPRy1D->SetMinimum(0.01);
448  hRECPPRy1D->Draw("hist");
449  cPPRy->cd();
450  cPPRy->SetLeftMargin(0.15);
451  cPPRy->SetRightMargin(0.05);
452  hRECPPRy1D->GetXaxis()->SetTitle("y, AOD - PPR cuts");
453  hRECPPRy1D->Draw("hist");
454  hRECPPRy1D->Draw("err same");
455  cPPRy->Update();
456 
457  cmc1->cd(17);
458  TH1D *hMCcTS1D = data->ShowProjection(icTS, stepGen);
459  Double_t maxcTS = hMCcTS1D->GetMaximum();
460  hMCcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
461  SetHistoProperties(hMCcTS1D,2,3,2,20,2,3005);
462  hMCcTS1D->GetXaxis()->SetTitleOffset(1.2);;
463  hMCcTS1D->Sumw2();
464  hMCcTS1D->SetMinimum(0.01);
465  hMCcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data");
466  hMCcTS1D->Draw("hist");
467  cmccTS->cd();
468  cmccTS->SetLeftMargin(0.15);
469  cmccTS->SetRightMargin(0.05);
470  hMCcTS1D->Draw("hist");
471  hMCcTS1D->Draw("err same");
472  cmccTS->Update();
473 
474  cmc1->cd(18);
475  TH1D *hMCAcccTS1D = data->ShowProjection(icTS, stepGenAcc);
476  maxcTS = hMCAcccTS1D->GetMaximum();
477  hMCAcccTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
478  SetHistoProperties(hMCAcccTS1D,8,3,8,20,8,3006);
479  hMCAcccTS1D->GetXaxis()->SetTitleOffset(1.2);;
480  hMCAcccTS1D->Sumw2();
481  hMCAcccTS1D->SetMinimum(0.01);
482  hMCAcccTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Acc cuts");
483  hMCAcccTS1D->Draw("hist");
484  cmcAcccTS->cd();
485  cmcAcccTS->SetLeftMargin(0.15);
486  cmcAcccTS->SetRightMargin(0.05);
487  hMCAcccTS1D->Draw("hist");
488  hMCAcccTS1D->Draw("err same");
489  cmcAcccTS->Update();
490 
491  cmc1->cd(19);
492  TH1D *hMCVertexcTS1D = data->ShowProjection(icTS, stepGenVertex);
493  maxcTS = hMCVertexcTS1D->GetMaximum();
494  hMCVertexcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
495  SetHistoProperties(hMCVertexcTS1D,49,3,49,20,49,3006);
496  hMCVertexcTS1D->GetXaxis()->SetTitleOffset(1.2);;
497  hMCVertexcTS1D->Sumw2();
498  hMCVertexcTS1D->SetMinimum(0.01);
499  hMCVertexcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Vertex cuts");
500  hMCVertexcTS1D->Draw("hist");
501  cmcVertexcTS->cd();
502  cmcVertexcTS->SetLeftMargin(0.15);
503  cmcVertexcTS->SetRightMargin(0.05);
504  hMCVertexcTS1D->Draw("hist");
505  hMCVertexcTS1D->Draw("err same");
506  cmcVertexcTS->Update();
507 
508  cmc1->cd(20);
509  TH1D *hMCRefitcTS1D = data->ShowProjection(icTS, stepGenRefit);
510  maxcTS = hMCRefitcTS1D->GetMaximum();
511  hMCRefitcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
512  SetHistoProperties(hMCRefitcTS1D,42,3,42,20,42,3006);
513  hMCRefitcTS1D->GetXaxis()->SetTitleOffset(1.2);;
514  hMCRefitcTS1D->Sumw2();
515  hMCRefitcTS1D->SetMinimum(0.01);
516  hMCRefitcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Refit cuts");
517  hMCRefitcTS1D->Draw("hist");
518  cmcRefitcTS->cd();
519  cmcRefitcTS->SetLeftMargin(0.15);
520  cmcRefitcTS->SetRightMargin(0.05);
521  hMCRefitcTS1D->Draw("hist");
522  hMCRefitcTS1D->Draw("err same");
523  cmcRefitcTS->Update();
524 
525  cmc1->cd(21);
526  TH1D *hRECcTS1D = data->ShowProjection(icTS, stepRec);
527  maxcTS = hRECcTS1D->GetMaximum();
528  hRECcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
529  SetHistoProperties(hRECcTS1D,4,3,4,20,4,3004);
530  hRECcTS1D->GetXaxis()->SetTitleOffset(1.2);;
531  hRECcTS1D->Sumw2();
532  hRECcTS1D->SetMinimum(0.01);
533  hRECcTS1D->Draw("hist");
534  ccTS->cd();
535  ccTS->SetLeftMargin(0.15);
536  ccTS->SetRightMargin(0.05);
537  hRECcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD");
538  hRECcTS1D->Draw("hist");
539  hRECcTS1D->Draw("err same");
540  ccTS->Update();
541 
542  cmc1->cd(22);
543  TH1D *hRECAcccTS1D = data->ShowProjection(icTS, stepRecAcc);
544  maxcTS = hRECAcccTS1D->GetMaximum();
545  hRECAcccTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
546  SetHistoProperties(hRECAcccTS1D,38,3,38,20,38,3007);
547  hRECAcccTS1D->GetXaxis()->SetTitleOffset(1.2);;
548  hRECAcccTS1D->Sumw2();
549  hRECAcccTS1D->SetMinimum(0.01);
550  hRECAcccTS1D->Draw("hist");
551  cAcccTS->cd();
552  cAcccTS->SetLeftMargin(0.15);
553  cAcccTS->SetRightMargin(0.05);
554  hRECAcccTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - Acc cuts");
555  hRECAcccTS1D->Draw("hist");
556  hRECAcccTS1D->Draw("err same");
557  cAcccTS->Update();
558 
559  cmc1->cd(23);
560  TH1D *hRECITSClcTS1D = data->ShowProjection(icTS, stepRecITSClus);
561  maxcTS = hRECITSClcTS1D->GetMaximum();
562  hRECITSClcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
563  SetHistoProperties(hRECITSClcTS1D,13,3,13,20,13,3019);
564  hRECITSClcTS1D->GetXaxis()->SetTitleOffset(1.2);;
565  hRECITSClcTS1D->Sumw2();
566  hRECITSClcTS1D->SetMinimum(0.01);
567  hRECITSClcTS1D->Draw("hist");
568  cITSClcTS->cd();
569  cITSClcTS->SetLeftMargin(0.15);
570  cITSClcTS->SetRightMargin(0.05);
571  hRECITSClcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - ITS Cl cuts");
572  hRECITSClcTS1D->Draw("hist");
573  hRECITSClcTS1D->Draw("err same");
574  cITSClcTS->Update();
575 
576  cmc1->cd(24);
577  TH1D *hRECPPRcTS1D = data->ShowProjection(icTS, stepRecPPR);
578  maxcTS = hRECPPRcTS1D->GetMaximum();
579  hRECPPRcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
580  SetHistoProperties(hRECPPRcTS1D,6,3,6,20,6,3019);
581  hRECPPRcTS1D->GetXaxis()->SetTitleOffset(1.2);;
582  hRECPPRcTS1D->Sumw2();
583  hRECPPRcTS1D->SetMinimum(0.01);
584  hRECPPRcTS1D->Draw("hist");
585  cPPRcTS->cd();
586  cPPRcTS->SetLeftMargin(0.15);
587  cPPRcTS->SetRightMargin(0.05);
588  hRECPPRcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - PPR cuts");
589  hRECPPRcTS1D->Draw("hist");
590  hRECPPRcTS1D->Draw("err same");
591  cPPRcTS->Update();
592 
593  // ptPi, ptK, cT
594  cmc2->Divide(8,3);
595 
596  cmc2->cd(1);
597  TH1D *hMCptPi1D = data->ShowProjection(iptPi, stepGen);
598  Double_t maxptPi = hMCptPi1D->GetMaximum();
599  hMCptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
600  hMCptPi1D->Sumw2();
601  SetHistoProperties(hMCptPi1D,2,3,2,20,2,3005);
602  hMCptPi1D->SetMinimum(0.01);
603  hMCptPi1D->GetXaxis()->SetTitleOffset(1.2);
604  hMCptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data");
605  hMCptPi1D->Draw("hist");
606  cmcptPi->cd();
607  cmcptPi->SetLeftMargin(0.15);
608  cmcptPi->SetRightMargin(0.05);
609  hMCptPi1D->Draw("hist");
610  hMCptPi1D->Draw("err same");
611  cmcptPi->Update();
612 
613  cmc2->cd(2);
614  TH1D *hMCAccptPi1D = data->ShowProjection(iptPi, stepGenAcc);
615  maxptPi = hMCAccptPi1D->GetMaximum();
616  hMCAccptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
617  hMCAccptPi1D->Sumw2();
618  hMCAccptPi1D->SetMinimum(0.01);
619  SetHistoProperties(hMCAccptPi1D,8,3,8,20,8,3006);
620  hMCAccptPi1D->GetXaxis()->SetTitleOffset(1.2);
621  hMCAccptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Acc cuts");
622  hMCAccptPi1D->Draw("hist");
623  cmcAccptPi->cd();
624  cmcAccptPi->SetLeftMargin(0.15);
625  cmcAccptPi->SetRightMargin(0.05);
626  hMCAccptPi1D->Draw("hist");
627  hMCAccptPi1D->Draw("err same");
628  cmcAccptPi->Update();
629 
630  cmc2->cd(3);
631  TH1D *hMCVertexptPi1D = data->ShowProjection(iptPi, stepGenVertex);
632  maxptPi = hMCVertexptPi1D->GetMaximum();
633  hMCVertexptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
634  hMCVertexptPi1D->Sumw2();
635  hMCVertexptPi1D->SetMinimum(0.01);
636  SetHistoProperties(hMCVertexptPi1D,49,3,49,20,49,3006);
637  hMCVertexptPi1D->GetXaxis()->SetTitleOffset(1.2);
638  hMCVertexptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Vertex cuts");
639  hMCVertexptPi1D->Draw("hist");
640  cmcVertexptPi->cd();
641  cmcVertexptPi->SetLeftMargin(0.15);
642  cmcVertexptPi->SetRightMargin(0.05);
643  hMCVertexptPi1D->Draw("hist");
644  hMCVertexptPi1D->Draw("err same");
645  cmcVertexptPi->Update();
646 
647  cmc2->cd(4);
648  TH1D *hMCRefitptPi1D = data->ShowProjection(iptPi, stepGenRefit);
649  maxptPi = hMCRefitptPi1D->GetMaximum();
650  hMCRefitptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
651  hMCRefitptPi1D->Sumw2();
652  hMCRefitptPi1D->SetMinimum(0.01);
653  SetHistoProperties(hMCRefitptPi1D,42,3,42,20,42,3006);
654  hMCRefitptPi1D->GetXaxis()->SetTitleOffset(1.2);
655  hMCRefitptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Refit cuts");
656  hMCRefitptPi1D->Draw("hist");
657  cmcRefitptPi->cd();
658  cmcRefitptPi->SetLeftMargin(0.15);
659  cmcRefitptPi->SetRightMargin(0.05);
660  hMCRefitptPi1D->Draw("hist");
661  hMCRefitptPi1D->Draw("err same");
662  cmcRefitptPi->Update();
663 
664  cmc2->cd(5);
665  TH1D *hRECptPi1D = data->ShowProjection(iptPi, stepRec);
666  maxptPi = hRECptPi1D->GetMaximum();
667  hRECptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
668  SetHistoProperties(hRECptPi1D,4,3,4,20,4,3004);
669  hRECptPi1D->GetXaxis()->SetTitleOffset(1.2);;
670  hRECptPi1D->Sumw2();
671  hRECptPi1D->SetMinimum(0.01);
672  hRECptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD");
673  hRECptPi1D->Draw("hist");
674  cptPi->cd();
675  cptPi->SetLeftMargin(0.15);
676  cptPi->SetRightMargin(0.05);
677  hRECptPi1D->Draw("hist");
678  hRECptPi1D->Draw("err same");
679  cptPi->Update();
680 
681  cmc2->cd(6);
682  TH1D *hRECAccptPi1D = data->ShowProjection(iptPi, stepRecAcc);
683  maxptPi = hRECAccptPi1D->GetMaximum();
684  hRECAccptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
685  SetHistoProperties(hRECAccptPi1D,38,3,38,20,38,3007);
686  hRECAccptPi1D->GetXaxis()->SetTitleOffset(1.2);;
687  hRECAccptPi1D->Sumw2();
688  hRECAccptPi1D->SetMinimum(0.01);
689  hRECAccptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - Acc cuts");
690  hRECAccptPi1D->Draw("hist");
691  cAccptPi->cd();
692  cAccptPi->SetLeftMargin(0.15);
693  cAccptPi->SetRightMargin(0.05);
694  hRECAccptPi1D->Draw("hist");
695  hRECAccptPi1D->Draw("err same");
696  cAccptPi->Update();
697 
698  cmc2->cd(7);
699  TH1D *hRECITSClptPi1D = data->ShowProjection(iptPi, stepRecITSClus);
700  maxptPi = hRECITSClptPi1D->GetMaximum();
701  hRECITSClptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
702  SetHistoProperties(hRECITSClptPi1D,13,3,13,20,13,3019);
703  hRECITSClptPi1D->GetXaxis()->SetTitleOffset(1.2);;
704  hRECITSClptPi1D->Sumw2();
705  hRECITSClptPi1D->SetMinimum(0.01);
706  hRECITSClptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - ITS Cl cuts");
707  hRECITSClptPi1D->Draw("hist");
708  cITSClptPi->cd();
709  cITSClptPi->SetLeftMargin(0.15);
710  cITSClptPi->SetRightMargin(0.05);
711  hRECITSClptPi1D->Draw("hist");
712  hRECITSClptPi1D->Draw("err same");
713  cITSClptPi->Update();
714 
715  cmc2->cd(8);
716  TH1D *hRECPPRptPi1D = data->ShowProjection(iptPi, stepRecPPR);
717  maxptPi = hRECPPRptPi1D->GetMaximum();
718  hRECPPRptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
719  SetHistoProperties(hRECPPRptPi1D,6,3,6,20,6,3019);
720  hRECPPRptPi1D->GetXaxis()->SetTitleOffset(1.2);;
721  hRECPPRptPi1D->Sumw2();
722  hRECPPRptPi1D->SetMinimum(0.01);
723  hRECPPRptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - PPR cuts");
724  hRECPPRptPi1D->Draw("hist");
725  cPPRptPi->cd();
726  cPPRptPi->SetLeftMargin(0.15);
727  cPPRptPi->SetRightMargin(0.05);
728  hRECPPRptPi1D->Draw("hist");
729  hRECPPRptPi1D->Draw("err same");
730  cPPRptPi->Update();
731 
732  cmc2->cd(9);
733  TH1D *hMCptK1D = data->ShowProjection(iptK, stepGen);
734  Double_t maxptK = hMCptK1D->GetMaximum();
735  hMCptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
736  SetHistoProperties(hMCptK1D,2,3,2,20,2,3005);
737  hMCptK1D->GetXaxis()->SetTitleOffset(1.2);;
738  hMCptK1D->Sumw2();
739  hMCptK1D->SetMinimum(0.01);
740  hMCptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data");
741  hMCptK1D->Draw("hist");
742  cmcptK->cd();
743  cmcptK->SetLeftMargin(0.15);
744  cmcptK->SetRightMargin(0.05);
745  hMCptK1D->Draw("hist");
746  hMCptK1D->Draw("err same");
747  cmcptK->Update();
748 
749  cmc2->cd(10);
750  TH1D *hMCAccptK1D = data->ShowProjection(iptK, stepGenAcc);
751  maxptK = hMCAccptK1D->GetMaximum();
752  hMCAccptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
753  SetHistoProperties(hMCAccptK1D,8,3,8,20,8,3006);
754  hMCAccptK1D->GetXaxis()->SetTitleOffset(1.2);;
755  hMCAccptK1D->Sumw2();
756  hMCAccptK1D->SetMinimum(0.01);
757  hMCAccptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Acc cuts");
758  hMCAccptK1D->Draw("hist");
759  cmcAccptK->cd();
760  cmcAccptK->SetLeftMargin(0.15);
761  cmcAccptK->SetRightMargin(0.05);
762  hMCAccptK1D->Draw("hist");
763  hMCAccptK1D->Draw("err same");
764  cmcAccptK->Update();
765 
766  cmc2->cd(11);
767  TH1D *hMCVertexptK1D = data->ShowProjection(iptK, stepGenVertex);
768  maxptK = hMCVertexptK1D->GetMaximum();
769  hMCVertexptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
770  SetHistoProperties(hMCVertexptK1D,49,3,49,20,49,3006);
771  hMCVertexptK1D->GetXaxis()->SetTitleOffset(1.2);;
772  hMCVertexptK1D->Sumw2();
773  hMCVertexptK1D->SetMinimum(0.01);
774  hMCVertexptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Vertex cuts");
775  hMCVertexptK1D->Draw("hist");
776  cmcVertexptK->cd();
777  cmcVertexptK->SetLeftMargin(0.15);
778  cmcVertexptK->SetRightMargin(0.05);
779  hMCVertexptK1D->Draw("hist");
780  hMCVertexptK1D->Draw("err same");
781  cmcVertexptK->Update();
782 
783  cmc2->cd(12);
784  TH1D *hMCRefitptK1D = data->ShowProjection(iptK, stepGenRefit);
785  maxptK = hMCRefitptK1D->GetMaximum();
786  hMCRefitptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
787  SetHistoProperties(hMCRefitptK1D,42,3,42,20,42,3006);
788  hMCRefitptK1D->GetXaxis()->SetTitleOffset(1.2);;
789  hMCRefitptK1D->Sumw2();
790  hMCRefitptK1D->SetMinimum(0.01);
791  hMCRefitptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Refit cuts");
792  hMCRefitptK1D->Draw("hist");
793  cmcRefitptK->cd();
794  cmcRefitptK->SetLeftMargin(0.15);
795  cmcRefitptK->SetRightMargin(0.05);
796  hMCRefitptK1D->Draw("hist");
797  hMCRefitptK1D->Draw("err same");
798  cmcRefitptK->Update();
799 
800  cmc2->cd(13);
801  TH1D *hRECptK1D = data->ShowProjection(iptK, stepRec);
802  maxptK = hRECptK1D->GetMaximum();
803  hRECptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
804  SetHistoProperties(hRECptK1D,4,3,4,20,4,3004);
805  hRECptK1D->GetXaxis()->SetTitleOffset(1.2);;
806  hRECptK1D->Sumw2();
807  hRECptK1D->SetMinimum(0.01);
808  hRECptK1D->Draw("hist");
809  cptK->cd();
810  cptK->SetLeftMargin(0.15);
811  cptK->SetRightMargin(0.05);
812  hRECptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD");
813  hRECptK1D->Draw("hist");
814  hRECptK1D->Draw("err same");
815  cptK->Update();
816 
817  cmc2->cd(14);
818  TH1D *hRECAccptK1D = data->ShowProjection(iptK, stepRecAcc);
819  maxptK = hRECAccptK1D->GetMaximum();
820  hRECAccptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
821  SetHistoProperties(hRECAccptK1D,38,3,38,20,38,3007);
822  hRECAccptK1D->GetXaxis()->SetTitleOffset(1.2);;
823  hRECAccptK1D->Sumw2();
824  hRECAccptK1D->SetMinimum(0.01);
825  hRECAccptK1D->Draw("hist");
826  cAccptK->cd();
827  cAccptK->SetLeftMargin(0.15);
828  cAccptK->SetRightMargin(0.05);
829  hRECAccptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - Acc cuts");
830  hRECAccptK1D->Draw("hist");
831  hRECAccptK1D->Draw("err same");
832  cAccptK->Update();
833 
834  cmc2->cd(15);
835  TH1D *hRECITSClptK1D = data->ShowProjection(iptK, stepRecITSClus);
836  maxptK = hRECITSClptK1D->GetMaximum();
837  hRECITSClptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
838  SetHistoProperties(hRECITSClptK1D,13,3,13,20,13,3019);
839  hRECITSClptK1D->GetXaxis()->SetTitleOffset(1.2);;
840  hRECITSClptK1D->Sumw2();
841  hRECITSClptK1D->SetMinimum(0.01);
842  hRECITSClptK1D->Draw("hist");
843  cITSClptK->cd();
844  cITSClptK->SetLeftMargin(0.15);
845  cITSClptK->SetRightMargin(0.05);
846  hRECITSClptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - ITS Cl cuts");
847  hRECITSClptK1D->Draw("hist");
848  hRECITSClptK1D->Draw("err same");
849  cITSClptK->Update();
850 
851  cmc2->cd(16);
852  TH1D *hRECPPRptK1D = data->ShowProjection(iptK, stepRecPPR);
853  maxptK = hRECPPRptK1D->GetMaximum();
854  hRECPPRptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
855  SetHistoProperties(hRECPPRptK1D,6,3,6,20,6,3019);
856  hRECPPRptK1D->GetXaxis()->SetTitleOffset(1.2);;
857  hRECPPRptK1D->Sumw2();
858  hRECPPRptK1D->SetMinimum(0.01);
859  hRECPPRptK1D->Draw("hist");
860  cPPRptK->cd();
861  cPPRptK->SetLeftMargin(0.15);
862  cPPRptK->SetRightMargin(0.05);
863  hRECPPRptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - PPR cuts");
864  hRECPPRptK1D->Draw("hist");
865  hRECPPRptK1D->Draw("err same");
866  cPPRptK->Update();
867 
868  cmc2->cd(17);
869  TH1D *hMCcT1D = data->ShowProjection(icT, stepGen);
870  Double_t maxcT = hMCcT1D->GetMaximum();
871  hMCcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
872  SetHistoProperties(hMCcT1D,2,3,2,20,2,3005);
873  hMCcT1D->GetXaxis()->SetTitleOffset(1.2);;
874  hMCcT1D->Sumw2();
875  hMCcT1D->SetMinimum(0.01);
876  hMCcT1D->GetXaxis()->SetTitle("ct (#mum), MC data");
877  hMCcT1D->Draw("hist");
878  cmccT->cd();
879  cmccT->SetLeftMargin(0.15);
880  cmccT->SetRightMargin(0.05);
881  hMCcT1D->Draw("hist");
882  hMCcT1D->Draw("err same");
883  cmccT->Update();
884 
885  cmc2->cd(18);
886  TH1D *hMCAcccT1D = data->ShowProjection(icT, stepGenAcc);
887  maxcT = hMCAcccT1D->GetMaximum();
888  hMCAcccT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
889  SetHistoProperties(hMCAcccT1D,8,3,8,20,8,3006);
890  hMCAcccT1D->GetXaxis()->SetTitleOffset(1.2);
891  hMCAcccT1D->Sumw2();
892  hMCAcccT1D->SetMinimum(0.01);
893  hMCAcccT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Acc cuts");
894  hMCAcccT1D->Draw("hist");
895  cmcAcccT->cd();
896  cmcAcccT->SetLeftMargin(0.15);
897  cmcAcccT->SetRightMargin(0.05);
898  hMCAcccT1D->Draw("hist");
899  hMCAcccT1D->Draw("err same");
900  cmcAcccT->Update();
901 
902  cmc2->cd(19);
903  TH1D *hMCVertexcT1D = data->ShowProjection(icT, stepGenVertex);
904  maxcT = hMCVertexcT1D->GetMaximum();
905  hMCVertexcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
906  SetHistoProperties(hMCVertexcT1D,49,3,49,20,49,3006);
907  hMCVertexcT1D->GetXaxis()->SetTitleOffset(1.2);
908  hMCVertexcT1D->Sumw2();
909  hMCVertexcT1D->SetMinimum(0.01);
910  hMCVertexcT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Vertex cuts");
911  hMCVertexcT1D->Draw("hist");
912  cmcVertexcT->cd();
913  cmcVertexcT->SetLeftMargin(0.15);
914  cmcVertexcT->SetRightMargin(0.05);
915  hMCVertexcT1D->Draw("hist");
916  hMCVertexcT1D->Draw("err same");
917  cmcVertexcT->Update();
918 
919  cmc2->cd(20);
920  TH1D *hMCRefitcT1D = data->ShowProjection(icT, stepGenRefit);
921  maxcT = hMCRefitcT1D->GetMaximum();
922  hMCRefitcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
923  SetHistoProperties(hMCRefitcT1D,42,3,42,20,42,3006);
924  hMCRefitcT1D->GetXaxis()->SetTitleOffset(1.2);
925  hMCRefitcT1D->Sumw2();
926  hMCRefitcT1D->SetMinimum(0.01);
927  hMCRefitcT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Refit cuts");
928  hMCRefitcT1D->Draw("hist");
929  cmcRefitcT->cd();
930  cmcRefitcT->SetLeftMargin(0.15);
931  cmcRefitcT->SetRightMargin(0.05);
932  hMCRefitcT1D->Draw("hist");
933  hMCRefitcT1D->Draw("err same");
934  cmcRefitcT->Update();
935 
936  cmc2->cd(21);
937  TH1D *hRECcT1D = data->ShowProjection(icT, stepRec);
938  maxcT = hRECcT1D->GetMaximum();
939  hRECcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
940  SetHistoProperties(hRECcT1D,4,3,4,20,4,3004);
941  hRECcT1D->GetXaxis()->SetTitleOffset(1.2);;
942  hRECcT1D->Sumw2();
943  hRECcT1D->SetMinimum(0.01);
944  hRECcT1D->Draw("hist");
945  ccT->cd();
946  ccT->SetLeftMargin(0.15);
947  ccT->SetRightMargin(0.05);
948  hRECcT1D->GetXaxis()->SetTitle("ct (#mum), AOD");
949  hRECcT1D->Draw("hist");
950  hRECcT1D->Draw("err same");
951  ccT->Update();
952 
953  cmc2->cd(22);
954  TH1D *hRECAcccT1D = data->ShowProjection(icT, stepRecAcc);
955  maxcT = hRECAcccT1D->GetMaximum();
956  hRECAcccT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
957  SetHistoProperties(hRECAcccT1D,38,3,38,20,38,3007);
958  hRECAcccT1D->GetXaxis()->SetTitleOffset(1.2);;
959  hRECAcccT1D->Sumw2();
960  hRECAcccT1D->SetMinimum(0.01);
961  hRECAcccT1D->Draw("hist");
962  cAcccT->cd();
963  cAcccT->SetLeftMargin(0.15);
964  cAcccT->SetRightMargin(0.05);
965  hRECAcccT1D->GetXaxis()->SetTitle("ct (#mum), AOD - Acc cuts");
966  hRECAcccT1D->Draw("hist");
967  hRECAcccT1D->Draw("err same");
968  cAcccT->Update();
969 
970  cmc2->cd(23);
971  TH1D *hRECITSClcT1D = data->ShowProjection(icT, stepRecITSClus);
972  maxcT = hRECITSClcT1D->GetMaximum();
973  hRECITSClcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
974  SetHistoProperties(hRECITSClcT1D,13,3,13,20,13,3019);
975  hRECITSClcT1D->GetXaxis()->SetTitleOffset(1.2);
976  hRECITSClcT1D->Sumw2();
977  hRECITSClcT1D->SetMinimum(0.01);
978  hRECITSClcT1D->Draw("hist");
979  cITSClcT->cd();
980  cITSClcT->SetLeftMargin(0.15);
981  cITSClcT->SetRightMargin(0.05);
982  hRECITSClcT1D->GetXaxis()->SetTitle("ct (#mum), AOD - ITS Cl cuts");
983  hRECITSClcT1D->Draw("hist");
984  hRECITSClcT1D->Draw("err same");
985  cITSClcT->Update();
986 
987  cmc2->cd(24);
988  TH1D *hRECPPRcT1D = data->ShowProjection(icT, stepRecPPR);
989  maxcT = hRECPPRcT1D->GetMaximum();
990  hRECPPRcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
991  SetHistoProperties(hRECPPRcT1D,6,3,6,20,6,3019);
992  hRECPPRcT1D->GetXaxis()->SetTitleOffset(1.2);
993  hRECPPRcT1D->Sumw2();
994  hRECPPRcT1D->SetMinimum(0.01);
995  hRECPPRcT1D->Draw("hist");
996  cPPRcT->cd();
997  cPPRcT->SetLeftMargin(0.15);
998  cPPRcT->SetRightMargin(0.05);
999  hRECPPRcT1D->GetXaxis()->SetTitle("ct (#mum), AOD - PPR cuts");
1000  hRECPPRcT1D->Draw("hist");
1001  hRECPPRcT1D->Draw("err same");
1002  cPPRcT->Update();
1003 
1004  //dca, d0pi, d0K
1005  cmc3->Divide(8,3);
1006  cmc3->cd(1);
1007  TH1D *hMCdca1D = data->ShowProjection(idca, stepGen);
1008  Double_t maxdca = hMCdca1D->GetMaximum();
1009  hMCdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1010  hMCdca1D->Sumw2();
1011  hMCdca1D->SetMinimum(0.01);
1012  SetHistoProperties(hMCdca1D,2,3,2,20,2,3005);
1013  hMCdca1D->GetXaxis()->SetTitleOffset(1.2);
1014  hMCdca1D->GetXaxis()->SetTitle("dca (#mum), MC data");
1015  hMCdca1D->Draw("hist");
1016  cmcdca->cd();
1017  cmcdca->SetLeftMargin(0.15);
1018  cmcdca->SetRightMargin(0.05);
1019  hMCdca1D->Draw("hist");
1020  hMCdca1D->Draw("err same");
1021  cmcdca->Update();
1022 
1023  cmc3->cd(2);
1024  TH1D *hMCAccdca1D = data->ShowProjection(idca, stepGenAcc);
1025  maxdca = hMCAccdca1D->GetMaximum();
1026  hMCAccdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1027  SetHistoProperties(hMCAccdca1D,8,3,8,20,8,3006);
1028  hMCAccdca1D->GetXaxis()->SetTitleOffset(1.2);;
1029  hMCAccdca1D->Sumw2();
1030  hMCAccdca1D->SetMinimum(0.01);
1031  hMCAccdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Acc cut");
1032  hMCAccdca1D->Draw("hist");
1033  cmcAccdca->cd();
1034  cmcAccdca->SetLeftMargin(0.15);
1035  cmcAccdca->SetRightMargin(0.05);
1036  hMCAccdca1D->Draw("hist");
1037  hMCAccdca1D->Draw("err same");
1038  cmcAccdca->Update();
1039 
1040  cmc3->cd(3);
1041  TH1D *hMCVertexdca1D = data->ShowProjection(idca, stepGenVertex);
1042  maxdca = hMCVertexdca1D->GetMaximum();
1043  hMCVertexdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1044  SetHistoProperties(hMCVertexdca1D,49,3,49,20,49,3006);
1045  hMCVertexdca1D->GetXaxis()->SetTitleOffset(1.2);;
1046  hMCVertexdca1D->Sumw2();
1047  hMCVertexdca1D->SetMinimum(0.01);
1048  hMCVertexdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Vertex cut");
1049  hMCVertexdca1D->Draw("hist");
1050  cmcVertexdca->cd();
1051  cmcVertexdca->SetLeftMargin(0.15);
1052  cmcVertexdca->SetRightMargin(0.05);
1053  hMCVertexdca1D->Draw("hist");
1054  hMCVertexdca1D->Draw("err same");
1055  cmcVertexdca->Update();
1056 
1057  cmc3->cd(4);
1058  TH1D *hMCRefitdca1D = data->ShowProjection(idca, stepGenRefit);
1059  maxdca = hMCRefitdca1D->GetMaximum();
1060  hMCRefitdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1061  SetHistoProperties(hMCRefitdca1D,42,3,42,20,42,3006);
1062  hMCRefitdca1D->GetXaxis()->SetTitleOffset(1.2);;
1063  hMCRefitdca1D->Sumw2();
1064  hMCRefitdca1D->SetMinimum(0.01);
1065  hMCRefitdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Refit cut");
1066  hMCRefitdca1D->Draw("hist");
1067  cmcRefitdca->cd();
1068  cmcRefitdca->SetLeftMargin(0.15);
1069  cmcRefitdca->SetRightMargin(0.05);
1070  hMCRefitdca1D->Draw("hist");
1071  hMCRefitdca1D->Draw("err same");
1072  cmcRefitdca->Update();
1073 
1074  cmc3->cd(5);
1075  TH1D *hRECdca1D = data->ShowProjection(idca, stepRec);
1076  maxdca = hRECdca1D->GetMaximum();
1077  hRECdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1078  SetHistoProperties(hRECdca1D,4,3,4,20,4,3004);
1079  hRECdca1D->GetXaxis()->SetTitleOffset(1.2);;
1080  hRECdca1D->Sumw2();
1081  hRECdca1D->SetMinimum(0.01);
1082  hRECdca1D->GetXaxis()->SetTitle("dca (#mum), AOD");
1083  hRECdca1D->Draw("hist");
1084  cdca->cd();
1085  cdca->SetLeftMargin(0.15);
1086  cdca->SetRightMargin(0.05);
1087  hRECdca1D->Draw("hist");
1088  hRECdca1D->Draw("err same");
1089  cdca->Update();
1090 
1091  cmc3->cd(6);
1092  TH1D *hRECAccdca1D = data->ShowProjection(idca, stepRecAcc);
1093  maxdca = hRECAccdca1D->GetMaximum();
1094  hRECAccdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1095  SetHistoProperties(hRECAccdca1D,38,3,38,20,38,3007);
1096  hRECAccdca1D->GetXaxis()->SetTitleOffset(1.2);;
1097  hRECAccdca1D->Sumw2();
1098  hRECAccdca1D->SetMinimum(0.01);
1099  hRECAccdca1D->GetXaxis()->SetTitle("dca (#mum), AOD - Acc cuts");
1100  hRECAccdca1D->Draw("hist");
1101  cAccdca->cd();
1102  cAccdca->SetLeftMargin(0.15);
1103  cAccdca->SetRightMargin(0.05);
1104  hRECAccdca1D->Draw("hist");
1105  hRECAccdca1D->Draw("err same");
1106  cAccdca->Update();
1107 
1108  cmc3->cd(6);
1109  TH1D *hRECITSCldca1D = data->ShowProjection(idca, stepRecITSClus);
1110  maxdca = hRECITSCldca1D->GetMaximum();
1111  hRECITSCldca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1112  SetHistoProperties(hRECITSCldca1D,13,3,13,20,13,3019);
1113  hRECITSCldca1D->GetXaxis()->SetTitleOffset(1.2);;
1114  hRECITSCldca1D->Sumw2();
1115  hRECITSCldca1D->SetMinimum(0.01);
1116  hRECITSCldca1D->GetXaxis()->SetTitle("dca (#mum), AOD - ITS Cl cuts");
1117  hRECITSCldca1D->Draw("hist");
1118  cITSCldca->cd();
1119  cITSCldca->SetLeftMargin(0.15);
1120  cITSCldca->SetRightMargin(0.05);
1121  hRECITSCldca1D->Draw("hist");
1122  hRECITSCldca1D->Draw("err same");
1123  cITSCldca->Update();
1124 
1125  cmc3->cd(7);
1126  TH1D *hRECPPRdca1D = data->ShowProjection(idca, stepRecPPR);
1127  maxdca = hRECPPRdca1D->GetMaximum();
1128  hRECPPRdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
1129  SetHistoProperties(hRECPPRdca1D,6,3,6,20,6,3019);
1130  hRECPPRdca1D->GetXaxis()->SetTitleOffset(1.2);;
1131  hRECPPRdca1D->Sumw2();
1132  hRECPPRdca1D->SetMinimum(0.01);
1133  hRECPPRdca1D->GetXaxis()->SetTitle("dca (#mum), AOD - PPR cuts");
1134  hRECPPRdca1D->Draw("hist");
1135  cPPRdca->cd();
1136  cPPRdca->SetLeftMargin(0.15);
1137  cPPRdca->SetRightMargin(0.05);
1138  hRECPPRdca1D->Draw("hist");
1139  hRECPPRdca1D->Draw("err same");
1140  cPPRdca->Update();
1141 
1142  cmc3->cd(8);
1143  TH1D *hMCd0pi1D = data->ShowProjection(id0pi, stepGen);
1144  Double_t maxd0pi = hMCd0pi1D->GetMaximum();
1145  hMCd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1146  hMCd0pi1D->Sumw2();
1147  hMCd0pi1D->SetMinimum(0.01);
1148  SetHistoProperties(hMCd0pi1D,2,3,2,20,2,3005);
1149  hMCd0pi1D->GetXaxis()->SetTitleOffset(1.2);
1150  hMCd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC data");
1151  hMCd0pi1D->Draw("hist");
1152  cmcd0pi->cd();
1153  cmcd0pi->SetLeftMargin(0.15);
1154  cmcd0pi->SetRightMargin(0.05);
1155  hMCd0pi1D->Draw("hist");
1156  hMCd0pi1D->Draw("err same");
1157  cmcd0pi->Update();
1158 
1159  cmc3->cd(9);
1160  TH1D *hMCAccd0pi1D = data->ShowProjection(id0pi, stepGenAcc);
1161  maxd0pi = hMCAccd0pi1D->GetMaximum();
1162  hMCAccd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1163  SetHistoProperties(hMCAccd0pi1D,8,3,8,20,8,3006);
1164  hMCAccd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
1165  hMCAccd0pi1D->Sumw2();
1166  hMCAccd0pi1D->SetMinimum(0.01);
1167  hMCAccd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Acc cut");
1168  hMCAccd0pi1D->Draw("hist");
1169  cmcAccd0pi->cd();
1170  cmcAccd0pi->SetLeftMargin(0.15);
1171  cmcAccd0pi->SetRightMargin(0.05);
1172  hMCAccd0pi1D->Draw("hist");
1173  hMCAccd0pi1D->Draw("err same");
1174  cmcAccd0pi->Update();
1175 
1176  cmc3->cd(10);
1177  TH1D *hMCVertexd0pi1D = data->ShowProjection(id0pi, stepGenVertex);
1178  maxd0pi = hMCVertexd0pi1D->GetMaximum();
1179  hMCVertexd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1180  SetHistoProperties(hMCVertexd0pi1D,49,3,49,20,49,3006);
1181  hMCVertexd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
1182  hMCVertexd0pi1D->Sumw2();
1183  hMCVertexd0pi1D->SetMinimum(0.01);
1184  hMCVertexd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Vertex cut");
1185  hMCVertexd0pi1D->Draw("hist");
1186  cmcVertexd0pi->cd();
1187  cmcVertexd0pi->SetLeftMargin(0.15);
1188  cmcVertexd0pi->SetRightMargin(0.05);
1189  hMCVertexd0pi1D->Draw("hist");
1190  hMCVertexd0pi1D->Draw("err same");
1191  cmcVertexd0pi->Update();
1192 
1193  cmc3->cd(11);
1194  TH1D *hMCRefitd0pi1D = data->ShowProjection(id0pi, stepGenRefit);
1195  maxd0pi = hMCRefitd0pi1D->GetMaximum();
1196  hMCRefitd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1197  SetHistoProperties(hMCRefitd0pi1D,42,3,42,20,42,3006);
1198  hMCRefitd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
1199  hMCRefitd0pi1D->Sumw2();
1200  hMCRefitd0pi1D->SetMinimum(0.01);
1201  hMCRefitd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Refit cut");
1202  hMCRefitd0pi1D->Draw("hist");
1203  cmcRefitd0pi->cd();
1204  cmcRefitd0pi->SetLeftMargin(0.15);
1205  cmcRefitd0pi->SetRightMargin(0.05);
1206  hMCRefitd0pi1D->Draw("hist");
1207  hMCRefitd0pi1D->Draw("err same");
1208  cmcRefitd0pi->Update();
1209 
1210  cmc3->cd(12);
1211  TH1D *hRECd0pi1D = data->ShowProjection(id0pi, stepRec);
1212  maxd0pi = hRECd0pi1D->GetMaximum();
1213  hRECd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1214  SetHistoProperties(hRECd0pi1D,4,3,4,20,4,3004);
1215  hRECd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
1216  hRECd0pi1D->Sumw2();
1217  hRECd0pi1D->SetMinimum(0.01);
1218  hRECd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD");
1219  hRECd0pi1D->Draw("hist");
1220  cd0pi->cd();
1221  cd0pi->SetLeftMargin(0.15);
1222  cd0pi->SetRightMargin(0.05);
1223  hRECd0pi1D->Draw("hist");
1224  hRECd0pi1D->Draw("err same");
1225  cd0pi->Update();
1226 
1227  cmc3->cd(13);
1228  TH1D *hRECAccd0pi1D = data->ShowProjection(id0pi, stepRecAcc);
1229  maxd0pi = hRECAccd0pi1D->GetMaximum();
1230  hRECAccd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1231  SetHistoProperties(hRECAccd0pi1D,38,3,38,20,38,3007);
1232  hRECAccd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
1233  hRECAccd0pi1D->Sumw2();
1234  hRECAccd0pi1D->SetMinimum(0.01);
1235  hRECAccd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - Acc cuts");
1236  hRECAccd0pi1D->Draw("hist");
1237  cAccd0pi->cd();
1238  cAccd0pi->SetLeftMargin(0.15);
1239  cAccd0pi->SetRightMargin(0.05);
1240  hRECAccd0pi1D->Draw("hist");
1241  hRECAccd0pi1D->Draw("err same");
1242  cAccd0pi->Update();
1243 
1244  cmc3->cd(14);
1245  TH1D *hRECITSCld0pi1D = data->ShowProjection(id0pi, stepRecITSClus);
1246  maxd0pi = hRECITSCld0pi1D->GetMaximum();
1247  hRECITSCld0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1248  SetHistoProperties(hRECITSCld0pi1D,13,3,13,20,13,3019);
1249  hRECITSCld0pi1D->GetXaxis()->SetTitleOffset(1.2);;
1250  hRECITSCld0pi1D->Sumw2();
1251  hRECITSCld0pi1D->SetMinimum(0.01);
1252  hRECITSCld0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - ITS Cl cuts");
1253  hRECITSCld0pi1D->Draw("hist");
1254  cITSCld0pi->cd();
1255  cITSCld0pi->SetLeftMargin(0.15);
1256  cITSCld0pi->SetRightMargin(0.05);
1257  hRECITSCld0pi1D->Draw("hist");
1258  hRECITSCld0pi1D->Draw("err same");
1259  cITSCld0pi->Update();
1260 
1261  cmc3->cd(15);
1262  TH1D *hRECPPRd0pi1D = data->ShowProjection(id0pi, stepRecPPR);
1263  maxd0pi = hRECPPRd0pi1D->GetMaximum();
1264  hRECPPRd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
1265  SetHistoProperties(hRECPPRd0pi1D,6,3,6,20,6,3019);
1266  hRECPPRd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
1267  hRECPPRd0pi1D->Sumw2();
1268  hRECPPRd0pi1D->SetMinimum(0.01);
1269  hRECPPRd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - PPR cuts");
1270  hRECPPRd0pi1D->Draw("hist");
1271  cPPRd0pi->cd();
1272  cPPRd0pi->SetLeftMargin(0.15);
1273  cPPRd0pi->SetRightMargin(0.05);
1274  hRECPPRd0pi1D->Draw("hist");
1275  hRECPPRd0pi1D->Draw("err same");
1276  cPPRd0pi->Update();
1277 
1278  cmc3->cd(16);
1279  TH1D *hMCd0K1D = data->ShowProjection(id0K, stepGen);
1280  Double_t maxd0K = hMCd0K1D->GetMaximum();
1281  hMCd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1282  hMCd0K1D->Sumw2();
1283  hMCd0K1D->SetMinimum(0.01);
1284  SetHistoProperties(hMCd0K1D,2,3,2,20,2,3005);
1285  hMCd0K1D->GetXaxis()->SetTitleOffset(1.2);
1286  hMCd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC data");
1287  hMCd0K1D->Draw("hist");
1288  cmcd0K->cd();
1289  cmcd0K->SetLeftMargin(0.15);
1290  cmcd0K->SetRightMargin(0.05);
1291  hMCd0K1D->Draw("hist");
1292  hMCd0K1D->Draw("err same");
1293  cmcd0K->Update();
1294 
1295  cmc3->cd(17);
1296  TH1D *hMCAccd0K1D = data->ShowProjection(id0K, stepGenAcc);
1297  maxd0K = hMCAccd0K1D->GetMaximum();
1298  hMCAccd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1299  SetHistoProperties(hMCAccd0K1D,8,3,8,20,8,3006);
1300  hMCAccd0K1D->GetXaxis()->SetTitleOffset(1.2);;
1301  hMCAccd0K1D->Sumw2();
1302  hMCAccd0K1D->SetMinimum(0.01);
1303  hMCAccd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Acc cut");
1304  hMCAccd0K1D->Draw("hist");
1305  cmcAccd0K->cd();
1306  cmcAccd0K->SetLeftMargin(0.15);
1307  cmcAccd0K->SetRightMargin(0.05);
1308  hMCAccd0K1D->Draw("hist");
1309  hMCAccd0K1D->Draw("err same");
1310  cmcAccd0K->Update();
1311 
1312  cmc3->cd(18);
1313  TH1D *hMCVertexd0K1D = data->ShowProjection(id0K, stepGenVertex);
1314  maxd0K = hMCVertexd0K1D->GetMaximum();
1315  hMCVertexd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1316  SetHistoProperties(hMCVertexd0K1D,49,3,49,20,49,3006);
1317  hMCVertexd0K1D->GetXaxis()->SetTitleOffset(1.2);;
1318  hMCVertexd0K1D->Sumw2();
1319  hMCVertexd0K1D->SetMinimum(0.01);
1320  hMCVertexd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Vertex cut");
1321  hMCVertexd0K1D->Draw("hist");
1322  cmcVertexd0K->cd();
1323  cmcVertexd0K->SetLeftMargin(0.15);
1324  cmcVertexd0K->SetRightMargin(0.05);
1325  hMCVertexd0K1D->Draw("hist");
1326  hMCVertexd0K1D->Draw("err same");
1327  cmcVertexd0K->Update();
1328 
1329  cmc3->cd(20);
1330  TH1D *hMCRefitd0K1D = data->ShowProjection(id0K, stepGenRefit);
1331  maxd0K = hMCRefitd0K1D->GetMaximum();
1332  hMCRefitd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1333  SetHistoProperties(hMCRefitd0K1D,42,3,42,20,42,3006);
1334  hMCRefitd0K1D->GetXaxis()->SetTitleOffset(1.2);;
1335  hMCRefitd0K1D->Sumw2();
1336  hMCRefitd0K1D->SetMinimum(0.01);
1337  hMCRefitd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Refit cut");
1338  hMCRefitd0K1D->Draw("hist");
1339  cmcRefitd0K->cd();
1340  cmcRefitd0K->SetLeftMargin(0.15);
1341  cmcRefitd0K->SetRightMargin(0.05);
1342  hMCRefitd0K1D->Draw("hist");
1343  hMCRefitd0K1D->Draw("err same");
1344  cmcRefitd0K->Update();
1345 
1346  cmc3->cd(21);
1347  TH1D *hRECd0K1D = data->ShowProjection(id0K, stepRec);
1348  maxd0K = hRECd0K1D->GetMaximum();
1349  hRECd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1350  SetHistoProperties(hRECd0K1D,4,3,4,20,4,3004);
1351  hRECd0K1D->GetXaxis()->SetTitleOffset(1.2);;
1352  hRECd0K1D->Sumw2();
1353  hRECd0K1D->SetMinimum(0.01);
1354  hRECd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD");
1355  hRECd0K1D->Draw("hist");
1356  cd0K->cd();
1357  cd0K->SetLeftMargin(0.15);
1358  cd0K->SetRightMargin(0.05);
1359  hRECd0K1D->Draw("hist");
1360  hRECd0K1D->Draw("err same");
1361  cd0K->Update();
1362 
1363  cmc3->cd(22);
1364  TH1D *hRECAccd0K1D = data->ShowProjection(id0K, stepRecAcc);
1365  maxd0K = hRECAccd0K1D->GetMaximum();
1366  hRECAccd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1367  SetHistoProperties(hRECAccd0K1D,38,3,38,20,38,3007);
1368  hRECAccd0K1D->GetXaxis()->SetTitleOffset(1.2);;
1369  hRECAccd0K1D->Sumw2();
1370  hRECAccd0K1D->SetMinimum(0.01);
1371  hRECAccd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - Acc cuts");
1372  hRECAccd0K1D->Draw("hist");
1373  cAccd0K->cd();
1374  cAccd0K->SetLeftMargin(0.15);
1375  cAccd0K->SetRightMargin(0.05);
1376  hRECAccd0K1D->Draw("hist");
1377  hRECAccd0K1D->Draw("err same");
1378  cAccd0K->Update();
1379 
1380  cmc3->cd(23);
1381  TH1D *hRECITSCld0K1D = data->ShowProjection(id0K, stepRecITSClus);
1382  maxd0K = hRECITSCld0K1D->GetMaximum();
1383  hRECITSCld0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1384  SetHistoProperties(hRECITSCld0K1D,13,3,13,20,13,3019);
1385  hRECITSCld0K1D->GetXaxis()->SetTitleOffset(1.2);;
1386  hRECITSCld0K1D->Sumw2();
1387  hRECITSCld0K1D->SetMinimum(0.01);
1388  hRECITSCld0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - ITS Cl cuts");
1389  hRECITSCld0K1D->Draw("hist");
1390  cITSCld0K->cd();
1391  cITSCld0K->SetLeftMargin(0.15);
1392  cITSCld0K->SetRightMargin(0.05);
1393  hRECITSCld0K1D->Draw("hist");
1394  hRECITSCld0K1D->Draw("err same");
1395  cITSCld0K->Update();
1396 
1397  cmc3->cd(24);
1398  TH1D *hRECPPRd0K1D = data->ShowProjection(id0K, stepRecPPR);
1399  maxd0K = hRECPPRd0K1D->GetMaximum();
1400  hRECPPRd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
1401  SetHistoProperties(hRECPPRd0K1D,6,3,6,20,6,3019);
1402  hRECPPRd0K1D->GetXaxis()->SetTitleOffset(1.2);;
1403  hRECPPRd0K1D->Sumw2();
1404  hRECPPRd0K1D->SetMinimum(0.01);
1405  hRECPPRd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - PPR cuts");
1406  hRECPPRd0K1D->Draw("hist");
1407  cPPRd0K->cd();
1408  cPPRd0K->SetLeftMargin(0.15);
1409  cPPRd0K->SetRightMargin(0.05);
1410  hRECPPRd0K1D->Draw("hist");
1411  hRECPPRd0K1D->Draw("err same");
1412  cPPRd0K->Update();
1413 
1414  //d0xd0, cosPointingAngle, Phi
1415  cmc4->Divide(8,3);
1416  cmc4->cd(1);
1417  TH1D *hMCd0xd01D = data->ShowProjection(id0xd0, stepGen);
1418  Double_t maxd0xd0 = hMCd0xd01D->GetMaximum();
1419  hMCd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1420  hMCd0xd01D->Sumw2();
1421  hMCd0xd01D->SetMinimum(0.01);
1422  SetHistoProperties(hMCd0xd01D,2,3,2,20,2,3005);
1423  hMCd0xd01D->GetXaxis()->SetTitleOffset(1.2);
1424  hMCd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC data");
1425  hMCd0xd01D->Draw("hist");
1426  cmcd0xd0->cd();
1427  cmcd0xd0->SetLeftMargin(0.15);
1428  cmcd0xd0->SetRightMargin(0.05);
1429  hMCd0xd01D->Draw("hist");
1430  hMCd0xd01D->Draw("err same");
1431  cmcd0xd0->Update();
1432 
1433  cmc4->cd(2);
1434  TH1D *hMCAccd0xd01D = data->ShowProjection(id0xd0, stepGenAcc);
1435  maxd0xd0 = hMCAccd0xd01D->GetMaximum();
1436  hMCAccd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1437  SetHistoProperties(hMCAccd0xd01D,8,3,8,20,8,3006);
1438  hMCAccd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
1439  hMCAccd0xd01D->Sumw2();
1440  hMCAccd0xd01D->SetMinimum(0.01);
1441  hMCAccd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Acc cut");
1442  hMCAccd0xd01D->Draw("hist");
1443  cmcAccd0xd0->cd();
1444  cmcAccd0xd0->SetLeftMargin(0.15);
1445  cmcAccd0xd0->SetRightMargin(0.05);
1446  hMCAccd0xd01D->Draw("hist");
1447  hMCAccd0xd01D->Draw("err same");
1448  cmcAccd0xd0->Update();
1449 
1450  cmc4->cd(3);
1451  TH1D *hMCVertexd0xd01D = data->ShowProjection(id0xd0, stepGenVertex);
1452  maxd0xd0 = hMCVertexd0xd01D->GetMaximum();
1453  hMCVertexd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1454  SetHistoProperties(hMCVertexd0xd01D,49,3,49,20,49,3006);
1455  hMCVertexd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
1456  hMCVertexd0xd01D->Sumw2();
1457  hMCVertexd0xd01D->SetMinimum(0.01);
1458  hMCVertexd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Vertex cut");
1459  hMCVertexd0xd01D->Draw("hist");
1460  cmcVertexd0xd0->cd();
1461  cmcVertexd0xd0->SetLeftMargin(0.15);
1462  cmcVertexd0xd0->SetRightMargin(0.05);
1463  hMCVertexd0xd01D->Draw("hist");
1464  hMCVertexd0xd01D->Draw("err same");
1465  cmcVertexd0xd0->Update();
1466 
1467  cmc4->cd(4);
1468  TH1D *hMCRefitd0xd01D = data->ShowProjection(id0xd0, stepGenRefit);
1469  maxd0xd0 = hMCRefitd0xd01D->GetMaximum();
1470  hMCRefitd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1471  SetHistoProperties(hMCRefitd0xd01D,42,3,42,20,42,3006);
1472  hMCRefitd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
1473  hMCRefitd0xd01D->Sumw2();
1474  hMCRefitd0xd01D->SetMinimum(0.01);
1475  hMCRefitd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Refit cut");
1476  hMCRefitd0xd01D->Draw("hist");
1477  cmcRefitd0xd0->cd();
1478  cmcRefitd0xd0->SetLeftMargin(0.15);
1479  cmcRefitd0xd0->SetRightMargin(0.05);
1480  hMCRefitd0xd01D->Draw("hist");
1481  hMCRefitd0xd01D->Draw("err same");
1482  cmcRefitd0xd0->Update();
1483 
1484  cmc4->cd(5);
1485  TH1D *hRECd0xd01D = data->ShowProjection(id0xd0, stepRec);
1486  maxd0xd0 = hRECd0xd01D->GetMaximum();
1487  hRECd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1488  SetHistoProperties(hRECd0xd01D,4,3,4,20,4,3004);
1489  hRECd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
1490  hRECd0xd01D->Sumw2();
1491  hRECd0xd01D->SetMinimum(0.01);
1492  hRECd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD");
1493  hRECd0xd01D->Draw("hist");
1494  cd0xd0->cd();
1495  cd0xd0->SetLeftMargin(0.15);
1496  cd0xd0->SetRightMargin(0.05);
1497  hRECd0xd01D->Draw("hist");
1498  hRECd0xd01D->Draw("err same");
1499  cd0xd0->Update();
1500 
1501  cmc4->cd(6);
1502  TH1D *hRECAccd0xd01D = data->ShowProjection(id0xd0, stepRecAcc);
1503  maxd0xd0 = hRECAccd0xd01D->GetMaximum();
1504  hRECAccd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1505  SetHistoProperties(hRECAccd0xd01D,38,3,38,20,38,3007);
1506  hRECAccd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
1507  hRECAccd0xd01D->Sumw2();
1508  hRECAccd0xd01D->SetMinimum(0.01);
1509  hRECAccd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - Acc cuts");
1510  hRECAccd0xd01D->Draw("hist");
1511  cAccd0xd0->cd();
1512  cAccd0xd0->SetLeftMargin(0.15);
1513  cAccd0xd0->SetRightMargin(0.05);
1514  hRECAccd0xd01D->Draw("hist");
1515  hRECAccd0xd01D->Draw("err same");
1516  cAccd0xd0->Update();
1517 
1518  cmc4->cd(7);
1519  TH1D *hRECITSCld0xd01D = data->ShowProjection(id0xd0, stepRecITSClus);
1520  maxd0xd0 = hRECITSCld0xd01D->GetMaximum();
1521  hRECITSCld0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1522  SetHistoProperties(hRECITSCld0xd01D,13,3,13,20,13,3019);
1523  hRECITSCld0xd01D->GetXaxis()->SetTitleOffset(1.2);;
1524  hRECITSCld0xd01D->Sumw2();
1525  hRECITSCld0xd01D->SetMinimum(0.01);
1526  hRECITSCld0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - ITS Cl cuts");
1527  hRECITSCld0xd01D->Draw("hist");
1528  cITSCld0xd0->cd();
1529  cITSCld0xd0->SetLeftMargin(0.15);
1530  cITSCld0xd0->SetRightMargin(0.05);
1531  hRECITSCld0xd01D->Draw("hist");
1532  hRECITSCld0xd01D->Draw("err same");
1533  cITSCld0xd0->Update();
1534 
1535  cmc4->cd(8);
1536  TH1D *hRECPPRd0xd01D = data->ShowProjection(id0xd0, stepRecPPR);
1537  maxd0xd0 = hRECPPRd0xd01D->GetMaximum();
1538  hRECPPRd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
1539  SetHistoProperties(hRECPPRd0xd01D,6,3,6,20,6,3019);
1540  hRECPPRd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
1541  hRECPPRd0xd01D->Sumw2();
1542  hRECPPRd0xd01D->SetMinimum(0.01);
1543  hRECPPRd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - PPR cuts");
1544  hRECPPRd0xd01D->Draw("hist");
1545  cPPRd0xd0->cd();
1546  cPPRd0xd0->SetLeftMargin(0.15);
1547  cPPRd0xd0->SetRightMargin(0.05);
1548  hRECPPRd0xd01D->Draw("hist");
1549  hRECPPRd0xd01D->Draw("err same");
1550  cPPRd0xd0->Update();
1551 
1552  cmc4->cd(9);
1553  TH1D *hMCpointing1D = data->ShowProjection(ipointing, stepGen);
1554  Double_t maxpointing = hMCpointing1D->GetMaximum();
1555  hMCpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1556  hMCpointing1D->Sumw2();
1557  hMCpointing1D->SetMinimum(0.01);
1558  SetHistoProperties(hMCpointing1D,2,3,2,20,2,3005);
1559  hMCpointing1D->GetXaxis()->SetTitleOffset(1.2);
1560  hMCpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC data");
1561  hMCpointing1D->Draw("hist");
1562  cmcpointing->cd();
1563  cmcpointing->SetLeftMargin(0.15);
1564  cmcpointing->SetRightMargin(0.05);
1565  hMCpointing1D->Draw("hist");
1566  hMCpointing1D->Draw("err same");
1567  cmcpointing->Update();
1568 
1569  cmc4->cd(10);
1570  TH1D *hMCAccpointing1D = data->ShowProjection(ipointing, stepGenAcc);
1571  maxpointing = hMCAccpointing1D->GetMaximum();
1572  hMCAccpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1573  SetHistoProperties(hMCAccpointing1D,8,3,8,20,8,3006);
1574  hMCAccpointing1D->GetXaxis()->SetTitleOffset(1.2);;
1575  hMCAccpointing1D->Sumw2();
1576  hMCAccpointing1D->SetMinimum(0.01);
1577  hMCAccpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Acc cut");
1578  hMCAccpointing1D->Draw("hist");
1579  cmcAccpointing->cd();
1580  cmcAccpointing->SetLeftMargin(0.15);
1581  cmcAccpointing->SetRightMargin(0.05);
1582  hMCAccpointing1D->Draw("hist");
1583  hMCAccpointing1D->Draw("err same");
1584  cmcAccpointing->Update();
1585 
1586  cmc4->cd(11);
1587  TH1D *hMCVertexpointing1D = data->ShowProjection(ipointing, stepGenVertex);
1588  maxpointing = hMCVertexpointing1D->GetMaximum();
1589  hMCVertexpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1590  SetHistoProperties(hMCVertexpointing1D,49,3,49,20,49,3006);
1591  hMCVertexpointing1D->GetXaxis()->SetTitleOffset(1.2);;
1592  hMCVertexpointing1D->Sumw2();
1593  hMCVertexpointing1D->SetMinimum(0.01);
1594  hMCVertexpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Vertex cut");
1595  hMCVertexpointing1D->Draw("hist");
1596  cmcVertexpointing->cd();
1597  cmcVertexpointing->SetLeftMargin(0.15);
1598  cmcVertexpointing->SetRightMargin(0.05);
1599  hMCVertexpointing1D->Draw("hist");
1600  hMCVertexpointing1D->Draw("err same");
1601  cmcVertexpointing->Update();
1602 
1603  cmc4->cd(12);
1604  TH1D *hMCRefitpointing1D = data->ShowProjection(ipointing, stepGenRefit);
1605  maxpointing = hMCRefitpointing1D->GetMaximum();
1606  hMCRefitpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1607  SetHistoProperties(hMCRefitpointing1D,42,3,42,20,42,3006);
1608  hMCRefitpointing1D->GetXaxis()->SetTitleOffset(1.2);;
1609  hMCRefitpointing1D->Sumw2();
1610  hMCRefitpointing1D->SetMinimum(0.01);
1611  hMCRefitpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Refit cut");
1612  hMCRefitpointing1D->Draw("hist");
1613  cmcRefitpointing->cd();
1614  cmcRefitpointing->SetLeftMargin(0.15);
1615  cmcRefitpointing->SetRightMargin(0.05);
1616  hMCRefitpointing1D->Draw("hist");
1617  hMCRefitpointing1D->Draw("err same");
1618  cmcRefitpointing->Update();
1619 
1620  cmc4->cd(13);
1621  TH1D *hRECpointing1D = data->ShowProjection(ipointing, stepRec);
1622  maxpointing = hRECpointing1D->GetMaximum();
1623  hRECpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1624  SetHistoProperties(hRECpointing1D,4,3,4,20,4,3004);
1625  hRECpointing1D->GetXaxis()->SetTitleOffset(1.2);;
1626  hRECpointing1D->Sumw2();
1627  hRECpointing1D->SetMinimum(0.01);
1628  hRECpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD");
1629  hRECpointing1D->Draw("hist");
1630  cpointing->cd();
1631  cpointing->SetLeftMargin(0.15);
1632  cpointing->SetRightMargin(0.05);
1633  hRECpointing1D->Draw("hist");
1634  hRECpointing1D->Draw("err same");
1635  cpointing->Update();
1636 
1637  cmc4->cd(14);
1638  TH1D *hRECAccpointing1D = data->ShowProjection(ipointing, stepRecAcc);
1639  maxpointing = hRECAccpointing1D->GetMaximum();
1640  hRECAccpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1641  SetHistoProperties(hRECAccpointing1D,38,3,38,20,38,3007);
1642  hRECAccpointing1D->GetXaxis()->SetTitleOffset(1.2);;
1643  hRECAccpointing1D->Sumw2();
1644  hRECAccpointing1D->SetMinimum(0.01);
1645  hRECAccpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - Acc cuts");
1646  hRECAccpointing1D->Draw("hist");
1647  cAccpointing->cd();
1648  cAccpointing->SetLeftMargin(0.15);
1649  cAccpointing->SetRightMargin(0.05);
1650  hRECAccpointing1D->Draw("hist");
1651  hRECAccpointing1D->Draw("err same");
1652  cAccpointing->Update();
1653 
1654  cmc4->cd(15);
1655  TH1D *hRECITSClpointing1D = data->ShowProjection(ipointing, stepRecITSClus);
1656  maxpointing = hRECITSClpointing1D->GetMaximum();
1657  hRECITSClpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1658  SetHistoProperties(hRECITSClpointing1D,13,3,13,20,13,3019);
1659  hRECITSClpointing1D->GetXaxis()->SetTitleOffset(1.2);;
1660  hRECITSClpointing1D->Sumw2();
1661  hRECITSClpointing1D->SetMinimum(0.01);
1662  hRECITSClpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - ITS Cl cuts");
1663  hRECITSClpointing1D->Draw("hist");
1664  cITSClpointing->cd();
1665  cITSClpointing->SetLeftMargin(0.15);
1666  cITSClpointing->SetRightMargin(0.05);
1667  hRECITSClpointing1D->Draw("hist");
1668  hRECITSClpointing1D->Draw("err same");
1669  cITSClpointing->Update();
1670 
1671  cmc4->cd(16);
1672  TH1D *hRECPPRpointing1D = data->ShowProjection(ipointing, stepRecPPR);
1673  maxpointing = hRECPPRpointing1D->GetMaximum();
1674  hRECPPRpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
1675  SetHistoProperties(hRECPPRpointing1D,6,3,6,20,6,3019);
1676  hRECPPRpointing1D->GetXaxis()->SetTitleOffset(1.2);;
1677  hRECPPRpointing1D->Sumw2();
1678  hRECPPRpointing1D->SetMinimum(0.01);
1679  hRECPPRpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - PPR cuts");
1680  hRECPPRpointing1D->Draw("hist");
1681  cPPRpointing->cd();
1682  cPPRpointing->SetLeftMargin(0.15);
1683  cPPRpointing->SetRightMargin(0.05);
1684  hRECPPRpointing1D->Draw("hist");
1685  hRECPPRpointing1D->Draw("err same");
1686  cPPRpointing->Update();
1687 
1688  cmc4->cd(17);
1689  TH1D *hMCphi1D = data->ShowProjection(iphi, stepGen);
1690  Double_t maxphi = hMCphi1D->GetMaximum();
1691  hMCphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1692  hMCphi1D->Sumw2();
1693  hMCphi1D->SetMinimum(0.01);
1694  SetHistoProperties(hMCphi1D,2,3,2,20,2,3005);
1695  hMCphi1D->GetXaxis()->SetTitleOffset(1.2);
1696  hMCphi1D->GetXaxis()->SetTitle("Phi, MC data");
1697  hMCphi1D->Draw("hist");
1698  cmcphi->cd();
1699  cmcphi->SetLeftMargin(0.15);
1700  cmcphi->SetRightMargin(0.05);
1701  hMCphi1D->Draw("hist");
1702  hMCphi1D->Draw("err same");
1703  cmcphi->Update();
1704 
1705  cmc4->cd(18);
1706  TH1D *hMCAccphi1D = data->ShowProjection(iphi, stepGenAcc);
1707  maxphi = hMCAccphi1D->GetMaximum();
1708  hMCAccphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1709  SetHistoProperties(hMCAccphi1D,8,3,8,20,8,3006);
1710  hMCAccphi1D->GetXaxis()->SetTitleOffset(1.2);;
1711  hMCAccphi1D->Sumw2();
1712  hMCAccphi1D->SetMinimum(0.01);
1713  hMCAccphi1D->GetXaxis()->SetTitle("Phi, MC - Acc cut");
1714  hMCAccphi1D->Draw("hist");
1715  cmcAccphi->cd();
1716  cmcAccphi->SetLeftMargin(0.15);
1717  cmcAccphi->SetRightMargin(0.05);
1718  hMCAccphi1D->Draw("hist");
1719  hMCAccphi1D->Draw("err same");
1720  cmcAccphi->Update();
1721 
1722  cmc4->cd(19);
1723  TH1D *hMCVertexphi1D = data->ShowProjection(iphi, stepGenVertex);
1724  maxphi = hMCVertexphi1D->GetMaximum();
1725  hMCVertexphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1726  SetHistoProperties(hMCVertexphi1D,49,3,49,20,49,3006);
1727  hMCVertexphi1D->GetXaxis()->SetTitleOffset(1.2);;
1728  hMCVertexphi1D->Sumw2();
1729  hMCVertexphi1D->SetMinimum(0.01);
1730  hMCVertexphi1D->GetXaxis()->SetTitle("Phi, MC - Vertex cut");
1731  hMCVertexphi1D->Draw("hist");
1732  cmcVertexphi->cd();
1733  cmcVertexphi->SetLeftMargin(0.15);
1734  cmcVertexphi->SetRightMargin(0.05);
1735  hMCVertexphi1D->Draw("hist");
1736  hMCVertexphi1D->Draw("err same");
1737  cmcVertexphi->Update();
1738 
1739  cmc4->cd(20);
1740  TH1D *hMCRefitphi1D = data->ShowProjection(iphi, stepGenRefit);
1741  maxphi = hMCRefitphi1D->GetMaximum();
1742  hMCRefitphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1743  SetHistoProperties(hMCRefitphi1D,42,3,42,20,42,3006);
1744  hMCRefitphi1D->GetXaxis()->SetTitleOffset(1.2);;
1745  hMCRefitphi1D->Sumw2();
1746  hMCRefitphi1D->SetMinimum(0.01);
1747  hMCRefitphi1D->GetXaxis()->SetTitle("Phi, MC - Refit cut");
1748  hMCRefitphi1D->Draw("hist");
1749  cmcRefitphi->cd();
1750  cmcRefitphi->SetLeftMargin(0.15);
1751  cmcRefitphi->SetRightMargin(0.05);
1752  hMCRefitphi1D->Draw("hist");
1753  hMCRefitphi1D->Draw("err same");
1754  cmcRefitphi->Update();
1755 
1756  cmc4->cd(21);
1757  TH1D *hRECphi1D = data->ShowProjection(iphi, stepRec);
1758  maxphi = hRECphi1D->GetMaximum();
1759  hRECphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1760  SetHistoProperties(hRECphi1D,4,3,4,20,4,3004);
1761  hRECphi1D->GetXaxis()->SetTitleOffset(1.2);;
1762  hRECphi1D->Sumw2();
1763  hRECphi1D->SetMinimum(0.01);
1764  hRECphi1D->GetXaxis()->SetTitle("Phi, AOD");
1765  hRECphi1D->Draw("hist");
1766  cphi->cd();
1767  cphi->SetLeftMargin(0.15);
1768  cphi->SetRightMargin(0.05);
1769  hRECphi1D->Draw("hist");
1770  hRECphi1D->Draw("err same");
1771  cphi->Update();
1772 
1773  cmc4->cd(22);
1774  TH1D *hRECAccphi1D = data->ShowProjection(iphi, stepRecAcc);
1775  maxphi = hRECAccphi1D->GetMaximum();
1776  hRECAccphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1777  SetHistoProperties(hRECAccphi1D,38,3,38,20,38,3007);
1778  hRECAccphi1D->GetXaxis()->SetTitleOffset(1.2);;
1779  hRECAccphi1D->Sumw2();
1780  hRECAccphi1D->SetMinimum(0.01);
1781  hRECAccphi1D->GetXaxis()->SetTitle("Phi, AOD - Acc cuts");
1782  hRECAccphi1D->Draw("hist");
1783  cAccphi->cd();
1784  cAccphi->SetLeftMargin(0.15);
1785  cAccphi->SetRightMargin(0.05);
1786  hRECAccphi1D->Draw("hist");
1787  hRECAccphi1D->Draw("err same");
1788  cAccphi->Update();
1789 
1790  cmc4->cd(23);
1791  TH1D *hRECITSClphi1D = data->ShowProjection(iphi, stepRecITSClus);
1792  maxphi = hRECITSClphi1D->GetMaximum();
1793  hRECITSClphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1794  SetHistoProperties(hRECITSClphi1D,13,3,13,20,13,3019);
1795  hRECITSClphi1D->GetXaxis()->SetTitleOffset(1.2);;
1796  hRECITSClphi1D->Sumw2();
1797  hRECITSClphi1D->SetMinimum(0.01);
1798  hRECITSClphi1D->GetXaxis()->SetTitle("Phi, AOD - ITS Cl cuts");
1799  hRECITSClphi1D->Draw("hist");
1800  cITSClphi->cd();
1801  cITSClphi->SetLeftMargin(0.15);
1802  cITSClphi->SetRightMargin(0.05);
1803  hRECITSClphi1D->Draw("hist");
1804  hRECITSClphi1D->Draw("err same");
1805  cITSClphi->Update();
1806 
1807  cmc4->cd(24);
1808  TH1D *hRECPPRphi1D = data->ShowProjection(iphi, stepRecPPR);
1809  maxphi = hRECPPRphi1D->GetMaximum();
1810  hRECPPRphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
1811  SetHistoProperties(hRECPPRphi1D,6,3,6,20,6,3019);
1812  hRECPPRphi1D->GetXaxis()->SetTitleOffset(1.2);;
1813  hRECPPRphi1D->Sumw2();
1814  hRECPPRphi1D->SetMinimum(0.01);
1815  hRECPPRphi1D->GetXaxis()->SetTitle("Phi, AOD - PPR cuts");
1816  hRECPPRphi1D->Draw("hist");
1817  cPPRphi->cd();
1818  cPPRphi->SetLeftMargin(0.15);
1819  cPPRphi->SetRightMargin(0.05);
1820  hRECPPRphi1D->Draw("hist");
1821  hRECPPRphi1D->Draw("err same");
1822  cPPRphi->Update();
1823 
1824  // printing on eps files
1825  cmc1->Print("Plots/dataMC_pt_y_cTS.eps");
1826  cmc2->Print("Plots/dataMC_ptPi_ptK_cT.eps");
1827  cmc3->Print("Plots/dataMC_dca_d0pi_d0K.eps");
1828  cmc4->Print("Plots/dataMC_d0xd0_pointing.eps");
1829  cmcpt->Print("Plots/pt_Gen.eps");
1830  cmcy->Print("Plots/y_Gen.eps");
1831  cmccTS->Print("Plots/cTS_Gen.eps");
1832  cmcptPi->Print("Plots/ptPi_Gen.eps");
1833  cmcptK->Print("Plots/ptK_Gen.eps");
1834  cmccT->Print("Plots/cT_Gen.eps");
1835  cmcdca->Print("Plots/dca_Gen.eps");
1836  cmcd0pi->Print("Plots/d0pi_Gen.eps");
1837  cmcd0K->Print("Plots/d0K_Gen.eps");
1838  cmcd0xd0->Print("Plots/d0xd0_Gen.eps");
1839  cmcpointing->Print("Plots/pointing_Gen.eps");
1840  cmcphi->Print("Plots/phi_Gen.eps");
1841 
1842  cmcAccpt->Print("Plots/pt_Gen_Acc.eps");
1843  cmcAccy->Print("Plots/y_Gen_Acc.eps");
1844  cmcAcccTS->Print("Plots/cTS_Gen_Acc.eps");
1845  cmcAccptPi->Print("Plots/ptPi_Gen_Acc.eps");
1846  cmcAccptK->Print("Plots/ptK_Gen_Acc.eps");
1847  cmcAcccT->Print("Plots/cT_Gen_Acc.eps");
1848  cmcAccdca->Print("Plots/dca_Gen_Acc.eps");
1849  cmcAccd0pi->Print("Plots/d0pi_Gen_Acc.eps");
1850  cmcAccd0K->Print("Plots/d0K_Gen_Acc.eps");
1851  cmcAccd0xd0->Print("Plots/d0xd0_Gen_Acc.eps");
1852  cmcAccpointing->Print("Plots/pointing_Gen_Acc.eps");
1853  cmcAccphi->Print("Plots/phi_Gen_Acc.eps");
1854 
1855  cmcVertexpt->Print("Plots/pt_Gen_Vertex.eps");
1856  cmcVertexy->Print("Plots/y_Gen_Vertex.eps");
1857  cmcVertexcTS->Print("Plots/cTS_Gen_Vertex.eps");
1858  cmcVertexptPi->Print("Plots/ptPi_Gen_Vertex.eps");
1859  cmcVertexptK->Print("Plots/ptK_Gen_Vertex.eps");
1860  cmcVertexcT->Print("Plots/cT_Gen_Vertex.eps");
1861  cmcVertexdca->Print("Plots/dca_Gen_Vertex.eps");
1862  cmcVertexd0pi->Print("Plots/d0pi_Gen_Vertex.eps");
1863  cmcVertexd0K->Print("Plots/d0K_Gen_Vertex.eps");
1864  cmcVertexd0xd0->Print("Plots/d0xd0_Gen_Vertex.eps");
1865  cmcVertexpointing->Print("Plots/pointing_Gen_Vertex.eps");
1866  cmcVertexphi->Print("Plots/phi_Gen_Vertex.eps");
1867 
1868  cmcRefitpt->Print("Plots/pt_Gen_Refit.eps");
1869  cmcRefity->Print("Plots/y_Gen_Refit.eps");
1870  cmcRefitcTS->Print("Plots/cTS_Gen_Refit.eps");
1871  cmcRefitptPi->Print("Plots/ptPi_Gen_Refit.eps");
1872  cmcRefitptK->Print("Plots/ptK_Gen_Refit.eps");
1873  cmcRefitcT->Print("Plots/cT_Gen_Refit.eps");
1874  cmcRefitdca->Print("Plots/dca_Gen_Refit.eps");
1875  cmcRefitd0pi->Print("Plots/d0pi_Gen_Refit.eps");
1876  cmcRefitd0K->Print("Plots/d0K_Gen_Refit.eps");
1877  cmcRefitd0xd0->Print("Plots/d0xd0_Gen_Refit.eps");
1878  cmcRefitpointing->Print("Plots/pointing_Gen_Refit.eps");
1879  cmcRefitphi->Print("Plots/phi_Gen_Refit.eps");
1880 
1881  cpt->Print("Plots/pt_Rec.eps");
1882  cy->Print("Plots/y_Rec.eps");
1883  ccTS->Print("Plots/cTS_Rec.eps");
1884  cptPi->Print("Plots/ptPi_Rec.eps");
1885  cptK->Print("Plots/ptK_Rec.eps");
1886  ccT->Print("Plots/cT_Rec.eps");
1887  cdca->Print("Plots/dca_Rec.eps");
1888  cd0pi->Print("Plots/d0pi_Rec.eps");
1889  cd0K->Print("Plots/d0K_Rec.eps");
1890  cd0xd0->Print("Plots/d0xd0_Rec.eps");
1891  cpointing->Print("Plots/pointing_Rec.eps");
1892  cphi->Print("Plots/phi_Rec.eps");
1893 
1894  cAccpt->Print("Plots/pt_Rec_Acc.eps");
1895  cAccy->Print("Plots/y_Rec_Acc.eps");
1896  cAcccTS->Print("Plots/cTS_Rec_Acc.eps");
1897  cAccptPi->Print("Plots/ptPi_Rec_Acc.eps");
1898  cAccptK->Print("Plots/ptK_Rec_Acc.eps");
1899  cAcccT->Print("Plots/cT_Rec_Acc.eps");
1900  cmcdca->Print("Plots/dca_Rec_Acc.eps");
1901  cAccd0pi->Print("Plots/d0pi_Rec_Acc.eps");
1902  cAccd0K->Print("Plots/d0K_Rec_Acc.eps");
1903  cAccd0xd0->Print("Plots/d0xd0_Rec_Acc.eps");
1904  cAccpointing->Print("Plots/pointing_Rec_Acc.eps");
1905  cAccphi->Print("Plots/phi_Rec_Acc.eps");
1906 
1907  cITSClpt->Print("Plots/pt_Rec_ITSCl.eps");
1908  cITSCly->Print("Plots/y_Rec_ITSCl.eps");
1909  cITSClcTS->Print("Plots/cTS_Rec_ITSCl.eps");
1910  cITSClptPi->Print("Plots/ptPi_Rec_ITSCl.eps");
1911  cITSClptK->Print("Plots/ptK_Rec_ITSCl.eps");
1912  cITSClcT->Print("Plots/cT_Rec_ITSCl.eps");
1913  cITSCldca->Print("Plots/dca_Rec_ITSCl.eps");
1914  cITSCld0pi->Print("Plots/d0pi_Rec_ITSCl.eps");
1915  cITSCld0K->Print("Plots/d0K_Rec_ITSCl.eps");
1916  cITSCld0xd0->Print("Plots/d0xd0_Rec_ITSCl.eps");
1917  cITSClpointing->Print("Plots/pointing_Rec_ITSCl.eps");
1918  cITSClphi->Print("Plots/phi_Rec_ITSCl.eps");
1919 
1920  cPPRpt->Print("Plots/pt_Rec_PPR.eps");
1921  cPPRy->Print("Plots/y_Rec_PPR.eps");
1922  cPPRcTS->Print("Plots/cTS_Rec_PPR.eps");
1923  cPPRptPi->Print("Plots/ptPi_Rec_PPR.eps");
1924  cPPRptK->Print("Plots/ptK_Rec_PPR.eps");
1925  cPPRcT->Print("Plots/cT_Rec_PPR.eps");
1926  cPPRdca->Print("Plots/dca_Rec_PPR.eps");
1927  cPPRd0pi->Print("Plots/d0pi_Rec_PPR.eps");
1928  cPPRd0K->Print("Plots/d0K_Rec_PPR.eps");
1929  cPPRd0xd0->Print("Plots/d0xd0_Rec_PPR.eps");
1930  cPPRpointing->Print("Plots/pointing_Rec_PPR.eps");
1931  cPPRphi->Print("Plots/phi_Rec_PPR.eps");
1932 
1933  // printing on gif files
1934  cmc1->Print("Plots/dataMC_pt_y_cTS.gif");
1935  cmc2->Print("Plots/dataMC_ptPi_ptK_cT.gif");
1936  cmc3->Print("Plots/dataMC_dca_d0pi_d0K.gif");
1937  cmc4->Print("Plots/dataMC_d0xd0_pointing.gif");
1938  cmcpt->Print("Plots/pt_Gen.gif");
1939  cmcy->Print("Plots/y_Gen.gif");
1940  cmccTS->Print("Plots/cTS_Gen.gif");
1941  cmcptPi->Print("Plots/ptPi_Gen.gif");
1942  cmcptK->Print("Plots/ptK_Gen.gif");
1943  cmccT->Print("Plots/cT_Gen.gif");
1944  cmcdca->Print("Plots/dca_Gen.gif");
1945  cmcd0pi->Print("Plots/d0pi_Gen.gif");
1946  cmcd0K->Print("Plots/d0K_Gen.gif");
1947  cmcd0xd0->Print("Plots/d0xd0_Gen.gif");
1948  cmcpointing->Print("Plots/pointing_Gen.gif");
1949  cmcphi->Print("Plots/phi_Gen.gif");
1950 
1951  cmcAccpt->Print("Plots/pt_Gen_Acc.gif");
1952  cmcAccy->Print("Plots/y_Gen_Acc.gif");
1953  cmcAcccTS->Print("Plots/cTS_Gen_Acc.gif");
1954  cmcAccptPi->Print("Plots/ptPi_Gen_Acc.gif");
1955  cmcAccptK->Print("Plots/ptK_Gen_Acc.gif");
1956  cmcAcccT->Print("Plots/cT_Gen_Acc.gif");
1957  cmcAccdca->Print("Plots/dca_Gen_Acc.gif");
1958  cmcAccd0pi->Print("Plots/d0pi_Gen_Acc.gif");
1959  cmcAccd0K->Print("Plots/d0K_Gen_Acc.gif");
1960  cmcAccd0xd0->Print("Plots/d0xd0_Gen_Acc.gif");
1961  cmcAccpointing->Print("Plots/pointing_Gen_Acc.gif");
1962  cmcAccphi->Print("Plots/phi_Gen_Acc.gif");
1963 
1964  cmcVertexpt->Print("Plots/pt_Gen_Vertex.gif");
1965  cmcVertexy->Print("Plots/y_Gen_Vertex.gif");
1966  cmcVertexcTS->Print("Plots/cTS_Gen_Vertex.gif");
1967  cmcVertexptPi->Print("Plots/ptPi_Gen_Vertex.gif");
1968  cmcVertexptK->Print("Plots/ptK_Gen_Vertex.gif");
1969  cmcVertexcT->Print("Plots/cT_Gen_Vertex.gif");
1970  cmcVertexdca->Print("Plots/dca_Gen_Vertex.gif");
1971  cmcVertexd0pi->Print("Plots/d0pi_Gen_Vertex.gif");
1972  cmcVertexd0K->Print("Plots/d0K_Gen_Vertex.gif");
1973  cmcVertexd0xd0->Print("Plots/d0xd0_Gen_Vertex.gif");
1974  cmcVertexpointing->Print("Plots/pointing_Gen_Vertex.gif");
1975  cmcVertexphi->Print("Plots/phi_Gen_Vertex.gif");
1976 
1977  cmcRefitpt->Print("Plots/pt_Gen_Refit.gif");
1978  cmcRefity->Print("Plots/y_Gen_Refit.gif");
1979  cmcRefitcTS->Print("Plots/cTS_Gen_Refit.gif");
1980  cmcRefitptPi->Print("Plots/ptPi_Gen_Refit.gif");
1981  cmcRefitptK->Print("Plots/ptK_Gen_Refit.gif");
1982  cmcRefitcT->Print("Plots/cT_Gen_Refit.gif");
1983  cmcRefitdca->Print("Plots/dca_Gen_Refit.gif");
1984  cmcRefitd0pi->Print("Plots/d0pi_Gen_Refit.gif");
1985  cmcRefitd0K->Print("Plots/d0K_Gen_Refit.gif");
1986  cmcRefitd0xd0->Print("Plots/d0xd0_Gen_Refit.gif");
1987  cmcRefitpointing->Print("Plots/pointing_Gen_Refit.gif");
1988  cmcRefitphi->Print("Plots/phi_Gen_Refit.gif");
1989 
1990  cpt->Print("Plots/pt_Rec.gif");
1991  cy->Print("Plots/y_Rec.gif");
1992  ccTS->Print("Plots/cTS_Rec.gif");
1993  cptPi->Print("Plots/ptPi_Rec.gif");
1994  cptK->Print("Plots/ptK_Rec.gif");
1995  ccT->Print("Plots/cT_Rec.gif");
1996  cdca->Print("Plots/dca_Rec.gif");
1997  cd0pi->Print("Plots/d0pi_Rec.gif");
1998  cd0K->Print("Plots/d0K_Rec.gif");
1999  cd0xd0->Print("Plots/d0xd0_Rec.gif");
2000  cpointing->Print("Plots/pointing_Rec.gif");
2001  cphi->Print("Plots/phi_Rec.gif");
2002 
2003  cAccpt->Print("Plots/pt_Rec_Acc.gif");
2004  cAccy->Print("Plots/y_Rec_Acc.gif");
2005  cAcccTS->Print("Plots/cTS_Rec_Acc.gif");
2006  cAccptPi->Print("Plots/ptPi_Rec_Acc.gif");
2007  cAccptK->Print("Plots/ptK_Rec_Acc.gif");
2008  cAcccT->Print("Plots/cT_Rec_Acc.gif");
2009  cmcdca->Print("Plots/dca_Rec_Acc.gif");
2010  cAccd0pi->Print("Plots/d0pi_Rec_Acc.gif");
2011  cAccd0K->Print("Plots/d0K_Rec_Acc.gif");
2012  cAccd0xd0->Print("Plots/d0xd0_Rec_Acc.gif");
2013  cAccpointing->Print("Plots/pointing_Rec_Acc.gif");
2014  cAccphi->Print("Plots/phi_Rec_Acc.gif");
2015 
2016  cITSClpt->Print("Plots/pt_Rec_ITSCl.gif");
2017  cITSCly->Print("Plots/y_Rec_ITSCl.gif");
2018  cITSClcTS->Print("Plots/cTS_Rec_ITSCl.gif");
2019  cITSClptPi->Print("Plots/ptPi_Rec_ITSCl.gif");
2020  cITSClptK->Print("Plots/ptK_Rec_ITSCl.gif");
2021  cITSClcT->Print("Plots/cT_Rec_ITSCl.gif");
2022  cITSCldca->Print("Plots/dca_Rec_ITSCl.gif");
2023  cITSCld0pi->Print("Plots/d0pi_Rec_ITSCl.gif");
2024  cITSCld0K->Print("Plots/d0K_Rec_ITSCl.gif");
2025  cITSCld0xd0->Print("Plots/d0xd0_Rec_ITSCl.gif");
2026  cITSClpointing->Print("Plots/pointing_Rec_ITSCl.gif");
2027  cITSClphi->Print("Plots/phi_Rec_ITSCl.gif");
2028 
2029  cPPRpt->Print("Plots/pt_Rec_PPR.gif");
2030  cPPRy->Print("Plots/y_Rec_PPR.gif");
2031  cPPRcTS->Print("Plots/cTS_Rec_PPR.gif");
2032  cPPRptPi->Print("Plots/ptPi_Rec_PPR.gif");
2033  cPPRptK->Print("Plots/ptK_Rec_PPR.gif");
2034  cPPRcT->Print("Plots/cT_Rec_PPR.gif");
2035  cPPRdca->Print("Plots/dca_Rec_PPR.gif");
2036  cPPRd0pi->Print("Plots/d0pi_Rec_PPR.gif");
2037  cPPRd0K->Print("Plots/d0K_Rec_PPR.gif");
2038  cPPRd0xd0->Print("Plots/d0xd0_Rec_PPR.gif");
2039  cPPRpointing->Print("Plots/pointing_Rec_PPR.gif");
2040  cPPRphi->Print("Plots/phi_Rec_PPR.gif");
2041 
2042  TFile* file_histo = new TFile("fileHisto.root","RECREATE");
2043  // MC
2044  hMCpt1D->Write("hMCpt1D");
2045  hMCy1D->Write("hMCy1D");
2046  hMCcTS1D->Write("hMCcTS1D");
2047  hMCptPi1D->Write("hMCptPi1D");
2048  hMCptK1D->Write("hMCptK1D");
2049  hMCcT1D->Write("hMCcT1D");
2050  hMCdca1D->Write("hMCdca1D");
2051  hMCd0pi1D->Write("hMCd0pi1D");
2052  hMCd0K1D->Write("hMCd0K1D");
2053  hMCd0xd01D->Write("hMCd0xd01D");
2054  hMCpointing1D->Write("hMCpointing1D");
2055  hMCphi1D->Write("hMCphi1D");
2056  // MC, Acc
2057  hMCAccpt1D->Write("hMCAccpt1D");
2058  hMCAccy1D->Write("hMCAccy1D");
2059  hMCAcccTS1D->Write("hMCAcccTS1D");
2060  hMCAccptPi1D->Write("hMCAccptPi1D");
2061  hMCAccptK1D->Write("hMCAccptK1D");
2062  hMCAcccT1D->Write("hMCAcccT1D");
2063  hMCAccdca1D->Write("hMCAccdca1D");
2064  hMCAccd0pi1D->Write("hMCAccd0pi1D");
2065  hMCAccd0K1D->Write("hMCAccd0K1D");
2066  hMCAccd0xd01D->Write("hMCAccd0xd01D");
2067  hMCAccpointing1D->Write("hMCAccpointing1D");
2068  hMCAccphi1D->Write("hMCAccphi1D");
2069  // REC
2070  hRECpt1D->Write("hRECpt1D");
2071  hRECy1D->Write("hRECy1D");
2072  hRECcTS1D->Write("hRECcTS1D");
2073  hRECptPi1D->Write("hRECptPi1D");
2074  hRECptK1D->Write("hRECptK1D");
2075  hRECcT1D->Write("hRECcT1D");
2076  hRECdca1D->Write("hRECdca1D");
2077  hRECd0pi1D->Write("hRECd0pi1D");
2078  hRECd0K1D->Write("hRECd0K1D");
2079  hRECd0xd01D->Write("hRECd0xd01D");
2080  hRECpointing1D->Write("hRECpointing1D");
2081  hRECphi1D->Write("hRECphi1D");
2082  // REC, Acc
2083  hRECAccpt1D->Write("hRECAccpt1D");
2084  hRECAccy1D->Write("hRECAccy1D");
2085  hRECAcccTS1D->Write("hRECAcccTS1D");
2086  hRECAccptPi1D->Write("hRECAccptPi1D");
2087  hRECAccptK1D->Write("hRECAccptK1D");
2088  hRECAcccT1D->Write("hRECAcccT1D");
2089  hRECAccdca1D->Write("hRECAccdca1D");
2090  hRECAccd0pi1D->Write("hRECAccd0pi1D");
2091  hRECAccd0K1D->Write("hRECAccd0K1D");
2092  hRECAccd0xd01D->Write("hRECAccd0xd01D");
2093  hRECAccpointing1D->Write("hRECAccpointing1D");
2094  hRECAccphi1D->Write("hRECAccphi1D");
2095  // REC, ITS Cl
2096  hRECITSClpt1D->Write("hRECITSClpt1D");
2097  hRECITSCly1D->Write("hRECITSCly1D");
2098  hRECITSClcTS1D->Write("hRECITSClcTS1D");
2099  hRECITSClptPi1D->Write("hRECITSClptPi1D");
2100  hRECITSClptK1D->Write("hRECITSClptK1D");
2101  hRECITSClcT1D->Write("hRECITSClcT1D");
2102  hRECITSCldca1D->Write("hRECITSCldca1D");
2103  hRECITSCld0pi1D->Write("hRECITSCld0pi1D");
2104  hRECITSCld0K1D->Write("hRECITSCld0K1D");
2105  hRECITSCld0xd01D->Write("hRECITSCld0xd01D");
2106  hRECITSClpointing1D->Write("hRECITSClpointing1D");
2107  hRECITSClphi1D->Write("hRECITSClphi1D");
2108  // REC, PPR
2109  hRECPPRpt1D->Write("hRECPPRpt1D");
2110  hRECPPRy1D->Write("hRECPPRy1D");
2111  hRECPPRcTS1D->Write("hRECPPRcTS1D");
2112  hRECPPRptPi1D->Write("hRECPPRptPi1D");
2113  hRECPPRptK1D->Write("hRECPPRptK1D");
2114  hRECPPRcT1D->Write("hRECPPRcT1D");
2115  hRECPPRdca1D->Write("hRECPPRdca1D");
2116  hRECPPRd0pi1D->Write("hRECPPRd0pi1D");
2117  hRECPPRd0K1D->Write("hRECPPRd0K1D");
2118  hRECPPRd0xd01D->Write("hRECPPRd0xd01D");
2119  hRECPPRpointing1D->Write("hRECPPRpointing1D");
2120  hRECPPRphi1D->Write("hRECPPRphi1D");
2121 
2122 }
2123 
2124 void SetHistoProperties(TH1D* h, Int_t lineColor, Int_t lineWidth, Int_t markerColor, Int_t markerStyle, Int_t fillColor, Int_t fillStyle){
2125 
2126  h->SetLineColor(lineColor);
2127  h->SetLineWidth(lineWidth);
2128  h->SetMarkerColor(markerColor);
2129  h->SetFillColor(fillColor);
2130  h->SetFillStyle(fillStyle);
2131  h->SetMarkerStyle(markerStyle);
2132 
2133  return;
2134 
2135 }
2136 
TBenchmark * gBenchmark
void ReadCFHeavyFlavourOutput()
TSystem * gSystem
TRandom * gRandom
void SetHistoProperties(TH1D *h, Int_t lineColor, Int_t lineWidth, Int_t markerColor, Int_t markerStyle, Int_t fillColor, Int_t fillStyle)