227 |
|
allStack.Add(histos[i]) |
228 |
|
|
229 |
|
d1 = ROOT.TH1F('noData','noData',nBins[v],xMin[v],xMax[v]) |
230 |
< |
datatitle='Data ' |
230 |
> |
datatitle='Data' |
231 |
> |
addFlag = '' |
232 |
> |
if 'Zee' in datanames and 'Zmm' in datanames: |
233 |
> |
addFlag = 'Z(l^{-}l^{+})H(b#bar{b})' |
234 |
> |
elif 'Zee' in datanames: |
235 |
> |
addFlag = 'Z(e^{-}e^{+})H(b#bar{b})' |
236 |
> |
elif 'Zmm' in datanames: |
237 |
> |
addFlag = 'Z(#mu^{-}#mu^{+})H(b#bar{b})' |
238 |
|
for i in range(0,len(datas)): |
239 |
|
d1.Add(datas[i],1) |
233 |
– |
if i ==0: |
234 |
– |
datatitle+=datanames[i] |
235 |
– |
else: |
236 |
– |
datatitle=datatitle+ ' + '+datanames[i] |
240 |
|
print "\033[1;32m\n\tDATA integral = %s\033[1;m"%d1.Integral() |
241 |
|
flow = d1.GetEntries()-d1.Integral() |
242 |
|
if flow > 0: |
276 |
|
Ymax = max(allStack.GetMaximum(),d1.GetMaximum())*1.7 |
277 |
|
if log: |
278 |
|
allStack.SetMinimum(0.05) |
279 |
< |
maxval = max( allStack.GetMaximum(),noStack.GetMaximum()) |
277 |
< |
Ymax = maxval*ROOT.TMath.Power(10,1.6*(ROOT.TMath.Log(1.6*(maxval/0.1))/ROOT.TMath.Log(10)))*(0.6*0.1) |
279 |
> |
Ymax = Ymax*ROOT.TMath.Power(10,1.6*(ROOT.TMath.Log(1.6*(maxval/0.1))/ROOT.TMath.Log(10)))*(0.6*0.1) |
280 |
|
ROOT.gPad.SetLogy() |
281 |
|
allStack.SetMaximum(Ymax) |
282 |
|
c.Update() |
291 |
|
|
292 |
|
tPrel = myText("CMS Preliminary",0.17,0.88,1.04) |
293 |
|
tLumi = myText("#sqrt{s} = %s, L = %s fb^{-1}"%(anaTag,(float(lumi_data)/1000.)),0.17,0.83) |
294 |
+ |
tAddFlag = myText(addFlag,0.17,0.78) |
295 |
|
|
296 |
|
unten.cd() |
297 |
|
ROOT.gPad.SetTicks(1,1) |