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],