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)
