cRest=ROOT.TChain('Events')
for s in['WW','ZZ','WZ','singleTop','DY','QCD20to600','QCD600to1000','QCD1000']:
cRest.Add('/data/schoef/convertedTuples_v22/copy/'+s+'/histo_'+s+'*.root')
cData=ROOT.TChain('Events')
cData.Add('/data/schoef/convertedTuples_v22/copy/data/histo_data*.root')
metb=[150,350]
htb=[400,750]
njetb=[2,3]
lpdg=''
mTCut='mT>20&&mT<120'
template_WJets_PosPdg=getCutPlotFromChain(cWJets,'nbtags',[0,1,2,3],mTCut+'&&'+nameAndCut(metb,htb,njetb,'pos',btagRequirement='None')[1],'weight',binningIsExplicit=True,addOverFlowBin='upper')
template_WJets_NegPdg=getCutPlotFromChain(cWJets,'nbtags',[0,1,2,3],mTCut+'&&'+nameAndCut(metb,htb,njetb,'neg',btagRequirement='None')[1],'weight',binningIsExplicit=True,addOverFlowBin='upper')
template_TTJets=getCutPlotFromChain(cTTJets,'nbtags',[0,1,2,3],mTCut+'&&'+nameAndCut(metb,htb,njetb,'',btagRequirement='None')[1],'weight',binningIsExplicit=True,addOverFlowBin='upper')
template_Rest_PosPdg=getCutPlotFromChain(cRest,'nbtags',[0,1,2,3],mTCut+'&&'+nameAndCut(metb,htb,njetb,'pos',btagRequirement='None')[1],'weight',binningIsExplicit=True,addOverFlowBin='upper')
template_Rest_NegPdg=getCutPlotFromChain(cRest,'nbtags',[0,1,2,3],mTCut+'&&'+nameAndCut(metb,htb,njetb,'neg',btagRequirement='None')[1],'weight',binningIsExplicit=True,addOverFlowBin='upper')
print "Nominal yields TT:",template_TTJets.Integral(),'WJets_PosPdg',template_WJets_PosPdg.Integral(),'WJets_NegPdg',template_WJets_NegPdg.Integral()
print "Nominal yields:",'Rest_PosPdg',template_Rest_PosPdg.Integral(),'Rest_NegPdg',template_Rest_NegPdg.Integral()
#template_WJets_PosPdg.Scale(1./template_WJets_PosPdg.GetBinContent(1))
#template_WJets_NegPdg.Scale(1./template_WJets_NegPdg.GetBinContent(1))
#template_TTJets.Scale(1./template_TTJets.GetBinContent(1))
hData_PosPdg=getCutPlotFromChain(cData,'nbtags',[0,1,2,3],mTCut+'&&'+nameAndCut(metb,htb,njetb,'pos',btagRequirement='None')[1],'weight',binningIsExplicit=True,addOverFlowBin='upper')
hData_NegPdg=getCutPlotFromChain(cData,'nbtags',[0,1,2,3],mTCut+'&&'+nameAndCut(metb,htb,njetb,'neg',btagRequirement='None')[1],'weight',binningIsExplicit=True,addOverFlowBin='upper')
#metBins = [[350, -1] ]
#njetCRBins = [[3,3], [2,3]]
#njetSRBins = [[4,4], [4,-1]]
binningWPt = range(0, 1000, 100) + [1100, 2000]
wGenPt_stacks = {}
wPlusGenPt_stacks = {}
wMinusGenPt_stacks = {}
wRecoPt_stacks = {}
wPlusRecoPt_stacks = {}
wMinusRecoPt_stacks = {}
for metb in metBins:
for htb in htBins:
for njetb in njetCRBins + njetSRBins:
for pdgSign in ['', 'pos', 'neg']:
name, cut = nameAndCut(metb, htb, njetb, pdgSign)
print name, cut
wGenPt_stacks[name] = getStack(
":XXX;gen. p_{T,W} (GeV);Number of Events",
binningWPt,
cut,
signals,
varfunc=wGenPt)
wGenPt_stacks[name][0].addOverFlowBin = "upper"
allStacks.append(wGenPt_stacks[name])
wRecoPt_stacks[name] = getStack(
":XXX;reco p_{T,W} (GeV);Number of Events",
binningWPt,
cut,
#metBins = [[350, -1] ]
#njetCRBins = [[3,3], [2,3]]
#njetSRBins = [[4,4], [4,-1]]
binningWPt = range(0,1000,100)+[1100,2000]
wGenPt_stacks={}
wPlusGenPt_stacks={}
wMinusGenPt_stacks={}
wRecoPt_stacks={}
wPlusRecoPt_stacks={}
wMinusRecoPt_stacks={}
for metb in metBins:
for htb in htBins:
for njetb in njetCRBins+njetSRBins:
for pdgSign in ['', 'pos','neg']:
name, cut = nameAndCut(metb, htb, njetb, pdgSign)
print name, cut
wGenPt_stacks[name] = getStack(":XXX;gen. p_{T,W} (GeV);Number of Events", binningWPt, cut, signals, varfunc = wGenPt)
wGenPt_stacks[name][0].addOverFlowBin = "upper"
allStacks.append(wGenPt_stacks[name])
wRecoPt_stacks[name] = getStack(":XXX;reco p_{T,W} (GeV);Number of Events", binningWPt, cut, signals, varfunc = wRecoPt)
wRecoPt_stacks[name][0].addOverFlowBin = "upper"
allStacks.append(wRecoPt_stacks[name])
execfile("../../RA4Analysis/plots/simplePlotsLoopKernel.py")
wGenPt_rwHisto={}
wRecoPt_rwHisto={}
'*.root')
cData = ROOT.TChain('Events')
cData.Add('/data/schoef/convertedTuples_v22/copy/data/histo_data*.root')
metb = [150, 350]
htb = [400, 750]
njetb = [2, 3]
lpdg = ''
mTCut = 'mT>20&&mT<120'
template_WJets_PosPdg = getCutPlotFromChain(
cWJets,
'nbtags', [0, 1, 2, 3],
mTCut + '&&' +
nameAndCut(metb, htb, njetb, 'pos', btagRequirement='None')[1],
'weight',
binningIsExplicit=True,
addOverFlowBin='upper')
template_WJets_NegPdg = getCutPlotFromChain(
cWJets,
'nbtags', [0, 1, 2, 3],
mTCut + '&&' +
nameAndCut(metb, htb, njetb, 'neg', btagRequirement='None')[1],
'weight',
binningIsExplicit=True,
addOverFlowBin='upper')
template_TTJets = getCutPlotFromChain(
cTTJets,
'nbtags', [0, 1, 2, 3],
mTCut + '&&' + nameAndCut(metb, htb, njetb, '', btagRequirement='None')[1],
