elif plotDir and not vardef.fullyBlinded(): canvas.addObs(obshist, title = plotConfig.obs.title) if vardef.name == 'count': printCounts(counters, plotConfig) elif vardef.name == args.bbb: print 'Bin-by-bin yield for variable', args.bbb printBinByBin(counters, plotConfig) elif vardef.name == args.chi2: print 'Chi2 for variable', args.chi2 printChi2(counters, plotConfig) if plotDir and vardef.name != 'count': canvas.xtitle = obshist.GetXaxis().GetTitle() canvas.ytitle = obshist.GetYaxis().GetTitle() canvas.selection = vardef.formSelection(plotConfig) if vardef.logy is None: logy = True addLinear = True else: logy = vardef.logy addLinear = False canvas.Update(logy = logy, ymax = vardef.ymax) if vardef.fullyBlinded(): # remove ratio pad. Hack to use SimpleCanvas interface
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) canvas.printWeb('monophoton/gjetsTFactor', 'distributions' + region) ########################################### ####### Transfer Factors ################## ########################################### methods = [('Data', gmets), ('MC', mcmets)] scanvas = SimpleCanvas(lumi=lumi) tfacts = []
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) canvas.printWeb('monophoton/gjetsTFactor', 'distributions'+region) ########################################### ####### Transfer Factors ################## ########################################### methods = [ ('Data', gmets), ('MC', mcmets) ] scanvas = SimpleCanvas(lumi = lumi) tfacts = []
canvas.Clear() canvas.legend.Clear() canvas.legend.setPosition(0.6, 0.7, 0.9, 0.9) dataName = 'dataHist'+skim dataHist = varDef.makeHist(dataName) if varDef.unit != '': dataHist.GetXaxis().SetTitle(varDef.title + ' (' + varDef.unit +')') else: dataHist.GetXaxis().SetTitle(varDef.title) dataCutString = cutString + ' && !(run > %s)' % config.runCutOff dataTree.Draw(varDef.expr+'>>'+varDef.name+'-'+dataName, 'weight * ('+cutString+')', 'goff') print dataHist.Integral() canvas.addObs(dataHist, title = 'Data '+skimm) for sample, color, mcTree in mcTrees: mcName = sample+skimm mcHist = varDef.makeHist(mcName) mcTree.Draw(varDef.expr+'>>'+varDef.name+'-'+mcName, str(lumi)+' * weight * ('+cutString+')', 'goff') canvas.addStacked(mcHist, title = sample, color = color) # canvas.ylimits = (0.001, -1.) canvas.xtitle = varDef.title canvas.ytitle = 'Events' canvas.printWeb('monophoton/phoMet/'+skim, skimm[:3]+'_'+jetsCut[0]+'_'+sign+'_'+varDef.name, logy = False)
title='Fake E_{T}^{miss}', color=ROOT.kRed, drawOpt='HIST') canvas.addStacked(sig, title='H(125)', color=ROOT.kBlue, drawOpt='HIST') bkgTrue.Scale(1., 'width') fakeTrue.Scale(fakeNorm / fakeTrue.GetSumOfWeights(), 'width') sigTrue.Scale(sigScale, 'width') canvas.addSignal(bkgTrue, title='True SM Bkgd.', color=ROOT.kGreen + 2, drawOpt='HIST') canvas.addSignal(fakeTrue, title='True Fake E_{T}^{miss}', color=ROOT.kRed + 2, drawOpt='HIST') canvas.addSignal(sigTrue, title='True H(125)', color=ROOT.kBlue + 2, drawOpt='HIST') canvas.addObs(data, drawOpt='EP') canvas.title = '#sigma#timesBR = %.2f, N_{fake} = %.0f' % (originalMu * sigScale, fakeNorm) canvas.xtitle = 'm_{T} (GeV)' canvas.ytitle = 'Events / GeV' canvas.printWeb('monophoton/fakemet', name, logy=False)