dasu/UltraFastSim/Plot.C

{
// Make plots
  TFile zh("ZHmmbbOut.root");
  TFile zj("ZmmOut.root");
  TFile zz("ZZOut.root");

  TH1F *zhmuonpt=zh.Get("Muons Pt");
  TH1F *zjmuonpt=zj.Get("Muons Pt");
  TH1F *zzmuonpt=zz.Get("Muons Pt");

  zjmuonpt->SetLineColor(kGreen);
  zzmuonpt->SetLineColor(kOrange+3);
  zhmuonpt->SetLineColor(kRed);

  zjmuonpt->Scale(1./zjmuonpt->Integral());
  zzmuonpt->Scale(1./zzmuonpt->Integral());
  zhmuonpt->Scale(1./zhmuonpt->Integral());

  TCanvas c1;

  zjmuonpt->Draw();
  zzmuonpt->Draw("same");
  zhmuonpt->Draw("same");

  TLegend *lg=new TLegend(0.5,0.5,0.7,0.7);
  lg->SetFillColor(kWhite);
  lg->AddEntry(zhmuonpt,"ZH, m(H)=120 GeV","l");
  lg->AddEntry(zjmuonpt,"Z+jets","l");
  lg->AddEntry(zzmuonpt,"ZZ","l");
  lg->Draw();

  /////////////////////////////

  TH1F *zhbjetpt=zh.Get("B-jets Et");
  TH1F *zjbjetpt=zj.Get("B-jets Et");
  TH1F *zzbjetpt=zz.Get("B-jets Et");

  zjbjetpt->SetLineColor(kGreen);
  zzbjetpt->SetLineColor(kOrange+3);
  zhbjetpt->SetLineColor(kRed);

  zjbjetpt->Scale(1./zjbjetpt->Integral());
  zzbjetpt->Scale(1./zzbjetpt->Integral());
  zhbjetpt->Scale(1./zhbjetpt->Integral());

  TCanvas c2;

  zjbjetpt->Draw();
  zzbjetpt->Draw("same");
  zhbjetpt->Draw("same");

  TLegend *lg2=new TLegend(0.5,0.5,0.7,0.7);
  lg2->SetFillColor(kWhite);
  lg2->AddEntry(zhbjetpt,"ZH, m(H)=120 GeV","l");
  lg2->AddEntry(zjbjetpt,"Z+jets","l");
  lg2->AddEntry(zzbjetpt,"ZZ","l");
  lg2->Draw();

  /////////////////////////////

  TH1F *zhzpt=zh.Get("Z Pt");
  TH1F *zjzpt=zj.Get("Z Pt");
  TH1F *zzzpt=zz.Get("Z Pt");

  zjzpt->SetLineColor(kGreen);
  zzzpt->SetLineColor(kOrange+3);
  zhzpt->SetLineColor(kRed);

  zjzpt->Scale(1./zjzpt->Integral());
  zzzpt->Scale(1./zzzpt->Integral());
  zhzpt->Scale(1./zhzpt->Integral());

  TCanvas c3;

  zjzpt->Draw();
  zzzpt->Draw("same");
  zhzpt->Draw("same");

  TLegend *lg3=new TLegend(0.5,0.5,0.7,0.7);
  lg3->SetFillColor(kWhite);
  lg3->AddEntry(zhzpt,"ZH, m(H)=120 GeV","l");
  lg3->AddEntry(zjzpt,"Z+jets","l");
  lg3->AddEntry(zzzpt,"ZZ","l");
  lg3->Draw();


  /////////////////////////////

  TH1F *zhhpt=zh.Get("H Pt");
  TH1F *zjhpt=zj.Get("H Pt");
  TH1F *zzhpt=zz.Get("H Pt");

  zjhpt->SetLineColor(kGreen);
  zzhpt->SetLineColor(kOrange+3);
  zhhpt->SetLineColor(kRed);

  zjhpt->Rebin(5);
  zzhpt->Rebin(5);
  zhhpt->Rebin(5);

  zjhpt->Scale(1./zjhpt->Integral());
  zzhpt->Scale(1./zzhpt->Integral());
  zhhpt->Scale(1./zhhpt->Integral());

  TCanvas c4;

  zjhpt->Draw();
  zzhpt->Draw("same");
  zhhpt->Draw("same");

  TLegend *lg4=new TLegend(0.5,0.5,0.7,0.7);
  lg4->SetFillColor(kWhite);
  lg4->AddEntry(zhhpt,"ZH, m(H)=120 GeV","l");
  lg4->AddEntry(zjhpt,"Z+jets","l");
  lg4->AddEntry(zzhpt,"ZZ","l");
  lg4->Draw();


  /////////////////////////////

  TH1F *zhzinvmass=zh.Get("Z invmass");
  TH1F *zjzinvmass=zj.Get("Z invmass");
  TH1F *zzzinvmass=zz.Get("Z invmass");

  zjzinvmass->SetLineColor(kGreen);
  zzzinvmass->SetLineColor(kOrange+3);
  zhzinvmass->SetLineColor(kRed);

  zjzinvmass->Scale(1./zjzinvmass->Integral());
  zzzinvmass->Scale(1./zzzinvmass->Integral());
  zhzinvmass->Scale(1./zhzinvmass->Integral());

  TCanvas c5;

  zjzinvmass->Draw();
  zzzinvmass->Draw("same");
  zhzinvmass->Draw("same");

  TLegend *lg5=new TLegend(0.5,0.5,0.7,0.7);
  lg5->SetFillColor(kWhite);
  lg5->AddEntry(zhzinvmass,"ZH, m(H)=120 GeV","l");
  lg5->AddEntry(zjzinvmass,"Z+jets","l");
  lg5->AddEntry(zzzinvmass,"ZZ","l");
  lg5->Draw();


  /////////////////////////////

  TH1F *zhhinvmass=zh.Get("H invmass preselection");
  TH1F *zjhinvmass=zj.Get("H invmass preselection");
  TH1F *zzhinvmass=zz.Get("H invmass preselection");

  zjhinvmass->SetLineColor(kGreen);
  zzhinvmass->SetLineColor(kOrange+3);
  zhhinvmass->SetLineColor(kRed);

  zjhinvmass->Rebin(5);
  zzhinvmass->Rebin(5);
  zhhinvmass->Rebin(5);

  zjhinvmass->Scale(1./zjhinvmass->Integral());
  zzhinvmass->Scale(1./zzhinvmass->Integral());
  zhhinvmass->Scale(1./zhhinvmass->Integral());

  TCanvas c6;

  zjhinvmass->Draw();
  zzhinvmass->Draw("same");
  zhhinvmass->Draw("same");

  TLegend *lg6=new TLegend(0.5,0.5,0.7,0.7);
  lg6->SetFillColor(kWhite);
  lg6->AddEntry(zhhinvmass,"ZH, m(H)=120 GeV","l");
  lg6->AddEntry(zjhinvmass,"Z+jets","l");
  lg6->AddEntry(zzhinvmass,"ZZ","l");
  lg6->Draw();

  /////////////////////////////

  TH1F *zhhinvmass2=zh.Get("H invmass");
  TH1F *zjhinvmass2=zj.Get("H invmass");
  TH1F *zzhinvmass2=zz.Get("H invmass");

  zjhinvmass2->SetLineColor(kGreen);
  zzhinvmass2->SetLineColor(kOrange+3);
  zhhinvmass2->SetLineColor(kRed);

  zjhinvmass2->Rebin(5);
  zzhinvmass2->Rebin(5);
  zhhinvmass2->Rebin(5);

  //zjhinvmass2->Scale(1./zjhinvmass2->Integral());
  zzhinvmass2->Scale(1./zzhinvmass2->Integral());
  zhhinvmass2->Scale(1./zhhinvmass2->Integral());

  TCanvas c7;

  zjhinvmass2->Draw();
  zzhinvmass2->Draw("same");
  zhhinvmass2->Draw("same");

  TLegend *lg7=new TLegend(0.5,0.5,0.7,0.7);
  lg7->SetFillColor(kWhite);
  lg7->AddEntry(zhhinvmass2,"ZH, m(H)=120 GeV","l");
  lg7->AddEntry(zjhinvmass2,"Z+jets","l");
  lg7->AddEntry(zzhinvmass2,"ZZ","l");
  lg7->Draw();

  /////////////////////////////

  TH1F *zhdphi=zh.Get("Dphi(Z,H)");
  TH1F *zjdphi=zj.Get("Dphi(Z,H)");
  TH1F *zzdphi=zz.Get("Dphi(Z,H)");

  zjdphi->SetLineColor(kGreen);
  zzdphi->SetLineColor(kOrange+3);
  zhdphi->SetLineColor(kRed);

  zjdphi->Scale(1./zjdphi->Integral());
  zzdphi->Scale(1./zzdphi->Integral());
  zhdphi->Scale(1./zhdphi->Integral());

  TCanvas c8;

  zjdphi->Draw();
  zzdphi->Draw("same");
  zhdphi->Draw("same");

  TLegend *lg8=new TLegend(0.5,0.5,0.7,0.7);
  lg8->SetFillColor(kWhite);
  lg8->AddEntry(zhdphi,"ZH, m(H)=120 GeV","l");
  lg8->AddEntry(zjdphi,"Z+jets","l");
  lg8->AddEntry(zzdphi,"ZZ","l");
  lg8->Draw();


}
Revision:	1.2
Committed:	Thu Feb 24 16:36:34 2011 UTC (14 years, 2 months ago) by aarond
Content type:	text/plain
Branch:	MAIN
CVS Tags:	v2
Changes since 1.1:	+1 -0 lines
Log Message:	test of commit
#	User	Rev	Content
1	dasu	1.1	{
2	aarond	1.2	// Make plots
3	dasu	1.1	TFile zh("ZHmmbbOut.root");
4			TFile zj("ZmmOut.root");
5			TFile zz("ZZOut.root");
6
7			TH1F *zhmuonpt=zh.Get("Muons Pt");
8			TH1F *zjmuonpt=zj.Get("Muons Pt");
9			TH1F *zzmuonpt=zz.Get("Muons Pt");
10
11			zjmuonpt->SetLineColor(kGreen);
12			zzmuonpt->SetLineColor(kOrange+3);
13			zhmuonpt->SetLineColor(kRed);
14
15			zjmuonpt->Scale(1./zjmuonpt->Integral());
16			zzmuonpt->Scale(1./zzmuonpt->Integral());
17			zhmuonpt->Scale(1./zhmuonpt->Integral());
18
19			TCanvas c1;
20
21			zjmuonpt->Draw();
22			zzmuonpt->Draw("same");
23			zhmuonpt->Draw("same");
24
25			TLegend *lg=new TLegend(0.5,0.5,0.7,0.7);
26			lg->SetFillColor(kWhite);
27			lg->AddEntry(zhmuonpt,"ZH, m(H)=120 GeV","l");
28			lg->AddEntry(zjmuonpt,"Z+jets","l");
29			lg->AddEntry(zzmuonpt,"ZZ","l");
30			lg->Draw();
31
32			/////////////////////////////
33
34			TH1F *zhbjetpt=zh.Get("B-jets Et");
35			TH1F *zjbjetpt=zj.Get("B-jets Et");
36			TH1F *zzbjetpt=zz.Get("B-jets Et");
37
38			zjbjetpt->SetLineColor(kGreen);
39			zzbjetpt->SetLineColor(kOrange+3);
40			zhbjetpt->SetLineColor(kRed);
41
42			zjbjetpt->Scale(1./zjbjetpt->Integral());
43			zzbjetpt->Scale(1./zzbjetpt->Integral());
44			zhbjetpt->Scale(1./zhbjetpt->Integral());
45
46			TCanvas c2;
47
48			zjbjetpt->Draw();
49			zzbjetpt->Draw("same");
50			zhbjetpt->Draw("same");
51
52			TLegend *lg2=new TLegend(0.5,0.5,0.7,0.7);
53			lg2->SetFillColor(kWhite);
54			lg2->AddEntry(zhbjetpt,"ZH, m(H)=120 GeV","l");
55			lg2->AddEntry(zjbjetpt,"Z+jets","l");
56			lg2->AddEntry(zzbjetpt,"ZZ","l");
57			lg2->Draw();
58
59			/////////////////////////////
60
61			TH1F *zhzpt=zh.Get("Z Pt");
62			TH1F *zjzpt=zj.Get("Z Pt");
63			TH1F *zzzpt=zz.Get("Z Pt");
64
65			zjzpt->SetLineColor(kGreen);
66			zzzpt->SetLineColor(kOrange+3);
67			zhzpt->SetLineColor(kRed);
68
69			zjzpt->Scale(1./zjzpt->Integral());
70			zzzpt->Scale(1./zzzpt->Integral());
71			zhzpt->Scale(1./zhzpt->Integral());
72
73			TCanvas c3;
74
75			zjzpt->Draw();
76			zzzpt->Draw("same");
77			zhzpt->Draw("same");
78
79			TLegend *lg3=new TLegend(0.5,0.5,0.7,0.7);
80			lg3->SetFillColor(kWhite);
81			lg3->AddEntry(zhzpt,"ZH, m(H)=120 GeV","l");
82			lg3->AddEntry(zjzpt,"Z+jets","l");
83			lg3->AddEntry(zzzpt,"ZZ","l");
84			lg3->Draw();
85
86
87			/////////////////////////////
88
89			TH1F *zhhpt=zh.Get("H Pt");
90			TH1F *zjhpt=zj.Get("H Pt");
91			TH1F *zzhpt=zz.Get("H Pt");
92
93			zjhpt->SetLineColor(kGreen);
94			zzhpt->SetLineColor(kOrange+3);
95			zhhpt->SetLineColor(kRed);
96
97			zjhpt->Rebin(5);
98			zzhpt->Rebin(5);
99			zhhpt->Rebin(5);
100
101			zjhpt->Scale(1./zjhpt->Integral());
102			zzhpt->Scale(1./zzhpt->Integral());
103			zhhpt->Scale(1./zhhpt->Integral());
104
105			TCanvas c4;
106
107			zjhpt->Draw();
108			zzhpt->Draw("same");
109			zhhpt->Draw("same");
110
111			TLegend *lg4=new TLegend(0.5,0.5,0.7,0.7);
112			lg4->SetFillColor(kWhite);
113			lg4->AddEntry(zhhpt,"ZH, m(H)=120 GeV","l");
114			lg4->AddEntry(zjhpt,"Z+jets","l");
115			lg4->AddEntry(zzhpt,"ZZ","l");
116			lg4->Draw();
117
118
119			/////////////////////////////
120
121			TH1F *zhzinvmass=zh.Get("Z invmass");
122			TH1F *zjzinvmass=zj.Get("Z invmass");
123			TH1F *zzzinvmass=zz.Get("Z invmass");
124
125			zjzinvmass->SetLineColor(kGreen);
126			zzzinvmass->SetLineColor(kOrange+3);
127			zhzinvmass->SetLineColor(kRed);
128
129			zjzinvmass->Scale(1./zjzinvmass->Integral());
130			zzzinvmass->Scale(1./zzzinvmass->Integral());
131			zhzinvmass->Scale(1./zhzinvmass->Integral());
132
133			TCanvas c5;
134
135			zjzinvmass->Draw();
136			zzzinvmass->Draw("same");
137			zhzinvmass->Draw("same");
138
139			TLegend *lg5=new TLegend(0.5,0.5,0.7,0.7);
140			lg5->SetFillColor(kWhite);
141			lg5->AddEntry(zhzinvmass,"ZH, m(H)=120 GeV","l");
142			lg5->AddEntry(zjzinvmass,"Z+jets","l");
143			lg5->AddEntry(zzzinvmass,"ZZ","l");
144			lg5->Draw();
145
146
147			/////////////////////////////
148
149			TH1F *zhhinvmass=zh.Get("H invmass preselection");
150			TH1F *zjhinvmass=zj.Get("H invmass preselection");
151			TH1F *zzhinvmass=zz.Get("H invmass preselection");
152
153			zjhinvmass->SetLineColor(kGreen);
154			zzhinvmass->SetLineColor(kOrange+3);
155			zhhinvmass->SetLineColor(kRed);
156
157			zjhinvmass->Rebin(5);
158			zzhinvmass->Rebin(5);
159			zhhinvmass->Rebin(5);
160
161			zjhinvmass->Scale(1./zjhinvmass->Integral());
162			zzhinvmass->Scale(1./zzhinvmass->Integral());
163			zhhinvmass->Scale(1./zhhinvmass->Integral());
164
165			TCanvas c6;
166
167			zjhinvmass->Draw();
168			zzhinvmass->Draw("same");
169			zhhinvmass->Draw("same");
170
171			TLegend *lg6=new TLegend(0.5,0.5,0.7,0.7);
172			lg6->SetFillColor(kWhite);
173			lg6->AddEntry(zhhinvmass,"ZH, m(H)=120 GeV","l");
174			lg6->AddEntry(zjhinvmass,"Z+jets","l");
175			lg6->AddEntry(zzhinvmass,"ZZ","l");
176			lg6->Draw();
177
178			/////////////////////////////
179
180			TH1F *zhhinvmass2=zh.Get("H invmass");
181			TH1F *zjhinvmass2=zj.Get("H invmass");
182			TH1F *zzhinvmass2=zz.Get("H invmass");
183
184			zjhinvmass2->SetLineColor(kGreen);
185			zzhinvmass2->SetLineColor(kOrange+3);
186			zhhinvmass2->SetLineColor(kRed);
187
188			zjhinvmass2->Rebin(5);
189			zzhinvmass2->Rebin(5);
190			zhhinvmass2->Rebin(5);
191
192			//zjhinvmass2->Scale(1./zjhinvmass2->Integral());
193			zzhinvmass2->Scale(1./zzhinvmass2->Integral());
194			zhhinvmass2->Scale(1./zhhinvmass2->Integral());
195
196			TCanvas c7;
197
198			zjhinvmass2->Draw();
199			zzhinvmass2->Draw("same");
200			zhhinvmass2->Draw("same");
201
202			TLegend *lg7=new TLegend(0.5,0.5,0.7,0.7);
203			lg7->SetFillColor(kWhite);
204			lg7->AddEntry(zhhinvmass2,"ZH, m(H)=120 GeV","l");
205			lg7->AddEntry(zjhinvmass2,"Z+jets","l");
206			lg7->AddEntry(zzhinvmass2,"ZZ","l");
207			lg7->Draw();
208
209			/////////////////////////////
210
211			TH1F *zhdphi=zh.Get("Dphi(Z,H)");
212			TH1F *zjdphi=zj.Get("Dphi(Z,H)");
213			TH1F *zzdphi=zz.Get("Dphi(Z,H)");
214
215			zjdphi->SetLineColor(kGreen);
216			zzdphi->SetLineColor(kOrange+3);
217			zhdphi->SetLineColor(kRed);
218
219			zjdphi->Scale(1./zjdphi->Integral());
220			zzdphi->Scale(1./zzdphi->Integral());
221			zhdphi->Scale(1./zhdphi->Integral());
222
223			TCanvas c8;
224
225			zjdphi->Draw();
226			zzdphi->Draw("same");
227			zhdphi->Draw("same");
228
229			TLegend *lg8=new TLegend(0.5,0.5,0.7,0.7);
230			lg8->SetFillColor(kWhite);
231			lg8->AddEntry(zhdphi,"ZH, m(H)=120 GeV","l");
232			lg8->AddEntry(zjdphi,"Z+jets","l");
233			lg8->AddEntry(zzdphi,"ZZ","l");
234			lg8->Draw();
235
236
237
238
239			}