yaxis.SetY1(simple.canvas.GetBottomMargin())
yaxis.SetY2(1. - simple.canvas.GetTopMargin())
simple.canvas.IsA().Destructor(simple.canvas)
simple.canvas = cnv
simple._needUpdate = False
simple.printWeb(plotDir, vardef.name)
if addLinear:
simple.ylimits = (0., -1.)
simple.minimum = -1.
vardef.ymax = -1.
simple._needUpdate = True
simple.printWeb(plotDir, vardef.name + 'Linear', logy = False)
# cleanup the mess
for obj in garbage:
obj.IsA().Destructor(obj)
cnv.IsA().Destructor(cnv)
else:
canvas.printWeb(plotDir, vardef.name)
if addLinear:
canvas.ylimits = (0., -1.)
canvas.minimum = -1.
vardef.ymax = -1.
canvas._needUpdate = True
canvas.printWeb(plotDir, vardef.name + 'Linear', logy = False)
outputFile.cd()
dmet.Write()
bmet.Write()
gmet.Write()
mcmet.Write()
dmets.append(dmet)
bmets.append(bmet)
gmets.append(gmet)
mcmets.append(mcmet)
canvas.cd()
canvas.Clear()
canvas.legend.Clear()
canvas.ylimits = (0.002, 2500)
canvas.SetLogy(True)
canvas.legend.setPosition(0.6, 0.7, 0.95, 0.9)
canvas.addStacked(bmet,
title='Background',
color=r.TColor.GetColor(0x55, 0x44, 0xff),
idx=-1)
canvas.addStacked(mcmet,
title='#gamma + jet MC',
color=r.TColor.GetColor(0xff, 0xaa, 0xcc),
idx=-1)
for probeEtaCut in probeEtaCuts:
finalCuts = allCuts + [probeEtaCut[0]]
cutString = ' && '.join(['(%s)' % c for c in finalCuts])
print cutString
label = recoilCut[1] + '_' + probeEtaCut[1]
canvas.Clear()
canvas.legend.Clear()
canvas.legend.setPosition(0.6, 0.7, 0.9, 0.9)
dataName = 'dataHist_' + label
dataHist = r.TH1D(dataName, '', 20, 0., 2.)
dataHist.Sumw2()
dataTree.Draw(varString + '>>' + dataName,
'weight * (' + cutString + ')', 'goff')
canvas.addObs(dataHist, title='Data')
for sample, color, mcTree in mcTrees:
mcName = sample + '_' + label
mcHist = r.TH1D(mcName, '', 20, 0., 2.)
mcHist.Sumw2()
mcTree.Draw(varString + '>>' + mcName,
str(lumi) + ' * weight * (' + cutString + ')', 'goff')
canvas.addStacked(mcHist, title=sample, color=color)
canvas.ylimits = (0.001, -1.)
canvas.printWeb('monophoton/phoMet', 'ratio_' + label, logy=True)
outputFile.cd()
dmet.Write()
bmet.Write()
gmet.Write()
mcmet.Write()
dmets.append(dmet)
bmets.append(bmet)
gmets.append(gmet)
mcmets.append(mcmet)
canvas.cd()
canvas.Clear()
canvas.legend.Clear()
canvas.ylimits = (0.002, 2500)
canvas.SetLogy(True)
canvas.legend.setPosition(0.6, 0.7, 0.95, 0.9)
canvas.addStacked(bmet, title = 'Background', color = r.TColor.GetColor(0x55, 0x44, 0xff), idx = -1)
canvas.addStacked(mcmet, title = '#gamma + jet MC', color = r.TColor.GetColor(0xff, 0xaa, 0xcc), idx = -1)
canvas.addObs(dmet, title = 'Data')
canvas.xtitle = canvas.obsHistogram().GetXaxis().GetTitle()
canvas.ytitle = canvas.obsHistogram().GetYaxis().GetTitle()
canvas.Update(logy = True)
gmet.Scale(1., 'width')
mcmet.Scale(1., 'width')
outputFile.cd()
dmet.Write()
bmet.Write()
gmet.Write()
mcmet.Write()
dmets.append(dmet)
bmets.append(bmet)
gmets.append(gmet)
mcmets.append(mcmet)
canvas.cd()
canvas.ylimits = (0.2, 2000.)
canvas.Clear()
canvas.legend.Clear()
canvas.ylimits = (0.2, 2000.)
canvas.SetLogy(True)
canvas.legend.setPosition(0.6, 0.7, 0.95, 0.9)
canvas.addStacked(bmet, title = 'Background', color = r.TColor.GetColor(0x55, 0x44, 0xff), idx = -1)
canvas.addStacked(mcmet, title = '#gamma + jet MC', color = r.TColor.GetColor(0xff, 0xaa, 0xcc), idx = -1)
canvas.addObs(dmet, title = 'Data')
canvas.xtitle = canvas.obsHistogram().GetXaxis().GetTitle()
for recoilCut in recoilCuts:
allCuts = cuts + [recoilCut[0]]
for probeEtaCut in probeEtaCuts:
finalCuts = allCuts + [probeEtaCut[0]]
cutString = " && ".join(["(%s)" % c for c in finalCuts])
print cutString
label = recoilCut[1] + "_" + probeEtaCut[1]
canvas.Clear()
canvas.legend.Clear()
canvas.legend.setPosition(0.6, 0.7, 0.9, 0.9)
dataName = "dataHist_" + label
dataHist = r.TH1D(dataName, "", 20, 0.0, 2.0)
dataHist.Sumw2()
dataTree.Draw(varString + ">>" + dataName, "weight * (" + cutString + ")", "goff")
canvas.addObs(dataHist, title="Data")
for sample, color, mcTree in mcTrees:
mcName = sample + "_" + label
mcHist = r.TH1D(mcName, "", 20, 0.0, 2.0)
mcHist.Sumw2()
mcTree.Draw(varString + ">>" + mcName, str(lumi) + " * weight * (" + cutString + ")", "goff")
canvas.addStacked(mcHist, title=sample, color=color)
canvas.ylimits = (0.001, -1.0)
canvas.printWeb("monophoton/phoMet", "ratio_" + label, logy=True)
