PWGHF/vertexingHF/macros/ReadCFHeavyFlavourOutput.C

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