[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / macros / PlotOutputMCCheck.C

#if !defined(__CINT__) || defined(__MAKECINT__)
#include <TCanvas.h>
#include <TH1.h>
#include <TH2.h>
#include <TStyle.h>
#include <TFile.h>
#endif

/* $Id$ */ 

// Macro to plot the output of AliAnalysisTaskCheckHFMCProd
// Author: F. Prino, prino@to.infn.it


void PlotOutputMCCheck(){
  TFile *fil=new TFile("AnalysisResultsMerged.root");
  TDirectoryFile* df=(TDirectoryFile*)fil->Get("HFMCCheck");
  TList* l=(TList*)df->Get("clistHFMCCheck");
  l->ls();

  TH1F* hNEvents=(TH1F*)l->FindObject("hNEvents");
  Int_t nAnalEv=hNEvents->GetBinContent(1);
  printf("Number of events= %d\n",nAnalEv);


  TCanvas* cv=new TCanvas("cv","Vertex");
  cv->Divide(3,3);
  cv->cd(1);
  TH1F* hSPD3DvX=(TH1F*)l->FindObject("hSPD3DvX");
  hSPD3DvX->Draw();
  cv->cd(2);
  TH1F* hSPD3DvY=(TH1F*)l->FindObject("hSPD3DvY");
  hSPD3DvY->Draw();
  cv->cd(3);
  TH1F* hSPD3DvZ=(TH1F*)l->FindObject("hSPD3DvZ");
  hSPD3DvZ->Draw();
  cv->cd(4);
  TH1F* hSPDZvX=(TH1F*)l->FindObject("hSPDZvX");
  hSPDZvX->Draw();
  cv->cd(5);
  TH1F* hSPDZvY=(TH1F*)l->FindObject("hSPDZvY");
  hSPDZvY->Draw();
  cv->cd(6);
  TH1F* hSPDZvZ=(TH1F*)l->FindObject("hSPDZvZ");
  hSPDZvZ->Draw();
  cv->cd(7);
  TH1F* hTRKvX=(TH1F*)l->FindObject("hTRKvX");
  hTRKvX->Draw();
  cv->cd(8);
  TH1F* hTRKvY=(TH1F*)l->FindObject("hTRKvY");
  hTRKvY->Draw();
  cv->cd(9);
  TH1F* hTRKvZ=(TH1F*)l->FindObject("hTRKvZ");
  hTRKvZ->Draw();

  TCanvas* c1=new TCanvas("c1","Multipl");
  c1->Divide(3,1);
  c1->cd(1);
  TH1F* hTracklets=(TH1F*)l->FindObject("hTracklets");
  hTracklets->Draw();
  c1->cd(2);
  TH1F* hTracks=(TH1F*)l->FindObject("hTracks");
  hTracks->Draw();
  c1->cd(3);
  TH1F* hSelTracks=(TH1F*)l->FindObject("hSelTracks");
  hSelTracks->Draw();

  TH1F* hncharmed=(TH1F*)l->FindObject("hncharmed");
  TCanvas* cn=new TCanvas("cn","ncharm");
  hncharmed->Draw("box");
  hncharmed->GetXaxis()->SetTitle("dNch/dy");
  hncharmed->GetYaxis()->SetTitle("N Charm hadrons in golden channels");
  cn->Update();

  TH1F* hnbvsnc=(TH1F*)l->FindObject("hnbvsnc");
  TCanvas* cnhf=new TCanvas("cnhf","nb/c");
  hnbvsnc->Draw("box");
  hnbvsnc->GetXaxis()->SetTitle("Nc");
  hnbvsnc->GetYaxis()->SetTitle("Nb");
  cnhf->Update();

  TH2F* hyptd0prompt=(TH2F*)l->FindObject("hyptd0prompt");
  TH2F* hyptd0feeddown=(TH2F*)l->FindObject("hyptd0feeddown");
  TH2F* hyptD02=(TH2F*)l->FindObject("hyptD02");
  TH2F* hyptD04=(TH2F*)l->FindObject("hyptD04");

  TCanvas* cd0=new TCanvas("cd0","D0");
  cd0->Divide(2,2);
  cd0->cd(1);
  hyptd0prompt->Draw("colz");
  cd0->cd(2);
  hyptd0feeddown->Draw("colz");
  cd0->cd(3);
  hyptD02->Draw("colz");
  cd0->cd(4);
  hyptD04->Draw("colz");

  TH2F* hyptdplusprompt=(TH2F*)l->FindObject("hyptdplusprompt");
  TH2F* hyptdplusfeeddown=(TH2F*)l->FindObject("hyptdplusfeeddown");
  TH2F* hyptDplusnonreson=(TH2F*)l->FindObject("hyptDplusnonreson");
  TH2F* hyptDplusreson=(TH2F*)l->FindObject("hyptDplusreson");

  TCanvas* cdplus=new TCanvas("cdplus","Dplus");
  cdplus->Divide(2,2);
  cdplus->cd(1);
  hyptdplusprompt->Draw("colz");
  cdplus->cd(2);
  hyptdplusfeeddown->Draw("colz");
  cdplus->cd(3);
  hyptDplusnonreson->Draw("colz");
  cdplus->cd(4);
  hyptDplusreson->Draw("colz");

  TH2F* hyptdsprompt=(TH2F*)l->FindObject("hyptdsprompt");
  TH2F* hyptdsfeeddown=(TH2F*)l->FindObject("hyptdsfeedown");
  TH2F* hyptdsphi=(TH2F*)l->FindObject("hyptdsphi");
  TH2F* hyptdsK0st=(TH2F*)l->FindObject("hyptdsk0st");
  
  TCanvas* cds=new TCanvas("cds","Ds");
  cds->Divide(2,2);
  cds->cd(1);
  hyptdsprompt->Draw("colz");
  cds->cd(2);
  hyptdsfeeddown->Draw("colz");
  cds->cd(3);
  hyptdsphi->Draw("colz");
  cds->cd(4);
  hyptdsK0st->Draw("colz");

  TH2F* hyptdstarprompt=(TH2F*)l->FindObject("hyptdstarprompt");
  TH2F* hyptdstarfeedown=(TH2F*)l->FindObject("hyptdstarfeedown");
  TH2F* hyptlcprompt=(TH2F*)l->FindObject("hyptlcprompt");
  TH2F* hyptlcfeedown=(TH2F*)l->FindObject("hyptlcfeedown");

  TCanvas* cdstlc=new TCanvas("cdstls","Dstar LambdaC");
  cdstlc->Divide(2,2);
  cdstlc->cd(1);
  hyptdstarprompt->Draw("colz");
  cdstlc->cd(2);
  hyptdstarfeedown->Draw("colz");
  cdstlc->cd(3);
  hyptlcprompt->Draw("colz");
  cdstlc->cd(4);
  hyptlcfeedown->Draw("colz");

}
Commit	Line	Data
7b6a4dcd	1	#if !defined(__CINT__) \|\| defined(__MAKECINT__)
	2	#include <TCanvas.h>
	3	#include <TH1.h>
	4	#include <TH2.h>
	5	#include <TStyle.h>
	6	#include <TFile.h>
	7	#endif
	8
	9	/* $Id$ */
	10
	11	// Macro to plot the output of AliAnalysisTaskCheckHFMCProd
	12	// Author: F. Prino, prino@to.infn.it
	13
	14
	15	void PlotOutputMCCheck(){
	16	TFile *fil=new TFile("AnalysisResultsMerged.root");
	17	TDirectoryFile* df=(TDirectoryFile*)fil->Get("HFMCCheck");
	18	TList* l=(TList*)df->Get("clistHFMCCheck");
	19	l->ls();
	20
	21	TH1F* hNEvents=(TH1F*)l->FindObject("hNEvents");
	22	Int_t nAnalEv=hNEvents->GetBinContent(1);
	23	printf("Number of events= %d\n",nAnalEv);
	24
	25
	26	TCanvas* cv=new TCanvas("cv","Vertex");
	27	cv->Divide(3,3);
	28	cv->cd(1);
	29	TH1F* hSPD3DvX=(TH1F*)l->FindObject("hSPD3DvX");
	30	hSPD3DvX->Draw();
	31	cv->cd(2);
	32	TH1F* hSPD3DvY=(TH1F*)l->FindObject("hSPD3DvY");
	33	hSPD3DvY->Draw();
	34	cv->cd(3);
	35	TH1F* hSPD3DvZ=(TH1F*)l->FindObject("hSPD3DvZ");
	36	hSPD3DvZ->Draw();
	37	cv->cd(4);
	38	TH1F* hSPDZvX=(TH1F*)l->FindObject("hSPDZvX");
	39	hSPDZvX->Draw();
	40	cv->cd(5);
	41	TH1F* hSPDZvY=(TH1F*)l->FindObject("hSPDZvY");
	42	hSPDZvY->Draw();
	43	cv->cd(6);
	44	TH1F* hSPDZvZ=(TH1F*)l->FindObject("hSPDZvZ");
	45	hSPDZvZ->Draw();
	46	cv->cd(7);
	47	TH1F* hTRKvX=(TH1F*)l->FindObject("hTRKvX");
	48	hTRKvX->Draw();
	49	cv->cd(8);
	50	TH1F* hTRKvY=(TH1F*)l->FindObject("hTRKvY");
	51	hTRKvY->Draw();
	52	cv->cd(9);
	53	TH1F* hTRKvZ=(TH1F*)l->FindObject("hTRKvZ");
	54	hTRKvZ->Draw();
	55
	56	TCanvas* c1=new TCanvas("c1","Multipl");
	57	c1->Divide(3,1);
	58	c1->cd(1);
	59	TH1F* hTracklets=(TH1F*)l->FindObject("hTracklets");
	60	hTracklets->Draw();
	61	c1->cd(2);
	62	TH1F* hTracks=(TH1F*)l->FindObject("hTracks");
	63	hTracks->Draw();
	64	c1->cd(3);
65	TH1F* hSelTracks=(TH1F*)l->FindObject("hSelTracks");
66	hSelTracks->Draw();
67
68	TH1F* hncharmed=(TH1F*)l->FindObject("hncharmed");
69	TCanvas* cn=new TCanvas("cn","ncharm");
70	hncharmed->Draw("box");
71	hncharmed->GetXaxis()->SetTitle("dNch/dy");
72	hncharmed->GetYaxis()->SetTitle("N Charm hadrons in golden channels");
73	cn->Update();
74
75	TH1F* hnbvsnc=(TH1F*)l->FindObject("hnbvsnc");
76	TCanvas* cnhf=new TCanvas("cnhf","nb/c");
77	hnbvsnc->Draw("box");
78	hnbvsnc->GetXaxis()->SetTitle("Nc");
79	hnbvsnc->GetYaxis()->SetTitle("Nb");
80	cnhf->Update();
81
82	TH2F* hyptd0prompt=(TH2F*)l->FindObject("hyptd0prompt");
83	TH2F* hyptd0feeddown=(TH2F*)l->FindObject("hyptd0feeddown");
84	TH2F* hyptD02=(TH2F*)l->FindObject("hyptD02");
85	TH2F* hyptD04=(TH2F*)l->FindObject("hyptD04");
86
87	TCanvas* cd0=new TCanvas("cd0","D0");
88	cd0->Divide(2,2);
89	cd0->cd(1);
90	hyptd0prompt->Draw("colz");
91	cd0->cd(2);
92	hyptd0feeddown->Draw("colz");
93	cd0->cd(3);
94	hyptD02->Draw("colz");
95	cd0->cd(4);
96	hyptD04->Draw("colz");
97
98	TH2F* hyptdplusprompt=(TH2F*)l->FindObject("hyptdplusprompt");
99	TH2F* hyptdplusfeeddown=(TH2F*)l->FindObject("hyptdplusfeeddown");
100	TH2F* hyptDplusnonreson=(TH2F*)l->FindObject("hyptDplusnonreson");
101	TH2F* hyptDplusreson=(TH2F*)l->FindObject("hyptDplusreson");
102
103	TCanvas* cdplus=new TCanvas("cdplus","Dplus");
104	cdplus->Divide(2,2);
105	cdplus->cd(1);
106	hyptdplusprompt->Draw("colz");
107	cdplus->cd(2);
108	hyptdplusfeeddown->Draw("colz");
109	cdplus->cd(3);
110	hyptDplusnonreson->Draw("colz");
111	cdplus->cd(4);
112	hyptDplusreson->Draw("colz");
113
114	TH2F* hyptdsprompt=(TH2F*)l->FindObject("hyptdsprompt");
115	TH2F* hyptdsfeeddown=(TH2F*)l->FindObject("hyptdsfeedown");
116	TH2F* hyptdsphi=(TH2F*)l->FindObject("hyptdsphi");
117	TH2F* hyptdsK0st=(TH2F*)l->FindObject("hyptdsk0st");
118
119	TCanvas* cds=new TCanvas("cds","Ds");
120	cds->Divide(2,2);
121	cds->cd(1);
122	hyptdsprompt->Draw("colz");
123	cds->cd(2);
124	hyptdsfeeddown->Draw("colz");
125	cds->cd(3);
126	hyptdsphi->Draw("colz");
127	cds->cd(4);
128	hyptdsK0st->Draw("colz");
129
130	TH2F* hyptdstarprompt=(TH2F*)l->FindObject("hyptdstarprompt");
131	TH2F* hyptdstarfeedown=(TH2F*)l->FindObject("hyptdstarfeedown");
132	TH2F* hyptlcprompt=(TH2F*)l->FindObject("hyptlcprompt");
133	TH2F* hyptlcfeedown=(TH2F*)l->FindObject("hyptlcfeedown");
134
135	TCanvas* cdstlc=new TCanvas("cdstls","Dstar LambdaC");
136	cdstlc->Divide(2,2);
137	cdstlc->cd(1);
138	hyptdstarprompt->Draw("colz");
139	cdstlc->cd(2);
140	hyptdstarfeedown->Draw("colz");
141	cdstlc->cd(3);
142	hyptlcprompt->Draw("colz");
143	cdstlc->cd(4);
144	hyptlcfeedown->Draw("colz");
145
146	}