5 |
|
|
6 |
|
char fName[1024]; |
7 |
|
|
8 |
< |
std::string zhName("ZH-plots-v5.root"); |
9 |
< |
std::string zzName("ZZ-plots-v5.root"); |
10 |
< |
std::string zjName("Zmm-v5-Plots.root"); |
8 |
> |
std::string zhName("ZHmmbb-v5-plots.root"); |
9 |
> |
std::string zzName("ZZmmbb-v5-plots.root"); |
10 |
> |
std::string zjName("Zmm-v5-plots.root"); |
11 |
|
|
12 |
|
// Make plots |
13 |
|
|
85 |
|
vector<string> hNames; |
86 |
|
hNames.push_back("B-jets Et "); |
87 |
|
hNames.push_back("B-jets Eta "); |
88 |
– |
hNames.push_back("B-jets Mult. "); |
88 |
|
hNames.push_back("B-jets Phi "); |
89 |
< |
hNames.push_back("Z Pt "); |
89 |
> |
hNames.push_back("B-jets Mult. "); |
90 |
> |
hNames.push_back("H invmass preselection "); |
91 |
|
hNames.push_back("H Phi "); |
92 |
+ |
hNames.push_back("Z Pt "); |
93 |
|
hNames.push_back("H Pt "); |
93 |
– |
hNames.push_back("H invmass "); |
94 |
– |
hNames.push_back("H invmass preselection "); |
94 |
|
hNames.push_back("Dphi(Z,H) "); |
95 |
+ |
hNames.push_back("H invmass "); |
96 |
|
|
97 |
|
for(int i = 0; i < geom.size(); i++) { |
98 |
|
for(int j = 0; j < bJetOP.size(); j++) { |
110 |
|
TH1F *zhHist=zh.Get(hName.c_str()); |
111 |
|
TH1F *zjHist=zj.Get(hName.c_str()); |
112 |
|
TH1F *zzHist=zz.Get(hName.c_str()); |
113 |
< |
|
113 |
> |
|
114 |
|
if(zhHist == 0 || zjHist == 0 || zzHist == 0) { |
115 |
|
cout << "Failed to find " << hName << endl; |
116 |
|
break; |
117 |
|
} |
118 |
+ |
|
119 |
+ |
if(i > 7) { |
120 |
+ |
zhHist.Rebin(10); |
121 |
+ |
zjHist.Rebin(10); |
122 |
+ |
zzHist.Rebin(10); |
123 |
+ |
} |
124 |
|
|
125 |
|
zjHist->SetLineColor(kGreen); |
126 |
< |
zzHist->SetLineColor(kOrange+3); |
126 |
> |
zzHist->SetLineColor(kOrange); |
127 |
|
zhHist->SetLineColor(kRed); |
128 |
< |
|
129 |
< |
zjHist->Scale(1./zjHist->Integral()); |
130 |
< |
zzHist->Scale(1./zzHist->Integral()); |
131 |
< |
zhHist->Scale(1./zhHist->Integral()); |
128 |
> |
//zjHist->SetFillColor(kYellow); |
129 |
> |
zzHist->SetFillColor(kBlue); |
130 |
> |
zhHist->SetFillColor(kBlack); |
131 |
> |
|
132 |
> |
// zjHist->Scale(1./zjHist->Integral()); |
133 |
> |
// zzHist->Scale(1./zzHist->Integral()); |
134 |
> |
// zhHist->Scale(1./zhHist->Integral()); |
135 |
|
|
136 |
|
TCanvas *c = new TCanvas(); |
137 |
|
|
153 |
|
strcat(fName, ".png"); |
154 |
|
c->SaveAs(fName); |
155 |
|
|
156 |
< |
// delete lg2; |
157 |
< |
// delete c; |
158 |
< |
|
156 |
> |
TCanvas *cStack = new TCanvas(); |
157 |
> |
THStack *stack = new THStack("Stack", hName.c_str()); |
158 |
> |
stack->Add(zjHist); |
159 |
> |
stack->Add(zzHist); |
160 |
> |
stack->Add(zhHist); |
161 |
> |
gPad->SetLogy(1); |
162 |
> |
stack->Draw(); |
163 |
> |
lg2->Draw(); |
164 |
> |
fName[strlen(fName)-4] = 0; |
165 |
> |
strcat(fName, "-stack.png"); |
166 |
> |
cStack->SaveAs(fName); |
167 |
> |
gPad->SetLogy(0); |
168 |
> |
|
169 |
|
} |
170 |
|
} |
171 |
|
} |