#include <Riostream.h>

extern TRandom *gRandom;
extern TBenchmark *gBenchmark;
extern TSystem *gSystem;

void ReadCFHeavyFlavourOutput() {

	// example macro to read the CF produced for D->Kpi
	// 8 steps considered
	// 12 variables
	
	gROOT->SetStyle("Plain");
	gStyle->SetPalette(1);
	gStyle->SetOptStat(0);
	gStyle->SetPalette(1);
	gStyle->SetCanvasColor(0);
	gStyle->SetFrameFillColor(0);
	gStyle->SetOptTitle(0);
	
	gSystem->SetIncludePath("-I. -I$ALICE_ROOT/include  -I$ROOTSYS/include");
	gSystem->Load("libANALYSIS.so");
	gSystem->Load("libANALYSISalice.so");
	gSystem->Load("$ALICE_ROOT/CORRFW/libCORRFW.so") ;
	
	// Setting up the container grid... 
	
	// Flag the sel steps 
	Int_t stepGen=0;
	Int_t stepGenAcc=1;
	Int_t stepGenVertex=2;
	Int_t stepGenRefit=3;
	Int_t stepRec=4;
	Int_t stepRecAcc=5;
	Int_t stepRecITSClus=6;
	Int_t stepRecPPR=7;
	
	// The sensitive variables, their indeces
	Int_t ipt =0;
	Int_t iy=1;
	UInt_t icTS  = 2;   //cTS stands for cosThetaStar
	UInt_t iptPi  = 3;
	UInt_t iptK  = 4;
	UInt_t icT  = 5;
	UInt_t idca  = 6;
	UInt_t id0pi  = 7;
	UInt_t id0K  = 8;
	UInt_t id0xd0  = 9;
	UInt_t ipointing  = 10;
	UInt_t iphi  = 11;
	
	// Get the container
	TFile* f = new TFile("AnalysisResults.root");	
	TDirectoryFile* d = (TDirectoryFile*)f->Get("PWG3_D2H_CFtaskD0toKpi");
	AliCFContainer *data = (AliCFContainer*) (d->Get("container"));
	
	// Make some 1 & 2-D projections.
	// MC level
	TCanvas *cmc1 =new TCanvas("cmc1","The MC distributions",0,0,1600,1200);  // canvas for all steps, 1st three variables
	TCanvas *cmc2 =new TCanvas("cmc2","The MC distributions",0,0,1600,1200);  // canvas for all steps, 2nd three variables 

	TCanvas *cmc3 =new TCanvas("cmc3","The MC distributions",0,0,1600,1200);  // canvas for all steps, 2nd three variables 
	TCanvas *cmc4 =new TCanvas("cmc4","The MC distributions",0,0,1600,765);  // canvas for all steps, 2nd three variables 
	TCanvas *cmcpt = new TCanvas("cmcpt","pt distribution from MC",50,50,550,550);
	TCanvas *cmcy = new TCanvas("cmcy","y distribution from MC",50,50,550,550);
	TCanvas *cmccTS = new TCanvas("cmcTS","cosThetaStar distribution from MC",50,50,550,550);  // cTS stands for cosThetaStar
	TCanvas *cmcptPi = new TCanvas("cmcptPi","pt_pi distribution from MC",50,50,550,550);
	TCanvas *cmcptK = new TCanvas("cmcptK","pt_K distribution from MC",50,50,550,550);
	TCanvas *cmccT = new TCanvas("cmccT","cT distribution from MC",50,50,550,550);
	TCanvas *cmcdca = new TCanvas("cmcdca","dca distribution from MC",50,50,550,550);
	TCanvas *cmcd0pi = new TCanvas("cmcd0pi","d0pi distribution from MC",50,50,550,550);
	TCanvas *cmcd0K = new TCanvas("cmcd0K","d0K distribution from MC",50,50,550,550);
	TCanvas *cmcd0xd0 = new TCanvas("cmcd0xd0","d0xd0 distribution from MC",50,50,550,550);
	TCanvas *cmcpointing = new TCanvas("cmcpointing","cosPointingAngle distribution from MC",50,50,550,550);
	TCanvas *cmcphi = new TCanvas("cmcphi","Phi distribution from MC",50,50,550,550);

	// MC Acceptance level
	TCanvas *cmcAccpt = new TCanvas("cmcAccpt","pt distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccy = new TCanvas("cmcAccy","y distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAcccTS = new TCanvas("cmcAccTS","cosThetaStar distribution from MC - Acceptance cuts",50,50,550,550);  // cTS stands for cosThetaStar
	TCanvas *cmcAccptPi = new TCanvas("cmcAccptPi","pt_pi distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccptK = new TCanvas("cmcAccptK","pt_K distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAcccT = new TCanvas("cmcAcccT","cT distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccdca = new TCanvas("cmcAccdca","dca distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccd0pi = new TCanvas("cmcAccd0pi","d0pi distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccd0K = new TCanvas("cmcAccd0K","d0K distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccd0xd0 = new TCanvas("cmcAccd0xd0","d0xd0 distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccpointing = new TCanvas("cmcAccpointing","cosPointingAngle distribution from MC - Acceptance cuts",50,50,550,550);
	TCanvas *cmcAccphi = new TCanvas("cmcAccphi","Phi distribution from MC - Acceptance cuts",50,50,550,550);

	// MC Vertex level
	TCanvas *cmcVertexpt = new TCanvas("cmcVertexpt","pt distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexy = new TCanvas("cmcVertexy","y distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexcTS = new TCanvas("cmcVertexTS","cosThetaStar distribution from MC - Vertex cuts",50,50,550,550);  // cTS stands for cosThetaStar
	TCanvas *cmcVertexptPi = new TCanvas("cmcVertexptPi","pt_pi distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexptK = new TCanvas("cmcVertexptK","pt_K distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexcT = new TCanvas("cmcVertexcT","cT distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexdca = new TCanvas("cmcVertexdca","dca distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexd0pi = new TCanvas("cmcVertexd0pi","d0pi distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexd0K = new TCanvas("cmcVertexd0K","d0K distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexd0xd0 = new TCanvas("cmcVertexd0xd0","d0xd0 distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexpointing = new TCanvas("cmcVertexpointing","cosPointingAngle distribution from MC - Vertex cuts",50,50,550,550);
	TCanvas *cmcVertexphi = new TCanvas("cmcVertexphi","Phi distribution from MC - Vertex cuts",50,50,550,550);

	// MC Refit level
	TCanvas *cmcRefitpt = new TCanvas("cmcRefitpt","pt distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefity = new TCanvas("cmcRefity","y distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitcTS = new TCanvas("cmcRefitTS","cosThetaStar distribution from MC - Refit cuts",50,50,550,550);  // cTS stands for cosThetaStar
	TCanvas *cmcRefitptPi = new TCanvas("cmcRefitptPi","pt_pi distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitptK = new TCanvas("cmcRefitptK","pt_K distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitcT = new TCanvas("cmcRefitcT","cT distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitdca = new TCanvas("cmcRefitdca","dca distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitd0pi = new TCanvas("cmcRefitd0pi","d0pi distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitd0K = new TCanvas("cmcRefitd0K","d0K distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitd0xd0 = new TCanvas("cmcRefitd0xd0","d0xd0 distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitpointing = new TCanvas("cmcRefitpointing","cosPointingAngle distribution from MC - Refit cuts",50,50,550,550);
	TCanvas *cmcRefitphi = new TCanvas("cmcRefitphi","Phi distribution from MC - Refit cuts",50,50,550,550);

	// Reco-aod level
	TCanvas *cpt = new TCanvas("cpt","pt distribution from reco aod",50,50,550,550);
	TCanvas *cy = new TCanvas("cy","y distribution from reco aod",50,50,550,550);
	TCanvas *ccTS = new TCanvas("cTS","cosThetaStar distribution from reco aod",50,50,550,550); 
	TCanvas *cptPi = new TCanvas("cptPi","pt_pi distribution from reco aod",50,50,550,550);
	TCanvas *cptK = new TCanvas("cptK","pt_K distribution from reco aod",50,50,550,550);
	TCanvas *ccT = new TCanvas("ccT","cT distribution from reco aod",50,50,550,550);
	TCanvas *cdca = new TCanvas("cdca","dca distribution from reco aod",50,50,550,550);
	TCanvas *cd0pi = new TCanvas("cd0pi","d0pi distribution from reco aod",50,50,550,550);
	TCanvas *cd0K = new TCanvas("cd0K","d0K distribution from reco aod",50,50,550,550);
	TCanvas *cd0xd0 = new TCanvas("cd0xd0","d0xd0 distribution from reco aod",50,50,550,550);
	TCanvas *cpointing = new TCanvas("cpointing","cosPointingAngle distribution from reco aod",50,50,550,550);
	TCanvas *cphi = new TCanvas("cphi","Phi distribution from reco aod",50,50,550,550);

	// Reco-aod Acceptance level
	TCanvas *cAccpt = new TCanvas("cAccpt","pt distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccy = new TCanvas("cAccy","y distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAcccTS = new TCanvas("cAccTS","cosThetaStar distribution from reco aod - Acceptance cuts",50,50,550,550); 
	TCanvas *cAccptPi = new TCanvas("cAccptPi","pt_pi distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccptK = new TCanvas("cAccptK","pt_K distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccdca = new TCanvas("cAccdca","dca distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccd0pi = new TCanvas("cAccd0pi","d0pi distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccd0K = new TCanvas("cAccd0K","d0K distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccd0xd0 = new TCanvas("cAccd0xd0","d0xd0 distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccpointing = new TCanvas("cAccpointing","cosPointingAngle distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAcccT = new TCanvas("cAcccT","cT distribution from reco aod - Acceptance cuts",50,50,550,550);
	TCanvas *cAccphi = new TCanvas("cAccphi","Phi distribution from reco aod - Acceptance cuts",50,50,550,550);

	// Reco-aod ITS Clusters level
	TCanvas *cITSClpt = new TCanvas("cITSClpt","pt distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSCly = new TCanvas("cITSCly","y distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSClcTS = new TCanvas("cITSClTS","cosThetaStar distribution from reco aod - ITS Clusters cuts",50,50,550,550); 
	TCanvas *cITSClptPi = new TCanvas("cITSClptPi","pt_pi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSClptK = new TCanvas("cITSClptK","pt_K distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSClcT = new TCanvas("cITSClcT","cT distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSCldca = new TCanvas("cITSCldca","dca distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSCld0pi = new TCanvas("cITSCld0pi","d0pi distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSCld0K = new TCanvas("cITSCld0K","d0K distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSCld0xd0 = new TCanvas("cITSCld0xd0","d0xd0 distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSClpointing = new TCanvas("cITSClpointing","cosPointingAngle distribution from reco aod - ITS Clusters cuts",50,50,550,550);
	TCanvas *cITSClphi = new TCanvas("cITSClphi","Phi distribution from reco aod - ITS Clusters cuts",50,50,550,550);

	// Reco-aod PPR level
	TCanvas *cPPRpt = new TCanvas("cPPRpt","pt distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRy = new TCanvas("cPPRy","y distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRcTS = new TCanvas("cPPRTS","cosThetaStar distribution from reco aod - PPR cuts",50,50,550,550); 
	TCanvas *cPPRptPi = new TCanvas("cPPRptPi","pt_pi distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRptK = new TCanvas("cPPRptK","pt_K distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRcT = new TCanvas("cPPRcT","cT distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRdca = new TCanvas("cPPRdca","dca distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRd0pi = new TCanvas("cPPRd0pi","d0pi distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRd0K = new TCanvas("cPPRd0K","d0K distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRd0xd0 = new TCanvas("cPPRd0xd0","d0xd0 distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRpointing = new TCanvas("cPPRpointing","cosPointingAngle distribution from reco aod - PPR cuts",50,50,550,550);
	TCanvas *cPPRphi = new TCanvas("cPPRphi","Phi distribution from reco aod - PPR cuts",50,50,550,550);

	// pt, y, cosThetaStar for all steps
	cmc1->Divide(8,3);   

	cmc1->cd(1);
	TH1D *hMCpt1D = data->ShowProjection(ipt, stepGen);
	Double_t maxpt = hMCpt1D->GetMaximum();
	hMCpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	hMCpt1D->Sumw2();
	hMCpt1D->SetMinimum(0.01);
	SetHistoProperties(hMCpt1D,2,3,2,20,2,3005);
	hMCpt1D->GetXaxis()->SetTitleOffset(1.2);
	hMCpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC data");
	hMCpt1D->Draw("hist");
	cmcpt->cd();
	cmcpt->SetLeftMargin(0.15);
	cmcpt->SetRightMargin(0.05);
	hMCpt1D->Draw("hist");
	hMCpt1D->Draw("err same");
	cmcpt->Update();
	
	cmc1->cd(2);
	TH1D *hMCAccpt1D = data->ShowProjection(ipt, stepGenAcc);
	maxpt = hMCAccpt1D->GetMaximum();
	hMCAccpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	SetHistoProperties(hMCAccpt1D,8,3,8,20,8,3006);
	hMCAccpt1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccpt1D->Sumw2();
	hMCAccpt1D->SetMinimum(0.01);
	hMCAccpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Acc cut");
	hMCAccpt1D->Draw("hist");
	cmcAccpt->cd();
	cmcAccpt->SetLeftMargin(0.15);
	cmcAccpt->SetRightMargin(0.05);
	hMCAccpt1D->Draw("hist");
	hMCAccpt1D->Draw("err same");
	cmcAccpt->Update();

	cmc1->cd(3);
	TH1D *hMCVertexpt1D = data->ShowProjection(ipt, stepGenVertex);
	maxpt = hMCVertexpt1D->GetMaximum();
	hMCVertexpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	SetHistoProperties(hMCVertexpt1D,49,3,49,20,49,3006);
	hMCVertexpt1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexpt1D->Sumw2();
	hMCVertexpt1D->SetMinimum(0.01);
	hMCVertexpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Vertex cut");
	hMCVertexpt1D->Draw("hist");
	cmcVertexpt->cd();
	cmcVertexpt->SetLeftMargin(0.15);
	cmcVertexpt->SetRightMargin(0.05);
	hMCVertexpt1D->Draw("hist");
	hMCVertexpt1D->Draw("err same");
	cmcVertexpt->Update();

	cmc1->cd(4);
	TH1D *hMCRefitpt1D = data->ShowProjection(ipt, stepGenRefit);
	maxpt = hMCRefitpt1D->GetMaximum();
	hMCRefitpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	SetHistoProperties(hMCRefitpt1D,42,3,42,20,42,3006);
	hMCRefitpt1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitpt1D->Sumw2();
	hMCRefitpt1D->SetMinimum(0.01);
	hMCRefitpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Refit cut");
	hMCRefitpt1D->Draw("hist");
	cmcRefitpt->cd();
	cmcRefitpt->SetLeftMargin(0.15);
	cmcRefitpt->SetRightMargin(0.05);
	hMCRefitpt1D->Draw("hist");
	hMCRefitpt1D->Draw("err same");
	cmcRefitpt->Update();

	cmc1->cd(5);
	TH1D *hRECpt1D = data->ShowProjection(ipt, stepRec);
	maxpt = hRECpt1D->GetMaximum();
	hRECpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	SetHistoProperties(hRECpt1D,4,3,4,20,4,3004);
	hRECpt1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECpt1D->Sumw2();
	hRECpt1D->SetMinimum(0.01);
	hRECpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD");
	hRECpt1D->Draw("hist");
	cpt->cd();
	cpt->SetLeftMargin(0.15);
	cpt->SetRightMargin(0.05);
	hRECpt1D->Draw("hist");
	hRECpt1D->Draw("err same");
	cpt->Update();

	cmc1->cd(6);
	TH1D *hRECAccpt1D = data->ShowProjection(ipt, stepRecAcc);
	maxpt = hRECAccpt1D->GetMaximum();
	hRECAccpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	SetHistoProperties(hRECAccpt1D,38,3,38,20,38,3007);
	hRECAccpt1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccpt1D->Sumw2();
	hRECAccpt1D->SetMinimum(0.01);
	hRECAccpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - Acc cuts");
	hRECAccpt1D->Draw("hist");
	cAccpt->cd();
	cAccpt->SetLeftMargin(0.15);
	cAccpt->SetRightMargin(0.05);
	hRECAccpt1D->Draw("hist");
	hRECAccpt1D->Draw("err same");
	cAccpt->Update();

	cmc1->cd(7);
	TH1D *hRECITSClpt1D = data->ShowProjection(ipt, stepRecITSClus);
	maxpt = hRECITSClpt1D->GetMaximum();
	hRECITSClpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	SetHistoProperties(hRECITSClpt1D,13,3,13,20,13,3019);
	hRECITSClpt1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSClpt1D->Sumw2();
	hRECITSClpt1D->SetMinimum(0.01);
	hRECITSClpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - ITS Cl cuts");
	hRECITSClpt1D->Draw("hist");
	cITSClpt->cd();
	cITSClpt->SetLeftMargin(0.15);
	cITSClpt->SetRightMargin(0.05);
	hRECITSClpt1D->Draw("hist");
	hRECITSClpt1D->Draw("err same");
	cITSClpt->Update();

	cmc1->cd(8);
	TH1D *hRECPPRpt1D = data->ShowProjection(ipt, stepRecPPR);
	maxpt = hRECPPRpt1D->GetMaximum();
	hRECPPRpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2);
	SetHistoProperties(hRECPPRpt1D,6,3,6,20,6,3019);
	hRECPPRpt1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRpt1D->Sumw2();
	hRECPPRpt1D->SetMinimum(0.01);
	hRECPPRpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD - PPR cuts");
	hRECPPRpt1D->Draw("hist");
	cPPRpt->cd();
	cPPRpt->SetLeftMargin(0.15);
	cPPRpt->SetRightMargin(0.05);
	hRECPPRpt1D->Draw("hist");
	hRECPPRpt1D->Draw("err same");
	cPPRpt->Update();

	cmc1->cd(9);
	TH1D *hMCy1D = data->ShowProjection(iy, stepGen);
	Double_t maxy = hMCy1D->GetMaximum();
	hMCy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hMCy1D,2,3,2,20,2,3005);
	hMCy1D->SetLineColor(2);
	hMCy1D->SetLineWidth(3);
	hMCy1D->SetMarkerColor(2);
	hMCy1D->SetFillColor(2);
	hMCy1D->SetFillStyle(3005);
	hMCy1D->SetMarkerStyle(20);
	hMCy1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCy1D->Sumw2();
	hMCy1D->SetMinimum(0.01);
	hMCy1D->GetXaxis()->SetTitle("y, MC data");
	hMCy1D->Draw("hist");
	cmcy->cd();
	cmcy->SetLeftMargin(0.15);
	cmcy->SetRightMargin(0.05);
	hMCy1D->Draw("hist");
	hMCy1D->Draw("err same");
	cmcy->Update();

	cmc1->cd(10);
	TH1D *hMCAccy1D = data->ShowProjection(iy, stepGenAcc);
	maxy = hMCAccy1D->GetMaximum();
	hMCAccy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hMCAccy1D,8,3,8,20,8,3006);
	hMCAccy1D->GetXaxis()->SetTitleOffset(1.2);
	hMCAccy1D->Sumw2();
	hMCAccy1D->SetMinimum(0.01);
	hMCAccy1D->GetXaxis()->SetTitle("y, MC data - Acc cuts");
	hMCAccy1D->Draw("hist");
	cmcAccy->cd();
	cmcAccy->SetLeftMargin(0.15);
	cmcAccy->SetRightMargin(0.05);
	hMCAccy1D->Draw("hist");
	hMCAccy1D->Draw("err same");
	cmcy->Update();

	cmc1->cd(11);
	TH1D *hMCVertexy1D = data->ShowProjection(iy, stepGenVertex);
	maxy = hMCVertexy1D->GetMaximum();
	hMCVertexy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hMCVertexy1D,49,3,49,20,49,3006);
	hMCVertexy1D->GetXaxis()->SetTitleOffset(1.2);
	hMCVertexy1D->Sumw2();
	hMCVertexy1D->SetMinimum(0.01);
	hMCVertexy1D->GetXaxis()->SetTitle("y, MC data - Vertex cuts");
	hMCVertexy1D->Draw("hist");
	cmcVertexy->cd();
	cmcVertexy->SetLeftMargin(0.15);
	cmcVertexy->SetRightMargin(0.05);
	hMCVertexy1D->Draw("hist");
	hMCVertexy1D->Draw("err same");
	cmcy->Update();

	cmc1->cd(12);
	TH1D *hMCRefity1D = data->ShowProjection(iy, stepGenRefit);
	maxy = hMCRefity1D->GetMaximum();
	hMCRefity1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hMCRefity1D,42,3,42,20,42,3006);
	hMCRefity1D->GetXaxis()->SetTitleOffset(1.2);
	hMCRefity1D->Sumw2();
	hMCRefity1D->SetMinimum(0.01);
	hMCRefity1D->GetXaxis()->SetTitle("y, MC data - Refit cuts");
	hMCRefity1D->Draw("hist");
	cmcRefity->cd();
	cmcRefity->SetLeftMargin(0.15);
	cmcRefity->SetRightMargin(0.05);
	hMCRefity1D->Draw("hist");
	hMCRefity1D->Draw("err same");
	cmcy->Update();

	cmc1->cd(13);
	TH1D *hRECy1D = data->ShowProjection(iy, stepRec);
	maxy = hRECy1D->GetMaximum();
	hRECy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hRECy1D,4,3,4,20,4,3004);
	hRECy1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECy1D->Sumw2();
	hRECy1D->SetMinimum(0.01);
	hRECy1D->Draw("hist");
	cy->cd();
	cy->SetLeftMargin(0.15);
	cy->SetRightMargin(0.05);
	hRECy1D->GetXaxis()->SetTitle("y, AOD");
	hRECy1D->Draw("hist");
	hRECy1D->Draw("err same");
	cy->Update();

	cmc1->cd(14);
	TH1D *hRECAccy1D = data->ShowProjection(iy, stepRecAcc);
	maxy = hRECAccy1D->GetMaximum();
	hRECAccy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hRECAccy1D,38,3,38,20,38,3007);
	hRECAccy1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccy1D->Sumw2();
	hRECAccy1D->SetMinimum(0.01);
	hRECAccy1D->Draw("hist");
	cAccy->cd();
	cAccy->SetLeftMargin(0.15);
	cAccy->SetRightMargin(0.05);
	hRECAccy1D->GetXaxis()->SetTitle("y, AOD - Acc cuts");
	hRECAccy1D->Draw("hist");
	hRECAccy1D->Draw("err same");
	cAccy->Update();

	cmc1->cd(15);
	TH1D *hRECITSCly1D = data->ShowProjection(iy, stepRecITSClus);
	maxy = hRECITSCly1D->GetMaximum();
	hRECITSCly1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hRECITSCly1D,13,3,13,20,13,3019);
	hRECITSCly1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSCly1D->Sumw2();
	hRECITSCly1D->SetMinimum(0.01);
	hRECITSCly1D->Draw("hist");
	cITSCly->cd();
	cITSCly->SetLeftMargin(0.15);
	cITSCly->SetRightMargin(0.05);
	hRECITSCly1D->GetXaxis()->SetTitle("y, AOD - ITS Cl cuts");
	hRECITSCly1D->Draw("hist");
	hRECITSCly1D->Draw("err same");
	cITSCly->Update();

	cmc1->cd(16);
	TH1D *hRECPPRy1D = data->ShowProjection(iy, stepRecPPR);
	maxy = hRECPPRy1D->GetMaximum();
	hRECPPRy1D->GetYaxis()->SetRangeUser(0,maxy*1.2);
	SetHistoProperties(hRECPPRy1D,6,3,6,20,6,3019);
	hRECPPRy1D->SetMarkerStyle(20);
	hRECPPRy1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRy1D->Sumw2();
	hRECPPRy1D->SetMinimum(0.01);
	hRECPPRy1D->Draw("hist");
	cPPRy->cd();
	cPPRy->SetLeftMargin(0.15);
	cPPRy->SetRightMargin(0.05);
	hRECPPRy1D->GetXaxis()->SetTitle("y, AOD - PPR cuts");
	hRECPPRy1D->Draw("hist");
	hRECPPRy1D->Draw("err same");
	cPPRy->Update();

	cmc1->cd(17);
	TH1D *hMCcTS1D = data->ShowProjection(icTS, stepGen);
	Double_t maxcTS = hMCcTS1D->GetMaximum();
	hMCcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hMCcTS1D,2,3,2,20,2,3005);
	hMCcTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCcTS1D->Sumw2();
	hMCcTS1D->SetMinimum(0.01);
	hMCcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data");
	hMCcTS1D->Draw("hist");
	cmccTS->cd();
	cmccTS->SetLeftMargin(0.15);
	cmccTS->SetRightMargin(0.05);
	hMCcTS1D->Draw("hist");
	hMCcTS1D->Draw("err same");
	cmccTS->Update();

	cmc1->cd(18);
	TH1D *hMCAcccTS1D = data->ShowProjection(icTS, stepGenAcc);
	maxcTS = hMCAcccTS1D->GetMaximum();
	hMCAcccTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hMCAcccTS1D,8,3,8,20,8,3006);
	hMCAcccTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAcccTS1D->Sumw2();
	hMCAcccTS1D->SetMinimum(0.01);
	hMCAcccTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Acc cuts");
	hMCAcccTS1D->Draw("hist");
	cmcAcccTS->cd();
	cmcAcccTS->SetLeftMargin(0.15);
	cmcAcccTS->SetRightMargin(0.05);
	hMCAcccTS1D->Draw("hist");
	hMCAcccTS1D->Draw("err same");
	cmcAcccTS->Update();

	cmc1->cd(19);
	TH1D *hMCVertexcTS1D = data->ShowProjection(icTS, stepGenVertex);
	maxcTS = hMCVertexcTS1D->GetMaximum();
	hMCVertexcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hMCVertexcTS1D,49,3,49,20,49,3006);
	hMCVertexcTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexcTS1D->Sumw2();
	hMCVertexcTS1D->SetMinimum(0.01);
	hMCVertexcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Vertex cuts");
	hMCVertexcTS1D->Draw("hist");
	cmcVertexcTS->cd();
	cmcVertexcTS->SetLeftMargin(0.15);
	cmcVertexcTS->SetRightMargin(0.05);
	hMCVertexcTS1D->Draw("hist");
	hMCVertexcTS1D->Draw("err same");
	cmcVertexcTS->Update();

	cmc1->cd(20);
	TH1D *hMCRefitcTS1D = data->ShowProjection(icTS, stepGenRefit);
	maxcTS = hMCRefitcTS1D->GetMaximum();
	hMCRefitcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hMCRefitcTS1D,42,3,42,20,42,3006);
	hMCRefitcTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitcTS1D->Sumw2();
	hMCRefitcTS1D->SetMinimum(0.01);
	hMCRefitcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data - Refit cuts");
	hMCRefitcTS1D->Draw("hist");
	cmcRefitcTS->cd();
	cmcRefitcTS->SetLeftMargin(0.15);
	cmcRefitcTS->SetRightMargin(0.05);
	hMCRefitcTS1D->Draw("hist");
	hMCRefitcTS1D->Draw("err same");
	cmcRefitcTS->Update();

	cmc1->cd(21);
	TH1D *hRECcTS1D = data->ShowProjection(icTS, stepRec);
	maxcTS = hRECcTS1D->GetMaximum();
	hRECcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hRECcTS1D,4,3,4,20,4,3004);
	hRECcTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECcTS1D->Sumw2();
	hRECcTS1D->SetMinimum(0.01);
	hRECcTS1D->Draw("hist");
	ccTS->cd();
	ccTS->SetLeftMargin(0.15);
	ccTS->SetRightMargin(0.05);
	hRECcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD");
	hRECcTS1D->Draw("hist");
	hRECcTS1D->Draw("err same");
	ccTS->Update();

	cmc1->cd(22);
	TH1D *hRECAcccTS1D = data->ShowProjection(icTS, stepRecAcc);
	maxcTS = hRECAcccTS1D->GetMaximum();
	hRECAcccTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hRECAcccTS1D,38,3,38,20,38,3007);
	hRECAcccTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAcccTS1D->Sumw2();
	hRECAcccTS1D->SetMinimum(0.01);
	hRECAcccTS1D->Draw("hist");
	cAcccTS->cd();
	cAcccTS->SetLeftMargin(0.15);
	cAcccTS->SetRightMargin(0.05);
	hRECAcccTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - Acc cuts");
	hRECAcccTS1D->Draw("hist");
	hRECAcccTS1D->Draw("err same");
	cAcccTS->Update();

	cmc1->cd(23);
	TH1D *hRECITSClcTS1D = data->ShowProjection(icTS, stepRecITSClus);
	maxcTS = hRECITSClcTS1D->GetMaximum();
	hRECITSClcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hRECITSClcTS1D,13,3,13,20,13,3019);
	hRECITSClcTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSClcTS1D->Sumw2();
	hRECITSClcTS1D->SetMinimum(0.01);
	hRECITSClcTS1D->Draw("hist");
	cITSClcTS->cd();
	cITSClcTS->SetLeftMargin(0.15);
	cITSClcTS->SetRightMargin(0.05);
	hRECITSClcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - ITS Cl cuts");
	hRECITSClcTS1D->Draw("hist");
	hRECITSClcTS1D->Draw("err same");
	cITSClcTS->Update();

	cmc1->cd(24);
	TH1D *hRECPPRcTS1D = data->ShowProjection(icTS, stepRecPPR);
	maxcTS = hRECPPRcTS1D->GetMaximum();
	hRECPPRcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2);
	SetHistoProperties(hRECPPRcTS1D,6,3,6,20,6,3019);
	hRECPPRcTS1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRcTS1D->Sumw2();
	hRECPPRcTS1D->SetMinimum(0.01);
	hRECPPRcTS1D->Draw("hist");
	cPPRcTS->cd();
	cPPRcTS->SetLeftMargin(0.15);
	cPPRcTS->SetRightMargin(0.05);
	hRECPPRcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD - PPR cuts");
	hRECPPRcTS1D->Draw("hist");
	hRECPPRcTS1D->Draw("err same");
	cPPRcTS->Update();

	// ptPi, ptK, cT
	cmc2->Divide(8,3); 

	cmc2->cd(1);
	TH1D *hMCptPi1D = data->ShowProjection(iptPi, stepGen);
	Double_t maxptPi = hMCptPi1D->GetMaximum();
	hMCptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	hMCptPi1D->Sumw2();
	SetHistoProperties(hMCptPi1D,2,3,2,20,2,3005);
	hMCptPi1D->SetMinimum(0.01);
	hMCptPi1D->GetXaxis()->SetTitleOffset(1.2);
	hMCptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data");
	hMCptPi1D->Draw("hist");
	cmcptPi->cd();
	cmcptPi->SetLeftMargin(0.15);
	cmcptPi->SetRightMargin(0.05);
	hMCptPi1D->Draw("hist");
	hMCptPi1D->Draw("err same");
	cmcptPi->Update();

	cmc2->cd(2);
	TH1D *hMCAccptPi1D = data->ShowProjection(iptPi, stepGenAcc);
	maxptPi = hMCAccptPi1D->GetMaximum();
	hMCAccptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	hMCAccptPi1D->Sumw2();
	hMCAccptPi1D->SetMinimum(0.01);
	SetHistoProperties(hMCAccptPi1D,8,3,8,20,8,3006);
	hMCAccptPi1D->GetXaxis()->SetTitleOffset(1.2);
	hMCAccptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Acc cuts");
	hMCAccptPi1D->Draw("hist");
	cmcAccptPi->cd();
	cmcAccptPi->SetLeftMargin(0.15);
	cmcAccptPi->SetRightMargin(0.05);
	hMCAccptPi1D->Draw("hist");
	hMCAccptPi1D->Draw("err same");
	cmcAccptPi->Update();

	cmc2->cd(3);
	TH1D *hMCVertexptPi1D = data->ShowProjection(iptPi, stepGenVertex);
	maxptPi = hMCVertexptPi1D->GetMaximum();
	hMCVertexptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	hMCVertexptPi1D->Sumw2();
	hMCVertexptPi1D->SetMinimum(0.01);
	SetHistoProperties(hMCVertexptPi1D,49,3,49,20,49,3006);
	hMCVertexptPi1D->GetXaxis()->SetTitleOffset(1.2);
	hMCVertexptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Vertex cuts");
	hMCVertexptPi1D->Draw("hist");
	cmcVertexptPi->cd();
	cmcVertexptPi->SetLeftMargin(0.15);
	cmcVertexptPi->SetRightMargin(0.05);
	hMCVertexptPi1D->Draw("hist");
	hMCVertexptPi1D->Draw("err same");
	cmcVertexptPi->Update();

	cmc2->cd(4);
	TH1D *hMCRefitptPi1D = data->ShowProjection(iptPi, stepGenRefit);
	maxptPi = hMCRefitptPi1D->GetMaximum();
	hMCRefitptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	hMCRefitptPi1D->Sumw2();
	hMCRefitptPi1D->SetMinimum(0.01);
	SetHistoProperties(hMCRefitptPi1D,42,3,42,20,42,3006);
	hMCRefitptPi1D->GetXaxis()->SetTitleOffset(1.2);
	hMCRefitptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data - Refit cuts");
	hMCRefitptPi1D->Draw("hist");
	cmcRefitptPi->cd();
	cmcRefitptPi->SetLeftMargin(0.15);
	cmcRefitptPi->SetRightMargin(0.05);
	hMCRefitptPi1D->Draw("hist");
	hMCRefitptPi1D->Draw("err same");
	cmcRefitptPi->Update();

	cmc2->cd(5);
	TH1D *hRECptPi1D = data->ShowProjection(iptPi, stepRec);
	maxptPi = hRECptPi1D->GetMaximum();
	hRECptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	SetHistoProperties(hRECptPi1D,4,3,4,20,4,3004);
	hRECptPi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECptPi1D->Sumw2();
	hRECptPi1D->SetMinimum(0.01);
	hRECptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD");
	hRECptPi1D->Draw("hist");
	cptPi->cd();
	cptPi->SetLeftMargin(0.15);
	cptPi->SetRightMargin(0.05);
	hRECptPi1D->Draw("hist");
	hRECptPi1D->Draw("err same");
	cptPi->Update();

	cmc2->cd(6);
	TH1D *hRECAccptPi1D = data->ShowProjection(iptPi, stepRecAcc);
	maxptPi = hRECAccptPi1D->GetMaximum();
	hRECAccptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	SetHistoProperties(hRECAccptPi1D,38,3,38,20,38,3007);
	hRECAccptPi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccptPi1D->Sumw2();
	hRECAccptPi1D->SetMinimum(0.01);
	hRECAccptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - Acc cuts");
	hRECAccptPi1D->Draw("hist");
	cAccptPi->cd();
	cAccptPi->SetLeftMargin(0.15);
	cAccptPi->SetRightMargin(0.05);
	hRECAccptPi1D->Draw("hist");
	hRECAccptPi1D->Draw("err same");
	cAccptPi->Update();

	cmc2->cd(7);
	TH1D *hRECITSClptPi1D = data->ShowProjection(iptPi, stepRecITSClus);
	maxptPi = hRECITSClptPi1D->GetMaximum();
	hRECITSClptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	SetHistoProperties(hRECITSClptPi1D,13,3,13,20,13,3019);
	hRECITSClptPi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSClptPi1D->Sumw2();
	hRECITSClptPi1D->SetMinimum(0.01);
	hRECITSClptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - ITS Cl cuts");
	hRECITSClptPi1D->Draw("hist");
	cITSClptPi->cd();
	cITSClptPi->SetLeftMargin(0.15);
	cITSClptPi->SetRightMargin(0.05);
	hRECITSClptPi1D->Draw("hist");
	hRECITSClptPi1D->Draw("err same");
	cITSClptPi->Update();

	cmc2->cd(8);
	TH1D *hRECPPRptPi1D = data->ShowProjection(iptPi, stepRecPPR);
	maxptPi = hRECPPRptPi1D->GetMaximum();
	hRECPPRptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2);
	SetHistoProperties(hRECPPRptPi1D,6,3,6,20,6,3019);
	hRECPPRptPi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRptPi1D->Sumw2();
	hRECPPRptPi1D->SetMinimum(0.01);
	hRECPPRptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD - PPR cuts");
	hRECPPRptPi1D->Draw("hist");
	cPPRptPi->cd();
	cPPRptPi->SetLeftMargin(0.15);
	cPPRptPi->SetRightMargin(0.05);
	hRECPPRptPi1D->Draw("hist");
	hRECPPRptPi1D->Draw("err same");
	cPPRptPi->Update();

	cmc2->cd(9);
	TH1D *hMCptK1D = data->ShowProjection(iptK, stepGen);
	Double_t maxptK = hMCptK1D->GetMaximum();
	hMCptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hMCptK1D,2,3,2,20,2,3005);
	hMCptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCptK1D->Sumw2();
	hMCptK1D->SetMinimum(0.01);
	hMCptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data");
	hMCptK1D->Draw("hist");
	cmcptK->cd();
	cmcptK->SetLeftMargin(0.15);
	cmcptK->SetRightMargin(0.05);
	hMCptK1D->Draw("hist");
	hMCptK1D->Draw("err same");
	cmcptK->Update();

	cmc2->cd(10);
	TH1D *hMCAccptK1D = data->ShowProjection(iptK, stepGenAcc);
	maxptK = hMCAccptK1D->GetMaximum();
	hMCAccptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hMCAccptK1D,8,3,8,20,8,3006);
	hMCAccptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccptK1D->Sumw2();
	hMCAccptK1D->SetMinimum(0.01);
	hMCAccptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Acc cuts");
	hMCAccptK1D->Draw("hist");
	cmcAccptK->cd();
	cmcAccptK->SetLeftMargin(0.15);
	cmcAccptK->SetRightMargin(0.05);
	hMCAccptK1D->Draw("hist");
	hMCAccptK1D->Draw("err same");
	cmcAccptK->Update();

	cmc2->cd(11);
	TH1D *hMCVertexptK1D = data->ShowProjection(iptK, stepGenVertex);
	maxptK = hMCVertexptK1D->GetMaximum();
	hMCVertexptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hMCVertexptK1D,49,3,49,20,49,3006);
	hMCVertexptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexptK1D->Sumw2();
	hMCVertexptK1D->SetMinimum(0.01);
	hMCVertexptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Vertex cuts");
	hMCVertexptK1D->Draw("hist");
	cmcVertexptK->cd();
	cmcVertexptK->SetLeftMargin(0.15);
	cmcVertexptK->SetRightMargin(0.05);
	hMCVertexptK1D->Draw("hist");
	hMCVertexptK1D->Draw("err same");
	cmcVertexptK->Update();

	cmc2->cd(12);
	TH1D *hMCRefitptK1D = data->ShowProjection(iptK, stepGenRefit);
	maxptK = hMCRefitptK1D->GetMaximum();
	hMCRefitptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hMCRefitptK1D,42,3,42,20,42,3006);
	hMCRefitptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitptK1D->Sumw2();
	hMCRefitptK1D->SetMinimum(0.01);
	hMCRefitptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data - Refit cuts");
	hMCRefitptK1D->Draw("hist");
	cmcRefitptK->cd();
	cmcRefitptK->SetLeftMargin(0.15);
	cmcRefitptK->SetRightMargin(0.05);
	hMCRefitptK1D->Draw("hist");
	hMCRefitptK1D->Draw("err same");
	cmcRefitptK->Update();

	cmc2->cd(13);
	TH1D *hRECptK1D = data->ShowProjection(iptK, stepRec);
	maxptK = hRECptK1D->GetMaximum();
	hRECptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hRECptK1D,4,3,4,20,4,3004);
	hRECptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECptK1D->Sumw2();
	hRECptK1D->SetMinimum(0.01);
	hRECptK1D->Draw("hist");
	cptK->cd();
	cptK->SetLeftMargin(0.15);
	cptK->SetRightMargin(0.05);
	hRECptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD");
	hRECptK1D->Draw("hist");
	hRECptK1D->Draw("err same");
	cptK->Update();

	cmc2->cd(14);
	TH1D *hRECAccptK1D = data->ShowProjection(iptK, stepRecAcc);
	maxptK = hRECAccptK1D->GetMaximum();
	hRECAccptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hRECAccptK1D,38,3,38,20,38,3007);
	hRECAccptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccptK1D->Sumw2();
	hRECAccptK1D->SetMinimum(0.01);
	hRECAccptK1D->Draw("hist");
	cAccptK->cd();
	cAccptK->SetLeftMargin(0.15);
	cAccptK->SetRightMargin(0.05);
	hRECAccptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - Acc cuts");
	hRECAccptK1D->Draw("hist");
	hRECAccptK1D->Draw("err same");
	cAccptK->Update();

	cmc2->cd(15);
	TH1D *hRECITSClptK1D = data->ShowProjection(iptK, stepRecITSClus);
	maxptK = hRECITSClptK1D->GetMaximum();
	hRECITSClptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hRECITSClptK1D,13,3,13,20,13,3019);
	hRECITSClptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSClptK1D->Sumw2();
	hRECITSClptK1D->SetMinimum(0.01);
	hRECITSClptK1D->Draw("hist");
	cITSClptK->cd();
	cITSClptK->SetLeftMargin(0.15);
	cITSClptK->SetRightMargin(0.05);
	hRECITSClptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - ITS Cl cuts");
	hRECITSClptK1D->Draw("hist");
	hRECITSClptK1D->Draw("err same");
	cITSClptK->Update();

	cmc2->cd(16);
	TH1D *hRECPPRptK1D = data->ShowProjection(iptK, stepRecPPR);
	maxptK = hRECPPRptK1D->GetMaximum();
	hRECPPRptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2);
	SetHistoProperties(hRECPPRptK1D,6,3,6,20,6,3019);
	hRECPPRptK1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRptK1D->Sumw2();
	hRECPPRptK1D->SetMinimum(0.01);
	hRECPPRptK1D->Draw("hist");
	cPPRptK->cd();
	cPPRptK->SetLeftMargin(0.15);
	cPPRptK->SetRightMargin(0.05);
	hRECPPRptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD - PPR cuts");
	hRECPPRptK1D->Draw("hist");
	hRECPPRptK1D->Draw("err same");
	cPPRptK->Update();

	cmc2->cd(17);
	TH1D *hMCcT1D = data->ShowProjection(icT, stepGen);
	Double_t maxcT = hMCcT1D->GetMaximum();
	hMCcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hMCcT1D,2,3,2,20,2,3005);
	hMCcT1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCcT1D->Sumw2();
	hMCcT1D->SetMinimum(0.01);
	hMCcT1D->GetXaxis()->SetTitle("ct (#mum), MC data");
	hMCcT1D->Draw("hist");
	cmccT->cd();
	cmccT->SetLeftMargin(0.15);
	cmccT->SetRightMargin(0.05);
	hMCcT1D->Draw("hist");
	hMCcT1D->Draw("err same");
	cmccT->Update();

	cmc2->cd(18);
	TH1D *hMCAcccT1D = data->ShowProjection(icT, stepGenAcc);
	maxcT = hMCAcccT1D->GetMaximum();
	hMCAcccT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hMCAcccT1D,8,3,8,20,8,3006);
	hMCAcccT1D->GetXaxis()->SetTitleOffset(1.2);
	hMCAcccT1D->Sumw2();
	hMCAcccT1D->SetMinimum(0.01);
	hMCAcccT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Acc cuts");
	hMCAcccT1D->Draw("hist");
	cmcAcccT->cd();
	cmcAcccT->SetLeftMargin(0.15);
	cmcAcccT->SetRightMargin(0.05);
	hMCAcccT1D->Draw("hist");
	hMCAcccT1D->Draw("err same");
	cmcAcccT->Update();

	cmc2->cd(19);
	TH1D *hMCVertexcT1D = data->ShowProjection(icT, stepGenVertex);
	maxcT = hMCVertexcT1D->GetMaximum();
	hMCVertexcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hMCVertexcT1D,49,3,49,20,49,3006);
	hMCVertexcT1D->GetXaxis()->SetTitleOffset(1.2);
	hMCVertexcT1D->Sumw2();
	hMCVertexcT1D->SetMinimum(0.01);
	hMCVertexcT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Vertex cuts");
	hMCVertexcT1D->Draw("hist");
	cmcVertexcT->cd();
	cmcVertexcT->SetLeftMargin(0.15);
	cmcVertexcT->SetRightMargin(0.05);
	hMCVertexcT1D->Draw("hist");
	hMCVertexcT1D->Draw("err same");
	cmcVertexcT->Update();

	cmc2->cd(20);
	TH1D *hMCRefitcT1D = data->ShowProjection(icT, stepGenRefit);
	maxcT = hMCRefitcT1D->GetMaximum();
	hMCRefitcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hMCRefitcT1D,42,3,42,20,42,3006);
	hMCRefitcT1D->GetXaxis()->SetTitleOffset(1.2);
	hMCRefitcT1D->Sumw2();
	hMCRefitcT1D->SetMinimum(0.01);
	hMCRefitcT1D->GetXaxis()->SetTitle("ct (#mum), MC data - Refit cuts");
	hMCRefitcT1D->Draw("hist");
	cmcRefitcT->cd();
	cmcRefitcT->SetLeftMargin(0.15);
	cmcRefitcT->SetRightMargin(0.05);
	hMCRefitcT1D->Draw("hist");
	hMCRefitcT1D->Draw("err same");
	cmcRefitcT->Update();

	cmc2->cd(21);
	TH1D *hRECcT1D = data->ShowProjection(icT, stepRec);
	maxcT = hRECcT1D->GetMaximum();
	hRECcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hRECcT1D,4,3,4,20,4,3004);
	hRECcT1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECcT1D->Sumw2();
	hRECcT1D->SetMinimum(0.01);
	hRECcT1D->Draw("hist");
	ccT->cd();
	ccT->SetLeftMargin(0.15);
	ccT->SetRightMargin(0.05);
	hRECcT1D->GetXaxis()->SetTitle("ct (#mum), AOD");
	hRECcT1D->Draw("hist");
	hRECcT1D->Draw("err same");
	ccT->Update();

	cmc2->cd(22);
	TH1D *hRECAcccT1D = data->ShowProjection(icT, stepRecAcc);
	maxcT = hRECAcccT1D->GetMaximum();
	hRECAcccT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hRECAcccT1D,38,3,38,20,38,3007);
	hRECAcccT1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAcccT1D->Sumw2();
	hRECAcccT1D->SetMinimum(0.01);
	hRECAcccT1D->Draw("hist");
	cAcccT->cd();
	cAcccT->SetLeftMargin(0.15);
	cAcccT->SetRightMargin(0.05);
	hRECAcccT1D->GetXaxis()->SetTitle("ct (#mum), AOD - Acc cuts");
	hRECAcccT1D->Draw("hist");
	hRECAcccT1D->Draw("err same");
	cAcccT->Update();

	cmc2->cd(23);
	TH1D *hRECITSClcT1D = data->ShowProjection(icT, stepRecITSClus);
	maxcT = hRECITSClcT1D->GetMaximum();
	hRECITSClcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hRECITSClcT1D,13,3,13,20,13,3019);
	hRECITSClcT1D->GetXaxis()->SetTitleOffset(1.2);
	hRECITSClcT1D->Sumw2();
	hRECITSClcT1D->SetMinimum(0.01);
	hRECITSClcT1D->Draw("hist");
	cITSClcT->cd();
	cITSClcT->SetLeftMargin(0.15);
	cITSClcT->SetRightMargin(0.05);
	hRECITSClcT1D->GetXaxis()->SetTitle("ct (#mum), AOD - ITS Cl cuts");
	hRECITSClcT1D->Draw("hist");
	hRECITSClcT1D->Draw("err same");
	cITSClcT->Update();

	cmc2->cd(24);
	TH1D *hRECPPRcT1D = data->ShowProjection(icT, stepRecPPR);
	maxcT = hRECPPRcT1D->GetMaximum();
	hRECPPRcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2);
	SetHistoProperties(hRECPPRcT1D,6,3,6,20,6,3019);
	hRECPPRcT1D->GetXaxis()->SetTitleOffset(1.2);
	hRECPPRcT1D->Sumw2();
	hRECPPRcT1D->SetMinimum(0.01);
	hRECPPRcT1D->Draw("hist");
	cPPRcT->cd();
	cPPRcT->SetLeftMargin(0.15);
	cPPRcT->SetRightMargin(0.05);
	hRECPPRcT1D->GetXaxis()->SetTitle("ct (#mum), AOD - PPR cuts");
	hRECPPRcT1D->Draw("hist");
	hRECPPRcT1D->Draw("err same");
	cPPRcT->Update();

	//dca, d0pi, d0K
	cmc3->Divide(8,3);
	cmc3->cd(1);
	TH1D *hMCdca1D = data->ShowProjection(idca, stepGen);
	Double_t maxdca = hMCdca1D->GetMaximum();
	hMCdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	hMCdca1D->Sumw2();
	hMCdca1D->SetMinimum(0.01);
	SetHistoProperties(hMCdca1D,2,3,2,20,2,3005);
	hMCdca1D->GetXaxis()->SetTitleOffset(1.2);
	hMCdca1D->GetXaxis()->SetTitle("dca (#mum), MC data");
	hMCdca1D->Draw("hist");
	cmcdca->cd();
	cmcdca->SetLeftMargin(0.15);
	cmcdca->SetRightMargin(0.05);
	hMCdca1D->Draw("hist");
	hMCdca1D->Draw("err same");
	cmcdca->Update();

	cmc3->cd(2);
	TH1D *hMCAccdca1D = data->ShowProjection(idca, stepGenAcc);
	maxdca = hMCAccdca1D->GetMaximum();
	hMCAccdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	SetHistoProperties(hMCAccdca1D,8,3,8,20,8,3006);
	hMCAccdca1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccdca1D->Sumw2();
	hMCAccdca1D->SetMinimum(0.01);
	hMCAccdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Acc cut");
	hMCAccdca1D->Draw("hist");
	cmcAccdca->cd();
	cmcAccdca->SetLeftMargin(0.15);
	cmcAccdca->SetRightMargin(0.05);
	hMCAccdca1D->Draw("hist");
	hMCAccdca1D->Draw("err same");
	cmcAccdca->Update();

	cmc3->cd(3);
	TH1D *hMCVertexdca1D = data->ShowProjection(idca, stepGenVertex);
	maxdca = hMCVertexdca1D->GetMaximum();
	hMCVertexdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	SetHistoProperties(hMCVertexdca1D,49,3,49,20,49,3006);
	hMCVertexdca1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexdca1D->Sumw2();
	hMCVertexdca1D->SetMinimum(0.01);
	hMCVertexdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Vertex cut");
	hMCVertexdca1D->Draw("hist");
	cmcVertexdca->cd();
	cmcVertexdca->SetLeftMargin(0.15);
	cmcVertexdca->SetRightMargin(0.05);
	hMCVertexdca1D->Draw("hist");
	hMCVertexdca1D->Draw("err same");
	cmcVertexdca->Update();

	cmc3->cd(4);
	TH1D *hMCRefitdca1D = data->ShowProjection(idca, stepGenRefit);
	maxdca = hMCRefitdca1D->GetMaximum();
	hMCRefitdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	SetHistoProperties(hMCRefitdca1D,42,3,42,20,42,3006);
	hMCRefitdca1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitdca1D->Sumw2();
	hMCRefitdca1D->SetMinimum(0.01);
	hMCRefitdca1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC - Refit cut");
	hMCRefitdca1D->Draw("hist");
	cmcRefitdca->cd();
	cmcRefitdca->SetLeftMargin(0.15);
	cmcRefitdca->SetRightMargin(0.05);
	hMCRefitdca1D->Draw("hist");
	hMCRefitdca1D->Draw("err same");
	cmcRefitdca->Update();

	cmc3->cd(5);
	TH1D *hRECdca1D = data->ShowProjection(idca, stepRec);
	maxdca = hRECdca1D->GetMaximum();
	hRECdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	SetHistoProperties(hRECdca1D,4,3,4,20,4,3004);
	hRECdca1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECdca1D->Sumw2();
	hRECdca1D->SetMinimum(0.01);
	hRECdca1D->GetXaxis()->SetTitle("dca (#mum), AOD");
	hRECdca1D->Draw("hist");
	cdca->cd();
	cdca->SetLeftMargin(0.15);
	cdca->SetRightMargin(0.05);
	hRECdca1D->Draw("hist");
	hRECdca1D->Draw("err same");
	cdca->Update();

	cmc3->cd(6);
	TH1D *hRECAccdca1D = data->ShowProjection(idca, stepRecAcc);
	maxdca = hRECAccdca1D->GetMaximum();
	hRECAccdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	SetHistoProperties(hRECAccdca1D,38,3,38,20,38,3007);
	hRECAccdca1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccdca1D->Sumw2();
	hRECAccdca1D->SetMinimum(0.01);
	hRECAccdca1D->GetXaxis()->SetTitle("dca (#mum), AOD - Acc cuts");
	hRECAccdca1D->Draw("hist");
	cAccdca->cd();
	cAccdca->SetLeftMargin(0.15);
	cAccdca->SetRightMargin(0.05);
	hRECAccdca1D->Draw("hist");
	hRECAccdca1D->Draw("err same");
	cAccdca->Update();

	cmc3->cd(6);
	TH1D *hRECITSCldca1D = data->ShowProjection(idca, stepRecITSClus);
	maxdca = hRECITSCldca1D->GetMaximum();
	hRECITSCldca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	SetHistoProperties(hRECITSCldca1D,13,3,13,20,13,3019);
	hRECITSCldca1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSCldca1D->Sumw2();
	hRECITSCldca1D->SetMinimum(0.01);
	hRECITSCldca1D->GetXaxis()->SetTitle("dca (#mum), AOD - ITS Cl cuts");
	hRECITSCldca1D->Draw("hist");
	cITSCldca->cd();
	cITSCldca->SetLeftMargin(0.15);
	cITSCldca->SetRightMargin(0.05);
	hRECITSCldca1D->Draw("hist");
	hRECITSCldca1D->Draw("err same");
	cITSCldca->Update();

	cmc3->cd(7);
	TH1D *hRECPPRdca1D = data->ShowProjection(idca, stepRecPPR);
	maxdca = hRECPPRdca1D->GetMaximum();
	hRECPPRdca1D->GetYaxis()->SetRangeUser(0,maxdca*1.2);
	SetHistoProperties(hRECPPRdca1D,6,3,6,20,6,3019);
	hRECPPRdca1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRdca1D->Sumw2();
	hRECPPRdca1D->SetMinimum(0.01);
	hRECPPRdca1D->GetXaxis()->SetTitle("dca (#mum), AOD - PPR cuts");
	hRECPPRdca1D->Draw("hist");
	cPPRdca->cd();
	cPPRdca->SetLeftMargin(0.15);
	cPPRdca->SetRightMargin(0.05);
	hRECPPRdca1D->Draw("hist");
	hRECPPRdca1D->Draw("err same");
	cPPRdca->Update();

	cmc3->cd(8);
	TH1D *hMCd0pi1D = data->ShowProjection(id0pi, stepGen);
	Double_t maxd0pi = hMCd0pi1D->GetMaximum();
	hMCd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	hMCd0pi1D->Sumw2();
	hMCd0pi1D->SetMinimum(0.01);
	SetHistoProperties(hMCd0pi1D,2,3,2,20,2,3005);
	hMCd0pi1D->GetXaxis()->SetTitleOffset(1.2);
	hMCd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC data");
	hMCd0pi1D->Draw("hist");
	cmcd0pi->cd();
	cmcd0pi->SetLeftMargin(0.15);
	cmcd0pi->SetRightMargin(0.05);
	hMCd0pi1D->Draw("hist");
	hMCd0pi1D->Draw("err same");
	cmcd0pi->Update();

	cmc3->cd(9);
	TH1D *hMCAccd0pi1D = data->ShowProjection(id0pi, stepGenAcc);
	maxd0pi = hMCAccd0pi1D->GetMaximum();
	hMCAccd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	SetHistoProperties(hMCAccd0pi1D,8,3,8,20,8,3006);
	hMCAccd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccd0pi1D->Sumw2();
	hMCAccd0pi1D->SetMinimum(0.01);
	hMCAccd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Acc cut");
	hMCAccd0pi1D->Draw("hist");
	cmcAccd0pi->cd();
	cmcAccd0pi->SetLeftMargin(0.15);
	cmcAccd0pi->SetRightMargin(0.05);
	hMCAccd0pi1D->Draw("hist");
	hMCAccd0pi1D->Draw("err same");
	cmcAccd0pi->Update();

	cmc3->cd(10);
	TH1D *hMCVertexd0pi1D = data->ShowProjection(id0pi, stepGenVertex);
	maxd0pi = hMCVertexd0pi1D->GetMaximum();
	hMCVertexd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	SetHistoProperties(hMCVertexd0pi1D,49,3,49,20,49,3006);
	hMCVertexd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexd0pi1D->Sumw2();
	hMCVertexd0pi1D->SetMinimum(0.01);
	hMCVertexd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Vertex cut");
	hMCVertexd0pi1D->Draw("hist");
	cmcVertexd0pi->cd();
	cmcVertexd0pi->SetLeftMargin(0.15);
	cmcVertexd0pi->SetRightMargin(0.05);
	hMCVertexd0pi1D->Draw("hist");
	hMCVertexd0pi1D->Draw("err same");
	cmcVertexd0pi->Update();

	cmc3->cd(11);
	TH1D *hMCRefitd0pi1D = data->ShowProjection(id0pi, stepGenRefit);
	maxd0pi = hMCRefitd0pi1D->GetMaximum();
	hMCRefitd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	SetHistoProperties(hMCRefitd0pi1D,42,3,42,20,42,3006);
	hMCRefitd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitd0pi1D->Sumw2();
	hMCRefitd0pi1D->SetMinimum(0.01);
	hMCRefitd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), MC - Refit cut");
	hMCRefitd0pi1D->Draw("hist");
	cmcRefitd0pi->cd();
	cmcRefitd0pi->SetLeftMargin(0.15);
	cmcRefitd0pi->SetRightMargin(0.05);
	hMCRefitd0pi1D->Draw("hist");
	hMCRefitd0pi1D->Draw("err same");
	cmcRefitd0pi->Update();

	cmc3->cd(12);
	TH1D *hRECd0pi1D = data->ShowProjection(id0pi, stepRec);
	maxd0pi = hRECd0pi1D->GetMaximum();
	hRECd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	SetHistoProperties(hRECd0pi1D,4,3,4,20,4,3004);
	hRECd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECd0pi1D->Sumw2();
	hRECd0pi1D->SetMinimum(0.01);
	hRECd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD");
	hRECd0pi1D->Draw("hist");
	cd0pi->cd();
	cd0pi->SetLeftMargin(0.15);
	cd0pi->SetRightMargin(0.05);
	hRECd0pi1D->Draw("hist");
	hRECd0pi1D->Draw("err same");
	cd0pi->Update();

	cmc3->cd(13);
	TH1D *hRECAccd0pi1D = data->ShowProjection(id0pi, stepRecAcc);
	maxd0pi = hRECAccd0pi1D->GetMaximum();
	hRECAccd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	SetHistoProperties(hRECAccd0pi1D,38,3,38,20,38,3007);
	hRECAccd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccd0pi1D->Sumw2();
	hRECAccd0pi1D->SetMinimum(0.01);
	hRECAccd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - Acc cuts");
	hRECAccd0pi1D->Draw("hist");
	cAccd0pi->cd();
	cAccd0pi->SetLeftMargin(0.15);
	cAccd0pi->SetRightMargin(0.05);
	hRECAccd0pi1D->Draw("hist");
	hRECAccd0pi1D->Draw("err same");
	cAccd0pi->Update();

	cmc3->cd(14);
	TH1D *hRECITSCld0pi1D = data->ShowProjection(id0pi, stepRecITSClus);
	maxd0pi = hRECITSCld0pi1D->GetMaximum();
	hRECITSCld0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	SetHistoProperties(hRECITSCld0pi1D,13,3,13,20,13,3019);
	hRECITSCld0pi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSCld0pi1D->Sumw2();
	hRECITSCld0pi1D->SetMinimum(0.01);
	hRECITSCld0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - ITS Cl cuts");
	hRECITSCld0pi1D->Draw("hist");
	cITSCld0pi->cd();
	cITSCld0pi->SetLeftMargin(0.15);
	cITSCld0pi->SetRightMargin(0.05);
	hRECITSCld0pi1D->Draw("hist");
	hRECITSCld0pi1D->Draw("err same");
	cITSCld0pi->Update();

	cmc3->cd(15);
	TH1D *hRECPPRd0pi1D = data->ShowProjection(id0pi, stepRecPPR);
	maxd0pi = hRECPPRd0pi1D->GetMaximum();
	hRECPPRd0pi1D->GetYaxis()->SetRangeUser(0,maxd0pi*1.2);
	SetHistoProperties(hRECPPRd0pi1D,6,3,6,20,6,3019);
	hRECPPRd0pi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRd0pi1D->Sumw2();
	hRECPPRd0pi1D->SetMinimum(0.01);
	hRECPPRd0pi1D->GetXaxis()->SetTitle("d0_{#pi} (#mum), AOD - PPR cuts");
	hRECPPRd0pi1D->Draw("hist");
	cPPRd0pi->cd();
	cPPRd0pi->SetLeftMargin(0.15);
	cPPRd0pi->SetRightMargin(0.05);
	hRECPPRd0pi1D->Draw("hist");
	hRECPPRd0pi1D->Draw("err same");
	cPPRd0pi->Update();

	cmc3->cd(16);
	TH1D *hMCd0K1D = data->ShowProjection(id0K, stepGen);
	Double_t maxd0K = hMCd0K1D->GetMaximum();
	hMCd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	hMCd0K1D->Sumw2();
	hMCd0K1D->SetMinimum(0.01);
	SetHistoProperties(hMCd0K1D,2,3,2,20,2,3005);
	hMCd0K1D->GetXaxis()->SetTitleOffset(1.2);
	hMCd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC data");
	hMCd0K1D->Draw("hist");
	cmcd0K->cd();
	cmcd0K->SetLeftMargin(0.15);
	cmcd0K->SetRightMargin(0.05);
	hMCd0K1D->Draw("hist");
	hMCd0K1D->Draw("err same");
	cmcd0K->Update();

	cmc3->cd(17);
	TH1D *hMCAccd0K1D = data->ShowProjection(id0K, stepGenAcc);
	maxd0K = hMCAccd0K1D->GetMaximum();
	hMCAccd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	SetHistoProperties(hMCAccd0K1D,8,3,8,20,8,3006);
	hMCAccd0K1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccd0K1D->Sumw2();
	hMCAccd0K1D->SetMinimum(0.01);
	hMCAccd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Acc cut");
	hMCAccd0K1D->Draw("hist");
	cmcAccd0K->cd();
	cmcAccd0K->SetLeftMargin(0.15);
	cmcAccd0K->SetRightMargin(0.05);
	hMCAccd0K1D->Draw("hist");
	hMCAccd0K1D->Draw("err same");
	cmcAccd0K->Update();

	cmc3->cd(18);
	TH1D *hMCVertexd0K1D = data->ShowProjection(id0K, stepGenVertex);
	maxd0K = hMCVertexd0K1D->GetMaximum();
	hMCVertexd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	SetHistoProperties(hMCVertexd0K1D,49,3,49,20,49,3006);
	hMCVertexd0K1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexd0K1D->Sumw2();
	hMCVertexd0K1D->SetMinimum(0.01);
	hMCVertexd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Vertex cut");
	hMCVertexd0K1D->Draw("hist");
	cmcVertexd0K->cd();
	cmcVertexd0K->SetLeftMargin(0.15);
	cmcVertexd0K->SetRightMargin(0.05);
	hMCVertexd0K1D->Draw("hist");
	hMCVertexd0K1D->Draw("err same");
	cmcVertexd0K->Update();

	cmc3->cd(20);
	TH1D *hMCRefitd0K1D = data->ShowProjection(id0K, stepGenRefit);
	maxd0K = hMCRefitd0K1D->GetMaximum();
	hMCRefitd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	SetHistoProperties(hMCRefitd0K1D,42,3,42,20,42,3006);
	hMCRefitd0K1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitd0K1D->Sumw2();
	hMCRefitd0K1D->SetMinimum(0.01);
	hMCRefitd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), MC - Refit cut");
	hMCRefitd0K1D->Draw("hist");
	cmcRefitd0K->cd();
	cmcRefitd0K->SetLeftMargin(0.15);
	cmcRefitd0K->SetRightMargin(0.05);
	hMCRefitd0K1D->Draw("hist");
	hMCRefitd0K1D->Draw("err same");
	cmcRefitd0K->Update();

	cmc3->cd(21);
	TH1D *hRECd0K1D = data->ShowProjection(id0K, stepRec);
	maxd0K = hRECd0K1D->GetMaximum();
	hRECd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	SetHistoProperties(hRECd0K1D,4,3,4,20,4,3004);
	hRECd0K1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECd0K1D->Sumw2();
	hRECd0K1D->SetMinimum(0.01);
	hRECd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD");
	hRECd0K1D->Draw("hist");
	cd0K->cd();
	cd0K->SetLeftMargin(0.15);
	cd0K->SetRightMargin(0.05);
	hRECd0K1D->Draw("hist");
	hRECd0K1D->Draw("err same");
	cd0K->Update();

	cmc3->cd(22);
	TH1D *hRECAccd0K1D = data->ShowProjection(id0K, stepRecAcc);
	maxd0K = hRECAccd0K1D->GetMaximum();
	hRECAccd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	SetHistoProperties(hRECAccd0K1D,38,3,38,20,38,3007);
	hRECAccd0K1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccd0K1D->Sumw2();
	hRECAccd0K1D->SetMinimum(0.01);
	hRECAccd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - Acc cuts");
	hRECAccd0K1D->Draw("hist");
	cAccd0K->cd();
	cAccd0K->SetLeftMargin(0.15);
	cAccd0K->SetRightMargin(0.05);
	hRECAccd0K1D->Draw("hist");
	hRECAccd0K1D->Draw("err same");
	cAccd0K->Update();

	cmc3->cd(23);
	TH1D *hRECITSCld0K1D = data->ShowProjection(id0K, stepRecITSClus);
	maxd0K = hRECITSCld0K1D->GetMaximum();
	hRECITSCld0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	SetHistoProperties(hRECITSCld0K1D,13,3,13,20,13,3019);
	hRECITSCld0K1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSCld0K1D->Sumw2();
	hRECITSCld0K1D->SetMinimum(0.01);
	hRECITSCld0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - ITS Cl cuts");
	hRECITSCld0K1D->Draw("hist");
	cITSCld0K->cd();
	cITSCld0K->SetLeftMargin(0.15);
	cITSCld0K->SetRightMargin(0.05);
	hRECITSCld0K1D->Draw("hist");
	hRECITSCld0K1D->Draw("err same");
	cITSCld0K->Update();

	cmc3->cd(24);
	TH1D *hRECPPRd0K1D = data->ShowProjection(id0K, stepRecPPR);
	maxd0K = hRECPPRd0K1D->GetMaximum();
	hRECPPRd0K1D->GetYaxis()->SetRangeUser(0,maxd0K*1.2);
	SetHistoProperties(hRECPPRd0K1D,6,3,6,20,6,3019);
	hRECPPRd0K1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRd0K1D->Sumw2();
	hRECPPRd0K1D->SetMinimum(0.01);
	hRECPPRd0K1D->GetXaxis()->SetTitle("d0_{K} (#mum), AOD - PPR cuts");
	hRECPPRd0K1D->Draw("hist");
	cPPRd0K->cd();
	cPPRd0K->SetLeftMargin(0.15);
	cPPRd0K->SetRightMargin(0.05);
	hRECPPRd0K1D->Draw("hist");
	hRECPPRd0K1D->Draw("err same");
	cPPRd0K->Update();

	//d0xd0, cosPointingAngle, Phi
	cmc4->Divide(8,3);
	cmc4->cd(1);
	TH1D *hMCd0xd01D = data->ShowProjection(id0xd0, stepGen);
	Double_t maxd0xd0 = hMCd0xd01D->GetMaximum();
	hMCd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	hMCd0xd01D->Sumw2();
	hMCd0xd01D->SetMinimum(0.01);
	SetHistoProperties(hMCd0xd01D,2,3,2,20,2,3005);
	hMCd0xd01D->GetXaxis()->SetTitleOffset(1.2);
	hMCd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC data");
	hMCd0xd01D->Draw("hist");
	cmcd0xd0->cd();
	cmcd0xd0->SetLeftMargin(0.15);
	cmcd0xd0->SetRightMargin(0.05);
	hMCd0xd01D->Draw("hist");
	hMCd0xd01D->Draw("err same");
	cmcd0xd0->Update();

	cmc4->cd(2);
	TH1D *hMCAccd0xd01D = data->ShowProjection(id0xd0, stepGenAcc);
	maxd0xd0 = hMCAccd0xd01D->GetMaximum();
	hMCAccd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	SetHistoProperties(hMCAccd0xd01D,8,3,8,20,8,3006);
	hMCAccd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccd0xd01D->Sumw2();
	hMCAccd0xd01D->SetMinimum(0.01);
	hMCAccd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Acc cut");
	hMCAccd0xd01D->Draw("hist");
	cmcAccd0xd0->cd();
	cmcAccd0xd0->SetLeftMargin(0.15);
	cmcAccd0xd0->SetRightMargin(0.05);
	hMCAccd0xd01D->Draw("hist");
	hMCAccd0xd01D->Draw("err same");
	cmcAccd0xd0->Update();

	cmc4->cd(3);
	TH1D *hMCVertexd0xd01D = data->ShowProjection(id0xd0, stepGenVertex);
	maxd0xd0 = hMCVertexd0xd01D->GetMaximum();
	hMCVertexd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	SetHistoProperties(hMCVertexd0xd01D,49,3,49,20,49,3006);
	hMCVertexd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexd0xd01D->Sumw2();
	hMCVertexd0xd01D->SetMinimum(0.01);
	hMCVertexd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Vertex cut");
	hMCVertexd0xd01D->Draw("hist");
	cmcVertexd0xd0->cd();
	cmcVertexd0xd0->SetLeftMargin(0.15);
	cmcVertexd0xd0->SetRightMargin(0.05);
	hMCVertexd0xd01D->Draw("hist");
	hMCVertexd0xd01D->Draw("err same");
	cmcVertexd0xd0->Update();

	cmc4->cd(4);
	TH1D *hMCRefitd0xd01D = data->ShowProjection(id0xd0, stepGenRefit);
	maxd0xd0 = hMCRefitd0xd01D->GetMaximum();
	hMCRefitd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	SetHistoProperties(hMCRefitd0xd01D,42,3,42,20,42,3006);
	hMCRefitd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitd0xd01D->Sumw2();
	hMCRefitd0xd01D->SetMinimum(0.01);
	hMCRefitd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), MC - Refit cut");
	hMCRefitd0xd01D->Draw("hist");
	cmcRefitd0xd0->cd();
	cmcRefitd0xd0->SetLeftMargin(0.15);
	cmcRefitd0xd0->SetRightMargin(0.05);
	hMCRefitd0xd01D->Draw("hist");
	hMCRefitd0xd01D->Draw("err same");
	cmcRefitd0xd0->Update();

	cmc4->cd(5);
	TH1D *hRECd0xd01D = data->ShowProjection(id0xd0, stepRec);
	maxd0xd0 = hRECd0xd01D->GetMaximum();
	hRECd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	SetHistoProperties(hRECd0xd01D,4,3,4,20,4,3004);
	hRECd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
	hRECd0xd01D->Sumw2();
	hRECd0xd01D->SetMinimum(0.01);
	hRECd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD");
	hRECd0xd01D->Draw("hist");
	cd0xd0->cd();
	cd0xd0->SetLeftMargin(0.15);
	cd0xd0->SetRightMargin(0.05);
	hRECd0xd01D->Draw("hist");
	hRECd0xd01D->Draw("err same");
	cd0xd0->Update();

	cmc4->cd(6);
	TH1D *hRECAccd0xd01D = data->ShowProjection(id0xd0, stepRecAcc);
	maxd0xd0 = hRECAccd0xd01D->GetMaximum();
	hRECAccd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	SetHistoProperties(hRECAccd0xd01D,38,3,38,20,38,3007);
	hRECAccd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccd0xd01D->Sumw2();
	hRECAccd0xd01D->SetMinimum(0.01);
	hRECAccd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - Acc cuts");
	hRECAccd0xd01D->Draw("hist");
	cAccd0xd0->cd();
	cAccd0xd0->SetLeftMargin(0.15);
	cAccd0xd0->SetRightMargin(0.05);
	hRECAccd0xd01D->Draw("hist");
	hRECAccd0xd01D->Draw("err same");
	cAccd0xd0->Update();

	cmc4->cd(7);
	TH1D *hRECITSCld0xd01D = data->ShowProjection(id0xd0, stepRecITSClus);
	maxd0xd0 = hRECITSCld0xd01D->GetMaximum();
	hRECITSCld0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	SetHistoProperties(hRECITSCld0xd01D,13,3,13,20,13,3019);
	hRECITSCld0xd01D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSCld0xd01D->Sumw2();
	hRECITSCld0xd01D->SetMinimum(0.01);
	hRECITSCld0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - ITS Cl cuts");
	hRECITSCld0xd01D->Draw("hist");
	cITSCld0xd0->cd();
	cITSCld0xd0->SetLeftMargin(0.15);
	cITSCld0xd0->SetRightMargin(0.05);
	hRECITSCld0xd01D->Draw("hist");
	hRECITSCld0xd01D->Draw("err same");
	cITSCld0xd0->Update();

	cmc4->cd(8);
	TH1D *hRECPPRd0xd01D = data->ShowProjection(id0xd0, stepRecPPR);
	maxd0xd0 = hRECPPRd0xd01D->GetMaximum();
	hRECPPRd0xd01D->GetYaxis()->SetRangeUser(0,maxd0xd0*1.2);
	SetHistoProperties(hRECPPRd0xd01D,6,3,6,20,6,3019);
	hRECPPRd0xd01D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRd0xd01D->Sumw2();
	hRECPPRd0xd01D->SetMinimum(0.01);
	hRECPPRd0xd01D->GetXaxis()->SetTitle("d0_{#pi}xd0_{K} (#mum^2), AOD - PPR cuts");
	hRECPPRd0xd01D->Draw("hist");
	cPPRd0xd0->cd();
	cPPRd0xd0->SetLeftMargin(0.15);
	cPPRd0xd0->SetRightMargin(0.05);
	hRECPPRd0xd01D->Draw("hist");
	hRECPPRd0xd01D->Draw("err same");
	cPPRd0xd0->Update();

	cmc4->cd(9);
	TH1D *hMCpointing1D = data->ShowProjection(ipointing, stepGen);
	Double_t maxpointing = hMCpointing1D->GetMaximum();
	hMCpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	hMCpointing1D->Sumw2();
	hMCpointing1D->SetMinimum(0.01);
	SetHistoProperties(hMCpointing1D,2,3,2,20,2,3005);
	hMCpointing1D->GetXaxis()->SetTitleOffset(1.2);
	hMCpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC data");
	hMCpointing1D->Draw("hist");
	cmcpointing->cd();
	cmcpointing->SetLeftMargin(0.15);
	cmcpointing->SetRightMargin(0.05);
	hMCpointing1D->Draw("hist");
	hMCpointing1D->Draw("err same");
	cmcpointing->Update();

	cmc4->cd(10);
	TH1D *hMCAccpointing1D = data->ShowProjection(ipointing, stepGenAcc);
	maxpointing = hMCAccpointing1D->GetMaximum();
	hMCAccpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	SetHistoProperties(hMCAccpointing1D,8,3,8,20,8,3006);
	hMCAccpointing1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccpointing1D->Sumw2();
	hMCAccpointing1D->SetMinimum(0.01);
	hMCAccpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Acc cut");
	hMCAccpointing1D->Draw("hist");
	cmcAccpointing->cd();
	cmcAccpointing->SetLeftMargin(0.15);
	cmcAccpointing->SetRightMargin(0.05);
	hMCAccpointing1D->Draw("hist");
	hMCAccpointing1D->Draw("err same");
	cmcAccpointing->Update();

	cmc4->cd(11);
	TH1D *hMCVertexpointing1D = data->ShowProjection(ipointing, stepGenVertex);
	maxpointing = hMCVertexpointing1D->GetMaximum();
	hMCVertexpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	SetHistoProperties(hMCVertexpointing1D,49,3,49,20,49,3006);
	hMCVertexpointing1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexpointing1D->Sumw2();
	hMCVertexpointing1D->SetMinimum(0.01);
	hMCVertexpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Vertex cut");
	hMCVertexpointing1D->Draw("hist");
	cmcVertexpointing->cd();
	cmcVertexpointing->SetLeftMargin(0.15);
	cmcVertexpointing->SetRightMargin(0.05);
	hMCVertexpointing1D->Draw("hist");
	hMCVertexpointing1D->Draw("err same");
	cmcVertexpointing->Update();

	cmc4->cd(12);
	TH1D *hMCRefitpointing1D = data->ShowProjection(ipointing, stepGenRefit);
	maxpointing = hMCRefitpointing1D->GetMaximum();
	hMCRefitpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	SetHistoProperties(hMCRefitpointing1D,42,3,42,20,42,3006);
	hMCRefitpointing1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitpointing1D->Sumw2();
	hMCRefitpointing1D->SetMinimum(0.01);
	hMCRefitpointing1D->GetXaxis()->SetTitle("cosPointingAngle, MC - Refit cut");
	hMCRefitpointing1D->Draw("hist");
	cmcRefitpointing->cd();
	cmcRefitpointing->SetLeftMargin(0.15);
	cmcRefitpointing->SetRightMargin(0.05);
	hMCRefitpointing1D->Draw("hist");
	hMCRefitpointing1D->Draw("err same");
	cmcRefitpointing->Update();

	cmc4->cd(13);
	TH1D *hRECpointing1D = data->ShowProjection(ipointing, stepRec);
	maxpointing = hRECpointing1D->GetMaximum();
	hRECpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	SetHistoProperties(hRECpointing1D,4,3,4,20,4,3004);
	hRECpointing1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECpointing1D->Sumw2();
	hRECpointing1D->SetMinimum(0.01);
	hRECpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD");
	hRECpointing1D->Draw("hist");
	cpointing->cd();
	cpointing->SetLeftMargin(0.15);
	cpointing->SetRightMargin(0.05);
	hRECpointing1D->Draw("hist");
	hRECpointing1D->Draw("err same");
	cpointing->Update();

	cmc4->cd(14);
	TH1D *hRECAccpointing1D = data->ShowProjection(ipointing, stepRecAcc);
	maxpointing = hRECAccpointing1D->GetMaximum();
	hRECAccpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	SetHistoProperties(hRECAccpointing1D,38,3,38,20,38,3007);
	hRECAccpointing1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccpointing1D->Sumw2();
	hRECAccpointing1D->SetMinimum(0.01);
	hRECAccpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - Acc cuts");
	hRECAccpointing1D->Draw("hist");
	cAccpointing->cd();
	cAccpointing->SetLeftMargin(0.15);
	cAccpointing->SetRightMargin(0.05);
	hRECAccpointing1D->Draw("hist");
	hRECAccpointing1D->Draw("err same");
	cAccpointing->Update();

	cmc4->cd(15);
	TH1D *hRECITSClpointing1D = data->ShowProjection(ipointing, stepRecITSClus);
	maxpointing = hRECITSClpointing1D->GetMaximum();
	hRECITSClpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	SetHistoProperties(hRECITSClpointing1D,13,3,13,20,13,3019);
	hRECITSClpointing1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSClpointing1D->Sumw2();
	hRECITSClpointing1D->SetMinimum(0.01);
	hRECITSClpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - ITS Cl cuts");
	hRECITSClpointing1D->Draw("hist");
	cITSClpointing->cd();
	cITSClpointing->SetLeftMargin(0.15);
	cITSClpointing->SetRightMargin(0.05);
	hRECITSClpointing1D->Draw("hist");
	hRECITSClpointing1D->Draw("err same");
	cITSClpointing->Update();

	cmc4->cd(16);
	TH1D *hRECPPRpointing1D = data->ShowProjection(ipointing, stepRecPPR);
	maxpointing = hRECPPRpointing1D->GetMaximum();
	hRECPPRpointing1D->GetYaxis()->SetRangeUser(0,maxpointing*1.2);
	SetHistoProperties(hRECPPRpointing1D,6,3,6,20,6,3019);
	hRECPPRpointing1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRpointing1D->Sumw2();
	hRECPPRpointing1D->SetMinimum(0.01);
	hRECPPRpointing1D->GetXaxis()->SetTitle("cosPointingAngle, AOD - PPR cuts");
	hRECPPRpointing1D->Draw("hist");
	cPPRpointing->cd();
	cPPRpointing->SetLeftMargin(0.15);
	cPPRpointing->SetRightMargin(0.05);
	hRECPPRpointing1D->Draw("hist");
	hRECPPRpointing1D->Draw("err same");
	cPPRpointing->Update();

	cmc4->cd(17);
	TH1D *hMCphi1D = data->ShowProjection(iphi, stepGen);
	Double_t maxphi = hMCphi1D->GetMaximum();
	hMCphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	hMCphi1D->Sumw2();
	hMCphi1D->SetMinimum(0.01);
	SetHistoProperties(hMCphi1D,2,3,2,20,2,3005);
	hMCphi1D->GetXaxis()->SetTitleOffset(1.2);
	hMCphi1D->GetXaxis()->SetTitle("Phi, MC data");
	hMCphi1D->Draw("hist");
	cmcphi->cd();
	cmcphi->SetLeftMargin(0.15);
	cmcphi->SetRightMargin(0.05);
	hMCphi1D->Draw("hist");
	hMCphi1D->Draw("err same");
	cmcphi->Update();

	cmc4->cd(18);
	TH1D *hMCAccphi1D = data->ShowProjection(iphi, stepGenAcc);
	maxphi = hMCAccphi1D->GetMaximum();
	hMCAccphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	SetHistoProperties(hMCAccphi1D,8,3,8,20,8,3006);
	hMCAccphi1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCAccphi1D->Sumw2();
	hMCAccphi1D->SetMinimum(0.01);
	hMCAccphi1D->GetXaxis()->SetTitle("Phi, MC - Acc cut");
	hMCAccphi1D->Draw("hist");
	cmcAccphi->cd();
	cmcAccphi->SetLeftMargin(0.15);
	cmcAccphi->SetRightMargin(0.05);
	hMCAccphi1D->Draw("hist");
	hMCAccphi1D->Draw("err same");
	cmcAccphi->Update();

	cmc4->cd(19);
	TH1D *hMCVertexphi1D = data->ShowProjection(iphi, stepGenVertex);
	maxphi = hMCVertexphi1D->GetMaximum();
	hMCVertexphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	SetHistoProperties(hMCVertexphi1D,49,3,49,20,49,3006);
	hMCVertexphi1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCVertexphi1D->Sumw2();
	hMCVertexphi1D->SetMinimum(0.01);
	hMCVertexphi1D->GetXaxis()->SetTitle("Phi, MC - Vertex cut");
	hMCVertexphi1D->Draw("hist");
	cmcVertexphi->cd();
	cmcVertexphi->SetLeftMargin(0.15);
	cmcVertexphi->SetRightMargin(0.05);
	hMCVertexphi1D->Draw("hist");
	hMCVertexphi1D->Draw("err same");
	cmcVertexphi->Update();

	cmc4->cd(20);
	TH1D *hMCRefitphi1D = data->ShowProjection(iphi, stepGenRefit);
	maxphi = hMCRefitphi1D->GetMaximum();
	hMCRefitphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	SetHistoProperties(hMCRefitphi1D,42,3,42,20,42,3006);
	hMCRefitphi1D->GetXaxis()->SetTitleOffset(1.2);;
	hMCRefitphi1D->Sumw2();
	hMCRefitphi1D->SetMinimum(0.01);
	hMCRefitphi1D->GetXaxis()->SetTitle("Phi, MC - Refit cut");
	hMCRefitphi1D->Draw("hist");
	cmcRefitphi->cd();
	cmcRefitphi->SetLeftMargin(0.15);
	cmcRefitphi->SetRightMargin(0.05);
	hMCRefitphi1D->Draw("hist");
	hMCRefitphi1D->Draw("err same");
	cmcRefitphi->Update();

	cmc4->cd(21);
	TH1D *hRECphi1D = data->ShowProjection(iphi, stepRec);
	maxphi = hRECphi1D->GetMaximum();
	hRECphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	SetHistoProperties(hRECphi1D,4,3,4,20,4,3004);
	hRECphi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECphi1D->Sumw2();
	hRECphi1D->SetMinimum(0.01);
	hRECphi1D->GetXaxis()->SetTitle("Phi, AOD");
	hRECphi1D->Draw("hist");
	cphi->cd();
	cphi->SetLeftMargin(0.15);
	cphi->SetRightMargin(0.05);
	hRECphi1D->Draw("hist");
	hRECphi1D->Draw("err same");
	cphi->Update();

	cmc4->cd(22);
	TH1D *hRECAccphi1D = data->ShowProjection(iphi, stepRecAcc);
	maxphi = hRECAccphi1D->GetMaximum();
	hRECAccphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	SetHistoProperties(hRECAccphi1D,38,3,38,20,38,3007);
	hRECAccphi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECAccphi1D->Sumw2();
	hRECAccphi1D->SetMinimum(0.01);
	hRECAccphi1D->GetXaxis()->SetTitle("Phi, AOD - Acc cuts");
	hRECAccphi1D->Draw("hist");
	cAccphi->cd();
	cAccphi->SetLeftMargin(0.15);
	cAccphi->SetRightMargin(0.05);
	hRECAccphi1D->Draw("hist");
	hRECAccphi1D->Draw("err same");
	cAccphi->Update();

	cmc4->cd(23);
	TH1D *hRECITSClphi1D = data->ShowProjection(iphi, stepRecITSClus);
	maxphi = hRECITSClphi1D->GetMaximum();
	hRECITSClphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	SetHistoProperties(hRECITSClphi1D,13,3,13,20,13,3019);
	hRECITSClphi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECITSClphi1D->Sumw2();
	hRECITSClphi1D->SetMinimum(0.01);
	hRECITSClphi1D->GetXaxis()->SetTitle("Phi, AOD - ITS Cl cuts");
	hRECITSClphi1D->Draw("hist");
	cITSClphi->cd();
	cITSClphi->SetLeftMargin(0.15);
	cITSClphi->SetRightMargin(0.05);
	hRECITSClphi1D->Draw("hist");
	hRECITSClphi1D->Draw("err same");
	cITSClphi->Update();

	cmc4->cd(24);
	TH1D *hRECPPRphi1D = data->ShowProjection(iphi, stepRecPPR);
	maxphi = hRECPPRphi1D->GetMaximum();
	hRECPPRphi1D->GetYaxis()->SetRangeUser(0,maxphi*1.2);
	SetHistoProperties(hRECPPRphi1D,6,3,6,20,6,3019);
	hRECPPRphi1D->GetXaxis()->SetTitleOffset(1.2);;
	hRECPPRphi1D->Sumw2();
	hRECPPRphi1D->SetMinimum(0.01);
	hRECPPRphi1D->GetXaxis()->SetTitle("Phi, AOD - PPR cuts");
	hRECPPRphi1D->Draw("hist");
	cPPRphi->cd();
	cPPRphi->SetLeftMargin(0.15);
	cPPRphi->SetRightMargin(0.05);
	hRECPPRphi1D->Draw("hist");
	hRECPPRphi1D->Draw("err same");
	cPPRphi->Update();

	// printing on eps files
	cmc1->Print("Plots/dataMC_pt_y_cTS.eps");
	cmc2->Print("Plots/dataMC_ptPi_ptK_cT.eps");
	cmc3->Print("Plots/dataMC_dca_d0pi_d0K.eps");
	cmc4->Print("Plots/dataMC_d0xd0_pointing.eps");
	cmcpt->Print("Plots/pt_Gen.eps");
	cmcy->Print("Plots/y_Gen.eps");
	cmccTS->Print("Plots/cTS_Gen.eps");
	cmcptPi->Print("Plots/ptPi_Gen.eps");
	cmcptK->Print("Plots/ptK_Gen.eps");
	cmccT->Print("Plots/cT_Gen.eps");
	cmcdca->Print("Plots/dca_Gen.eps");
	cmcd0pi->Print("Plots/d0pi_Gen.eps");
	cmcd0K->Print("Plots/d0K_Gen.eps");
	cmcd0xd0->Print("Plots/d0xd0_Gen.eps");
	cmcpointing->Print("Plots/pointing_Gen.eps");
	cmcphi->Print("Plots/phi_Gen.eps");

	cmcAccpt->Print("Plots/pt_Gen_Acc.eps");
	cmcAccy->Print("Plots/y_Gen_Acc.eps");
	cmcAcccTS->Print("Plots/cTS_Gen_Acc.eps");
	cmcAccptPi->Print("Plots/ptPi_Gen_Acc.eps");
	cmcAccptK->Print("Plots/ptK_Gen_Acc.eps");
	cmcAcccT->Print("Plots/cT_Gen_Acc.eps");
	cmcAccdca->Print("Plots/dca_Gen_Acc.eps");
	cmcAccd0pi->Print("Plots/d0pi_Gen_Acc.eps");
	cmcAccd0K->Print("Plots/d0K_Gen_Acc.eps");
	cmcAccd0xd0->Print("Plots/d0xd0_Gen_Acc.eps");
	cmcAccpointing->Print("Plots/pointing_Gen_Acc.eps");
	cmcAccphi->Print("Plots/phi_Gen_Acc.eps");

	cmcVertexpt->Print("Plots/pt_Gen_Vertex.eps");
	cmcVertexy->Print("Plots/y_Gen_Vertex.eps");
	cmcVertexcTS->Print("Plots/cTS_Gen_Vertex.eps");
	cmcVertexptPi->Print("Plots/ptPi_Gen_Vertex.eps");
	cmcVertexptK->Print("Plots/ptK_Gen_Vertex.eps");
	cmcVertexcT->Print("Plots/cT_Gen_Vertex.eps");
	cmcVertexdca->Print("Plots/dca_Gen_Vertex.eps");
	cmcVertexd0pi->Print("Plots/d0pi_Gen_Vertex.eps");
	cmcVertexd0K->Print("Plots/d0K_Gen_Vertex.eps");
	cmcVertexd0xd0->Print("Plots/d0xd0_Gen_Vertex.eps");
	cmcVertexpointing->Print("Plots/pointing_Gen_Vertex.eps");
	cmcVertexphi->Print("Plots/phi_Gen_Vertex.eps");

	cmcRefitpt->Print("Plots/pt_Gen_Refit.eps");
	cmcRefity->Print("Plots/y_Gen_Refit.eps");
	cmcRefitcTS->Print("Plots/cTS_Gen_Refit.eps");
	cmcRefitptPi->Print("Plots/ptPi_Gen_Refit.eps");
	cmcRefitptK->Print("Plots/ptK_Gen_Refit.eps");
	cmcRefitcT->Print("Plots/cT_Gen_Refit.eps");
	cmcRefitdca->Print("Plots/dca_Gen_Refit.eps");
	cmcRefitd0pi->Print("Plots/d0pi_Gen_Refit.eps");
	cmcRefitd0K->Print("Plots/d0K_Gen_Refit.eps");
	cmcRefitd0xd0->Print("Plots/d0xd0_Gen_Refit.eps");
	cmcRefitpointing->Print("Plots/pointing_Gen_Refit.eps");
	cmcRefitphi->Print("Plots/phi_Gen_Refit.eps");

	cpt->Print("Plots/pt_Rec.eps");
	cy->Print("Plots/y_Rec.eps");
	ccTS->Print("Plots/cTS_Rec.eps");
	cptPi->Print("Plots/ptPi_Rec.eps");
	cptK->Print("Plots/ptK_Rec.eps");
	ccT->Print("Plots/cT_Rec.eps");
	cdca->Print("Plots/dca_Rec.eps");
	cd0pi->Print("Plots/d0pi_Rec.eps");
	cd0K->Print("Plots/d0K_Rec.eps");
	cd0xd0->Print("Plots/d0xd0_Rec.eps");
	cpointing->Print("Plots/pointing_Rec.eps");
	cphi->Print("Plots/phi_Rec.eps");

	cAccpt->Print("Plots/pt_Rec_Acc.eps");
	cAccy->Print("Plots/y_Rec_Acc.eps");
	cAcccTS->Print("Plots/cTS_Rec_Acc.eps");
	cAccptPi->Print("Plots/ptPi_Rec_Acc.eps");
	cAccptK->Print("Plots/ptK_Rec_Acc.eps");
	cAcccT->Print("Plots/cT_Rec_Acc.eps");
	cmcdca->Print("Plots/dca_Rec_Acc.eps");
	cAccd0pi->Print("Plots/d0pi_Rec_Acc.eps");
	cAccd0K->Print("Plots/d0K_Rec_Acc.eps");
	cAccd0xd0->Print("Plots/d0xd0_Rec_Acc.eps");
	cAccpointing->Print("Plots/pointing_Rec_Acc.eps");
	cAccphi->Print("Plots/phi_Rec_Acc.eps");

	cITSClpt->Print("Plots/pt_Rec_ITSCl.eps");
	cITSCly->Print("Plots/y_Rec_ITSCl.eps");
	cITSClcTS->Print("Plots/cTS_Rec_ITSCl.eps");
	cITSClptPi->Print("Plots/ptPi_Rec_ITSCl.eps");
	cITSClptK->Print("Plots/ptK_Rec_ITSCl.eps");
	cITSClcT->Print("Plots/cT_Rec_ITSCl.eps");
	cITSCldca->Print("Plots/dca_Rec_ITSCl.eps");
	cITSCld0pi->Print("Plots/d0pi_Rec_ITSCl.eps");
	cITSCld0K->Print("Plots/d0K_Rec_ITSCl.eps");
	cITSCld0xd0->Print("Plots/d0xd0_Rec_ITSCl.eps");
	cITSClpointing->Print("Plots/pointing_Rec_ITSCl.eps");
	cITSClphi->Print("Plots/phi_Rec_ITSCl.eps");

	cPPRpt->Print("Plots/pt_Rec_PPR.eps");
	cPPRy->Print("Plots/y_Rec_PPR.eps");
	cPPRcTS->Print("Plots/cTS_Rec_PPR.eps");
	cPPRptPi->Print("Plots/ptPi_Rec_PPR.eps");
	cPPRptK->Print("Plots/ptK_Rec_PPR.eps");
	cPPRcT->Print("Plots/cT_Rec_PPR.eps");
	cPPRdca->Print("Plots/dca_Rec_PPR.eps");
	cPPRd0pi->Print("Plots/d0pi_Rec_PPR.eps");
	cPPRd0K->Print("Plots/d0K_Rec_PPR.eps");
	cPPRd0xd0->Print("Plots/d0xd0_Rec_PPR.eps");
	cPPRpointing->Print("Plots/pointing_Rec_PPR.eps");
	cPPRphi->Print("Plots/phi_Rec_PPR.eps");

	// printing on gif files
	cmc1->Print("Plots/dataMC_pt_y_cTS.gif");
	cmc2->Print("Plots/dataMC_ptPi_ptK_cT.gif");
	cmc3->Print("Plots/dataMC_dca_d0pi_d0K.gif");
	cmc4->Print("Plots/dataMC_d0xd0_pointing.gif");
	cmcpt->Print("Plots/pt_Gen.gif");
	cmcy->Print("Plots/y_Gen.gif");
	cmccTS->Print("Plots/cTS_Gen.gif");
	cmcptPi->Print("Plots/ptPi_Gen.gif");
	cmcptK->Print("Plots/ptK_Gen.gif");
	cmccT->Print("Plots/cT_Gen.gif");
	cmcdca->Print("Plots/dca_Gen.gif");
	cmcd0pi->Print("Plots/d0pi_Gen.gif");
	cmcd0K->Print("Plots/d0K_Gen.gif");
	cmcd0xd0->Print("Plots/d0xd0_Gen.gif");
	cmcpointing->Print("Plots/pointing_Gen.gif");
	cmcphi->Print("Plots/phi_Gen.gif");

	cmcAccpt->Print("Plots/pt_Gen_Acc.gif");
	cmcAccy->Print("Plots/y_Gen_Acc.gif");
	cmcAcccTS->Print("Plots/cTS_Gen_Acc.gif");
	cmcAccptPi->Print("Plots/ptPi_Gen_Acc.gif");
	cmcAccptK->Print("Plots/ptK_Gen_Acc.gif");
	cmcAcccT->Print("Plots/cT_Gen_Acc.gif");
	cmcAccdca->Print("Plots/dca_Gen_Acc.gif");
	cmcAccd0pi->Print("Plots/d0pi_Gen_Acc.gif");
	cmcAccd0K->Print("Plots/d0K_Gen_Acc.gif");
	cmcAccd0xd0->Print("Plots/d0xd0_Gen_Acc.gif");
	cmcAccpointing->Print("Plots/pointing_Gen_Acc.gif");
	cmcAccphi->Print("Plots/phi_Gen_Acc.gif");

	cmcVertexpt->Print("Plots/pt_Gen_Vertex.gif");
	cmcVertexy->Print("Plots/y_Gen_Vertex.gif");
	cmcVertexcTS->Print("Plots/cTS_Gen_Vertex.gif");
	cmcVertexptPi->Print("Plots/ptPi_Gen_Vertex.gif");
	cmcVertexptK->Print("Plots/ptK_Gen_Vertex.gif");
	cmcVertexcT->Print("Plots/cT_Gen_Vertex.gif");
	cmcVertexdca->Print("Plots/dca_Gen_Vertex.gif");
	cmcVertexd0pi->Print("Plots/d0pi_Gen_Vertex.gif");
	cmcVertexd0K->Print("Plots/d0K_Gen_Vertex.gif");
	cmcVertexd0xd0->Print("Plots/d0xd0_Gen_Vertex.gif");
	cmcVertexpointing->Print("Plots/pointing_Gen_Vertex.gif");
	cmcVertexphi->Print("Plots/phi_Gen_Vertex.gif");

	cmcRefitpt->Print("Plots/pt_Gen_Refit.gif");
	cmcRefity->Print("Plots/y_Gen_Refit.gif");
	cmcRefitcTS->Print("Plots/cTS_Gen_Refit.gif");
	cmcRefitptPi->Print("Plots/ptPi_Gen_Refit.gif");
	cmcRefitptK->Print("Plots/ptK_Gen_Refit.gif");
	cmcRefitcT->Print("Plots/cT_Gen_Refit.gif");
	cmcRefitdca->Print("Plots/dca_Gen_Refit.gif");
	cmcRefitd0pi->Print("Plots/d0pi_Gen_Refit.gif");
	cmcRefitd0K->Print("Plots/d0K_Gen_Refit.gif");
	cmcRefitd0xd0->Print("Plots/d0xd0_Gen_Refit.gif");
	cmcRefitpointing->Print("Plots/pointing_Gen_Refit.gif");
	cmcRefitphi->Print("Plots/phi_Gen_Refit.gif");

	cpt->Print("Plots/pt_Rec.gif");
	cy->Print("Plots/y_Rec.gif");
	ccTS->Print("Plots/cTS_Rec.gif");
	cptPi->Print("Plots/ptPi_Rec.gif");
	cptK->Print("Plots/ptK_Rec.gif");
	ccT->Print("Plots/cT_Rec.gif");
	cdca->Print("Plots/dca_Rec.gif");
	cd0pi->Print("Plots/d0pi_Rec.gif");
	cd0K->Print("Plots/d0K_Rec.gif");
	cd0xd0->Print("Plots/d0xd0_Rec.gif");
	cpointing->Print("Plots/pointing_Rec.gif");
	cphi->Print("Plots/phi_Rec.gif");

	cAccpt->Print("Plots/pt_Rec_Acc.gif");
	cAccy->Print("Plots/y_Rec_Acc.gif");
	cAcccTS->Print("Plots/cTS_Rec_Acc.gif");
	cAccptPi->Print("Plots/ptPi_Rec_Acc.gif");
	cAccptK->Print("Plots/ptK_Rec_Acc.gif");
	cAcccT->Print("Plots/cT_Rec_Acc.gif");
	cmcdca->Print("Plots/dca_Rec_Acc.gif");
	cAccd0pi->Print("Plots/d0pi_Rec_Acc.gif");
	cAccd0K->Print("Plots/d0K_Rec_Acc.gif");
	cAccd0xd0->Print("Plots/d0xd0_Rec_Acc.gif");
	cAccpointing->Print("Plots/pointing_Rec_Acc.gif");
	cAccphi->Print("Plots/phi_Rec_Acc.gif");

	cITSClpt->Print("Plots/pt_Rec_ITSCl.gif");
	cITSCly->Print("Plots/y_Rec_ITSCl.gif");
	cITSClcTS->Print("Plots/cTS_Rec_ITSCl.gif");
	cITSClptPi->Print("Plots/ptPi_Rec_ITSCl.gif");
	cITSClptK->Print("Plots/ptK_Rec_ITSCl.gif");
	cITSClcT->Print("Plots/cT_Rec_ITSCl.gif");
	cITSCldca->Print("Plots/dca_Rec_ITSCl.gif");
	cITSCld0pi->Print("Plots/d0pi_Rec_ITSCl.gif");
	cITSCld0K->Print("Plots/d0K_Rec_ITSCl.gif");
	cITSCld0xd0->Print("Plots/d0xd0_Rec_ITSCl.gif");
	cITSClpointing->Print("Plots/pointing_Rec_ITSCl.gif");
	cITSClphi->Print("Plots/phi_Rec_ITSCl.gif");

	cPPRpt->Print("Plots/pt_Rec_PPR.gif");
	cPPRy->Print("Plots/y_Rec_PPR.gif");
	cPPRcTS->Print("Plots/cTS_Rec_PPR.gif");
	cPPRptPi->Print("Plots/ptPi_Rec_PPR.gif");
	cPPRptK->Print("Plots/ptK_Rec_PPR.gif");
	cPPRcT->Print("Plots/cT_Rec_PPR.gif");
	cPPRdca->Print("Plots/dca_Rec_PPR.gif");
	cPPRd0pi->Print("Plots/d0pi_Rec_PPR.gif");
	cPPRd0K->Print("Plots/d0K_Rec_PPR.gif");
	cPPRd0xd0->Print("Plots/d0xd0_Rec_PPR.gif");
	cPPRpointing->Print("Plots/pointing_Rec_PPR.gif");
	cPPRphi->Print("Plots/phi_Rec_PPR.gif");

	TFile* file_histo = new TFile("fileHisto.root","RECREATE");
	// MC
	hMCpt1D->Write("hMCpt1D");
	hMCy1D->Write("hMCy1D");
	hMCcTS1D->Write("hMCcTS1D");
	hMCptPi1D->Write("hMCptPi1D");
	hMCptK1D->Write("hMCptK1D");
	hMCcT1D->Write("hMCcT1D");
	hMCdca1D->Write("hMCdca1D");
	hMCd0pi1D->Write("hMCd0pi1D");
	hMCd0K1D->Write("hMCd0K1D");
	hMCd0xd01D->Write("hMCd0xd01D");
	hMCpointing1D->Write("hMCpointing1D");
	hMCphi1D->Write("hMCphi1D");
	// MC, Acc
	hMCAccpt1D->Write("hMCAccpt1D");
	hMCAccy1D->Write("hMCAccy1D");
	hMCAcccTS1D->Write("hMCAcccTS1D");
	hMCAccptPi1D->Write("hMCAccptPi1D");
	hMCAccptK1D->Write("hMCAccptK1D");
	hMCAcccT1D->Write("hMCAcccT1D");
	hMCAccdca1D->Write("hMCAccdca1D");
	hMCAccd0pi1D->Write("hMCAccd0pi1D");
	hMCAccd0K1D->Write("hMCAccd0K1D");
	hMCAccd0xd01D->Write("hMCAccd0xd01D");
	hMCAccpointing1D->Write("hMCAccpointing1D");
	hMCAccphi1D->Write("hMCAccphi1D");
	// REC
	hRECpt1D->Write("hRECpt1D");
	hRECy1D->Write("hRECy1D");
	hRECcTS1D->Write("hRECcTS1D");
	hRECptPi1D->Write("hRECptPi1D");
	hRECptK1D->Write("hRECptK1D");
	hRECcT1D->Write("hRECcT1D");
	hRECdca1D->Write("hRECdca1D");
	hRECd0pi1D->Write("hRECd0pi1D");
	hRECd0K1D->Write("hRECd0K1D");
	hRECd0xd01D->Write("hRECd0xd01D");
	hRECpointing1D->Write("hRECpointing1D");
	hRECphi1D->Write("hRECphi1D");
	// REC, Acc
	hRECAccpt1D->Write("hRECAccpt1D");
	hRECAccy1D->Write("hRECAccy1D");
	hRECAcccTS1D->Write("hRECAcccTS1D");
	hRECAccptPi1D->Write("hRECAccptPi1D");
	hRECAccptK1D->Write("hRECAccptK1D");
	hRECAcccT1D->Write("hRECAcccT1D");
	hRECAccdca1D->Write("hRECAccdca1D");
	hRECAccd0pi1D->Write("hRECAccd0pi1D");
	hRECAccd0K1D->Write("hRECAccd0K1D");
	hRECAccd0xd01D->Write("hRECAccd0xd01D");
	hRECAccpointing1D->Write("hRECAccpointing1D");
	hRECAccphi1D->Write("hRECAccphi1D");
	// REC, ITS Cl
	hRECITSClpt1D->Write("hRECITSClpt1D");
	hRECITSCly1D->Write("hRECITSCly1D");
	hRECITSClcTS1D->Write("hRECITSClcTS1D");
	hRECITSClptPi1D->Write("hRECITSClptPi1D");
	hRECITSClptK1D->Write("hRECITSClptK1D");
	hRECITSClcT1D->Write("hRECITSClcT1D");
	hRECITSCldca1D->Write("hRECITSCldca1D");
	hRECITSCld0pi1D->Write("hRECITSCld0pi1D");
	hRECITSCld0K1D->Write("hRECITSCld0K1D");
	hRECITSCld0xd01D->Write("hRECITSCld0xd01D");
	hRECITSClpointing1D->Write("hRECITSClpointing1D");
	hRECITSClphi1D->Write("hRECITSClphi1D");
	// REC, PPR
	hRECPPRpt1D->Write("hRECPPRpt1D");
	hRECPPRy1D->Write("hRECPPRy1D");
	hRECPPRcTS1D->Write("hRECPPRcTS1D");
	hRECPPRptPi1D->Write("hRECPPRptPi1D");
	hRECPPRptK1D->Write("hRECPPRptK1D");
	hRECPPRcT1D->Write("hRECPPRcT1D");
	hRECPPRdca1D->Write("hRECPPRdca1D");
	hRECPPRd0pi1D->Write("hRECPPRd0pi1D");
	hRECPPRd0K1D->Write("hRECPPRd0K1D");
	hRECPPRd0xd01D->Write("hRECPPRd0xd01D");
	hRECPPRpointing1D->Write("hRECPPRpointing1D");
	hRECPPRphi1D->Write("hRECPPRphi1D");

}

void SetHistoProperties(TH1D* h, Int_t lineColor, Int_t lineWidth, Int_t markerColor, Int_t markerStyle, Int_t fillColor, Int_t fillStyle){

	h->SetLineColor(lineColor);
	h->SetLineWidth(lineWidth);
	h->SetMarkerColor(markerColor);
	h->SetFillColor(fillColor);
	h->SetFillStyle(fillStyle);
	h->SetMarkerStyle(markerStyle);

	return;

}