--- /dev/null
+#include <Riostream.h>
+
+extern TRandom *gRandom;
+extern TBenchmark *gBenchmark;
+extern TSystem *gSystem;
+
+void ReadCFHeavyFlavourOutput() {
+
+ // example macro to read the CF produced for D->Kpi
+ // 8 steps considered
+ // 12 variables
+
+ gROOT->SetStyle("Plain");
+ gStyle->SetPalette(1);
+ gStyle->SetOptStat(0);
+ gStyle->SetPalette(1);
+ gStyle->SetCanvasColor(0);
+ gStyle->SetFrameFillColor(0);
+ gStyle->SetOptTitle(0);
+
+ gSystem->SetIncludePath("-I. -I$ALICE_ROOT/include -I$ROOTSYS/include");
+ gSystem->Load("libANALYSIS.so");
+ gSystem->Load("libANALYSISalice.so");
+ gSystem->Load("$ALICE_ROOT/CORRFW/libCORRFW.so") ;
+
+ // Setting up the container grid...
+
+ // Flag the sel steps
+ Int_t stepGen=0;
+ Int_t stepGenAcc=1;
+ Int_t stepGenVertex=2;
+ Int_t stepGenRefit=3;
+ Int_t stepRec=4;
+ Int_t stepRecAcc=5;
+ Int_t stepRecITSClus=6;
+ Int_t stepRecPPR=7;
+
+ // The sensitive variables, their indeces
+ Int_t ipt =0;
+ Int_t iy=1;
+ UInt_t icTS = 2; //cTS stands for cosThetaStar
+ UInt_t iptPi = 3;
+ UInt_t iptK = 4;
+ UInt_t icT = 5;
+ UInt_t idca = 6;
+ UInt_t id0pi = 7;
+ UInt_t id0K = 8;
+ UInt_t id0xd0 = 9;
+ UInt_t ipointing = 10;
+ UInt_t iphi = 11;
+
+ // Get the container
+ TFile* f = new TFile("AnalysisResults.root");
+ TDirectoryFile* d = (TDirectoryFile*)f->Get("PWG3_D2H_CFtaskD0toKpi");
+ AliCFContainer *data = (AliCFContainer*) (d->Get("container"));
+
+ // Make some 1 & 2-D projections.
+ // MC level
+ TCanvas *cmc1 =new TCanvas("cmc1","The MC distributions",0,0,1600,1200); // canvas for all steps, 1st three variables
+ TCanvas *cmc2 =new TCanvas("cmc2","The MC distributions",0,0,1600,1200); // canvas for all steps, 2nd three variables
+
+ TCanvas *cmc3 =new TCanvas("cmc3","The MC distributions",0,0,1600,1200); // canvas for all steps, 2nd three variables
+ TCanvas *cmc4 =new TCanvas("cmc4","The MC distributions",0,0,1600,765); // canvas for all steps, 2nd three variables
+ TCanvas *cmcpt = new TCanvas("cmcpt","pt distribution from MC",50,50,550,550);
+ TCanvas *cmcy = new TCanvas("cmcy","y distribution from MC",50,50,550,550);
+ TCanvas *cmccTS = new TCanvas("cmcTS","cosThetaStar distribution from MC",50,50,550,550); // cTS stands for cosThetaStar
+ TCanvas *cmcptPi = new TCanvas("cmcptPi","pt_pi distribution from MC",50,50,550,550);
+ TCanvas *cmcptK = new TCanvas("cmcptK","pt_K distribution from MC",50,50,550,550);
+ TCanvas *cmccT = new TCanvas("cmccT","cT distribution from MC",50,50,550,550);
+ TCanvas *cmcdca = new TCanvas("cmcdca","dca distribution from MC",50,50,550,550);
+ TCanvas *cmcd0pi = new TCanvas("cmcd0pi","d0pi distribution from MC",50,50,550,550);
+ TCanvas *cmcd0K = new TCanvas("cmcd0K","d0K distribution from MC",50,50,550,550);
+ TCanvas *cmcd0xd0 = new TCanvas("cmcd0xd0","d0xd0 distribution from MC",50,50,550,550);
+ TCanvas *cmcpointing = new TCanvas("cmcpointing","cosPointingAngle distribution from MC",50,50,550,550);
+ TCanvas *cmcphi = new TCanvas("cmcphi","Phi distribution from MC",50,50,550,550);
+
+ // MC Acceptance level
+ TCanvas *cmcAccpt = new TCanvas("cmcAccpt","pt distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccy = new TCanvas("cmcAccy","y distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAcccTS = new TCanvas("cmcAccTS","cosThetaStar distribution from MC - Acceptance cuts",50,50,550,550); // cTS stands for cosThetaStar
+ TCanvas *cmcAccptPi = new TCanvas("cmcAccptPi","pt_pi distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccptK = new TCanvas("cmcAccptK","pt_K distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAcccT = new TCanvas("cmcAcccT","cT distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccdca = new TCanvas("cmcAccdca","dca distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccd0pi = new TCanvas("cmcAccd0pi","d0pi distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccd0K = new TCanvas("cmcAccd0K","d0K distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccd0xd0 = new TCanvas("cmcAccd0xd0","d0xd0 distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccpointing = new TCanvas("cmcAccpointing","cosPointingAngle distribution from MC - Acceptance cuts",50,50,550,550);
+ TCanvas *cmcAccphi = new TCanvas("cmcAccphi","Phi distribution from MC - Acceptance cuts",50,50,550,550);
+
+ // MC Vertex level
+ TCanvas *cmcVertexpt = new TCanvas("cmcVertexpt","pt distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexy = new TCanvas("cmcVertexy","y distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexcTS = new TCanvas("cmcVertexTS","cosThetaStar distribution from MC - Vertex cuts",50,50,550,550); // cTS stands for cosThetaStar
+ TCanvas *cmcVertexptPi = new TCanvas("cmcVertexptPi","pt_pi distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexptK = new TCanvas("cmcVertexptK","pt_K distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexcT = new TCanvas("cmcVertexcT","cT distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexdca = new TCanvas("cmcVertexdca","dca distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexd0pi = new TCanvas("cmcVertexd0pi","d0pi distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexd0K = new TCanvas("cmcVertexd0K","d0K distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexd0xd0 = new TCanvas("cmcVertexd0xd0","d0xd0 distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexpointing = new TCanvas("cmcVertexpointing","cosPointingAngle distribution from MC - Vertex cuts",50,50,550,550);
+ TCanvas *cmcVertexphi = new TCanvas("cmcVertexphi","Phi distribution from MC - Vertex cuts",50,50,550,550);
+
+ // MC Refit level
+ TCanvas *cmcRefitpt = new TCanvas("cmcRefitpt","pt distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefity = new TCanvas("cmcRefity","y distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitcTS = new TCanvas("cmcRefitTS","cosThetaStar distribution from MC - Refit cuts",50,50,550,550); // cTS stands for cosThetaStar
+ TCanvas *cmcRefitptPi = new TCanvas("cmcRefitptPi","pt_pi distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitptK = new TCanvas("cmcRefitptK","pt_K distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitcT = new TCanvas("cmcRefitcT","cT distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitdca = new TCanvas("cmcRefitdca","dca distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitd0pi = new TCanvas("cmcRefitd0pi","d0pi distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitd0K = new TCanvas("cmcRefitd0K","d0K distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitd0xd0 = new TCanvas("cmcRefitd0xd0","d0xd0 distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitpointing = new TCanvas("cmcRefitpointing","cosPointingAngle distribution from MC - Refit cuts",50,50,550,550);
+ TCanvas *cmcRefitphi = new TCanvas("cmcRefitphi","Phi distribution from MC - Refit cuts",50,50,550,550);
+
+ // Reco-aod level
+ TCanvas *cpt = new TCanvas("cpt","pt distribution from reco aod",50,50,550,550);
+ TCanvas *cy = new TCanvas("cy","y distribution from reco aod",50,50,550,550);
+ TCanvas *ccTS = new TCanvas("cTS","cosThetaStar distribution from reco aod",50,50,550,550);
+ TCanvas *cptPi = new TCanvas("cptPi","pt_pi distribution from reco aod",50,50,550,550);
+ TCanvas *cptK = new TCanvas("cptK","pt_K distribution from reco aod",50,50,550,550);
+ TCanvas *ccT = new TCanvas("ccT","cT distribution from reco aod",50,50,550,550);
+ TCanvas *cdca = new TCanvas("cdca","dca distribution from reco aod",50,50,550,550);
+ TCanvas *cd0pi = new TCanvas("cd0pi","d0pi distribution from reco aod",50,50,550,550);
+ TCanvas *cd0K = new TCanvas("cd0K","d0K distribution from reco aod",50,50,550,550);
+ TCanvas *cd0xd0 = new TCanvas("cd0xd0","d0xd0 distribution from reco aod",50,50,550,550);
+ TCanvas *cpointing = new TCanvas("cpointing","cosPointingAngle distribution from reco aod",50,50,550,550);
+ TCanvas *cphi = new TCanvas("cphi","Phi distribution from reco aod",50,50,550,550);
+
+ // Reco-aod Acceptance level
+ TCanvas *cAccpt = new TCanvas("cAccpt","pt distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccy = new TCanvas("cAccy","y distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAcccTS = new TCanvas("cAccTS","cosThetaStar distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccptPi = new TCanvas("cAccptPi","pt_pi distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccptK = new TCanvas("cAccptK","pt_K distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccdca = new TCanvas("cAccdca","dca distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccd0pi = new TCanvas("cAccd0pi","d0pi distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccd0K = new TCanvas("cAccd0K","d0K distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccd0xd0 = new TCanvas("cAccd0xd0","d0xd0 distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccpointing = new TCanvas("cAccpointing","cosPointingAngle distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAcccT = new TCanvas("cAcccT","cT distribution from reco aod - Acceptance cuts",50,50,550,550);
+ TCanvas *cAccphi = new TCanvas("cAccphi","Phi distribution from reco aod - Acceptance cuts",50,50,550,550);
+
+ // Reco-aod ITS Clusters level
+ TCanvas *cITSClpt = new TCanvas("cITSClpt","pt distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSCly = new TCanvas("cITSCly","y distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSClcTS = new TCanvas("cITSClTS","cosThetaStar distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSClptPi = new TCanvas("cITSClptPi","pt_pi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSClptK = new TCanvas("cITSClptK","pt_K distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSClcT = new TCanvas("cITSClcT","cT distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSCldca = new TCanvas("cITSCldca","dca distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSCld0pi = new TCanvas("cITSCld0pi","d0pi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSCld0K = new TCanvas("cITSCld0K","d0K distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSCld0xd0 = new TCanvas("cITSCld0xd0","d0xd0 distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSClpointing = new TCanvas("cITSClpointing","cosPointingAngle distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+ TCanvas *cITSClphi = new TCanvas("cITSClphi","Phi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
+
+ // Reco-aod PPR level
+ TCanvas *cPPRpt = new TCanvas("cPPRpt","pt distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRy = new TCanvas("cPPRy","y distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRcTS = new TCanvas("cPPRTS","cosThetaStar distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRptPi = new TCanvas("cPPRptPi","pt_pi distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRptK = new TCanvas("cPPRptK","pt_K distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRcT = new TCanvas("cPPRcT","cT distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRdca = new TCanvas("cPPRdca","dca distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRd0pi = new TCanvas("cPPRd0pi","d0pi distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRd0K = new TCanvas("cPPRd0K","d0K distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRd0xd0 = new TCanvas("cPPRd0xd0","d0xd0 distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRpointing = new TCanvas("cPPRpointing","cosPointingAngle distribution from reco aod - PPR cuts",50,50,550,550);
+ TCanvas *cPPRphi = new TCanvas("cPPRphi","Phi distribution from reco aod - PPR cuts",50,50,550,550);
+
+ // pt, y, cosThetaStar for all steps
+ cmc1->Divide(8,3);
+
+ cmc1->cd(1);
+ TH1D *hMCpt1D = data->ShowProjection(ipt, stepGen);
+ Double_t maxpt = hMCpt1D->GetMaximum();
+ hMCpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ hMCpt1D->Sumw2();
+ hMCpt1D->SetMinimum(0.01);
+ SetHistoProperties(hMCpt1D,2,3,2,20,2,3005);
+ hMCpt1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC data");
+ hMCpt1D->Draw("hist");
+ cmcpt->cd();
+ cmcpt->SetLeftMargin(0.15);
+ cmcpt->SetRightMargin(0.05);
+ hMCpt1D->Draw("hist");
+ hMCpt1D->Draw("err same");
+ cmcpt->Update();
+
+ cmc1->cd(2);
+ TH1D *hMCAccpt1D = data->ShowProjection(ipt, stepGenAcc);
+ maxpt = hMCAccpt1D->GetMaximum();
+ hMCAccpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ SetHistoProperties(hMCAccpt1D,8,3,8,20,8,3006);
+ hMCAccpt1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccpt1D->Sumw2();
+ hMCAccpt1D->SetMinimum(0.01);
+ hMCAccpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Acc cut");
+ hMCAccpt1D->Draw("hist");
+ cmcAccpt->cd();
+ cmcAccpt->SetLeftMargin(0.15);
+ cmcAccpt->SetRightMargin(0.05);
+ hMCAccpt1D->Draw("hist");
+ hMCAccpt1D->Draw("err same");
+ cmcAccpt->Update();
+
+ cmc1->cd(3);
+ TH1D *hMCVertexpt1D = data->ShowProjection(ipt, stepGenVertex);
+ maxpt = hMCVertexpt1D->GetMaximum();
+ hMCVertexpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ SetHistoProperties(hMCVertexpt1D,49,3,49,20,49,3006);
+ hMCVertexpt1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexpt1D->Sumw2();
+ hMCVertexpt1D->SetMinimum(0.01);
+ hMCVertexpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Vertex cut");
+ hMCVertexpt1D->Draw("hist");
+ cmcVertexpt->cd();
+ cmcVertexpt->SetLeftMargin(0.15);
+ cmcVertexpt->SetRightMargin(0.05);
+ hMCVertexpt1D->Draw("hist");
+ hMCVertexpt1D->Draw("err same");
+ cmcVertexpt->Update();
+
+ cmc1->cd(4);
+ TH1D *hMCRefitpt1D = data->ShowProjection(ipt, stepGenRefit);
+ maxpt = hMCRefitpt1D->GetMaximum();
+ hMCRefitpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ SetHistoProperties(hMCRefitpt1D,42,3,42,20,42,3006);
+ hMCRefitpt1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitpt1D->Sumw2();
+ hMCRefitpt1D->SetMinimum(0.01);
+ hMCRefitpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Refit cut");
+ hMCRefitpt1D->Draw("hist");
+ cmcRefitpt->cd();
+ cmcRefitpt->SetLeftMargin(0.15);
+ cmcRefitpt->SetRightMargin(0.05);
+ hMCRefitpt1D->Draw("hist");
+ hMCRefitpt1D->Draw("err same");
+ cmcRefitpt->Update();
+
+ cmc1->cd(5);
+ TH1D *hRECpt1D = data->ShowProjection(ipt, stepRec);
+ maxpt = hRECpt1D->GetMaximum();
+ hRECpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ SetHistoProperties(hRECpt1D,4,3,4,20,4,3004);
+ hRECpt1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECpt1D->Sumw2();
+ hRECpt1D->SetMinimum(0.01);
+ hRECpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD");
+ hRECpt1D->Draw("hist");
+ cpt->cd();
+ cpt->SetLeftMargin(0.15);
+ cpt->SetRightMargin(0.05);
+ hRECpt1D->Draw("hist");
+ hRECpt1D->Draw("err same");
+ cpt->Update();
+
+ cmc1->cd(6);
+ TH1D *hRECAccpt1D = data->ShowProjection(ipt, stepRecAcc);
+ maxpt = hRECAccpt1D->GetMaximum();
+ hRECAccpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ SetHistoProperties(hRECAccpt1D,38,3,38,20,38,3007);
+ hRECAccpt1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccpt1D->Sumw2();
+ hRECAccpt1D->SetMinimum(0.01);
+ hRECAccpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - Acc cuts");
+ hRECAccpt1D->Draw("hist");
+ cAccpt->cd();
+ cAccpt->SetLeftMargin(0.15);
+ cAccpt->SetRightMargin(0.05);
+ hRECAccpt1D->Draw("hist");
+ hRECAccpt1D->Draw("err same");
+ cAccpt->Update();
+
+ cmc1->cd(7);
+ TH1D *hRECITSClpt1D = data->ShowProjection(ipt, stepRecITSClus);
+ maxpt = hRECITSClpt1D->GetMaximum();
+ hRECITSClpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ SetHistoProperties(hRECITSClpt1D,13,3,13,20,13,3019);
+ hRECITSClpt1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSClpt1D->Sumw2();
+ hRECITSClpt1D->SetMinimum(0.01);
+ hRECITSClpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - ITS Cl cuts");
+ hRECITSClpt1D->Draw("hist");
+ cITSClpt->cd();
+ cITSClpt->SetLeftMargin(0.15);
+ cITSClpt->SetRightMargin(0.05);
+ hRECITSClpt1D->Draw("hist");
+ hRECITSClpt1D->Draw("err same");
+ cITSClpt->Update();
+
+ cmc1->cd(8);
+ TH1D *hRECPPRpt1D = data->ShowProjection(ipt, stepRecPPR);
+ maxpt = hRECPPRpt1D->GetMaximum();
+ hRECPPRpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
+ SetHistoProperties(hRECPPRpt1D,6,3,6,20,6,3019);
+ hRECPPRpt1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRpt1D->Sumw2();
+ hRECPPRpt1D->SetMinimum(0.01);
+ hRECPPRpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - PPR cuts");
+ hRECPPRpt1D->Draw("hist");
+ cPPRpt->cd();
+ cPPRpt->SetLeftMargin(0.15);
+ cPPRpt->SetRightMargin(0.05);
+ hRECPPRpt1D->Draw("hist");
+ hRECPPRpt1D->Draw("err same");
+ cPPRpt->Update();
+
+ cmc1->cd(9);
+ TH1D *hMCy1D = data->ShowProjection(iy, stepGen);
+ Double_t maxy = hMCy1D->GetMaximum();
+ hMCy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hMCy1D,2,3,2,20,2,3005);
+ hMCy1D->SetLineColor(2);
+ hMCy1D->SetLineWidth(3);
+ hMCy1D->SetMarkerColor(2);
+ hMCy1D->SetFillColor(2);
+ hMCy1D->SetFillStyle(3005);
+ hMCy1D->SetMarkerStyle(20);
+ hMCy1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCy1D->Sumw2();
+ hMCy1D->SetMinimum(0.01);
+ hMCy1D->GetXaxis()->SetTitle("y, MC data");
+ hMCy1D->Draw("hist");
+ cmcy->cd();
+ cmcy->SetLeftMargin(0.15);
+ cmcy->SetRightMargin(0.05);
+ hMCy1D->Draw("hist");
+ hMCy1D->Draw("err same");
+ cmcy->Update();
+
+ cmc1->cd(10);
+ TH1D *hMCAccy1D = data->ShowProjection(iy, stepGenAcc);
+ maxy = hMCAccy1D->GetMaximum();
+ hMCAccy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hMCAccy1D,8,3,8,20,8,3006);
+ hMCAccy1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCAccy1D->Sumw2();
+ hMCAccy1D->SetMinimum(0.01);
+ hMCAccy1D->GetXaxis()->SetTitle("y, MC data - Acc cuts");
+ hMCAccy1D->Draw("hist");
+ cmcAccy->cd();
+ cmcAccy->SetLeftMargin(0.15);
+ cmcAccy->SetRightMargin(0.05);
+ hMCAccy1D->Draw("hist");
+ hMCAccy1D->Draw("err same");
+ cmcy->Update();
+
+ cmc1->cd(11);
+ TH1D *hMCVertexy1D = data->ShowProjection(iy, stepGenVertex);
+ maxy = hMCVertexy1D->GetMaximum();
+ hMCVertexy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hMCVertexy1D,49,3,49,20,49,3006);
+ hMCVertexy1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCVertexy1D->Sumw2();
+ hMCVertexy1D->SetMinimum(0.01);
+ hMCVertexy1D->GetXaxis()->SetTitle("y, MC data - Vertex cuts");
+ hMCVertexy1D->Draw("hist");
+ cmcVertexy->cd();
+ cmcVertexy->SetLeftMargin(0.15);
+ cmcVertexy->SetRightMargin(0.05);
+ hMCVertexy1D->Draw("hist");
+ hMCVertexy1D->Draw("err same");
+ cmcy->Update();
+
+ cmc1->cd(12);
+ TH1D *hMCRefity1D = data->ShowProjection(iy, stepGenRefit);
+ maxy = hMCRefity1D->GetMaximum();
+ hMCRefity1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hMCRefity1D,42,3,42,20,42,3006);
+ hMCRefity1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCRefity1D->Sumw2();
+ hMCRefity1D->SetMinimum(0.01);
+ hMCRefity1D->GetXaxis()->SetTitle("y, MC data - Refit cuts");
+ hMCRefity1D->Draw("hist");
+ cmcRefity->cd();
+ cmcRefity->SetLeftMargin(0.15);
+ cmcRefity->SetRightMargin(0.05);
+ hMCRefity1D->Draw("hist");
+ hMCRefity1D->Draw("err same");
+ cmcy->Update();
+
+ cmc1->cd(13);
+ TH1D *hRECy1D = data->ShowProjection(iy, stepRec);
+ maxy = hRECy1D->GetMaximum();
+ hRECy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hRECy1D,4,3,4,20,4,3004);
+ hRECy1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECy1D->Sumw2();
+ hRECy1D->SetMinimum(0.01);
+ hRECy1D->Draw("hist");
+ cy->cd();
+ cy->SetLeftMargin(0.15);
+ cy->SetRightMargin(0.05);
+ hRECy1D->GetXaxis()->SetTitle("y, AOD");
+ hRECy1D->Draw("hist");
+ hRECy1D->Draw("err same");
+ cy->Update();
+
+ cmc1->cd(14);
+ TH1D *hRECAccy1D = data->ShowProjection(iy, stepRecAcc);
+ maxy = hRECAccy1D->GetMaximum();
+ hRECAccy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hRECAccy1D,38,3,38,20,38,3007);
+ hRECAccy1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccy1D->Sumw2();
+ hRECAccy1D->SetMinimum(0.01);
+ hRECAccy1D->Draw("hist");
+ cAccy->cd();
+ cAccy->SetLeftMargin(0.15);
+ cAccy->SetRightMargin(0.05);
+ hRECAccy1D->GetXaxis()->SetTitle("y, AOD - Acc cuts");
+ hRECAccy1D->Draw("hist");
+ hRECAccy1D->Draw("err same");
+ cAccy->Update();
+
+ cmc1->cd(15);
+ TH1D *hRECITSCly1D = data->ShowProjection(iy, stepRecITSClus);
+ maxy = hRECITSCly1D->GetMaximum();
+ hRECITSCly1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hRECITSCly1D,13,3,13,20,13,3019);
+ hRECITSCly1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSCly1D->Sumw2();
+ hRECITSCly1D->SetMinimum(0.01);
+ hRECITSCly1D->Draw("hist");
+ cITSCly->cd();
+ cITSCly->SetLeftMargin(0.15);
+ cITSCly->SetRightMargin(0.05);
+ hRECITSCly1D->GetXaxis()->SetTitle("y, AOD - ITS Cl cuts");
+ hRECITSCly1D->Draw("hist");
+ hRECITSCly1D->Draw("err same");
+ cITSCly->Update();
+
+ cmc1->cd(16);
+ TH1D *hRECPPRy1D = data->ShowProjection(iy, stepRecPPR);
+ maxy = hRECPPRy1D->GetMaximum();
+ hRECPPRy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
+ SetHistoProperties(hRECPPRy1D,6,3,6,20,6,3019);
+ hRECPPRy1D->SetMarkerStyle(20);
+ hRECPPRy1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRy1D->Sumw2();
+ hRECPPRy1D->SetMinimum(0.01);
+ hRECPPRy1D->Draw("hist");
+ cPPRy->cd();
+ cPPRy->SetLeftMargin(0.15);
+ cPPRy->SetRightMargin(0.05);
+ hRECPPRy1D->GetXaxis()->SetTitle("y, AOD - PPR cuts");
+ hRECPPRy1D->Draw("hist");
+ hRECPPRy1D->Draw("err same");
+ cPPRy->Update();
+
+ cmc1->cd(17);
+ TH1D *hMCcTS1D = data->ShowProjection(icTS, stepGen);
+ Double_t maxcTS = hMCcTS1D->GetMaximum();
+ hMCcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hMCcTS1D,2,3,2,20,2,3005);
+ hMCcTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCcTS1D->Sumw2();
+ hMCcTS1D->SetMinimum(0.01);
+ hMCcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data");
+ hMCcTS1D->Draw("hist");
+ cmccTS->cd();
+ cmccTS->SetLeftMargin(0.15);
+ cmccTS->SetRightMargin(0.05);
+ hMCcTS1D->Draw("hist");
+ hMCcTS1D->Draw("err same");
+ cmccTS->Update();
+
+ cmc1->cd(18);
+ TH1D *hMCAcccTS1D = data->ShowProjection(icTS, stepGenAcc);
+ maxcTS = hMCAcccTS1D->GetMaximum();
+ hMCAcccTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hMCAcccTS1D,8,3,8,20,8,3006);
+ hMCAcccTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAcccTS1D->Sumw2();
+ hMCAcccTS1D->SetMinimum(0.01);
+ hMCAcccTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Acc cuts");
+ hMCAcccTS1D->Draw("hist");
+ cmcAcccTS->cd();
+ cmcAcccTS->SetLeftMargin(0.15);
+ cmcAcccTS->SetRightMargin(0.05);
+ hMCAcccTS1D->Draw("hist");
+ hMCAcccTS1D->Draw("err same");
+ cmcAcccTS->Update();
+
+ cmc1->cd(19);
+ TH1D *hMCVertexcTS1D = data->ShowProjection(icTS, stepGenVertex);
+ maxcTS = hMCVertexcTS1D->GetMaximum();
+ hMCVertexcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hMCVertexcTS1D,49,3,49,20,49,3006);
+ hMCVertexcTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexcTS1D->Sumw2();
+ hMCVertexcTS1D->SetMinimum(0.01);
+ hMCVertexcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Vertex cuts");
+ hMCVertexcTS1D->Draw("hist");
+ cmcVertexcTS->cd();
+ cmcVertexcTS->SetLeftMargin(0.15);
+ cmcVertexcTS->SetRightMargin(0.05);
+ hMCVertexcTS1D->Draw("hist");
+ hMCVertexcTS1D->Draw("err same");
+ cmcVertexcTS->Update();
+
+ cmc1->cd(20);
+ TH1D *hMCRefitcTS1D = data->ShowProjection(icTS, stepGenRefit);
+ maxcTS = hMCRefitcTS1D->GetMaximum();
+ hMCRefitcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hMCRefitcTS1D,42,3,42,20,42,3006);
+ hMCRefitcTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitcTS1D->Sumw2();
+ hMCRefitcTS1D->SetMinimum(0.01);
+ hMCRefitcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Refit cuts");
+ hMCRefitcTS1D->Draw("hist");
+ cmcRefitcTS->cd();
+ cmcRefitcTS->SetLeftMargin(0.15);
+ cmcRefitcTS->SetRightMargin(0.05);
+ hMCRefitcTS1D->Draw("hist");
+ hMCRefitcTS1D->Draw("err same");
+ cmcRefitcTS->Update();
+
+ cmc1->cd(21);
+ TH1D *hRECcTS1D = data->ShowProjection(icTS, stepRec);
+ maxcTS = hRECcTS1D->GetMaximum();
+ hRECcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hRECcTS1D,4,3,4,20,4,3004);
+ hRECcTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECcTS1D->Sumw2();
+ hRECcTS1D->SetMinimum(0.01);
+ hRECcTS1D->Draw("hist");
+ ccTS->cd();
+ ccTS->SetLeftMargin(0.15);
+ ccTS->SetRightMargin(0.05);
+ hRECcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD");
+ hRECcTS1D->Draw("hist");
+ hRECcTS1D->Draw("err same");
+ ccTS->Update();
+
+ cmc1->cd(22);
+ TH1D *hRECAcccTS1D = data->ShowProjection(icTS, stepRecAcc);
+ maxcTS = hRECAcccTS1D->GetMaximum();
+ hRECAcccTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hRECAcccTS1D,38,3,38,20,38,3007);
+ hRECAcccTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAcccTS1D->Sumw2();
+ hRECAcccTS1D->SetMinimum(0.01);
+ hRECAcccTS1D->Draw("hist");
+ cAcccTS->cd();
+ cAcccTS->SetLeftMargin(0.15);
+ cAcccTS->SetRightMargin(0.05);
+ hRECAcccTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - Acc cuts");
+ hRECAcccTS1D->Draw("hist");
+ hRECAcccTS1D->Draw("err same");
+ cAcccTS->Update();
+
+ cmc1->cd(23);
+ TH1D *hRECITSClcTS1D = data->ShowProjection(icTS, stepRecITSClus);
+ maxcTS = hRECITSClcTS1D->GetMaximum();
+ hRECITSClcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hRECITSClcTS1D,13,3,13,20,13,3019);
+ hRECITSClcTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSClcTS1D->Sumw2();
+ hRECITSClcTS1D->SetMinimum(0.01);
+ hRECITSClcTS1D->Draw("hist");
+ cITSClcTS->cd();
+ cITSClcTS->SetLeftMargin(0.15);
+ cITSClcTS->SetRightMargin(0.05);
+ hRECITSClcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - ITS Cl cuts");
+ hRECITSClcTS1D->Draw("hist");
+ hRECITSClcTS1D->Draw("err same");
+ cITSClcTS->Update();
+
+ cmc1->cd(24);
+ TH1D *hRECPPRcTS1D = data->ShowProjection(icTS, stepRecPPR);
+ maxcTS = hRECPPRcTS1D->GetMaximum();
+ hRECPPRcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
+ SetHistoProperties(hRECPPRcTS1D,6,3,6,20,6,3019);
+ hRECPPRcTS1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRcTS1D->Sumw2();
+ hRECPPRcTS1D->SetMinimum(0.01);
+ hRECPPRcTS1D->Draw("hist");
+ cPPRcTS->cd();
+ cPPRcTS->SetLeftMargin(0.15);
+ cPPRcTS->SetRightMargin(0.05);
+ hRECPPRcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - PPR cuts");
+ hRECPPRcTS1D->Draw("hist");
+ hRECPPRcTS1D->Draw("err same");
+ cPPRcTS->Update();
+
+ // ptPi, ptK, cT
+ cmc2->Divide(8,3);
+
+ cmc2->cd(1);
+ TH1D *hMCptPi1D = data->ShowProjection(iptPi, stepGen);
+ Double_t maxptPi = hMCptPi1D->GetMaximum();
+ hMCptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ hMCptPi1D->Sumw2();
+ SetHistoProperties(hMCptPi1D,2,3,2,20,2,3005);
+ hMCptPi1D->SetMinimum(0.01);
+ hMCptPi1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data");
+ hMCptPi1D->Draw("hist");
+ cmcptPi->cd();
+ cmcptPi->SetLeftMargin(0.15);
+ cmcptPi->SetRightMargin(0.05);
+ hMCptPi1D->Draw("hist");
+ hMCptPi1D->Draw("err same");
+ cmcptPi->Update();
+
+ cmc2->cd(2);
+ TH1D *hMCAccptPi1D = data->ShowProjection(iptPi, stepGenAcc);
+ maxptPi = hMCAccptPi1D->GetMaximum();
+ hMCAccptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ hMCAccptPi1D->Sumw2();
+ hMCAccptPi1D->SetMinimum(0.01);
+ SetHistoProperties(hMCAccptPi1D,8,3,8,20,8,3006);
+ hMCAccptPi1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCAccptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Acc cuts");
+ hMCAccptPi1D->Draw("hist");
+ cmcAccptPi->cd();
+ cmcAccptPi->SetLeftMargin(0.15);
+ cmcAccptPi->SetRightMargin(0.05);
+ hMCAccptPi1D->Draw("hist");
+ hMCAccptPi1D->Draw("err same");
+ cmcAccptPi->Update();
+
+ cmc2->cd(3);
+ TH1D *hMCVertexptPi1D = data->ShowProjection(iptPi, stepGenVertex);
+ maxptPi = hMCVertexptPi1D->GetMaximum();
+ hMCVertexptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ hMCVertexptPi1D->Sumw2();
+ hMCVertexptPi1D->SetMinimum(0.01);
+ SetHistoProperties(hMCVertexptPi1D,49,3,49,20,49,3006);
+ hMCVertexptPi1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCVertexptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Vertex cuts");
+ hMCVertexptPi1D->Draw("hist");
+ cmcVertexptPi->cd();
+ cmcVertexptPi->SetLeftMargin(0.15);
+ cmcVertexptPi->SetRightMargin(0.05);
+ hMCVertexptPi1D->Draw("hist");
+ hMCVertexptPi1D->Draw("err same");
+ cmcVertexptPi->Update();
+
+ cmc2->cd(4);
+ TH1D *hMCRefitptPi1D = data->ShowProjection(iptPi, stepGenRefit);
+ maxptPi = hMCRefitptPi1D->GetMaximum();
+ hMCRefitptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ hMCRefitptPi1D->Sumw2();
+ hMCRefitptPi1D->SetMinimum(0.01);
+ SetHistoProperties(hMCRefitptPi1D,42,3,42,20,42,3006);
+ hMCRefitptPi1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCRefitptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Refit cuts");
+ hMCRefitptPi1D->Draw("hist");
+ cmcRefitptPi->cd();
+ cmcRefitptPi->SetLeftMargin(0.15);
+ cmcRefitptPi->SetRightMargin(0.05);
+ hMCRefitptPi1D->Draw("hist");
+ hMCRefitptPi1D->Draw("err same");
+ cmcRefitptPi->Update();
+
+ cmc2->cd(5);
+ TH1D *hRECptPi1D = data->ShowProjection(iptPi, stepRec);
+ maxptPi = hRECptPi1D->GetMaximum();
+ hRECptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ SetHistoProperties(hRECptPi1D,4,3,4,20,4,3004);
+ hRECptPi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECptPi1D->Sumw2();
+ hRECptPi1D->SetMinimum(0.01);
+ hRECptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD");
+ hRECptPi1D->Draw("hist");
+ cptPi->cd();
+ cptPi->SetLeftMargin(0.15);
+ cptPi->SetRightMargin(0.05);
+ hRECptPi1D->Draw("hist");
+ hRECptPi1D->Draw("err same");
+ cptPi->Update();
+
+ cmc2->cd(6);
+ TH1D *hRECAccptPi1D = data->ShowProjection(iptPi, stepRecAcc);
+ maxptPi = hRECAccptPi1D->GetMaximum();
+ hRECAccptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ SetHistoProperties(hRECAccptPi1D,38,3,38,20,38,3007);
+ hRECAccptPi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccptPi1D->Sumw2();
+ hRECAccptPi1D->SetMinimum(0.01);
+ hRECAccptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - Acc cuts");
+ hRECAccptPi1D->Draw("hist");
+ cAccptPi->cd();
+ cAccptPi->SetLeftMargin(0.15);
+ cAccptPi->SetRightMargin(0.05);
+ hRECAccptPi1D->Draw("hist");
+ hRECAccptPi1D->Draw("err same");
+ cAccptPi->Update();
+
+ cmc2->cd(7);
+ TH1D *hRECITSClptPi1D = data->ShowProjection(iptPi, stepRecITSClus);
+ maxptPi = hRECITSClptPi1D->GetMaximum();
+ hRECITSClptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ SetHistoProperties(hRECITSClptPi1D,13,3,13,20,13,3019);
+ hRECITSClptPi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSClptPi1D->Sumw2();
+ hRECITSClptPi1D->SetMinimum(0.01);
+ hRECITSClptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - ITS Cl cuts");
+ hRECITSClptPi1D->Draw("hist");
+ cITSClptPi->cd();
+ cITSClptPi->SetLeftMargin(0.15);
+ cITSClptPi->SetRightMargin(0.05);
+ hRECITSClptPi1D->Draw("hist");
+ hRECITSClptPi1D->Draw("err same");
+ cITSClptPi->Update();
+
+ cmc2->cd(8);
+ TH1D *hRECPPRptPi1D = data->ShowProjection(iptPi, stepRecPPR);
+ maxptPi = hRECPPRptPi1D->GetMaximum();
+ hRECPPRptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
+ SetHistoProperties(hRECPPRptPi1D,6,3,6,20,6,3019);
+ hRECPPRptPi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRptPi1D->Sumw2();
+ hRECPPRptPi1D->SetMinimum(0.01);
+ hRECPPRptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - PPR cuts");
+ hRECPPRptPi1D->Draw("hist");
+ cPPRptPi->cd();
+ cPPRptPi->SetLeftMargin(0.15);
+ cPPRptPi->SetRightMargin(0.05);
+ hRECPPRptPi1D->Draw("hist");
+ hRECPPRptPi1D->Draw("err same");
+ cPPRptPi->Update();
+
+ cmc2->cd(9);
+ TH1D *hMCptK1D = data->ShowProjection(iptK, stepGen);
+ Double_t maxptK = hMCptK1D->GetMaximum();
+ hMCptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hMCptK1D,2,3,2,20,2,3005);
+ hMCptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCptK1D->Sumw2();
+ hMCptK1D->SetMinimum(0.01);
+ hMCptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data");
+ hMCptK1D->Draw("hist");
+ cmcptK->cd();
+ cmcptK->SetLeftMargin(0.15);
+ cmcptK->SetRightMargin(0.05);
+ hMCptK1D->Draw("hist");
+ hMCptK1D->Draw("err same");
+ cmcptK->Update();
+
+ cmc2->cd(10);
+ TH1D *hMCAccptK1D = data->ShowProjection(iptK, stepGenAcc);
+ maxptK = hMCAccptK1D->GetMaximum();
+ hMCAccptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hMCAccptK1D,8,3,8,20,8,3006);
+ hMCAccptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccptK1D->Sumw2();
+ hMCAccptK1D->SetMinimum(0.01);
+ hMCAccptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Acc cuts");
+ hMCAccptK1D->Draw("hist");
+ cmcAccptK->cd();
+ cmcAccptK->SetLeftMargin(0.15);
+ cmcAccptK->SetRightMargin(0.05);
+ hMCAccptK1D->Draw("hist");
+ hMCAccptK1D->Draw("err same");
+ cmcAccptK->Update();
+
+ cmc2->cd(11);
+ TH1D *hMCVertexptK1D = data->ShowProjection(iptK, stepGenVertex);
+ maxptK = hMCVertexptK1D->GetMaximum();
+ hMCVertexptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hMCVertexptK1D,49,3,49,20,49,3006);
+ hMCVertexptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexptK1D->Sumw2();
+ hMCVertexptK1D->SetMinimum(0.01);
+ hMCVertexptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Vertex cuts");
+ hMCVertexptK1D->Draw("hist");
+ cmcVertexptK->cd();
+ cmcVertexptK->SetLeftMargin(0.15);
+ cmcVertexptK->SetRightMargin(0.05);
+ hMCVertexptK1D->Draw("hist");
+ hMCVertexptK1D->Draw("err same");
+ cmcVertexptK->Update();
+
+ cmc2->cd(12);
+ TH1D *hMCRefitptK1D = data->ShowProjection(iptK, stepGenRefit);
+ maxptK = hMCRefitptK1D->GetMaximum();
+ hMCRefitptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hMCRefitptK1D,42,3,42,20,42,3006);
+ hMCRefitptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitptK1D->Sumw2();
+ hMCRefitptK1D->SetMinimum(0.01);
+ hMCRefitptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Refit cuts");
+ hMCRefitptK1D->Draw("hist");
+ cmcRefitptK->cd();
+ cmcRefitptK->SetLeftMargin(0.15);
+ cmcRefitptK->SetRightMargin(0.05);
+ hMCRefitptK1D->Draw("hist");
+ hMCRefitptK1D->Draw("err same");
+ cmcRefitptK->Update();
+
+ cmc2->cd(13);
+ TH1D *hRECptK1D = data->ShowProjection(iptK, stepRec);
+ maxptK = hRECptK1D->GetMaximum();
+ hRECptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hRECptK1D,4,3,4,20,4,3004);
+ hRECptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECptK1D->Sumw2();
+ hRECptK1D->SetMinimum(0.01);
+ hRECptK1D->Draw("hist");
+ cptK->cd();
+ cptK->SetLeftMargin(0.15);
+ cptK->SetRightMargin(0.05);
+ hRECptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD");
+ hRECptK1D->Draw("hist");
+ hRECptK1D->Draw("err same");
+ cptK->Update();
+
+ cmc2->cd(14);
+ TH1D *hRECAccptK1D = data->ShowProjection(iptK, stepRecAcc);
+ maxptK = hRECAccptK1D->GetMaximum();
+ hRECAccptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hRECAccptK1D,38,3,38,20,38,3007);
+ hRECAccptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccptK1D->Sumw2();
+ hRECAccptK1D->SetMinimum(0.01);
+ hRECAccptK1D->Draw("hist");
+ cAccptK->cd();
+ cAccptK->SetLeftMargin(0.15);
+ cAccptK->SetRightMargin(0.05);
+ hRECAccptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - Acc cuts");
+ hRECAccptK1D->Draw("hist");
+ hRECAccptK1D->Draw("err same");
+ cAccptK->Update();
+
+ cmc2->cd(15);
+ TH1D *hRECITSClptK1D = data->ShowProjection(iptK, stepRecITSClus);
+ maxptK = hRECITSClptK1D->GetMaximum();
+ hRECITSClptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hRECITSClptK1D,13,3,13,20,13,3019);
+ hRECITSClptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSClptK1D->Sumw2();
+ hRECITSClptK1D->SetMinimum(0.01);
+ hRECITSClptK1D->Draw("hist");
+ cITSClptK->cd();
+ cITSClptK->SetLeftMargin(0.15);
+ cITSClptK->SetRightMargin(0.05);
+ hRECITSClptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - ITS Cl cuts");
+ hRECITSClptK1D->Draw("hist");
+ hRECITSClptK1D->Draw("err same");
+ cITSClptK->Update();
+
+ cmc2->cd(16);
+ TH1D *hRECPPRptK1D = data->ShowProjection(iptK, stepRecPPR);
+ maxptK = hRECPPRptK1D->GetMaximum();
+ hRECPPRptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
+ SetHistoProperties(hRECPPRptK1D,6,3,6,20,6,3019);
+ hRECPPRptK1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRptK1D->Sumw2();
+ hRECPPRptK1D->SetMinimum(0.01);
+ hRECPPRptK1D->Draw("hist");
+ cPPRptK->cd();
+ cPPRptK->SetLeftMargin(0.15);
+ cPPRptK->SetRightMargin(0.05);
+ hRECPPRptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - PPR cuts");
+ hRECPPRptK1D->Draw("hist");
+ hRECPPRptK1D->Draw("err same");
+ cPPRptK->Update();
+
+ cmc2->cd(17);
+ TH1D *hMCcT1D = data->ShowProjection(icT, stepGen);
+ Double_t maxcT = hMCcT1D->GetMaximum();
+ hMCcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hMCcT1D,2,3,2,20,2,3005);
+ hMCcT1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCcT1D->Sumw2();
+ hMCcT1D->SetMinimum(0.01);
+ hMCcT1D->GetXaxis()->SetTitle("ct (#mum), MC data");
+ hMCcT1D->Draw("hist");
+ cmccT->cd();
+ cmccT->SetLeftMargin(0.15);
+ cmccT->SetRightMargin(0.05);
+ hMCcT1D->Draw("hist");
+ hMCcT1D->Draw("err same");
+ cmccT->Update();
+
+ cmc2->cd(18);
+ TH1D *hMCAcccT1D = data->ShowProjection(icT, stepGenAcc);
+ maxcT = hMCAcccT1D->GetMaximum();
+ hMCAcccT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hMCAcccT1D,8,3,8,20,8,3006);
+ hMCAcccT1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCAcccT1D->Sumw2();
+ hMCAcccT1D->SetMinimum(0.01);
+ hMCAcccT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Acc cuts");
+ hMCAcccT1D->Draw("hist");
+ cmcAcccT->cd();
+ cmcAcccT->SetLeftMargin(0.15);
+ cmcAcccT->SetRightMargin(0.05);
+ hMCAcccT1D->Draw("hist");
+ hMCAcccT1D->Draw("err same");
+ cmcAcccT->Update();
+
+ cmc2->cd(19);
+ TH1D *hMCVertexcT1D = data->ShowProjection(icT, stepGenVertex);
+ maxcT = hMCVertexcT1D->GetMaximum();
+ hMCVertexcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hMCVertexcT1D,49,3,49,20,49,3006);
+ hMCVertexcT1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCVertexcT1D->Sumw2();
+ hMCVertexcT1D->SetMinimum(0.01);
+ hMCVertexcT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Vertex cuts");
+ hMCVertexcT1D->Draw("hist");
+ cmcVertexcT->cd();
+ cmcVertexcT->SetLeftMargin(0.15);
+ cmcVertexcT->SetRightMargin(0.05);
+ hMCVertexcT1D->Draw("hist");
+ hMCVertexcT1D->Draw("err same");
+ cmcVertexcT->Update();
+
+ cmc2->cd(20);
+ TH1D *hMCRefitcT1D = data->ShowProjection(icT, stepGenRefit);
+ maxcT = hMCRefitcT1D->GetMaximum();
+ hMCRefitcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hMCRefitcT1D,42,3,42,20,42,3006);
+ hMCRefitcT1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCRefitcT1D->Sumw2();
+ hMCRefitcT1D->SetMinimum(0.01);
+ hMCRefitcT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Refit cuts");
+ hMCRefitcT1D->Draw("hist");
+ cmcRefitcT->cd();
+ cmcRefitcT->SetLeftMargin(0.15);
+ cmcRefitcT->SetRightMargin(0.05);
+ hMCRefitcT1D->Draw("hist");
+ hMCRefitcT1D->Draw("err same");
+ cmcRefitcT->Update();
+
+ cmc2->cd(21);
+ TH1D *hRECcT1D = data->ShowProjection(icT, stepRec);
+ maxcT = hRECcT1D->GetMaximum();
+ hRECcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hRECcT1D,4,3,4,20,4,3004);
+ hRECcT1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECcT1D->Sumw2();
+ hRECcT1D->SetMinimum(0.01);
+ hRECcT1D->Draw("hist");
+ ccT->cd();
+ ccT->SetLeftMargin(0.15);
+ ccT->SetRightMargin(0.05);
+ hRECcT1D->GetXaxis()->SetTitle("ct (#mum), AOD");
+ hRECcT1D->Draw("hist");
+ hRECcT1D->Draw("err same");
+ ccT->Update();
+
+ cmc2->cd(22);
+ TH1D *hRECAcccT1D = data->ShowProjection(icT, stepRecAcc);
+ maxcT = hRECAcccT1D->GetMaximum();
+ hRECAcccT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hRECAcccT1D,38,3,38,20,38,3007);
+ hRECAcccT1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAcccT1D->Sumw2();
+ hRECAcccT1D->SetMinimum(0.01);
+ hRECAcccT1D->Draw("hist");
+ cAcccT->cd();
+ cAcccT->SetLeftMargin(0.15);
+ cAcccT->SetRightMargin(0.05);
+ hRECAcccT1D->GetXaxis()->SetTitle("ct (#mum), AOD - Acc cuts");
+ hRECAcccT1D->Draw("hist");
+ hRECAcccT1D->Draw("err same");
+ cAcccT->Update();
+
+ cmc2->cd(23);
+ TH1D *hRECITSClcT1D = data->ShowProjection(icT, stepRecITSClus);
+ maxcT = hRECITSClcT1D->GetMaximum();
+ hRECITSClcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hRECITSClcT1D,13,3,13,20,13,3019);
+ hRECITSClcT1D->GetXaxis()->SetTitleOffset(1.2);
+ hRECITSClcT1D->Sumw2();
+ hRECITSClcT1D->SetMinimum(0.01);
+ hRECITSClcT1D->Draw("hist");
+ cITSClcT->cd();
+ cITSClcT->SetLeftMargin(0.15);
+ cITSClcT->SetRightMargin(0.05);
+ hRECITSClcT1D->GetXaxis()->SetTitle("ct (#mum), AOD - ITS Cl cuts");
+ hRECITSClcT1D->Draw("hist");
+ hRECITSClcT1D->Draw("err same");
+ cITSClcT->Update();
+
+ cmc2->cd(24);
+ TH1D *hRECPPRcT1D = data->ShowProjection(icT, stepRecPPR);
+ maxcT = hRECPPRcT1D->GetMaximum();
+ hRECPPRcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
+ SetHistoProperties(hRECPPRcT1D,6,3,6,20,6,3019);
+ hRECPPRcT1D->GetXaxis()->SetTitleOffset(1.2);
+ hRECPPRcT1D->Sumw2();
+ hRECPPRcT1D->SetMinimum(0.01);
+ hRECPPRcT1D->Draw("hist");
+ cPPRcT->cd();
+ cPPRcT->SetLeftMargin(0.15);
+ cPPRcT->SetRightMargin(0.05);
+ hRECPPRcT1D->GetXaxis()->SetTitle("ct (#mum), AOD - PPR cuts");
+ hRECPPRcT1D->Draw("hist");
+ hRECPPRcT1D->Draw("err same");
+ cPPRcT->Update();
+
+ //dca, d0pi, d0K
+ cmc3->Divide(8,3);
+ cmc3->cd(1);
+ TH1D *hMCdca1D = data->ShowProjection(idca, stepGen);
+ Double_t maxdca = hMCdca1D->GetMaximum();
+ hMCdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ hMCdca1D->Sumw2();
+ hMCdca1D->SetMinimum(0.01);
+ SetHistoProperties(hMCdca1D,2,3,2,20,2,3005);
+ hMCdca1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCdca1D->GetXaxis()->SetTitle("dca (#mum), MC data");
+ hMCdca1D->Draw("hist");
+ cmcdca->cd();
+ cmcdca->SetLeftMargin(0.15);
+ cmcdca->SetRightMargin(0.05);
+ hMCdca1D->Draw("hist");
+ hMCdca1D->Draw("err same");
+ cmcdca->Update();
+
+ cmc3->cd(2);
+ TH1D *hMCAccdca1D = data->ShowProjection(idca, stepGenAcc);
+ maxdca = hMCAccdca1D->GetMaximum();
+ hMCAccdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ SetHistoProperties(hMCAccdca1D,8,3,8,20,8,3006);
+ hMCAccdca1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccdca1D->Sumw2();
+ hMCAccdca1D->SetMinimum(0.01);
+ hMCAccdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Acc cut");
+ hMCAccdca1D->Draw("hist");
+ cmcAccdca->cd();
+ cmcAccdca->SetLeftMargin(0.15);
+ cmcAccdca->SetRightMargin(0.05);
+ hMCAccdca1D->Draw("hist");
+ hMCAccdca1D->Draw("err same");
+ cmcAccdca->Update();
+
+ cmc3->cd(3);
+ TH1D *hMCVertexdca1D = data->ShowProjection(idca, stepGenVertex);
+ maxdca = hMCVertexdca1D->GetMaximum();
+ hMCVertexdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ SetHistoProperties(hMCVertexdca1D,49,3,49,20,49,3006);
+ hMCVertexdca1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexdca1D->Sumw2();
+ hMCVertexdca1D->SetMinimum(0.01);
+ hMCVertexdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Vertex cut");
+ hMCVertexdca1D->Draw("hist");
+ cmcVertexdca->cd();
+ cmcVertexdca->SetLeftMargin(0.15);
+ cmcVertexdca->SetRightMargin(0.05);
+ hMCVertexdca1D->Draw("hist");
+ hMCVertexdca1D->Draw("err same");
+ cmcVertexdca->Update();
+
+ cmc3->cd(4);
+ TH1D *hMCRefitdca1D = data->ShowProjection(idca, stepGenRefit);
+ maxdca = hMCRefitdca1D->GetMaximum();
+ hMCRefitdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ SetHistoProperties(hMCRefitdca1D,42,3,42,20,42,3006);
+ hMCRefitdca1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitdca1D->Sumw2();
+ hMCRefitdca1D->SetMinimum(0.01);
+ hMCRefitdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Refit cut");
+ hMCRefitdca1D->Draw("hist");
+ cmcRefitdca->cd();
+ cmcRefitdca->SetLeftMargin(0.15);
+ cmcRefitdca->SetRightMargin(0.05);
+ hMCRefitdca1D->Draw("hist");
+ hMCRefitdca1D->Draw("err same");
+ cmcRefitdca->Update();
+
+ cmc3->cd(5);
+ TH1D *hRECdca1D = data->ShowProjection(idca, stepRec);
+ maxdca = hRECdca1D->GetMaximum();
+ hRECdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ SetHistoProperties(hRECdca1D,4,3,4,20,4,3004);
+ hRECdca1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECdca1D->Sumw2();
+ hRECdca1D->SetMinimum(0.01);
+ hRECdca1D->GetXaxis()->SetTitle("dca (#mum), AOD");
+ hRECdca1D->Draw("hist");
+ cdca->cd();
+ cdca->SetLeftMargin(0.15);
+ cdca->SetRightMargin(0.05);
+ hRECdca1D->Draw("hist");
+ hRECdca1D->Draw("err same");
+ cdca->Update();
+
+ cmc3->cd(6);
+ TH1D *hRECAccdca1D = data->ShowProjection(idca, stepRecAcc);
+ maxdca = hRECAccdca1D->GetMaximum();
+ hRECAccdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ SetHistoProperties(hRECAccdca1D,38,3,38,20,38,3007);
+ hRECAccdca1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccdca1D->Sumw2();
+ hRECAccdca1D->SetMinimum(0.01);
+ hRECAccdca1D->GetXaxis()->SetTitle("dca (#mum), AOD - Acc cuts");
+ hRECAccdca1D->Draw("hist");
+ cAccdca->cd();
+ cAccdca->SetLeftMargin(0.15);
+ cAccdca->SetRightMargin(0.05);
+ hRECAccdca1D->Draw("hist");
+ hRECAccdca1D->Draw("err same");
+ cAccdca->Update();
+
+ cmc3->cd(6);
+ TH1D *hRECITSCldca1D = data->ShowProjection(idca, stepRecITSClus);
+ maxdca = hRECITSCldca1D->GetMaximum();
+ hRECITSCldca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ SetHistoProperties(hRECITSCldca1D,13,3,13,20,13,3019);
+ hRECITSCldca1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSCldca1D->Sumw2();
+ hRECITSCldca1D->SetMinimum(0.01);
+ hRECITSCldca1D->GetXaxis()->SetTitle("dca (#mum), AOD - ITS Cl cuts");
+ hRECITSCldca1D->Draw("hist");
+ cITSCldca->cd();
+ cITSCldca->SetLeftMargin(0.15);
+ cITSCldca->SetRightMargin(0.05);
+ hRECITSCldca1D->Draw("hist");
+ hRECITSCldca1D->Draw("err same");
+ cITSCldca->Update();
+
+ cmc3->cd(7);
+ TH1D *hRECPPRdca1D = data->ShowProjection(idca, stepRecPPR);
+ maxdca = hRECPPRdca1D->GetMaximum();
+ hRECPPRdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
+ SetHistoProperties(hRECPPRdca1D,6,3,6,20,6,3019);
+ hRECPPRdca1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRdca1D->Sumw2();
+ hRECPPRdca1D->SetMinimum(0.01);
+ hRECPPRdca1D->GetXaxis()->SetTitle("dca (#mum), AOD - PPR cuts");
+ hRECPPRdca1D->Draw("hist");
+ cPPRdca->cd();
+ cPPRdca->SetLeftMargin(0.15);
+ cPPRdca->SetRightMargin(0.05);
+ hRECPPRdca1D->Draw("hist");
+ hRECPPRdca1D->Draw("err same");
+ cPPRdca->Update();
+
+ cmc3->cd(8);
+ TH1D *hMCd0pi1D = data->ShowProjection(id0pi, stepGen);
+ Double_t maxd0pi = hMCd0pi1D->GetMaximum();
+ hMCd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ hMCd0pi1D->Sumw2();
+ hMCd0pi1D->SetMinimum(0.01);
+ SetHistoProperties(hMCd0pi1D,2,3,2,20,2,3005);
+ hMCd0pi1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC data");
+ hMCd0pi1D->Draw("hist");
+ cmcd0pi->cd();
+ cmcd0pi->SetLeftMargin(0.15);
+ cmcd0pi->SetRightMargin(0.05);
+ hMCd0pi1D->Draw("hist");
+ hMCd0pi1D->Draw("err same");
+ cmcd0pi->Update();
+
+ cmc3->cd(9);
+ TH1D *hMCAccd0pi1D = data->ShowProjection(id0pi, stepGenAcc);
+ maxd0pi = hMCAccd0pi1D->GetMaximum();
+ hMCAccd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ SetHistoProperties(hMCAccd0pi1D,8,3,8,20,8,3006);
+ hMCAccd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccd0pi1D->Sumw2();
+ hMCAccd0pi1D->SetMinimum(0.01);
+ hMCAccd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Acc cut");
+ hMCAccd0pi1D->Draw("hist");
+ cmcAccd0pi->cd();
+ cmcAccd0pi->SetLeftMargin(0.15);
+ cmcAccd0pi->SetRightMargin(0.05);
+ hMCAccd0pi1D->Draw("hist");
+ hMCAccd0pi1D->Draw("err same");
+ cmcAccd0pi->Update();
+
+ cmc3->cd(10);
+ TH1D *hMCVertexd0pi1D = data->ShowProjection(id0pi, stepGenVertex);
+ maxd0pi = hMCVertexd0pi1D->GetMaximum();
+ hMCVertexd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ SetHistoProperties(hMCVertexd0pi1D,49,3,49,20,49,3006);
+ hMCVertexd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexd0pi1D->Sumw2();
+ hMCVertexd0pi1D->SetMinimum(0.01);
+ hMCVertexd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Vertex cut");
+ hMCVertexd0pi1D->Draw("hist");
+ cmcVertexd0pi->cd();
+ cmcVertexd0pi->SetLeftMargin(0.15);
+ cmcVertexd0pi->SetRightMargin(0.05);
+ hMCVertexd0pi1D->Draw("hist");
+ hMCVertexd0pi1D->Draw("err same");
+ cmcVertexd0pi->Update();
+
+ cmc3->cd(11);
+ TH1D *hMCRefitd0pi1D = data->ShowProjection(id0pi, stepGenRefit);
+ maxd0pi = hMCRefitd0pi1D->GetMaximum();
+ hMCRefitd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ SetHistoProperties(hMCRefitd0pi1D,42,3,42,20,42,3006);
+ hMCRefitd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitd0pi1D->Sumw2();
+ hMCRefitd0pi1D->SetMinimum(0.01);
+ hMCRefitd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Refit cut");
+ hMCRefitd0pi1D->Draw("hist");
+ cmcRefitd0pi->cd();
+ cmcRefitd0pi->SetLeftMargin(0.15);
+ cmcRefitd0pi->SetRightMargin(0.05);
+ hMCRefitd0pi1D->Draw("hist");
+ hMCRefitd0pi1D->Draw("err same");
+ cmcRefitd0pi->Update();
+
+ cmc3->cd(12);
+ TH1D *hRECd0pi1D = data->ShowProjection(id0pi, stepRec);
+ maxd0pi = hRECd0pi1D->GetMaximum();
+ hRECd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ SetHistoProperties(hRECd0pi1D,4,3,4,20,4,3004);
+ hRECd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECd0pi1D->Sumw2();
+ hRECd0pi1D->SetMinimum(0.01);
+ hRECd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD");
+ hRECd0pi1D->Draw("hist");
+ cd0pi->cd();
+ cd0pi->SetLeftMargin(0.15);
+ cd0pi->SetRightMargin(0.05);
+ hRECd0pi1D->Draw("hist");
+ hRECd0pi1D->Draw("err same");
+ cd0pi->Update();
+
+ cmc3->cd(13);
+ TH1D *hRECAccd0pi1D = data->ShowProjection(id0pi, stepRecAcc);
+ maxd0pi = hRECAccd0pi1D->GetMaximum();
+ hRECAccd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ SetHistoProperties(hRECAccd0pi1D,38,3,38,20,38,3007);
+ hRECAccd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccd0pi1D->Sumw2();
+ hRECAccd0pi1D->SetMinimum(0.01);
+ hRECAccd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - Acc cuts");
+ hRECAccd0pi1D->Draw("hist");
+ cAccd0pi->cd();
+ cAccd0pi->SetLeftMargin(0.15);
+ cAccd0pi->SetRightMargin(0.05);
+ hRECAccd0pi1D->Draw("hist");
+ hRECAccd0pi1D->Draw("err same");
+ cAccd0pi->Update();
+
+ cmc3->cd(14);
+ TH1D *hRECITSCld0pi1D = data->ShowProjection(id0pi, stepRecITSClus);
+ maxd0pi = hRECITSCld0pi1D->GetMaximum();
+ hRECITSCld0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ SetHistoProperties(hRECITSCld0pi1D,13,3,13,20,13,3019);
+ hRECITSCld0pi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSCld0pi1D->Sumw2();
+ hRECITSCld0pi1D->SetMinimum(0.01);
+ hRECITSCld0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - ITS Cl cuts");
+ hRECITSCld0pi1D->Draw("hist");
+ cITSCld0pi->cd();
+ cITSCld0pi->SetLeftMargin(0.15);
+ cITSCld0pi->SetRightMargin(0.05);
+ hRECITSCld0pi1D->Draw("hist");
+ hRECITSCld0pi1D->Draw("err same");
+ cITSCld0pi->Update();
+
+ cmc3->cd(15);
+ TH1D *hRECPPRd0pi1D = data->ShowProjection(id0pi, stepRecPPR);
+ maxd0pi = hRECPPRd0pi1D->GetMaximum();
+ hRECPPRd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
+ SetHistoProperties(hRECPPRd0pi1D,6,3,6,20,6,3019);
+ hRECPPRd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRd0pi1D->Sumw2();
+ hRECPPRd0pi1D->SetMinimum(0.01);
+ hRECPPRd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - PPR cuts");
+ hRECPPRd0pi1D->Draw("hist");
+ cPPRd0pi->cd();
+ cPPRd0pi->SetLeftMargin(0.15);
+ cPPRd0pi->SetRightMargin(0.05);
+ hRECPPRd0pi1D->Draw("hist");
+ hRECPPRd0pi1D->Draw("err same");
+ cPPRd0pi->Update();
+
+ cmc3->cd(16);
+ TH1D *hMCd0K1D = data->ShowProjection(id0K, stepGen);
+ Double_t maxd0K = hMCd0K1D->GetMaximum();
+ hMCd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ hMCd0K1D->Sumw2();
+ hMCd0K1D->SetMinimum(0.01);
+ SetHistoProperties(hMCd0K1D,2,3,2,20,2,3005);
+ hMCd0K1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC data");
+ hMCd0K1D->Draw("hist");
+ cmcd0K->cd();
+ cmcd0K->SetLeftMargin(0.15);
+ cmcd0K->SetRightMargin(0.05);
+ hMCd0K1D->Draw("hist");
+ hMCd0K1D->Draw("err same");
+ cmcd0K->Update();
+
+ cmc3->cd(17);
+ TH1D *hMCAccd0K1D = data->ShowProjection(id0K, stepGenAcc);
+ maxd0K = hMCAccd0K1D->GetMaximum();
+ hMCAccd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ SetHistoProperties(hMCAccd0K1D,8,3,8,20,8,3006);
+ hMCAccd0K1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccd0K1D->Sumw2();
+ hMCAccd0K1D->SetMinimum(0.01);
+ hMCAccd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Acc cut");
+ hMCAccd0K1D->Draw("hist");
+ cmcAccd0K->cd();
+ cmcAccd0K->SetLeftMargin(0.15);
+ cmcAccd0K->SetRightMargin(0.05);
+ hMCAccd0K1D->Draw("hist");
+ hMCAccd0K1D->Draw("err same");
+ cmcAccd0K->Update();
+
+ cmc3->cd(18);
+ TH1D *hMCVertexd0K1D = data->ShowProjection(id0K, stepGenVertex);
+ maxd0K = hMCVertexd0K1D->GetMaximum();
+ hMCVertexd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ SetHistoProperties(hMCVertexd0K1D,49,3,49,20,49,3006);
+ hMCVertexd0K1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexd0K1D->Sumw2();
+ hMCVertexd0K1D->SetMinimum(0.01);
+ hMCVertexd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Vertex cut");
+ hMCVertexd0K1D->Draw("hist");
+ cmcVertexd0K->cd();
+ cmcVertexd0K->SetLeftMargin(0.15);
+ cmcVertexd0K->SetRightMargin(0.05);
+ hMCVertexd0K1D->Draw("hist");
+ hMCVertexd0K1D->Draw("err same");
+ cmcVertexd0K->Update();
+
+ cmc3->cd(20);
+ TH1D *hMCRefitd0K1D = data->ShowProjection(id0K, stepGenRefit);
+ maxd0K = hMCRefitd0K1D->GetMaximum();
+ hMCRefitd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ SetHistoProperties(hMCRefitd0K1D,42,3,42,20,42,3006);
+ hMCRefitd0K1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitd0K1D->Sumw2();
+ hMCRefitd0K1D->SetMinimum(0.01);
+ hMCRefitd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Refit cut");
+ hMCRefitd0K1D->Draw("hist");
+ cmcRefitd0K->cd();
+ cmcRefitd0K->SetLeftMargin(0.15);
+ cmcRefitd0K->SetRightMargin(0.05);
+ hMCRefitd0K1D->Draw("hist");
+ hMCRefitd0K1D->Draw("err same");
+ cmcRefitd0K->Update();
+
+ cmc3->cd(21);
+ TH1D *hRECd0K1D = data->ShowProjection(id0K, stepRec);
+ maxd0K = hRECd0K1D->GetMaximum();
+ hRECd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ SetHistoProperties(hRECd0K1D,4,3,4,20,4,3004);
+ hRECd0K1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECd0K1D->Sumw2();
+ hRECd0K1D->SetMinimum(0.01);
+ hRECd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD");
+ hRECd0K1D->Draw("hist");
+ cd0K->cd();
+ cd0K->SetLeftMargin(0.15);
+ cd0K->SetRightMargin(0.05);
+ hRECd0K1D->Draw("hist");
+ hRECd0K1D->Draw("err same");
+ cd0K->Update();
+
+ cmc3->cd(22);
+ TH1D *hRECAccd0K1D = data->ShowProjection(id0K, stepRecAcc);
+ maxd0K = hRECAccd0K1D->GetMaximum();
+ hRECAccd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ SetHistoProperties(hRECAccd0K1D,38,3,38,20,38,3007);
+ hRECAccd0K1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccd0K1D->Sumw2();
+ hRECAccd0K1D->SetMinimum(0.01);
+ hRECAccd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - Acc cuts");
+ hRECAccd0K1D->Draw("hist");
+ cAccd0K->cd();
+ cAccd0K->SetLeftMargin(0.15);
+ cAccd0K->SetRightMargin(0.05);
+ hRECAccd0K1D->Draw("hist");
+ hRECAccd0K1D->Draw("err same");
+ cAccd0K->Update();
+
+ cmc3->cd(23);
+ TH1D *hRECITSCld0K1D = data->ShowProjection(id0K, stepRecITSClus);
+ maxd0K = hRECITSCld0K1D->GetMaximum();
+ hRECITSCld0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ SetHistoProperties(hRECITSCld0K1D,13,3,13,20,13,3019);
+ hRECITSCld0K1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSCld0K1D->Sumw2();
+ hRECITSCld0K1D->SetMinimum(0.01);
+ hRECITSCld0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - ITS Cl cuts");
+ hRECITSCld0K1D->Draw("hist");
+ cITSCld0K->cd();
+ cITSCld0K->SetLeftMargin(0.15);
+ cITSCld0K->SetRightMargin(0.05);
+ hRECITSCld0K1D->Draw("hist");
+ hRECITSCld0K1D->Draw("err same");
+ cITSCld0K->Update();
+
+ cmc3->cd(24);
+ TH1D *hRECPPRd0K1D = data->ShowProjection(id0K, stepRecPPR);
+ maxd0K = hRECPPRd0K1D->GetMaximum();
+ hRECPPRd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
+ SetHistoProperties(hRECPPRd0K1D,6,3,6,20,6,3019);
+ hRECPPRd0K1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRd0K1D->Sumw2();
+ hRECPPRd0K1D->SetMinimum(0.01);
+ hRECPPRd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - PPR cuts");
+ hRECPPRd0K1D->Draw("hist");
+ cPPRd0K->cd();
+ cPPRd0K->SetLeftMargin(0.15);
+ cPPRd0K->SetRightMargin(0.05);
+ hRECPPRd0K1D->Draw("hist");
+ hRECPPRd0K1D->Draw("err same");
+ cPPRd0K->Update();
+
+ //d0xd0, cosPointingAngle, Phi
+ cmc4->Divide(8,3);
+ cmc4->cd(1);
+ TH1D *hMCd0xd01D = data->ShowProjection(id0xd0, stepGen);
+ Double_t maxd0xd0 = hMCd0xd01D->GetMaximum();
+ hMCd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ hMCd0xd01D->Sumw2();
+ hMCd0xd01D->SetMinimum(0.01);
+ SetHistoProperties(hMCd0xd01D,2,3,2,20,2,3005);
+ hMCd0xd01D->GetXaxis()->SetTitleOffset(1.2);
+ hMCd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC data");
+ hMCd0xd01D->Draw("hist");
+ cmcd0xd0->cd();
+ cmcd0xd0->SetLeftMargin(0.15);
+ cmcd0xd0->SetRightMargin(0.05);
+ hMCd0xd01D->Draw("hist");
+ hMCd0xd01D->Draw("err same");
+ cmcd0xd0->Update();
+
+ cmc4->cd(2);
+ TH1D *hMCAccd0xd01D = data->ShowProjection(id0xd0, stepGenAcc);
+ maxd0xd0 = hMCAccd0xd01D->GetMaximum();
+ hMCAccd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ SetHistoProperties(hMCAccd0xd01D,8,3,8,20,8,3006);
+ hMCAccd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccd0xd01D->Sumw2();
+ hMCAccd0xd01D->SetMinimum(0.01);
+ hMCAccd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Acc cut");
+ hMCAccd0xd01D->Draw("hist");
+ cmcAccd0xd0->cd();
+ cmcAccd0xd0->SetLeftMargin(0.15);
+ cmcAccd0xd0->SetRightMargin(0.05);
+ hMCAccd0xd01D->Draw("hist");
+ hMCAccd0xd01D->Draw("err same");
+ cmcAccd0xd0->Update();
+
+ cmc4->cd(3);
+ TH1D *hMCVertexd0xd01D = data->ShowProjection(id0xd0, stepGenVertex);
+ maxd0xd0 = hMCVertexd0xd01D->GetMaximum();
+ hMCVertexd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ SetHistoProperties(hMCVertexd0xd01D,49,3,49,20,49,3006);
+ hMCVertexd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexd0xd01D->Sumw2();
+ hMCVertexd0xd01D->SetMinimum(0.01);
+ hMCVertexd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Vertex cut");
+ hMCVertexd0xd01D->Draw("hist");
+ cmcVertexd0xd0->cd();
+ cmcVertexd0xd0->SetLeftMargin(0.15);
+ cmcVertexd0xd0->SetRightMargin(0.05);
+ hMCVertexd0xd01D->Draw("hist");
+ hMCVertexd0xd01D->Draw("err same");
+ cmcVertexd0xd0->Update();
+
+ cmc4->cd(4);
+ TH1D *hMCRefitd0xd01D = data->ShowProjection(id0xd0, stepGenRefit);
+ maxd0xd0 = hMCRefitd0xd01D->GetMaximum();
+ hMCRefitd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ SetHistoProperties(hMCRefitd0xd01D,42,3,42,20,42,3006);
+ hMCRefitd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitd0xd01D->Sumw2();
+ hMCRefitd0xd01D->SetMinimum(0.01);
+ hMCRefitd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Refit cut");
+ hMCRefitd0xd01D->Draw("hist");
+ cmcRefitd0xd0->cd();
+ cmcRefitd0xd0->SetLeftMargin(0.15);
+ cmcRefitd0xd0->SetRightMargin(0.05);
+ hMCRefitd0xd01D->Draw("hist");
+ hMCRefitd0xd01D->Draw("err same");
+ cmcRefitd0xd0->Update();
+
+ cmc4->cd(5);
+ TH1D *hRECd0xd01D = data->ShowProjection(id0xd0, stepRec);
+ maxd0xd0 = hRECd0xd01D->GetMaximum();
+ hRECd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ SetHistoProperties(hRECd0xd01D,4,3,4,20,4,3004);
+ hRECd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECd0xd01D->Sumw2();
+ hRECd0xd01D->SetMinimum(0.01);
+ hRECd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD");
+ hRECd0xd01D->Draw("hist");
+ cd0xd0->cd();
+ cd0xd0->SetLeftMargin(0.15);
+ cd0xd0->SetRightMargin(0.05);
+ hRECd0xd01D->Draw("hist");
+ hRECd0xd01D->Draw("err same");
+ cd0xd0->Update();
+
+ cmc4->cd(6);
+ TH1D *hRECAccd0xd01D = data->ShowProjection(id0xd0, stepRecAcc);
+ maxd0xd0 = hRECAccd0xd01D->GetMaximum();
+ hRECAccd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ SetHistoProperties(hRECAccd0xd01D,38,3,38,20,38,3007);
+ hRECAccd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccd0xd01D->Sumw2();
+ hRECAccd0xd01D->SetMinimum(0.01);
+ hRECAccd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - Acc cuts");
+ hRECAccd0xd01D->Draw("hist");
+ cAccd0xd0->cd();
+ cAccd0xd0->SetLeftMargin(0.15);
+ cAccd0xd0->SetRightMargin(0.05);
+ hRECAccd0xd01D->Draw("hist");
+ hRECAccd0xd01D->Draw("err same");
+ cAccd0xd0->Update();
+
+ cmc4->cd(7);
+ TH1D *hRECITSCld0xd01D = data->ShowProjection(id0xd0, stepRecITSClus);
+ maxd0xd0 = hRECITSCld0xd01D->GetMaximum();
+ hRECITSCld0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ SetHistoProperties(hRECITSCld0xd01D,13,3,13,20,13,3019);
+ hRECITSCld0xd01D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSCld0xd01D->Sumw2();
+ hRECITSCld0xd01D->SetMinimum(0.01);
+ hRECITSCld0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - ITS Cl cuts");
+ hRECITSCld0xd01D->Draw("hist");
+ cITSCld0xd0->cd();
+ cITSCld0xd0->SetLeftMargin(0.15);
+ cITSCld0xd0->SetRightMargin(0.05);
+ hRECITSCld0xd01D->Draw("hist");
+ hRECITSCld0xd01D->Draw("err same");
+ cITSCld0xd0->Update();
+
+ cmc4->cd(8);
+ TH1D *hRECPPRd0xd01D = data->ShowProjection(id0xd0, stepRecPPR);
+ maxd0xd0 = hRECPPRd0xd01D->GetMaximum();
+ hRECPPRd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
+ SetHistoProperties(hRECPPRd0xd01D,6,3,6,20,6,3019);
+ hRECPPRd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRd0xd01D->Sumw2();
+ hRECPPRd0xd01D->SetMinimum(0.01);
+ hRECPPRd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - PPR cuts");
+ hRECPPRd0xd01D->Draw("hist");
+ cPPRd0xd0->cd();
+ cPPRd0xd0->SetLeftMargin(0.15);
+ cPPRd0xd0->SetRightMargin(0.05);
+ hRECPPRd0xd01D->Draw("hist");
+ hRECPPRd0xd01D->Draw("err same");
+ cPPRd0xd0->Update();
+
+ cmc4->cd(9);
+ TH1D *hMCpointing1D = data->ShowProjection(ipointing, stepGen);
+ Double_t maxpointing = hMCpointing1D->GetMaximum();
+ hMCpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ hMCpointing1D->Sumw2();
+ hMCpointing1D->SetMinimum(0.01);
+ SetHistoProperties(hMCpointing1D,2,3,2,20,2,3005);
+ hMCpointing1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC data");
+ hMCpointing1D->Draw("hist");
+ cmcpointing->cd();
+ cmcpointing->SetLeftMargin(0.15);
+ cmcpointing->SetRightMargin(0.05);
+ hMCpointing1D->Draw("hist");
+ hMCpointing1D->Draw("err same");
+ cmcpointing->Update();
+
+ cmc4->cd(10);
+ TH1D *hMCAccpointing1D = data->ShowProjection(ipointing, stepGenAcc);
+ maxpointing = hMCAccpointing1D->GetMaximum();
+ hMCAccpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ SetHistoProperties(hMCAccpointing1D,8,3,8,20,8,3006);
+ hMCAccpointing1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccpointing1D->Sumw2();
+ hMCAccpointing1D->SetMinimum(0.01);
+ hMCAccpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Acc cut");
+ hMCAccpointing1D->Draw("hist");
+ cmcAccpointing->cd();
+ cmcAccpointing->SetLeftMargin(0.15);
+ cmcAccpointing->SetRightMargin(0.05);
+ hMCAccpointing1D->Draw("hist");
+ hMCAccpointing1D->Draw("err same");
+ cmcAccpointing->Update();
+
+ cmc4->cd(11);
+ TH1D *hMCVertexpointing1D = data->ShowProjection(ipointing, stepGenVertex);
+ maxpointing = hMCVertexpointing1D->GetMaximum();
+ hMCVertexpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ SetHistoProperties(hMCVertexpointing1D,49,3,49,20,49,3006);
+ hMCVertexpointing1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexpointing1D->Sumw2();
+ hMCVertexpointing1D->SetMinimum(0.01);
+ hMCVertexpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Vertex cut");
+ hMCVertexpointing1D->Draw("hist");
+ cmcVertexpointing->cd();
+ cmcVertexpointing->SetLeftMargin(0.15);
+ cmcVertexpointing->SetRightMargin(0.05);
+ hMCVertexpointing1D->Draw("hist");
+ hMCVertexpointing1D->Draw("err same");
+ cmcVertexpointing->Update();
+
+ cmc4->cd(12);
+ TH1D *hMCRefitpointing1D = data->ShowProjection(ipointing, stepGenRefit);
+ maxpointing = hMCRefitpointing1D->GetMaximum();
+ hMCRefitpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ SetHistoProperties(hMCRefitpointing1D,42,3,42,20,42,3006);
+ hMCRefitpointing1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitpointing1D->Sumw2();
+ hMCRefitpointing1D->SetMinimum(0.01);
+ hMCRefitpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Refit cut");
+ hMCRefitpointing1D->Draw("hist");
+ cmcRefitpointing->cd();
+ cmcRefitpointing->SetLeftMargin(0.15);
+ cmcRefitpointing->SetRightMargin(0.05);
+ hMCRefitpointing1D->Draw("hist");
+ hMCRefitpointing1D->Draw("err same");
+ cmcRefitpointing->Update();
+
+ cmc4->cd(13);
+ TH1D *hRECpointing1D = data->ShowProjection(ipointing, stepRec);
+ maxpointing = hRECpointing1D->GetMaximum();
+ hRECpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ SetHistoProperties(hRECpointing1D,4,3,4,20,4,3004);
+ hRECpointing1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECpointing1D->Sumw2();
+ hRECpointing1D->SetMinimum(0.01);
+ hRECpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD");
+ hRECpointing1D->Draw("hist");
+ cpointing->cd();
+ cpointing->SetLeftMargin(0.15);
+ cpointing->SetRightMargin(0.05);
+ hRECpointing1D->Draw("hist");
+ hRECpointing1D->Draw("err same");
+ cpointing->Update();
+
+ cmc4->cd(14);
+ TH1D *hRECAccpointing1D = data->ShowProjection(ipointing, stepRecAcc);
+ maxpointing = hRECAccpointing1D->GetMaximum();
+ hRECAccpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ SetHistoProperties(hRECAccpointing1D,38,3,38,20,38,3007);
+ hRECAccpointing1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccpointing1D->Sumw2();
+ hRECAccpointing1D->SetMinimum(0.01);
+ hRECAccpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - Acc cuts");
+ hRECAccpointing1D->Draw("hist");
+ cAccpointing->cd();
+ cAccpointing->SetLeftMargin(0.15);
+ cAccpointing->SetRightMargin(0.05);
+ hRECAccpointing1D->Draw("hist");
+ hRECAccpointing1D->Draw("err same");
+ cAccpointing->Update();
+
+ cmc4->cd(15);
+ TH1D *hRECITSClpointing1D = data->ShowProjection(ipointing, stepRecITSClus);
+ maxpointing = hRECITSClpointing1D->GetMaximum();
+ hRECITSClpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ SetHistoProperties(hRECITSClpointing1D,13,3,13,20,13,3019);
+ hRECITSClpointing1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSClpointing1D->Sumw2();
+ hRECITSClpointing1D->SetMinimum(0.01);
+ hRECITSClpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - ITS Cl cuts");
+ hRECITSClpointing1D->Draw("hist");
+ cITSClpointing->cd();
+ cITSClpointing->SetLeftMargin(0.15);
+ cITSClpointing->SetRightMargin(0.05);
+ hRECITSClpointing1D->Draw("hist");
+ hRECITSClpointing1D->Draw("err same");
+ cITSClpointing->Update();
+
+ cmc4->cd(16);
+ TH1D *hRECPPRpointing1D = data->ShowProjection(ipointing, stepRecPPR);
+ maxpointing = hRECPPRpointing1D->GetMaximum();
+ hRECPPRpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
+ SetHistoProperties(hRECPPRpointing1D,6,3,6,20,6,3019);
+ hRECPPRpointing1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRpointing1D->Sumw2();
+ hRECPPRpointing1D->SetMinimum(0.01);
+ hRECPPRpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - PPR cuts");
+ hRECPPRpointing1D->Draw("hist");
+ cPPRpointing->cd();
+ cPPRpointing->SetLeftMargin(0.15);
+ cPPRpointing->SetRightMargin(0.05);
+ hRECPPRpointing1D->Draw("hist");
+ hRECPPRpointing1D->Draw("err same");
+ cPPRpointing->Update();
+
+ cmc4->cd(17);
+ TH1D *hMCphi1D = data->ShowProjection(iphi, stepGen);
+ Double_t maxphi = hMCphi1D->GetMaximum();
+ hMCphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ hMCphi1D->Sumw2();
+ hMCphi1D->SetMinimum(0.01);
+ SetHistoProperties(hMCphi1D,2,3,2,20,2,3005);
+ hMCphi1D->GetXaxis()->SetTitleOffset(1.2);
+ hMCphi1D->GetXaxis()->SetTitle("Phi, MC data");
+ hMCphi1D->Draw("hist");
+ cmcphi->cd();
+ cmcphi->SetLeftMargin(0.15);
+ cmcphi->SetRightMargin(0.05);
+ hMCphi1D->Draw("hist");
+ hMCphi1D->Draw("err same");
+ cmcphi->Update();
+
+ cmc4->cd(18);
+ TH1D *hMCAccphi1D = data->ShowProjection(iphi, stepGenAcc);
+ maxphi = hMCAccphi1D->GetMaximum();
+ hMCAccphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ SetHistoProperties(hMCAccphi1D,8,3,8,20,8,3006);
+ hMCAccphi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCAccphi1D->Sumw2();
+ hMCAccphi1D->SetMinimum(0.01);
+ hMCAccphi1D->GetXaxis()->SetTitle("Phi, MC - Acc cut");
+ hMCAccphi1D->Draw("hist");
+ cmcAccphi->cd();
+ cmcAccphi->SetLeftMargin(0.15);
+ cmcAccphi->SetRightMargin(0.05);
+ hMCAccphi1D->Draw("hist");
+ hMCAccphi1D->Draw("err same");
+ cmcAccphi->Update();
+
+ cmc4->cd(19);
+ TH1D *hMCVertexphi1D = data->ShowProjection(iphi, stepGenVertex);
+ maxphi = hMCVertexphi1D->GetMaximum();
+ hMCVertexphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ SetHistoProperties(hMCVertexphi1D,49,3,49,20,49,3006);
+ hMCVertexphi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCVertexphi1D->Sumw2();
+ hMCVertexphi1D->SetMinimum(0.01);
+ hMCVertexphi1D->GetXaxis()->SetTitle("Phi, MC - Vertex cut");
+ hMCVertexphi1D->Draw("hist");
+ cmcVertexphi->cd();
+ cmcVertexphi->SetLeftMargin(0.15);
+ cmcVertexphi->SetRightMargin(0.05);
+ hMCVertexphi1D->Draw("hist");
+ hMCVertexphi1D->Draw("err same");
+ cmcVertexphi->Update();
+
+ cmc4->cd(20);
+ TH1D *hMCRefitphi1D = data->ShowProjection(iphi, stepGenRefit);
+ maxphi = hMCRefitphi1D->GetMaximum();
+ hMCRefitphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ SetHistoProperties(hMCRefitphi1D,42,3,42,20,42,3006);
+ hMCRefitphi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hMCRefitphi1D->Sumw2();
+ hMCRefitphi1D->SetMinimum(0.01);
+ hMCRefitphi1D->GetXaxis()->SetTitle("Phi, MC - Refit cut");
+ hMCRefitphi1D->Draw("hist");
+ cmcRefitphi->cd();
+ cmcRefitphi->SetLeftMargin(0.15);
+ cmcRefitphi->SetRightMargin(0.05);
+ hMCRefitphi1D->Draw("hist");
+ hMCRefitphi1D->Draw("err same");
+ cmcRefitphi->Update();
+
+ cmc4->cd(21);
+ TH1D *hRECphi1D = data->ShowProjection(iphi, stepRec);
+ maxphi = hRECphi1D->GetMaximum();
+ hRECphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ SetHistoProperties(hRECphi1D,4,3,4,20,4,3004);
+ hRECphi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECphi1D->Sumw2();
+ hRECphi1D->SetMinimum(0.01);
+ hRECphi1D->GetXaxis()->SetTitle("Phi, AOD");
+ hRECphi1D->Draw("hist");
+ cphi->cd();
+ cphi->SetLeftMargin(0.15);
+ cphi->SetRightMargin(0.05);
+ hRECphi1D->Draw("hist");
+ hRECphi1D->Draw("err same");
+ cphi->Update();
+
+ cmc4->cd(22);
+ TH1D *hRECAccphi1D = data->ShowProjection(iphi, stepRecAcc);
+ maxphi = hRECAccphi1D->GetMaximum();
+ hRECAccphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ SetHistoProperties(hRECAccphi1D,38,3,38,20,38,3007);
+ hRECAccphi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECAccphi1D->Sumw2();
+ hRECAccphi1D->SetMinimum(0.01);
+ hRECAccphi1D->GetXaxis()->SetTitle("Phi, AOD - Acc cuts");
+ hRECAccphi1D->Draw("hist");
+ cAccphi->cd();
+ cAccphi->SetLeftMargin(0.15);
+ cAccphi->SetRightMargin(0.05);
+ hRECAccphi1D->Draw("hist");
+ hRECAccphi1D->Draw("err same");
+ cAccphi->Update();
+
+ cmc4->cd(23);
+ TH1D *hRECITSClphi1D = data->ShowProjection(iphi, stepRecITSClus);
+ maxphi = hRECITSClphi1D->GetMaximum();
+ hRECITSClphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ SetHistoProperties(hRECITSClphi1D,13,3,13,20,13,3019);
+ hRECITSClphi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECITSClphi1D->Sumw2();
+ hRECITSClphi1D->SetMinimum(0.01);
+ hRECITSClphi1D->GetXaxis()->SetTitle("Phi, AOD - ITS Cl cuts");
+ hRECITSClphi1D->Draw("hist");
+ cITSClphi->cd();
+ cITSClphi->SetLeftMargin(0.15);
+ cITSClphi->SetRightMargin(0.05);
+ hRECITSClphi1D->Draw("hist");
+ hRECITSClphi1D->Draw("err same");
+ cITSClphi->Update();
+
+ cmc4->cd(24);
+ TH1D *hRECPPRphi1D = data->ShowProjection(iphi, stepRecPPR);
+ maxphi = hRECPPRphi1D->GetMaximum();
+ hRECPPRphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
+ SetHistoProperties(hRECPPRphi1D,6,3,6,20,6,3019);
+ hRECPPRphi1D->GetXaxis()->SetTitleOffset(1.2);;
+ hRECPPRphi1D->Sumw2();
+ hRECPPRphi1D->SetMinimum(0.01);
+ hRECPPRphi1D->GetXaxis()->SetTitle("Phi, AOD - PPR cuts");
+ hRECPPRphi1D->Draw("hist");
+ cPPRphi->cd();
+ cPPRphi->SetLeftMargin(0.15);
+ cPPRphi->SetRightMargin(0.05);
+ hRECPPRphi1D->Draw("hist");
+ hRECPPRphi1D->Draw("err same");
+ cPPRphi->Update();
+
+ // printing on eps files
+ cmc1->Print("Plots/dataMC_pt_y_cTS.eps");
+ cmc2->Print("Plots/dataMC_ptPi_ptK_cT.eps");
+ cmc3->Print("Plots/dataMC_dca_d0pi_d0K.eps");
+ cmc4->Print("Plots/dataMC_d0xd0_pointing.eps");
+ cmcpt->Print("Plots/pt_Gen.eps");
+ cmcy->Print("Plots/y_Gen.eps");
+ cmccTS->Print("Plots/cTS_Gen.eps");
+ cmcptPi->Print("Plots/ptPi_Gen.eps");
+ cmcptK->Print("Plots/ptK_Gen.eps");
+ cmccT->Print("Plots/cT_Gen.eps");
+ cmcdca->Print("Plots/dca_Gen.eps");
+ cmcd0pi->Print("Plots/d0pi_Gen.eps");
+ cmcd0K->Print("Plots/d0K_Gen.eps");
+ cmcd0xd0->Print("Plots/d0xd0_Gen.eps");
+ cmcpointing->Print("Plots/pointing_Gen.eps");
+ cmcphi->Print("Plots/phi_Gen.eps");
+
+ cmcAccpt->Print("Plots/pt_Gen_Acc.eps");
+ cmcAccy->Print("Plots/y_Gen_Acc.eps");
+ cmcAcccTS->Print("Plots/cTS_Gen_Acc.eps");
+ cmcAccptPi->Print("Plots/ptPi_Gen_Acc.eps");
+ cmcAccptK->Print("Plots/ptK_Gen_Acc.eps");
+ cmcAcccT->Print("Plots/cT_Gen_Acc.eps");
+ cmcAccdca->Print("Plots/dca_Gen_Acc.eps");
+ cmcAccd0pi->Print("Plots/d0pi_Gen_Acc.eps");
+ cmcAccd0K->Print("Plots/d0K_Gen_Acc.eps");
+ cmcAccd0xd0->Print("Plots/d0xd0_Gen_Acc.eps");
+ cmcAccpointing->Print("Plots/pointing_Gen_Acc.eps");
+ cmcAccphi->Print("Plots/phi_Gen_Acc.eps");
+
+ cmcVertexpt->Print("Plots/pt_Gen_Vertex.eps");
+ cmcVertexy->Print("Plots/y_Gen_Vertex.eps");
+ cmcVertexcTS->Print("Plots/cTS_Gen_Vertex.eps");
+ cmcVertexptPi->Print("Plots/ptPi_Gen_Vertex.eps");
+ cmcVertexptK->Print("Plots/ptK_Gen_Vertex.eps");
+ cmcVertexcT->Print("Plots/cT_Gen_Vertex.eps");
+ cmcVertexdca->Print("Plots/dca_Gen_Vertex.eps");
+ cmcVertexd0pi->Print("Plots/d0pi_Gen_Vertex.eps");
+ cmcVertexd0K->Print("Plots/d0K_Gen_Vertex.eps");
+ cmcVertexd0xd0->Print("Plots/d0xd0_Gen_Vertex.eps");
+ cmcVertexpointing->Print("Plots/pointing_Gen_Vertex.eps");
+ cmcVertexphi->Print("Plots/phi_Gen_Vertex.eps");
+
+ cmcRefitpt->Print("Plots/pt_Gen_Refit.eps");
+ cmcRefity->Print("Plots/y_Gen_Refit.eps");
+ cmcRefitcTS->Print("Plots/cTS_Gen_Refit.eps");
+ cmcRefitptPi->Print("Plots/ptPi_Gen_Refit.eps");
+ cmcRefitptK->Print("Plots/ptK_Gen_Refit.eps");
+ cmcRefitcT->Print("Plots/cT_Gen_Refit.eps");
+ cmcRefitdca->Print("Plots/dca_Gen_Refit.eps");
+ cmcRefitd0pi->Print("Plots/d0pi_Gen_Refit.eps");
+ cmcRefitd0K->Print("Plots/d0K_Gen_Refit.eps");
+ cmcRefitd0xd0->Print("Plots/d0xd0_Gen_Refit.eps");
+ cmcRefitpointing->Print("Plots/pointing_Gen_Refit.eps");
+ cmcRefitphi->Print("Plots/phi_Gen_Refit.eps");
+
+ cpt->Print("Plots/pt_Rec.eps");
+ cy->Print("Plots/y_Rec.eps");
+ ccTS->Print("Plots/cTS_Rec.eps");
+ cptPi->Print("Plots/ptPi_Rec.eps");
+ cptK->Print("Plots/ptK_Rec.eps");
+ ccT->Print("Plots/cT_Rec.eps");
+ cdca->Print("Plots/dca_Rec.eps");
+ cd0pi->Print("Plots/d0pi_Rec.eps");
+ cd0K->Print("Plots/d0K_Rec.eps");
+ cd0xd0->Print("Plots/d0xd0_Rec.eps");
+ cpointing->Print("Plots/pointing_Rec.eps");
+ cphi->Print("Plots/phi_Rec.eps");
+
+ cAccpt->Print("Plots/pt_Rec_Acc.eps");
+ cAccy->Print("Plots/y_Rec_Acc.eps");
+ cAcccTS->Print("Plots/cTS_Rec_Acc.eps");
+ cAccptPi->Print("Plots/ptPi_Rec_Acc.eps");
+ cAccptK->Print("Plots/ptK_Rec_Acc.eps");
+ cAcccT->Print("Plots/cT_Rec_Acc.eps");
+ cmcdca->Print("Plots/dca_Rec_Acc.eps");
+ cAccd0pi->Print("Plots/d0pi_Rec_Acc.eps");
+ cAccd0K->Print("Plots/d0K_Rec_Acc.eps");
+ cAccd0xd0->Print("Plots/d0xd0_Rec_Acc.eps");
+ cAccpointing->Print("Plots/pointing_Rec_Acc.eps");
+ cAccphi->Print("Plots/phi_Rec_Acc.eps");
+
+ cITSClpt->Print("Plots/pt_Rec_ITSCl.eps");
+ cITSCly->Print("Plots/y_Rec_ITSCl.eps");
+ cITSClcTS->Print("Plots/cTS_Rec_ITSCl.eps");
+ cITSClptPi->Print("Plots/ptPi_Rec_ITSCl.eps");
+ cITSClptK->Print("Plots/ptK_Rec_ITSCl.eps");
+ cITSClcT->Print("Plots/cT_Rec_ITSCl.eps");
+ cITSCldca->Print("Plots/dca_Rec_ITSCl.eps");
+ cITSCld0pi->Print("Plots/d0pi_Rec_ITSCl.eps");
+ cITSCld0K->Print("Plots/d0K_Rec_ITSCl.eps");
+ cITSCld0xd0->Print("Plots/d0xd0_Rec_ITSCl.eps");
+ cITSClpointing->Print("Plots/pointing_Rec_ITSCl.eps");
+ cITSClphi->Print("Plots/phi_Rec_ITSCl.eps");
+
+ cPPRpt->Print("Plots/pt_Rec_PPR.eps");
+ cPPRy->Print("Plots/y_Rec_PPR.eps");
+ cPPRcTS->Print("Plots/cTS_Rec_PPR.eps");
+ cPPRptPi->Print("Plots/ptPi_Rec_PPR.eps");
+ cPPRptK->Print("Plots/ptK_Rec_PPR.eps");
+ cPPRcT->Print("Plots/cT_Rec_PPR.eps");
+ cPPRdca->Print("Plots/dca_Rec_PPR.eps");
+ cPPRd0pi->Print("Plots/d0pi_Rec_PPR.eps");
+ cPPRd0K->Print("Plots/d0K_Rec_PPR.eps");
+ cPPRd0xd0->Print("Plots/d0xd0_Rec_PPR.eps");
+ cPPRpointing->Print("Plots/pointing_Rec_PPR.eps");
+ cPPRphi->Print("Plots/phi_Rec_PPR.eps");
+
+ // printing on gif files
+ cmc1->Print("Plots/dataMC_pt_y_cTS.gif");
+ cmc2->Print("Plots/dataMC_ptPi_ptK_cT.gif");
+ cmc3->Print("Plots/dataMC_dca_d0pi_d0K.gif");
+ cmc4->Print("Plots/dataMC_d0xd0_pointing.gif");
+ cmcpt->Print("Plots/pt_Gen.gif");
+ cmcy->Print("Plots/y_Gen.gif");
+ cmccTS->Print("Plots/cTS_Gen.gif");
+ cmcptPi->Print("Plots/ptPi_Gen.gif");
+ cmcptK->Print("Plots/ptK_Gen.gif");
+ cmccT->Print("Plots/cT_Gen.gif");
+ cmcdca->Print("Plots/dca_Gen.gif");
+ cmcd0pi->Print("Plots/d0pi_Gen.gif");
+ cmcd0K->Print("Plots/d0K_Gen.gif");
+ cmcd0xd0->Print("Plots/d0xd0_Gen.gif");
+ cmcpointing->Print("Plots/pointing_Gen.gif");
+ cmcphi->Print("Plots/phi_Gen.gif");
+
+ cmcAccpt->Print("Plots/pt_Gen_Acc.gif");
+ cmcAccy->Print("Plots/y_Gen_Acc.gif");
+ cmcAcccTS->Print("Plots/cTS_Gen_Acc.gif");
+ cmcAccptPi->Print("Plots/ptPi_Gen_Acc.gif");
+ cmcAccptK->Print("Plots/ptK_Gen_Acc.gif");
+ cmcAcccT->Print("Plots/cT_Gen_Acc.gif");
+ cmcAccdca->Print("Plots/dca_Gen_Acc.gif");
+ cmcAccd0pi->Print("Plots/d0pi_Gen_Acc.gif");
+ cmcAccd0K->Print("Plots/d0K_Gen_Acc.gif");
+ cmcAccd0xd0->Print("Plots/d0xd0_Gen_Acc.gif");
+ cmcAccpointing->Print("Plots/pointing_Gen_Acc.gif");
+ cmcAccphi->Print("Plots/phi_Gen_Acc.gif");
+
+ cmcVertexpt->Print("Plots/pt_Gen_Vertex.gif");
+ cmcVertexy->Print("Plots/y_Gen_Vertex.gif");
+ cmcVertexcTS->Print("Plots/cTS_Gen_Vertex.gif");
+ cmcVertexptPi->Print("Plots/ptPi_Gen_Vertex.gif");
+ cmcVertexptK->Print("Plots/ptK_Gen_Vertex.gif");
+ cmcVertexcT->Print("Plots/cT_Gen_Vertex.gif");
+ cmcVertexdca->Print("Plots/dca_Gen_Vertex.gif");
+ cmcVertexd0pi->Print("Plots/d0pi_Gen_Vertex.gif");
+ cmcVertexd0K->Print("Plots/d0K_Gen_Vertex.gif");
+ cmcVertexd0xd0->Print("Plots/d0xd0_Gen_Vertex.gif");
+ cmcVertexpointing->Print("Plots/pointing_Gen_Vertex.gif");
+ cmcVertexphi->Print("Plots/phi_Gen_Vertex.gif");
+
+ cmcRefitpt->Print("Plots/pt_Gen_Refit.gif");
+ cmcRefity->Print("Plots/y_Gen_Refit.gif");
+ cmcRefitcTS->Print("Plots/cTS_Gen_Refit.gif");
+ cmcRefitptPi->Print("Plots/ptPi_Gen_Refit.gif");
+ cmcRefitptK->Print("Plots/ptK_Gen_Refit.gif");
+ cmcRefitcT->Print("Plots/cT_Gen_Refit.gif");
+ cmcRefitdca->Print("Plots/dca_Gen_Refit.gif");
+ cmcRefitd0pi->Print("Plots/d0pi_Gen_Refit.gif");
+ cmcRefitd0K->Print("Plots/d0K_Gen_Refit.gif");
+ cmcRefitd0xd0->Print("Plots/d0xd0_Gen_Refit.gif");
+ cmcRefitpointing->Print("Plots/pointing_Gen_Refit.gif");
+ cmcRefitphi->Print("Plots/phi_Gen_Refit.gif");
+
+ cpt->Print("Plots/pt_Rec.gif");
+ cy->Print("Plots/y_Rec.gif");
+ ccTS->Print("Plots/cTS_Rec.gif");
+ cptPi->Print("Plots/ptPi_Rec.gif");
+ cptK->Print("Plots/ptK_Rec.gif");
+ ccT->Print("Plots/cT_Rec.gif");
+ cdca->Print("Plots/dca_Rec.gif");
+ cd0pi->Print("Plots/d0pi_Rec.gif");
+ cd0K->Print("Plots/d0K_Rec.gif");
+ cd0xd0->Print("Plots/d0xd0_Rec.gif");
+ cpointing->Print("Plots/pointing_Rec.gif");
+ cphi->Print("Plots/phi_Rec.gif");
+
+ cAccpt->Print("Plots/pt_Rec_Acc.gif");
+ cAccy->Print("Plots/y_Rec_Acc.gif");
+ cAcccTS->Print("Plots/cTS_Rec_Acc.gif");
+ cAccptPi->Print("Plots/ptPi_Rec_Acc.gif");
+ cAccptK->Print("Plots/ptK_Rec_Acc.gif");
+ cAcccT->Print("Plots/cT_Rec_Acc.gif");
+ cmcdca->Print("Plots/dca_Rec_Acc.gif");
+ cAccd0pi->Print("Plots/d0pi_Rec_Acc.gif");
+ cAccd0K->Print("Plots/d0K_Rec_Acc.gif");
+ cAccd0xd0->Print("Plots/d0xd0_Rec_Acc.gif");
+ cAccpointing->Print("Plots/pointing_Rec_Acc.gif");
+ cAccphi->Print("Plots/phi_Rec_Acc.gif");
+
+ cITSClpt->Print("Plots/pt_Rec_ITSCl.gif");
+ cITSCly->Print("Plots/y_Rec_ITSCl.gif");
+ cITSClcTS->Print("Plots/cTS_Rec_ITSCl.gif");
+ cITSClptPi->Print("Plots/ptPi_Rec_ITSCl.gif");
+ cITSClptK->Print("Plots/ptK_Rec_ITSCl.gif");
+ cITSClcT->Print("Plots/cT_Rec_ITSCl.gif");
+ cITSCldca->Print("Plots/dca_Rec_ITSCl.gif");
+ cITSCld0pi->Print("Plots/d0pi_Rec_ITSCl.gif");
+ cITSCld0K->Print("Plots/d0K_Rec_ITSCl.gif");
+ cITSCld0xd0->Print("Plots/d0xd0_Rec_ITSCl.gif");
+ cITSClpointing->Print("Plots/pointing_Rec_ITSCl.gif");
+ cITSClphi->Print("Plots/phi_Rec_ITSCl.gif");
+
+ cPPRpt->Print("Plots/pt_Rec_PPR.gif");
+ cPPRy->Print("Plots/y_Rec_PPR.gif");
+ cPPRcTS->Print("Plots/cTS_Rec_PPR.gif");
+ cPPRptPi->Print("Plots/ptPi_Rec_PPR.gif");
+ cPPRptK->Print("Plots/ptK_Rec_PPR.gif");
+ cPPRcT->Print("Plots/cT_Rec_PPR.gif");
+ cPPRdca->Print("Plots/dca_Rec_PPR.gif");
+ cPPRd0pi->Print("Plots/d0pi_Rec_PPR.gif");
+ cPPRd0K->Print("Plots/d0K_Rec_PPR.gif");
+ cPPRd0xd0->Print("Plots/d0xd0_Rec_PPR.gif");
+ cPPRpointing->Print("Plots/pointing_Rec_PPR.gif");
+ cPPRphi->Print("Plots/phi_Rec_PPR.gif");
+
+ TFile* file_histo = new TFile("fileHisto.root","RECREATE");
+ // MC
+ hMCpt1D->Write("hMCpt1D");
+ hMCy1D->Write("hMCy1D");
+ hMCcTS1D->Write("hMCcTS1D");
+ hMCptPi1D->Write("hMCptPi1D");
+ hMCptK1D->Write("hMCptK1D");
+ hMCcT1D->Write("hMCcT1D");
+ hMCdca1D->Write("hMCdca1D");
+ hMCd0pi1D->Write("hMCd0pi1D");
+ hMCd0K1D->Write("hMCd0K1D");
+ hMCd0xd01D->Write("hMCd0xd01D");
+ hMCpointing1D->Write("hMCpointing1D");
+ hMCphi1D->Write("hMCphi1D");
+ // MC, Acc
+ hMCAccpt1D->Write("hMCAccpt1D");
+ hMCAccy1D->Write("hMCAccy1D");
+ hMCAcccTS1D->Write("hMCAcccTS1D");
+ hMCAccptPi1D->Write("hMCAccptPi1D");
+ hMCAccptK1D->Write("hMCAccptK1D");
+ hMCAcccT1D->Write("hMCAcccT1D");
+ hMCAccdca1D->Write("hMCAccdca1D");
+ hMCAccd0pi1D->Write("hMCAccd0pi1D");
+ hMCAccd0K1D->Write("hMCAccd0K1D");
+ hMCAccd0xd01D->Write("hMCAccd0xd01D");
+ hMCAccpointing1D->Write("hMCAccpointing1D");
+ hMCAccphi1D->Write("hMCAccphi1D");
+ // REC
+ hRECpt1D->Write("hRECpt1D");
+ hRECy1D->Write("hRECy1D");
+ hRECcTS1D->Write("hRECcTS1D");
+ hRECptPi1D->Write("hRECptPi1D");
+ hRECptK1D->Write("hRECptK1D");
+ hRECcT1D->Write("hRECcT1D");
+ hRECdca1D->Write("hRECdca1D");
+ hRECd0pi1D->Write("hRECd0pi1D");
+ hRECd0K1D->Write("hRECd0K1D");
+ hRECd0xd01D->Write("hRECd0xd01D");
+ hRECpointing1D->Write("hRECpointing1D");
+ hRECphi1D->Write("hRECphi1D");
+ // REC, Acc
+ hRECAccpt1D->Write("hRECAccpt1D");
+ hRECAccy1D->Write("hRECAccy1D");
+ hRECAcccTS1D->Write("hRECAcccTS1D");
+ hRECAccptPi1D->Write("hRECAccptPi1D");
+ hRECAccptK1D->Write("hRECAccptK1D");
+ hRECAcccT1D->Write("hRECAcccT1D");
+ hRECAccdca1D->Write("hRECAccdca1D");
+ hRECAccd0pi1D->Write("hRECAccd0pi1D");
+ hRECAccd0K1D->Write("hRECAccd0K1D");
+ hRECAccd0xd01D->Write("hRECAccd0xd01D");
+ hRECAccpointing1D->Write("hRECAccpointing1D");
+ hRECAccphi1D->Write("hRECAccphi1D");
+ // REC, ITS Cl
+ hRECITSClpt1D->Write("hRECITSClpt1D");
+ hRECITSCly1D->Write("hRECITSCly1D");
+ hRECITSClcTS1D->Write("hRECITSClcTS1D");
+ hRECITSClptPi1D->Write("hRECITSClptPi1D");
+ hRECITSClptK1D->Write("hRECITSClptK1D");
+ hRECITSClcT1D->Write("hRECITSClcT1D");
+ hRECITSCldca1D->Write("hRECITSCldca1D");
+ hRECITSCld0pi1D->Write("hRECITSCld0pi1D");
+ hRECITSCld0K1D->Write("hRECITSCld0K1D");
+ hRECITSCld0xd01D->Write("hRECITSCld0xd01D");
+ hRECITSClpointing1D->Write("hRECITSClpointing1D");
+ hRECITSClphi1D->Write("hRECITSClphi1D");
+ // REC, PPR
+ hRECPPRpt1D->Write("hRECPPRpt1D");
+ hRECPPRy1D->Write("hRECPPRy1D");
+ hRECPPRcTS1D->Write("hRECPPRcTS1D");
+ hRECPPRptPi1D->Write("hRECPPRptPi1D");
+ hRECPPRptK1D->Write("hRECPPRptK1D");
+ hRECPPRcT1D->Write("hRECPPRcT1D");
+ hRECPPRdca1D->Write("hRECPPRdca1D");
+ hRECPPRd0pi1D->Write("hRECPPRd0pi1D");
+ hRECPPRd0K1D->Write("hRECPPRd0K1D");
+ hRECPPRd0xd01D->Write("hRECPPRd0xd01D");
+ hRECPPRpointing1D->Write("hRECPPRpointing1D");
+ hRECPPRphi1D->Write("hRECPPRphi1D");
+
+}
+
+void SetHistoProperties(TH1D* h, Int_t lineColor, Int_t lineWidth, Int_t markerColor, Int_t markerStyle, Int_t fillColor, Int_t fillStyle){
+
+ h->SetLineColor(lineColor);
+ h->SetLineWidth(lineWidth);
+ h->SetMarkerColor(markerColor);
+ h->SetFillColor(fillColor);
+ h->SetFillStyle(fillStyle);
+ h->SetMarkerStyle(markerStyle);
+
+ return;
+
+}
+