[u/mrichter/AliRoot.git] / HLT / QA / tasks / macros / drawGlobalESDHistograms.C

// $Id$
/*
 * Plotting macro for comparing offline- and HLT- ESD trees from
 * HLT-OFFLINE-GLOBAL-comparison.root produced using $ALICE_ROOT/HLT/QA/tasks/AliAnalysisTaskHLT.*
 *
 * Usage: aliroot drawGlobalESDHistograms.C'("HLT-OFFLINE-GLOBAL-comparison.root")'
 *
 * It saves the canvas with the output histograms in a png file.
 *
 * @ingroup alihlt_qa
 * @author Camilla.Stokkevag@student.uib.no, Kalliopi.Kanaki@ift.uib.no 
 */

void drawGlobalESDHistograms(const char* filename="HLT-OFFLINE-GLOBAL-comparison.root"){
 
 gROOT->SetStyle("Plain");
 gStyle->SetPalette(1);
 gStyle->SetOptStat(10);

 TFile *f1 = TFile::Open(filename);
 TList *l1 = (TList*)f1->Get("global_histograms");
 

//=================================================================//
//--------------- READ HISTOGRAMS ---------------------------------//
//=================================================================//

 TH1F *h1 = l1->FindObject("fNcluster_hlt");
 TH1F *h2 = l1->FindObject("fNcluster_off"); 
 h1->SetTitle("cluster distribution");
 h1->GetXaxis()->SetTitle("# of clusters per track");
 h2->SetLineColor(2);
 if(h1->GetMaximum() >= h2->GetMaximum())  h2->SetMaximum(1.1*h1->GetMaximum());
 else h1->SetMaximum(1.1*h2->GetMaximum());
  

 TH1F *h3  = l1->FindObject("fDCA_hlt");
 TH1F *h4  = l1->FindObject("fDCA_off"); 
 TH1F *hSG = l1->FindObject("fDCA_hltSG"); 
 h3->SetTitle("DCA between track and vertex on XY plane");
 h3->GetXaxis()->SetTitle("DCAr (cm)");
 h4->SetLineColor(2);
 hSG->SetLineColor(kBlue);
 if(h3->GetMaximum() >= h4->GetMaximum()) h4->SetMaximum(1.1*h3->GetMaximum());
 else h3->SetMaximum(1.1*h4->GetMaximum());
 

 TH1F *h5 = l1->FindObject("fMult_hlt");
 TH1F *h6 = l1->FindObject("fMult_off"); 
 h5->SetTitle("track multiplicity");
 h6->SetLineColor(2);
 if(h5->GetMaximum() > h6->GetMaximum()) h6->SetMaximum(1.1*h5->GetMaximum());
 else h5->SetMaximum(1.1*h6->GetMaximum());


 TH1F *h7 = l1->FindObject("fCharge_hlt");
 TH1F *h8 = l1->FindObject("fCharge_off");
 h7->GetXaxis()->SetTitle("polarity"); 
 h7->SetTitle("charge distribution");
 h8->SetLineColor(2);
 if(h7->GetMaximum() > h8->GetMaximum()) h8->SetMaximum(1.1*h7->GetMaximum());
 else h7->SetMaximum(1.1*h8->GetMaximum());
 

 TH1F *h9  = l1->FindObject("fMomentum_hlt");
 TH1F *h10 = l1->FindObject("fMomentum_off"); 
 h9->GetXaxis()->SetTitle("p_{t} (GeV/c)"); 
 h9->SetTitle("transverse momentum");
 h10->SetLineColor(2);
 if(h9->GetMaximum() > h10->GetMaximum()) h10->SetMaximum(1.1*h9->GetMaximum());
 else h9->SetMaximum(1.1*h10->GetMaximum());


 TH1F *h11 = l1->FindObject("fEta_hlt");
 TH1F *h12 = l1->FindObject("fEta_off"); 
 h11->SetTitle("pseudorapidity");
 h12->SetLineColor(2);
 if(h11->GetMaximum() > h12->GetMaximum()) h12->SetMaximum(1.1*h11->GetMaximum());
 else h11->SetMaximum(1.1*h12->GetMaximum());


 TH1F *h13 = l1->FindObject("fXvertex_hlt");
 TH1F *h14 = l1->FindObject("fXvertex_off"); 
 h13->GetXaxis()->SetTitle("X (cm)");
 h13->SetTitle("X of primary vertex");
 h14->SetLineColor(2);
 if(h13->GetMaximum() > h14->GetMaximum()) h14->SetMaximum(1.1*h13->GetMaximum());
 else h13->SetMaximum(1.1*h14->GetMaximum());

 TH1F *h15 = l1->FindObject("fYvertex_hlt");
 TH1F *h16 = l1->FindObject("fYvertex_off"); 
 h15->GetXaxis()->SetTitle("Y (cm)");
 h15->SetTitle("Y of primary vertex");
 h16->SetLineColor(2);
 if(h15->GetMaximum() > h16->GetMaximum()) h16->SetMaximum(1.1*h15->GetMaximum());
 else h15->SetMaximum(1.1*h16->GetMaximum());

 TH1F *h17 = l1->FindObject("fZvertex_hlt");
 TH1F *h18 = l1->FindObject("fZvertex_off"); 
 h17->GetXaxis()->SetTitle("Z (cm)");
 h17->SetTitle("Z of primary vertex");
 h18->SetLineColor(2);
 if(h17->GetMaximum() > h18->GetMaximum()) h18->SetMaximum(1.1*h17->GetMaximum());
 else h17->SetMaximum(1.1*h18->GetMaximum());


//  TH2F *h15 = l1->FindObject("fXYvertex_off");
//  h15->GetXaxis()->SetTitle("X (cm)");
//  h15->GetYaxis()->SetTitle("Y (cm)");
//  h15->SetTitle("XY primary vertex Offline");
// 
//  TH2F *h16 = l1->FindObject("fXYvertex_hlt"); 
//  h16->GetXaxis()->SetTitle("X (cm)");
//  h16->GetYaxis()->SetTitle("Y (cm)");
//  h16->SetTitle("XY of primary vertex HLT");

 TLegend *leg1 = new TLegend(0.6,0.6,0.8,0.8);
 leg1->SetFillColor(10);
 leg1->SetLineColor(10);
 leg1->AddEntry(h1,"HLT", "l");
 leg1->AddEntry(h2,"OFF", "l");


//=================================================================//
//--------------------- DRAW HISTOGRAMS ---------------------------//
//=================================================================//


 TCanvas *c1 = new TCanvas("c1","HLT vs offline",1300,800);
 c1->Divide(3,3);

 c1->cd(1);
 h1->Draw();
 h2->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st1 = (TPaveStats*)h1->FindObject("stats");
 st1->SetLineColor(0);

 gPad->Update();
 TPaveStats *st2 = (TPaveStats*)h2->FindObject("stats");
 st2->SetY1NDC(st1->GetY1NDC()-0.05);
 st2->SetY2NDC(st1->GetY2NDC()-0.05);
 st2->SetLineColor(0);
 st2->SetTextColor(h2->GetLineColor());
 st2->SetFillStyle(0);
 st2->Draw();

 //======================

 c1->cd(2)->SetLogy();
 h3->Draw();
 h4->Draw("sames");
 hSG->Draw("sames");
 leg1->Draw("same");

 gPad->Update();
 TPaveStats *st3 = (TPaveStats*)h3->FindObject("stats");
 st3->SetLineColor(0);

 gPad->Update();
 TPaveStats *st4 = (TPaveStats*)h4->FindObject("stats");
 st4->SetY1NDC(st3->GetY1NDC()-0.05);
 st4->SetY2NDC(st3->GetY2NDC()-0.05);
 st4->SetLineColor(0);
 st4->SetTextColor(h4->GetLineColor());
 st4->SetFillStyle(0);
 st4->Draw();

 gPad->Update();
 TPaveStats *stSG = (TPaveStats*)hSG->FindObject("stats");
 stSG->SetY1NDC(st4->GetY1NDC()-0.05);
 stSG->SetY2NDC(st4->GetY2NDC()-0.05);
 stSG->SetLineColor(0);
 stSG->SetTextColor(hSG->GetLineColor());
 stSG->SetFillStyle(0);
 stSG->Draw();
 
//======================

 c1->cd(3);
 h5->Draw();
 h6->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st5 = (TPaveStats*)h5->FindObject("stats");
 st5->SetLineColor(0);

 gPad->Update();
 TPaveStats *st6 = (TPaveStats*)h6->FindObject("stats");
 st6->SetY1NDC(st5->GetY1NDC()-0.05);
 st6->SetY2NDC(st5->GetY2NDC()-0.05);
 st6->SetLineColor(0);
 st6->SetTextColor(h6->GetLineColor());
 st6->SetFillStyle(0);
 st6->Draw();
 
//======================

 c1->cd(4);
 h7->Draw();
 h8->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st7 = (TPaveStats*)h7->FindObject("stats");
 st7->SetLineColor(0);

 gPad->Update();
 TPaveStats *st8 = (TPaveStats*)h8->FindObject("stats");
 st8->SetY1NDC(st7->GetY1NDC()-0.05);
 st8->SetY2NDC(st7->GetY2NDC()-0.05);
 st8->SetLineColor(0);
 st8->SetTextColor(h8->GetLineColor());
 st8->SetFillStyle(0);
 st8->Draw();
 
//======================
 
 c1->cd(5);
 h9->Draw();
 h10->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st9 = (TPaveStats*)h9->FindObject("stats");
 st9->SetLineColor(0);

 gPad->Update();
 TPaveStats *st10 = (TPaveStats*)h10->FindObject("stats");
 st10->SetY1NDC(st9->GetY1NDC()-0.05);
 st10->SetY2NDC(st9->GetY2NDC()-0.05);
 st10->SetLineColor(0);
 st10->SetTextColor(h10->GetLineColor());
 st10->SetFillStyle(0);
 st10->Draw();
 
//======================
 
 c1->cd(6);
 h11->Draw();
 h12->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st11 = (TPaveStats*)h11->FindObject("stats");
 st11->SetLineColor(0);

 gPad->Update();
 TPaveStats *st12 = (TPaveStats*)h12->FindObject("stats");
 st12->SetY1NDC(st11->GetY1NDC()-0.05);
 st12->SetY2NDC(st11->GetY2NDC()-0.05);
 st12->SetLineColor(0);
 st12->SetTextColor(h12->GetLineColor());
 st12->SetFillStyle(0);
 st12->Draw();
 
//======================

 c1->cd(7);
 h13->Draw();
 h14->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st13 = (TPaveStats*)h13->FindObject("stats");
 st13->SetLineColor(0);

 gPad->Update();
 TPaveStats *st14 = (TPaveStats*)h14->FindObject("stats");
 st14->SetY1NDC(st13->GetY1NDC()-0.05);
 st14->SetY2NDC(st13->GetY2NDC()-0.05);
 st14->SetLineColor(0);
 st14->SetTextColor(h14->GetLineColor());
 st14->SetFillStyle(0);
 st14->Draw();

//======================

 c1->cd(8);
 h15->Draw();
 h16->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st15 = (TPaveStats*)h15->FindObject("stats");
 st15->SetLineColor(0);

 gPad->Update();
 TPaveStats *st16 = (TPaveStats*)h16->FindObject("stats");
 st16->SetY1NDC(st15->GetY1NDC()-0.05);
 st16->SetY2NDC(st15->GetY2NDC()-0.05);
 st16->SetLineColor(0);
 st16->SetTextColor(h16->GetLineColor());
 st16->SetFillStyle(0);
 st16->Draw();

//======================

 c1->cd(9);
 h17->Draw();
 h18->Draw("sames");
 leg1->Draw("same");
 
 gPad->Update();
 TPaveStats *st17 = (TPaveStats*)h17->FindObject("stats");
 st17->SetLineColor(0);

 gPad->Update();
 TPaveStats *st18 = (TPaveStats*)h18->FindObject("stats");
 st18->SetY1NDC(st17->GetY1NDC()-0.05);
 st18->SetY2NDC(st17->GetY2NDC()-0.05);
 st18->SetLineColor(0);
 st18->SetTextColor(h18->GetLineColor());
 st18->SetFillStyle(0);
 st18->Draw();

//======================

 c1->SaveAs("HLT-offline.png");  

}
Commit	Line	Data
0740807a	1	// $Id$
	2	/*
	3	* Plotting macro for comparing offline- and HLT- ESD trees from
	4	* HLT-OFFLINE-GLOBAL-comparison.root produced using $ALICE_ROOT/HLT/QA/tasks/AliAnalysisTaskHLT.*
	5	*
	6	* Usage: aliroot drawGlobalESDHistograms.C'("HLT-OFFLINE-GLOBAL-comparison.root")'
	7	*
	8	* It saves the canvas with the output histograms in a png file.
	9	*
	10	* @ingroup alihlt_qa
	11	* @author Camilla.Stokkevag@student.uib.no, Kalliopi.Kanaki@ift.uib.no
	12	*/
	13
	14	void drawGlobalESDHistograms(const char* filename="HLT-OFFLINE-GLOBAL-comparison.root"){
	15
	16	gROOT->SetStyle("Plain");
	17	gStyle->SetPalette(1);
	18	gStyle->SetOptStat(10);
	19
	20	TFile *f1 = TFile::Open(filename);
	21	TList l1 = (TList)f1->Get("global_histograms");
	22
	23
	24	//=================================================================//
	25	//--------------- READ HISTOGRAMS ---------------------------------//
	26	//=================================================================//
	27
	28	TH1F *h1 = l1->FindObject("fNcluster_hlt");
	29	TH1F *h2 = l1->FindObject("fNcluster_off");
	30	h1->SetTitle("cluster distribution");
	31	h1->GetXaxis()->SetTitle("# of clusters per track");
	32	h2->SetLineColor(2);
	33	if(h1->GetMaximum() >= h2->GetMaximum()) h2->SetMaximum(1.1*h1->GetMaximum());
	34	else h1->SetMaximum(1.1*h2->GetMaximum());
	35
	36
	37
3a510b95	38	TH1F *h3 = l1->FindObject("fDCA_hlt");
	39	TH1F *h4 = l1->FindObject("fDCA_off");
	40	TH1F *hSG = l1->FindObject("fDCA_hltSG");
0740807a	41	h3->SetTitle("DCA between track and vertex on XY plane");
	42	h3->GetXaxis()->SetTitle("DCAr (cm)");
	43	h4->SetLineColor(2);
3a510b95	44	hSG->SetLineColor(kBlue);
0740807a	45	if(h3->GetMaximum() >= h4->GetMaximum()) h4->SetMaximum(1.1*h3->GetMaximum());
	46	else h3->SetMaximum(1.1*h4->GetMaximum());
	47
	48
	49	TH1F *h5 = l1->FindObject("fMult_hlt");
	50	TH1F *h6 = l1->FindObject("fMult_off");
	51	h5->SetTitle("track multiplicity");
	52	h6->SetLineColor(2);
	53	if(h5->GetMaximum() > h6->GetMaximum()) h6->SetMaximum(1.1*h5->GetMaximum());
	54	else h5->SetMaximum(1.1*h6->GetMaximum());
	55
	56
	57	TH1F *h7 = l1->FindObject("fCharge_hlt");
	58	TH1F *h8 = l1->FindObject("fCharge_off");
	59	h7->GetXaxis()->SetTitle("polarity");
	60	h7->SetTitle("charge distribution");
	61	h8->SetLineColor(2);
	62	if(h7->GetMaximum() > h8->GetMaximum()) h8->SetMaximum(1.1*h7->GetMaximum());
	63	else h7->SetMaximum(1.1*h8->GetMaximum());
	64
	65
	66	TH1F *h9 = l1->FindObject("fMomentum_hlt");
	67	TH1F *h10 = l1->FindObject("fMomentum_off");
	68	h9->GetXaxis()->SetTitle("p_{t} (GeV/c)");
	69	h9->SetTitle("transverse momentum");
	70	h10->SetLineColor(2);
	71	if(h9->GetMaximum() > h10->GetMaximum()) h10->SetMaximum(1.1*h9->GetMaximum());
	72	else h9->SetMaximum(1.1*h10->GetMaximum());
	73
	74
	75	TH1F *h11 = l1->FindObject("fEta_hlt");
	76	TH1F *h12 = l1->FindObject("fEta_off");
	77	h11->SetTitle("pseudorapidity");
	78	h12->SetLineColor(2);
	79	if(h11->GetMaximum() > h12->GetMaximum()) h12->SetMaximum(1.1*h11->GetMaximum());
	80	else h11->SetMaximum(1.1*h12->GetMaximum());
	81
	82
	83	TH1F *h13 = l1->FindObject("fXvertex_hlt");
	84	TH1F *h14 = l1->FindObject("fXvertex_off");
	85	h13->GetXaxis()->SetTitle("X (cm)");
	86	h13->SetTitle("X of primary vertex");
	87	h14->SetLineColor(2);
	88	if(h13->GetMaximum() > h14->GetMaximum()) h14->SetMaximum(1.1*h13->GetMaximum());
	89	else h13->SetMaximum(1.1*h14->GetMaximum());
	90
	91	TH1F *h15 = l1->FindObject("fYvertex_hlt");
	92	TH1F *h16 = l1->FindObject("fYvertex_off");
	93	h15->GetXaxis()->SetTitle("Y (cm)");
	94	h15->SetTitle("Y of primary vertex");
	95	h16->SetLineColor(2);
	96	if(h15->GetMaximum() > h16->GetMaximum()) h16->SetMaximum(1.1*h15->GetMaximum());
	97	else h15->SetMaximum(1.1*h16->GetMaximum());
	98
	99	TH1F *h17 = l1->FindObject("fZvertex_hlt");
	100	TH1F *h18 = l1->FindObject("fZvertex_off");
	101	h17->GetXaxis()->SetTitle("Z (cm)");
	102	h17->SetTitle("Z of primary vertex");
	103	h18->SetLineColor(2);
	104	if(h17->GetMaximum() > h18->GetMaximum()) h18->SetMaximum(1.1*h17->GetMaximum());
	105	else h17->SetMaximum(1.1*h18->GetMaximum());
	106
	107
	108	// TH2F *h15 = l1->FindObject("fXYvertex_off");
109	// h15->GetXaxis()->SetTitle("X (cm)");
110	// h15->GetYaxis()->SetTitle("Y (cm)");
111	// h15->SetTitle("XY primary vertex Offline");
112	//
113	// TH2F *h16 = l1->FindObject("fXYvertex_hlt");
114	// h16->GetXaxis()->SetTitle("X (cm)");
115	// h16->GetYaxis()->SetTitle("Y (cm)");
116	// h16->SetTitle("XY of primary vertex HLT");
117
118	TLegend *leg1 = new TLegend(0.6,0.6,0.8,0.8);
119	leg1->SetFillColor(10);
120	leg1->SetLineColor(10);
121	leg1->AddEntry(h1,"HLT", "l");
122	leg1->AddEntry(h2,"OFF", "l");
123
124
125
126	//=================================================================//
127	//--------------------- DRAW HISTOGRAMS ---------------------------//
128	//=================================================================//
129
130
131
132	TCanvas *c1 = new TCanvas("c1","HLT vs offline",1300,800);
133	c1->Divide(3,3);
134
135	c1->cd(1);
136	h1->Draw();
137	h2->Draw("sames");
138	leg1->Draw("same");
139
140	gPad->Update();
141	TPaveStats st1 = (TPaveStats)h1->FindObject("stats");
142	st1->SetLineColor(0);
143
144	gPad->Update();
145	TPaveStats st2 = (TPaveStats)h2->FindObject("stats");
146	st2->SetY1NDC(st1->GetY1NDC()-0.05);
147	st2->SetY2NDC(st1->GetY2NDC()-0.05);
148	st2->SetLineColor(0);
149	st2->SetTextColor(h2->GetLineColor());
150	st2->SetFillStyle(0);
151	st2->Draw();
152
153	//======================
154
155	c1->cd(2)->SetLogy();
156	h3->Draw();
157	h4->Draw("sames");
3a510b95	158	hSG->Draw("sames");
0740807a	159	leg1->Draw("same");
	160
	161	gPad->Update();
	162	TPaveStats st3 = (TPaveStats)h3->FindObject("stats");
	163	st3->SetLineColor(0);
	164
	165	gPad->Update();
	166	TPaveStats st4 = (TPaveStats)h4->FindObject("stats");
	167	st4->SetY1NDC(st3->GetY1NDC()-0.05);
	168	st4->SetY2NDC(st3->GetY2NDC()-0.05);
	169	st4->SetLineColor(0);
	170	st4->SetTextColor(h4->GetLineColor());
	171	st4->SetFillStyle(0);
	172	st4->Draw();
3a510b95	173
	174	gPad->Update();
	175	TPaveStats stSG = (TPaveStats)hSG->FindObject("stats");
	176	stSG->SetY1NDC(st4->GetY1NDC()-0.05);
	177	stSG->SetY2NDC(st4->GetY2NDC()-0.05);
	178	stSG->SetLineColor(0);
	179	stSG->SetTextColor(hSG->GetLineColor());
	180	stSG->SetFillStyle(0);
	181	stSG->Draw();
0740807a	182
	183	//======================
	184
	185	c1->cd(3);
	186	h5->Draw();
	187	h6->Draw("sames");
	188	leg1->Draw("same");
	189
	190	gPad->Update();
	191	TPaveStats st5 = (TPaveStats)h5->FindObject("stats");
	192	st5->SetLineColor(0);
	193
	194	gPad->Update();
	195	TPaveStats st6 = (TPaveStats)h6->FindObject("stats");
	196	st6->SetY1NDC(st5->GetY1NDC()-0.05);
	197	st6->SetY2NDC(st5->GetY2NDC()-0.05);
	198	st6->SetLineColor(0);
	199	st6->SetTextColor(h6->GetLineColor());
	200	st6->SetFillStyle(0);
	201	st6->Draw();
	202
	203	//======================
	204
	205	c1->cd(4);
	206	h7->Draw();
	207	h8->Draw("sames");
	208	leg1->Draw("same");
	209
	210	gPad->Update();
	211	TPaveStats st7 = (TPaveStats)h7->FindObject("stats");
	212	st7->SetLineColor(0);
	213
	214	gPad->Update();
	215	TPaveStats st8 = (TPaveStats)h8->FindObject("stats");
	216	st8->SetY1NDC(st7->GetY1NDC()-0.05);
	217	st8->SetY2NDC(st7->GetY2NDC()-0.05);
	218	st8->SetLineColor(0);
	219	st8->SetTextColor(h8->GetLineColor());
	220	st8->SetFillStyle(0);
	221	st8->Draw();
	222
	223	//======================
	224
	225	c1->cd(5);
	226	h9->Draw();
	227	h10->Draw("sames");
	228	leg1->Draw("same");
	229
	230	gPad->Update();
	231	TPaveStats st9 = (TPaveStats)h9->FindObject("stats");
	232	st9->SetLineColor(0);
	233
	234	gPad->Update();
	235	TPaveStats st10 = (TPaveStats)h10->FindObject("stats");
	236	st10->SetY1NDC(st9->GetY1NDC()-0.05);
	237	st10->SetY2NDC(st9->GetY2NDC()-0.05);
	238	st10->SetLineColor(0);
	239	st10->SetTextColor(h10->GetLineColor());
	240	st10->SetFillStyle(0);
	241	st10->Draw();
	242
	243	//======================
	244
	245	c1->cd(6);
246	h11->Draw();
247	h12->Draw("sames");
248	leg1->Draw("same");
249
250	gPad->Update();
251	TPaveStats st11 = (TPaveStats)h11->FindObject("stats");
252	st11->SetLineColor(0);
253
254	gPad->Update();
255	TPaveStats st12 = (TPaveStats)h12->FindObject("stats");
256	st12->SetY1NDC(st11->GetY1NDC()-0.05);
257	st12->SetY2NDC(st11->GetY2NDC()-0.05);
258	st12->SetLineColor(0);
259	st12->SetTextColor(h12->GetLineColor());
260	st12->SetFillStyle(0);
261	st12->Draw();
262
263	//======================
264
265	c1->cd(7);
266	h13->Draw();
267	h14->Draw("sames");
268	leg1->Draw("same");
269
270	gPad->Update();
271	TPaveStats st13 = (TPaveStats)h13->FindObject("stats");
272	st13->SetLineColor(0);
273
274	gPad->Update();
275	TPaveStats st14 = (TPaveStats)h14->FindObject("stats");
276	st14->SetY1NDC(st13->GetY1NDC()-0.05);
277	st14->SetY2NDC(st13->GetY2NDC()-0.05);
278	st14->SetLineColor(0);
279	st14->SetTextColor(h14->GetLineColor());
280	st14->SetFillStyle(0);
281	st14->Draw();
282
283	//======================
284
285	c1->cd(8);
286	h15->Draw();
287	h16->Draw("sames");
288	leg1->Draw("same");
289
290	gPad->Update();
291	TPaveStats st15 = (TPaveStats)h15->FindObject("stats");
292	st15->SetLineColor(0);
293
294	gPad->Update();
295	TPaveStats st16 = (TPaveStats)h16->FindObject("stats");
296	st16->SetY1NDC(st15->GetY1NDC()-0.05);
297	st16->SetY2NDC(st15->GetY2NDC()-0.05);
298	st16->SetLineColor(0);
299	st16->SetTextColor(h16->GetLineColor());
300	st16->SetFillStyle(0);
301	st16->Draw();
302
303	//======================
304
305	c1->cd(9);
306	h17->Draw();
307	h18->Draw("sames");
308	leg1->Draw("same");
309
310	gPad->Update();
311	TPaveStats st17 = (TPaveStats)h17->FindObject("stats");
312	st17->SetLineColor(0);
313
314	gPad->Update();
315	TPaveStats st18 = (TPaveStats)h18->FindObject("stats");
316	st18->SetY1NDC(st17->GetY1NDC()-0.05);
317	st18->SetY2NDC(st17->GetY2NDC()-0.05);
318	st18->SetLineColor(0);
319	st18->SetTextColor(h18->GetLineColor());
320	st18->SetFillStyle(0);
321	st18->Draw();
322
323	//======================
324
325	c1->SaveAs("HLT-offline.png");
326
327	}