TTJets.setFillProperties(1001,ROOT.kGreen-5) qcdPlotters=[] for sample in [ "QCD_Pt1000to1400_50ns","QCD_Pt120to170_50ns","QCD_Pt1400to1800_50ns","QCD_Pt170to300_50ns","QCD_Pt1800to2400_50ns","QCD_Pt2400to3200_50ns","QCD_Pt300to470_50ns","QCD_Pt3200_50ns","QCD_Pt470to600_50ns","QCD_Pt600to800_50ns","QCD_Pt800to1000_50ns","QCD_Pt80to120_50ns" ]: qcdPlotters.append(TreePlotter('samples/'+sample+'.root','tree')) qcdPlotters[-1].setupFromFile('samples/'+sample+'.pck') qcdPlotters[-1].addCorrectionFactor('xsec','xsec',0.0,'lnN') qcdPlotters[-1].addCorrectionFactor('genWeight','genWeight',0.0,'lnN') QCD = MergedPlotter(qcdPlotters) QCD.setFillProperties(1001,ROOT.kAzure-9) #QCD.setFillProperties(1001,ROOT.kGreen-5) singleMu = TreePlotter('samples/SingleMuon_Run2015B.root','tree') singleMu.addCorrectionFactor('hlt','(HLT_ISOMU_50ns==1||(HLT_MU==1))',0.0,'lnN') singleEle = TreePlotter('samples/SingleElectron_Run2015B.root','tree') singleEle.addCorrectionFactor('hlt','((HLT_ISOELE_50ns==1||(HLT_ELE==1))&&(!(HLT_ISOMU_50ns==1||(HLT_MU==1))))',0.0,'lnN') jetHT = TreePlotter('samples/JetHT_Run2015B.root','tree') #jetHT.addCorrectionFactor('hlt','((HLT_JJ==1||HLT_HT==1)&&(!(HLT_ISOMU_50ns==1||(HLT_MU==1)||(HLT_ISOELE_50ns==1||(HLT_ELE==1)))))',0.0,'lnN')
QCDPtPlotters = getPlotters(nonResTemplate, False) QCDhtPlotters = getPlotters('QCD_HT', False) QCDherwigPlotters = getPlotters('QCD_Pt-', False) QCD = MergedPlotter(QCDPtPlotters) QCDht = MergedPlotter(QCDhtPlotters) QCDherwig = MergedPlotter(QCDherwigPlotters) # # # DATAPlotters = getPlotters('SingleJet',True) # # data=MergedPlotter(DATAPlotters) # # SigPlotters = getPlotters('QstarToQW_2000',False) # sig = MergedPlotter(SigPlotters) # # QCD.setFillProperties(1001, 921) QCDht.setLineProperties(1, 633, 2) QCDht.setFillProperties(1001, 0) # QCDherwig.setLineProperties(1, 434, 2) QCDherwig.setFillProperties(1001, 0) # # sig.setFillProperties(3001,ROOT.kGreen+2) #Stack # jjStack = StackPlotter() # jjStack.addPlotter(QCD,"QCD","QCD multijet","background") # jjStack.addPlotter(data,"data_obs","Data","data") for p in purities:
vvPlotters[-1].addCorrectionFactor('genWeight', 'tree') vvPlotters[-1].addCorrectionFactor('puWeight', 'tree') vv = MergedPlotter(vvPlotters) WWLNUJJ = TreePlotter('samples/VBF_RadionToWW_narrow_3000.root', 'tree') WWLNUJJ.setupFromFile('samples/VBF_RadionToWW_narrow_3000.pck') WWLNUJJ.addCorrectionFactor('xsec', 'tree') WWLNUJJ.addCorrectionFactor(100, 'flat') WWLNUJJ.addCorrectionFactor('genWeight', 'tree') WWLNUJJ.addCorrectionFactor('puWeight', 'tree') WWLNUJJ.setFillProperties(0, ROOT.kWhite) WWLNUJJ.setLineProperties(1, ROOT.kRed, 3) #Fill properties WJets.setFillProperties(1001, ROOT.kAzure - 9) vv.setFillProperties(1001, ROOT.kOrange) ttO.setFillProperties(1001, ROOT.kSpring - 5) ttM.setFillProperties(1001, ROOT.kTeal - 1) #ZJets.setFillProperties(1001,ROOT.kAzure+5) #GJets.setFillProperties(1001,ROOT.kYellow) QCD.setFillProperties(1001, ROOT.kGray) data = MergedPlotter(dataPlotters) #Stack for lnu+J lnujjStack = StackPlotter() #lnujjStack.addPlotter(QCD,"QCD","QCD multijet","background") lnujjStack.addPlotter(ttO, "tt", "t#bar{t} (other)", "background") lnujjStack.addPlotter(WJets, "WJets", "W+Jets", "background") lnujjStack.addPlotter(ttM, "tt", "t#bar{t} (W)", "background")
from CMGTools.VVResonances.plotting.StackPlotter import StackPlotter #create the W+jets plotters wjPlotters=[] for sample in ['WJetsToLNu_HT100to200','WJetsToLNu_HT200to400','WJetsToLNu_HT400to600','WJetsToLNu_HT600toInf']: wjPlotters.append(TreePlotter('samples/'+sample+'.root','tree')) wjPlotters[-1].setupFromFile('samples/'+sample+'.pck') wjPlotters[-1].addCorrectionFactor('xsec','xsec',0.0,'lnN') wjPlotters[-1].addCorrectionFactor('genWeight','genWeight',0.0,'lnN') WJets = MergedPlotter(wjPlotters) WJets.setFillProperties(1001,ROOT.kAzure-9) RSGWWLNuQQ = TreePlotter('samples/RSGravToWWToLNQQ_kMpl01_2500.root','tree') RSGWWLNuQQ.setupFromFile('samples/RSGravToWWToLNQQ_kMpl01_2500.pck') RSGWWLNuQQ.setFillProperties(0,ROOT.kWhite) RSGWWLNuQQ.setLineProperties(1,ROOT.kOrange+10,3) RSGWWLNuQQ.addCorrectionFactor('xsec',0.001,0.0,'lnN') #Stack vvStack = StackPlotter() vvStack.addPlotter(WJets,"W+jets","W+Jets","background") vvStack.addPlotter(RSGWWLNuQQ,"RSG2000","RSGWW #rightarrow l#nu QQ","signal")