[u/mrichter/AliRoot.git] / PWG2 / SPECTRA / macros / drawProtonResults.C

void drawProtonResults(const char* esdFileName) {
  //Macro to visualize the proton ratio results
  //It also visualizes the QA plots
  drawQAPlots(esdFileName);
}

//___________________________________________________//
void drawQAPlots(const char* esdFileName) {
  //Draws the QA plots from the output of the analysis
  //=========================================================//
  //List of cuts
  TFile *fCutFile = TFile::Open("ListOfCuts.root");
  TCanvas *cListOfCuts = dynamic_cast<TCanvas *>(fCutFile->Get("cListOfCuts"));
  cListOfCuts->Draw();

  //=========================================================//
  //QA plots
  TFile *f = TFile::Open(esdFileName);
  TList *listQA = dynamic_cast<TList *>(f->Get("outputQAList"));
  TList *gListGlobalQA = dynamic_cast<TList *>(listQA->At(0));

  //================QA plots================//
  TList *fQA2DList = dynamic_cast<TList *>(gListGlobalQA->At(0));

  //Accepted protons
  TList *fQAProtonsAcceptedList = dynamic_cast<TList *>(gListGlobalQA->At(1));
  TH1F *gProtonsITSClustersPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(0));
  TH1F *gProtonsChi2PerClusterITSPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(1));
  TH1F *gProtonsTPCClustersPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(2));
  TH1F *gProtonsChi2PerClusterTPCPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(3));
  TH1F *gProtonsExtCov11Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(4));
  TH1F *gProtonsExtCov22Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(5));
  TH1F *gProtonsExtCov33Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(6));
  TH1F *gProtonsExtCov44Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(7));
  TH1F *gProtonsExtCov55Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(8));
  TH1F *gProtonsSigmaToVertexPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(9));
  TH1F *gProtonsSigmaToVertexTPCPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(10));
  TH1F *gProtonsDCAXYPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(11));
  TH1F *gProtonsDCAXYTPCPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(12));
  TH1F *gProtonsDCAZPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(13));
  TH1F *gProtonsDCAZTPCPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(14));
  TH1F *gProtonsConstrainChi2Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(15));
  TH1F *gProtonsITSRefitPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(16));
  TH1F *gProtonsTPCRefitPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(17));
  TH1F *gProtonsESDpidPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(18));
  TH1F *gProtonsTPCpidPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(19));
  TH1F *gProtonsPointOnITSLayer1Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(20));
  TH1F *gProtonsPointOnITSLayer2Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(21));
  TH1F *gProtonsPointOnITSLayer3Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(22));
  TH1F *gProtonsPointOnITSLayer4Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(23));
  TH1F *gProtonsPointOnITSLayer5Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(24));
  TH1F *gProtonsPointOnITSLayer6Pass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(25));
  TH1F *gProtonsNumberOfTPCdEdxPointsPass = dynamic_cast<TH1F *>(fQAProtonsAcceptedList->At(26));

  //Rejected protons
  TList *fQAProtonsRejectedList = dynamic_cast<TList *>(gListGlobalQA->At(2));
  TH1F *gProtonsITSClustersReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(0));
  TH1F *gProtonsChi2PerClusterITSReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(1));
  TH1F *gProtonsTPCClustersReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(2));
  TH1F *gProtonsChi2PerClusterTPCReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(3));
  TH1F *gProtonsExtCov11Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(4));
  TH1F *gProtonsExtCov22Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(5));
  TH1F *gProtonsExtCov33Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(6));
  TH1F *gProtonsExtCov44Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(7));
  TH1F *gProtonsExtCov55Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(8));
  TH1F *gProtonsSigmaToVertexReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(9));
  TH1F *gProtonsSigmaToVertexTPCReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(10));
  TH1F *gProtonsDCAXYReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(11));
  TH1F *gProtonsDCAXYTPCReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(12));
  TH1F *gProtonsDCAZReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(13));
  TH1F *gProtonsDCAZTPCReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(14));
  TH1F *gProtonsConstrainChi2Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(15));
  TH1F *gProtonsITSRefitReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(16));
  TH1F *gProtonsTPCRefitReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(17));
  TH1F *gProtonsESDpidReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(18));
  TH1F *gProtonsTPCpidReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(19));
  TH1F *gProtonsPointOnITSLayer1Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(20));
  TH1F *gProtonsPointOnITSLayer2Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(21));
  TH1F *gProtonsPointOnITSLayer3Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(22));
  TH1F *gProtonsPointOnITSLayer4Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(23));
  TH1F *gProtonsPointOnITSLayer5Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(24));
  TH1F *gProtonsPointOnITSLayer6Reject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(25));
  TH1F *gProtonsNumberOfTPCdEdxPointsReject = dynamic_cast<TH1F *>(fQAProtonsRejectedList->At(26));

  //Accepted antiprotons
  TList *fQAAntiProtonsAcceptedList = dynamic_cast<TList *>(gListGlobalQA->At(3));
  TH1F *gAntiProtonsITSClustersPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(0));
  TH1F *gAntiProtonsChi2PerClusterITSPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(1));
  TH1F *gAntiProtonsTPCClustersPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(2));
  TH1F *gAntiProtonsChi2PerClusterTPCPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(3));
  TH1F *gAntiProtonsExtCov11Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(4));
  TH1F *gAntiProtonsExtCov22Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(5));
  TH1F *gAntiProtonsExtCov33Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(6));
  TH1F *gAntiProtonsExtCov44Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(7));
  TH1F *gAntiProtonsExtCov55Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(8));
  TH1F *gAntiProtonsSigmaToVertexPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(9));
  TH1F *gAntiProtonsSigmaToVertexTPCPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(10));
  TH1F *gAntiProtonsDCAXYPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(11));
  TH1F *gAntiProtonsDCAXYTPCPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(12));
  TH1F *gAntiProtonsDCAZPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(13));
  TH1F *gAntiProtonsDCAZTPCPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(14));
  TH1F *gAntiProtonsConstrainChi2Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(15));
  TH1F *gAntiProtonsITSRefitPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(16));
  TH1F *gAntiProtonsTPCRefitPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(17));
  TH1F *gAntiProtonsESDpidPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(18));
  TH1F *gAntiProtonsTPCpidPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(19));
  TH1F *gAntiProtonsPointOnITSLayer1Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(20));
  TH1F *gAntiProtonsPointOnITSLayer2Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(21));
  TH1F *gAntiProtonsPointOnITSLayer3Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(22));
  TH1F *gAntiProtonsPointOnITSLayer4Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(23));
  TH1F *gAntiProtonsPointOnITSLayer5Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(24));
  TH1F *gAntiProtonsPointOnITSLayer6Pass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(25));
  TH1F *gAntiProtonsNumberOfTPCdEdxPointsPass = dynamic_cast<TH1F *>(fQAAntiProtonsAcceptedList->At(26));

  //Rejected antiprotons
  TList *fQAAntiProtonsRejectedList = dynamic_cast<TList *>(gListGlobalQA->At(4));
  TH1F *gAntiProtonsITSClustersReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(0));
  TH1F *gAntiProtonsChi2PerClusterITSReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(1));
  TH1F *gAntiProtonsTPCClustersReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(2));
  TH1F *gAntiProtonsChi2PerClusterTPCReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(3));
  TH1F *gAntiProtonsExtCov11Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(4));
  TH1F *gAntiProtonsExtCov22Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(5));
  TH1F *gAntiProtonsExtCov33Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(6));
  TH1F *gAntiProtonsExtCov44Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(7));
  TH1F *gAntiProtonsExtCov55Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(8));
  TH1F *gAntiProtonsSigmaToVertexReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(9));
  TH1F *gAntiProtonsSigmaToVertexTPCReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(10));
  TH1F *gAntiProtonsDCAXYReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(11));
  TH1F *gAntiProtonsDCAXYTPCReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(12));
  TH1F *gAntiProtonsDCAZReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(13));
  TH1F *gAntiProtonsDCAZTPCReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(14));
  TH1F *gAntiProtonsConstrainChi2Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(15));
  TH1F *gAntiProtonsITSRefitReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(16));
  TH1F *gAntiProtonsTPCRefitReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(17));
  TH1F *gAntiProtonsESDpidReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(18));
  TH1F *gAntiProtonsTPCpidReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(19));
  TH1F *gAntiProtonsPointOnITSLayer1Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(20));
  TH1F *gAntiProtonsPointOnITSLayer2Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(21));
  TH1F *gAntiProtonsPointOnITSLayer3Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(22));
  TH1F *gAntiProtonsPointOnITSLayer4Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(23));
  TH1F *gAntiProtonsPointOnITSLayer5Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(24));
  TH1F *gAntiProtonsPointOnITSLayer6Reject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(25));
  TH1F *gAntiProtonsNumberOfTPCdEdxPointsReject = dynamic_cast<TH1F *>(fQAAntiProtonsRejectedList->At(26));

  //__________________________________________________//
  TCanvas *c1 = new TCanvas("c1","",0,0,600,400);
  c1->SetFillColor(10); c1->SetHighLightColor(10);
  c1->Divide(2,1);
  c1->cd(1); gProtonsITSClustersPass->Draw(); 
  gProtonsITSClustersReject->Draw("same");
  c1->cd(2); gAntiProtonsITSClustersPass->Draw(); 
  gAntiProtonsITSClustersReject->Draw("same");

  TCanvas *c2 = new TCanvas("c2","",0,100,600,400);
  c2->SetFillColor(10); c2->SetHighLightColor(10);
  c2->Divide(2,1);
  c2->cd(1); gProtonsChi2PerClusterITSPass->Draw(); 
  gProtonsChi2PerClusterITSReject->Draw("same");
  c2->cd(2); gAntiProtonsChi2PerClusterITSPass->Draw(); 
  gAntiProtonsChi2PerClusterITSReject->Draw("same");

  TCanvas *c3 = new TCanvas("c3","",0,200,600,400);
  c3->SetFillColor(10); c3->SetHighLightColor(10);
  c3->Divide(2,1);
  c3->cd(1); gProtonsTPCClustersPass->Draw();
  gProtonsTPCClustersReject->Draw("same");
  c3->cd(2); gAntiProtonsTPCClustersPass->Draw();
  gAntiProtonsTPCClustersReject->Draw("same");

  TCanvas *c4 = new TCanvas("c4","",0,300,600,400);
  c4->SetFillColor(10); c4->SetHighLightColor(10);
  c4->Divide(2,1);
  c4->cd(1); gProtonsChi2PerClusterTPCPass->Draw(); 
  gProtonsChi2PerClusterTPCReject->Draw("same");
  c4->cd(2); gAntiProtonsChi2PerClusterTPCPass->Draw(); 
  gAntiProtonsChi2PerClusterTPCReject->Draw("same");

  TCanvas *c5 = new TCanvas("c5","",0,400,600,400);
  c5->SetFillColor(10); c5->SetHighLightColor(10);
  c5->Divide(2,1);
  c5->cd(1); gProtonsExtCov11Pass->Draw(); 
  gProtonsExtCov11Reject->Draw("same");
  c5->cd(2); gAntiProtonsExtCov11Pass->Draw(); 
  gAntiProtonsExtCov11Reject->Draw("same");

  TCanvas *c6 = new TCanvas("c6","",0,500,600,400);
  c6->SetFillColor(10); c6->SetHighLightColor(10);
  c6->Divide(2,1);
  c6->cd(1); gProtonsExtCov22Pass->Draw(); 
  gProtonsExtCov22Reject->Draw("same");
  c6->cd(2); gAntiProtonsExtCov22Pass->Draw(); 
  gAntiProtonsExtCov22Reject->Draw("same");

  TCanvas *c7 = new TCanvas("c7","",600,0,600,400);
  c7->SetFillColor(10); c7->SetHighLightColor(10);
  c7->Divide(2,1);
  c7->cd(1); gProtonsExtCov33Pass->Draw(); 
  gProtonsExtCov33Reject->Draw("same");
  c7->cd(2); gAntiProtonsExtCov33Pass->Draw(); 
  gAntiProtonsExtCov33Reject->Draw("same");

  TCanvas *c8 = new TCanvas("c8","",600,100,600,400);
  c8->SetFillColor(10); c8->SetHighLightColor(10);
  c8->Divide(2,1);
  c8->cd(1); gProtonsExtCov44Pass->Draw(); 
  gProtonsExtCov44Reject->Draw("same");
  c8->cd(2); gAntiProtonsExtCov44Pass->Draw(); 
  gAntiProtonsExtCov44Reject->Draw("same");

  TCanvas *c9 = new TCanvas("c9","",600,200,600,400);
  c9->SetFillColor(10); c9->SetHighLightColor(10);
  c9->Divide(2,1);
  c9->cd(1); gProtonsExtCov55Pass->Draw(); 
  gProtonsExtCov55Reject->Draw("same");
  c9->cd(2); gAntiProtonsExtCov55Pass->Draw(); 
  gAntiProtonsExtCov55Reject->Draw("same");

  TCanvas *c10 = new TCanvas("c10","",600,300,600,400);
  c10->SetFillColor(10); c10->SetHighLightColor(10);
  c10->Divide(2,1);
  c10->cd(1); gProtonsSigmaToVertexPass->Draw(); 
  gProtonsSigmaToVertexReject->Draw("same");
  c10->cd(2); gAntiProtonsSigmaToVertexPass->Draw(); 
  gAntiProtonsSigmaToVertexReject->Draw("same");

  TCanvas *c11 = new TCanvas("c11","",600,400,600,400);
  c11->SetFillColor(10); c11->SetHighLightColor(10);
  c11->Divide(2,1);
  c11->cd(1); gProtonsSigmaToVertexTPCPass->Draw(); 
  gProtonsSigmaToVertexTPCReject->Draw("same");
  c11->cd(2); gAntiProtonsSigmaToVertexTPCPass->Draw(); 
  gAntiProtonsSigmaToVertexTPCReject->Draw("same");

  TCanvas *c12 = new TCanvas("c12","",600,500,600,400);
  c12->SetFillColor(10); c12->SetHighLightColor(10);
  c12->Divide(2,1);
  c12->cd(1); gProtonsDCAXYPass->Draw(); 
  gProtonsDCAXYReject->Draw("same");
  c12->cd(2); gAntiProtonsDCAXYPass->Draw(); 
  gAntiProtonsDCAXYReject->Draw("same");

  TCanvas *c13 = new TCanvas("c13","",1200,0,600,400);
  c13->SetFillColor(10); c13->SetHighLightColor(10);
  c13->Divide(2,1);
  c13->cd(1); gProtonsDCAXYTPCPass->Draw(); 
  gProtonsDCAXYTPCReject->Draw("same");
  c13->cd(2); gAntiProtonsDCAXYTPCPass->Draw(); 
  gAntiProtonsDCAXYTPCReject->Draw("same");

  TCanvas *c14 = new TCanvas("c14","",1200,100,600,400);
  c14->SetFillColor(10); c14->SetHighLightColor(10);
  c14->Divide(2,1);
  c14->cd(1); gProtonsDCAZPass->Draw(); 
  gProtonsDCAZReject->Draw("same");
  c14->cd(2); gAntiProtonsDCAZPass->Draw(); 
  gAntiProtonsDCAZReject->Draw("same");

  TCanvas *c15 = new TCanvas("c15","",1200,200,600,400);
  c15->SetFillColor(10); c15->SetHighLightColor(10);
  c15->Divide(2,1);
  c15->cd(1); gProtonsDCAZTPCPass->Draw(); 
  gProtonsDCAZTPCReject->Draw("same");
  c15->cd(2); gAntiProtonsDCAZTPCPass->Draw(); 
  gAntiProtonsDCAZTPCReject->Draw("same");

  TCanvas *c16 = new TCanvas("c16","",1200,300,600,400);
  c16->SetFillColor(10); c16->SetHighLightColor(10);
  c16->Divide(2,1);
  c16->cd(1); gProtonsNumberOfTPCdEdxPointsPass->Draw(); 
  gProtonsNumberOfTPCdEdxPointsReject->Draw("same");
  c16->cd(2); gAntiProtonsNumberOfTPCdEdxPointsPass->Draw(); 
  gAntiProtonsNumberOfTPCdEdxPointsReject->Draw("same");

  //================Vertex QA================//
  TList *gListVertexQA = dynamic_cast<TList *>(listQA->At(1));
  TH1F *gHistVx = dynamic_cast<TH1F *>(gListVertexQA->At(0));
  TH1F *gHistVxAccepted = dynamic_cast<TH1F *>(gListVertexQA->At(1));
  gHistVxAccepted->SetFillColor(10);
  TH1F *gHistVy = dynamic_cast<TH1F *>(gListVertexQA->At(2));
  TH1F *gHistVyAccepted = dynamic_cast<TH1F *>(gListVertexQA->At(3));
  gHistVyAccepted->SetFillColor(10);
  TH1F *gHistVz = dynamic_cast<TH1F *>(gListVertexQA->At(4));
  TH1F *gHistVzAccepted = dynamic_cast<TH1F *>(gListVertexQA->At(5));
  gHistVzAccepted->SetFillColor(10);

  TCanvas *cVertex = new TCanvas("cVertex","Veretx QA",0,0,900,400);
  cVertex->SetFillColor(10); cVertex->SetHighLightColor(10);
  cVertex->Divide(3,1);
  cVertex->cd(1); gHistVx->Draw(); gHistVxAccepted->Draw("same");
  cVertex->cd(2); gHistVy->Draw(); gHistVyAccepted->Draw("same");
  cVertex->cd(3); gHistVz->Draw(); gHistVzAccepted->Draw("same");
}
Commit	Line	Data
735cc63d	1	void drawProtonResults(const char* esdFileName) {
	2	//Macro to visualize the proton ratio results
	3	//It also visualizes the QA plots
	4	drawQAPlots(esdFileName);
	5	}
	6
	7	//___________________________________________________//
	8	void drawQAPlots(const char* esdFileName) {
	9	//Draws the QA plots from the output of the analysis
	10	//=========================================================//
	11	//List of cuts
	12	TFile *fCutFile = TFile::Open("ListOfCuts.root");
	13	TCanvas cListOfCuts = dynamic_cast<TCanvas >(fCutFile->Get("cListOfCuts"));
	14	cListOfCuts->Draw();
	15
	16	//=========================================================//
	17	//QA plots
	18	TFile *f = TFile::Open(esdFileName);
	19	TList listQA = dynamic_cast<TList >(f->Get("outputQAList"));
	20	TList gListGlobalQA = dynamic_cast<TList >(listQA->At(0));
	21
	22	//================QA plots================//
	23	TList fQA2DList = dynamic_cast<TList >(gListGlobalQA->At(0));
	24
	25	//Accepted protons
	26	TList fQAProtonsAcceptedList = dynamic_cast<TList >(gListGlobalQA->At(1));
	27	TH1F gProtonsITSClustersPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(0));
	28	TH1F gProtonsChi2PerClusterITSPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(1));
	29	TH1F gProtonsTPCClustersPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(2));
	30	TH1F gProtonsChi2PerClusterTPCPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(3));
	31	TH1F gProtonsExtCov11Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(4));
	32	TH1F gProtonsExtCov22Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(5));
	33	TH1F gProtonsExtCov33Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(6));
	34	TH1F gProtonsExtCov44Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(7));
	35	TH1F gProtonsExtCov55Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(8));
	36	TH1F gProtonsSigmaToVertexPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(9));
	37	TH1F gProtonsSigmaToVertexTPCPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(10));
	38	TH1F gProtonsDCAXYPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(11));
	39	TH1F gProtonsDCAXYTPCPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(12));
	40	TH1F gProtonsDCAZPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(13));
	41	TH1F gProtonsDCAZTPCPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(14));
	42	TH1F gProtonsConstrainChi2Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(15));
	43	TH1F gProtonsITSRefitPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(16));
	44	TH1F gProtonsTPCRefitPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(17));
	45	TH1F gProtonsESDpidPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(18));
	46	TH1F gProtonsTPCpidPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(19));
	47	TH1F gProtonsPointOnITSLayer1Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(20));
	48	TH1F gProtonsPointOnITSLayer2Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(21));
	49	TH1F gProtonsPointOnITSLayer3Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(22));
	50	TH1F gProtonsPointOnITSLayer4Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(23));
	51	TH1F gProtonsPointOnITSLayer5Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(24));
	52	TH1F gProtonsPointOnITSLayer6Pass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(25));
	53	TH1F gProtonsNumberOfTPCdEdxPointsPass = dynamic_cast<TH1F >(fQAProtonsAcceptedList->At(26));
	54
	55	//Rejected protons
	56	TList fQAProtonsRejectedList = dynamic_cast<TList >(gListGlobalQA->At(2));
	57	TH1F gProtonsITSClustersReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(0));
	58	TH1F gProtonsChi2PerClusterITSReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(1));
	59	TH1F gProtonsTPCClustersReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(2));
	60	TH1F gProtonsChi2PerClusterTPCReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(3));
	61	TH1F gProtonsExtCov11Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(4));
	62	TH1F gProtonsExtCov22Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(5));
	63	TH1F gProtonsExtCov33Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(6));
	64	TH1F gProtonsExtCov44Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(7));
