dasu/UltraFastSim/Plot.C

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