Esempio n. 1
0
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')
Esempio n. 2
0
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:
Esempio n. 3
0
    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")
Esempio n. 4
0
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")