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)