#include 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; }