65	TH1F gProtonsExtCov55Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(8));
66	TH1F gProtonsSigmaToVertexReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(9));
67	TH1F gProtonsSigmaToVertexTPCReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(10));
68	TH1F gProtonsDCAXYReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(11));
69	TH1F gProtonsDCAXYTPCReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(12));
70	TH1F gProtonsDCAZReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(13));
71	TH1F gProtonsDCAZTPCReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(14));
72	TH1F gProtonsConstrainChi2Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(15));
73	TH1F gProtonsITSRefitReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(16));
74	TH1F gProtonsTPCRefitReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(17));
75	TH1F gProtonsESDpidReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(18));
76	TH1F gProtonsTPCpidReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(19));
77	TH1F gProtonsPointOnITSLayer1Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(20));
78	TH1F gProtonsPointOnITSLayer2Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(21));
79	TH1F gProtonsPointOnITSLayer3Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(22));
80	TH1F gProtonsPointOnITSLayer4Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(23));
81	TH1F gProtonsPointOnITSLayer5Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(24));
82	TH1F gProtonsPointOnITSLayer6Reject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(25));
83	TH1F gProtonsNumberOfTPCdEdxPointsReject = dynamic_cast<TH1F >(fQAProtonsRejectedList->At(26));
84
85	//Accepted antiprotons
86	TList fQAAntiProtonsAcceptedList = dynamic_cast<TList >(gListGlobalQA->At(3));
87	TH1F gAntiProtonsITSClustersPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(0));
88	TH1F gAntiProtonsChi2PerClusterITSPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(1));
89	TH1F gAntiProtonsTPCClustersPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(2));
90	TH1F gAntiProtonsChi2PerClusterTPCPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(3));
91	TH1F gAntiProtonsExtCov11Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(4));
92	TH1F gAntiProtonsExtCov22Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(5));
93	TH1F gAntiProtonsExtCov33Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(6));
94	TH1F gAntiProtonsExtCov44Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(7));
95	TH1F gAntiProtonsExtCov55Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(8));
96	TH1F gAntiProtonsSigmaToVertexPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(9));
97	TH1F gAntiProtonsSigmaToVertexTPCPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(10));
98	TH1F gAntiProtonsDCAXYPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(11));
99	TH1F gAntiProtonsDCAXYTPCPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(12));
100	TH1F gAntiProtonsDCAZPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(13));
101	TH1F gAntiProtonsDCAZTPCPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(14));
102	TH1F gAntiProtonsConstrainChi2Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(15));
103	TH1F gAntiProtonsITSRefitPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(16));
104	TH1F gAntiProtonsTPCRefitPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(17));
105	TH1F gAntiProtonsESDpidPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(18));
106	TH1F gAntiProtonsTPCpidPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(19));
107	TH1F gAntiProtonsPointOnITSLayer1Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(20));
108	TH1F gAntiProtonsPointOnITSLayer2Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(21));
109	TH1F gAntiProtonsPointOnITSLayer3Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(22));
110	TH1F gAntiProtonsPointOnITSLayer4Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(23));
111	TH1F gAntiProtonsPointOnITSLayer5Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(24));
112	TH1F gAntiProtonsPointOnITSLayer6Pass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(25));
113	TH1F gAntiProtonsNumberOfTPCdEdxPointsPass = dynamic_cast<TH1F >(fQAAntiProtonsAcceptedList->At(26));
114
115	//Rejected antiprotons
116	TList fQAAntiProtonsRejectedList = dynamic_cast<TList >(gListGlobalQA->At(4));
117	TH1F gAntiProtonsITSClustersReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(0));
118	TH1F gAntiProtonsChi2PerClusterITSReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(1));
119	TH1F gAntiProtonsTPCClustersReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(2));
120	TH1F gAntiProtonsChi2PerClusterTPCReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(3));
121	TH1F gAntiProtonsExtCov11Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(4));
122	TH1F gAntiProtonsExtCov22Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(5));
123	TH1F gAntiProtonsExtCov33Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(6));
124	TH1F gAntiProtonsExtCov44Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(7));
125	TH1F gAntiProtonsExtCov55Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(8));
126	TH1F gAntiProtonsSigmaToVertexReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(9));
127	TH1F gAntiProtonsSigmaToVertexTPCReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(10));
128	TH1F gAntiProtonsDCAXYReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(11));
129	TH1F gAntiProtonsDCAXYTPCReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(12));
130	TH1F gAntiProtonsDCAZReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(13));
131	TH1F gAntiProtonsDCAZTPCReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(14));
132	TH1F gAntiProtonsConstrainChi2Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(15));
133	TH1F gAntiProtonsITSRefitReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(16));
134	TH1F gAntiProtonsTPCRefitReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(17));
135	TH1F gAntiProtonsESDpidReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(18));
136	TH1F gAntiProtonsTPCpidReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(19));
137	TH1F gAntiProtonsPointOnITSLayer1Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(20));
138	TH1F gAntiProtonsPointOnITSLayer2Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(21));
139	TH1F gAntiProtonsPointOnITSLayer3Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(22));
140	TH1F gAntiProtonsPointOnITSLayer4Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(23));
141	TH1F gAntiProtonsPointOnITSLayer5Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(24));
142	TH1F gAntiProtonsPointOnITSLayer6Reject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(25));
143	TH1F gAntiProtonsNumberOfTPCdEdxPointsReject = dynamic_cast<TH1F >(fQAAntiProtonsRejectedList->At(26));
144
145	//__________________________________________________//
146	TCanvas *c1 = new TCanvas("c1","",0,0,600,400);
147	c1->SetFillColor(10); c1->SetHighLightColor(10);
148	c1->Divide(2,1);
149	c1->cd(1); gProtonsITSClustersPass->Draw();
150	gProtonsITSClustersReject->Draw("same");
151	c1->cd(2); gAntiProtonsITSClustersPass->Draw();
152	gAntiProtonsITSClustersReject->Draw("same");
153
154	TCanvas *c2 = new TCanvas("c2","",0,100,600,400);
155	c2->SetFillColor(10); c2->SetHighLightColor(10);
156	c2->Divide(2,1);
157	c2->cd(1); gProtonsChi2PerClusterITSPass->Draw();
158	gProtonsChi2PerClusterITSReject->Draw("same");
159	c2->cd(2); gAntiProtonsChi2PerClusterITSPass->Draw();
160	gAntiProtonsChi2PerClusterITSReject->Draw("same");
161
162	TCanvas *c3 = new TCanvas("c3","",0,200,600,400);
163	c3->SetFillColor(10); c3->SetHighLightColor(10);
164	c3->Divide(2,1);
165	c3->cd(1); gProtonsTPCClustersPass->Draw();
166	gProtonsTPCClustersReject->Draw("same");
167	c3->cd(2); gAntiProtonsTPCClustersPass->Draw();
168	gAntiProtonsTPCClustersReject->Draw("same");
169
170	TCanvas *c4 = new TCanvas("c4","",0,300,600,400);
171	c4->SetFillColor(10); c4->SetHighLightColor(10);
172	c4->Divide(2,1);
173	c4->cd(1); gProtonsChi2PerClusterTPCPass->Draw();
174	gProtonsChi2PerClusterTPCReject->Draw("same");
175	c4->cd(2); gAntiProtonsChi2PerClusterTPCPass->Draw();
176	gAntiProtonsChi2PerClusterTPCReject->Draw("same");
177
178	TCanvas *c5 = new TCanvas("c5","",0,400,600,400);
179	c5->SetFillColor(10); c5->SetHighLightColor(10);
180	c5->Divide(2,1);
181	c5->cd(1); gProtonsExtCov11Pass->Draw();
182	gProtonsExtCov11Reject->Draw("same");
183	c5->cd(2); gAntiProtonsExtCov11Pass->Draw();
184	gAntiProtonsExtCov11Reject->Draw("same");
185
186	TCanvas *c6 = new TCanvas("c6","",0,500,600,400);
187	c6->SetFillColor(10); c6->SetHighLightColor(10);
188	c6->Divide(2,1);
189	c6->cd(1); gProtonsExtCov22Pass->Draw();
190	gProtonsExtCov22Reject->Draw("same");
191	c6->cd(2); gAntiProtonsExtCov22Pass->Draw();
192	gAntiProtonsExtCov22Reject->Draw("same");
193
194	TCanvas *c7 = new TCanvas("c7","",600,0,600,400);
195	c7->SetFillColor(10); c7->SetHighLightColor(10);
196	c7->Divide(2,1);
197	c7->cd(1); gProtonsExtCov33Pass->Draw();
198	gProtonsExtCov33Reject->Draw("same");
199	c7->cd(2); gAntiProtonsExtCov33Pass->Draw();
200	gAntiProtonsExtCov33Reject->Draw("same");
201
202	TCanvas *c8 = new TCanvas("c8","",600,100,600,400);
203	c8->SetFillColor(10); c8->SetHighLightColor(10);
204	c8->Divide(2,1);
205	c8->cd(1); gProtonsExtCov44Pass->Draw();
206	gProtonsExtCov44Reject->Draw("same");
207	c8->cd(2); gAntiProtonsExtCov44Pass->Draw();
208	gAntiProtonsExtCov44Reject->Draw("same");
209
210	TCanvas *c9 = new TCanvas("c9","",600,200,600,400);
211	c9->SetFillColor(10); c9->SetHighLightColor(10);
212	c9->Divide(2,1);
213	c9->cd(1); gProtonsExtCov55Pass->Draw();
214	gProtonsExtCov55Reject->Draw("same");
215	c9->cd(2); gAntiProtonsExtCov55Pass->Draw();
216	gAntiProtonsExtCov55Reject->Draw("same");
217
218	TCanvas *c10 = new TCanvas("c10","",600,300,600,400);
219	c10->SetFillColor(10); c10->SetHighLightColor(10);
220	c10->Divide(2,1);
221	c10->cd(1); gProtonsSigmaToVertexPass->Draw();
222	gProtonsSigmaToVertexReject->Draw("same");
223	c10->cd(2); gAntiProtonsSigmaToVertexPass->Draw();
224	gAntiProtonsSigmaToVertexReject->Draw("same");
225
226	TCanvas *c11 = new TCanvas("c11","",600,400,600,400);
227	c11->SetFillColor(10); c11->SetHighLightColor(10);
228	c11->Divide(2,1);
229	c11->cd(1); gProtonsSigmaToVertexTPCPass->Draw();
230	gProtonsSigmaToVertexTPCReject->Draw("same");
231	c11->cd(2); gAntiProtonsSigmaToVertexTPCPass->Draw();
232	gAntiProtonsSigmaToVertexTPCReject->Draw("same");
233
234	TCanvas *c12 = new TCanvas("c12","",600,500,600,400);
235	c12->SetFillColor(10); c12->SetHighLightColor(10);
236	c12->Divide(2,1);
237	c12->cd(1); gProtonsDCAXYPass->Draw();
238	gProtonsDCAXYReject->Draw("same");
239	c12->cd(2); gAntiProtonsDCAXYPass->Draw();
240	gAntiProtonsDCAXYReject->Draw("same");
241
242	TCanvas *c13 = new TCanvas("c13","",1200,0,600,400);
243	c13->SetFillColor(10); c13->SetHighLightColor(10);
244	c13->Divide(2,1);
245	c13->cd(1); gProtonsDCAXYTPCPass->Draw();
246	gProtonsDCAXYTPCReject->Draw("same");
247	c13->cd(2); gAntiProtonsDCAXYTPCPass->Draw();
248	gAntiProtonsDCAXYTPCReject->Draw("same");
249
250	TCanvas *c14 = new TCanvas("c14","",1200,100,600,400);
251	c14->SetFillColor(10); c14->SetHighLightColor(10);
252	c14->Divide(2,1);
253	c14->cd(1); gProtonsDCAZPass->Draw();
254	gProtonsDCAZReject->Draw("same");
255	c14->cd(2); gAntiProtonsDCAZPass->Draw();
256	gAntiProtonsDCAZReject->Draw("same");
257
258	TCanvas *c15 = new TCanvas("c15","",1200,200,600,400);
259	c15->SetFillColor(10); c15->SetHighLightColor(10);
260	c15->Divide(2,1);
261	c15->cd(1); gProtonsDCAZTPCPass->Draw();
262	gProtonsDCAZTPCReject->Draw("same");
263	c15->cd(2); gAntiProtonsDCAZTPCPass->Draw();
264	gAntiProtonsDCAZTPCReject->Draw("same");
265
266	TCanvas *c16 = new TCanvas("c16","",1200,300,600,400);
267	c16->SetFillColor(10); c16->SetHighLightColor(10);
268	c16->Divide(2,1);
269	c16->cd(1); gProtonsNumberOfTPCdEdxPointsPass->Draw();
270	gProtonsNumberOfTPCdEdxPointsReject->Draw("same");
271	c16->cd(2); gAntiProtonsNumberOfTPCdEdxPointsPass->Draw();
272	gAntiProtonsNumberOfTPCdEdxPointsReject->Draw("same");
273
274	//================Vertex QA================//
275	TList gListVertexQA = dynamic_cast<TList >(listQA->At(1));
276	TH1F gHistVx = dynamic_cast<TH1F >(gListVertexQA->At(0));
277	TH1F gHistVxAccepted = dynamic_cast<TH1F >(gListVertexQA->At(1));
278	gHistVxAccepted->SetFillColor(10);
279	TH1F gHistVy = dynamic_cast<TH1F >(gListVertexQA->At(2));
280	TH1F gHistVyAccepted = dynamic_cast<TH1F >(gListVertexQA->At(3));
281	gHistVyAccepted->SetFillColor(10);
282	TH1F gHistVz = dynamic_cast<TH1F >(gListVertexQA->At(4));
283	TH1F gHistVzAccepted = dynamic_cast<TH1F >(gListVertexQA->At(5));
284	gHistVzAccepted->SetFillColor(10);
285
286	TCanvas *cVertex = new TCanvas("cVertex","Veretx QA",0,0,900,400);
287	cVertex->SetFillColor(10); cVertex->SetHighLightColor(10);
288	cVertex->Divide(3,1);
289	cVertex->cd(1); gHistVx->Draw(); gHistVxAccepted->Draw("same");
290	cVertex->cd(2); gHistVy->Draw(); gHistVyAccepted->Draw("same");
291	cVertex->cd(3); gHistVz->Draw(); gHistVzAccepted->Draw("same");
292	}
293