zpt_cut, met_cut= '100', '100'
if not os.path.exists(outdir): os.system('mkdir '+outdir)
tag = tag0+'_'+'test'
outTag=outdir+'/'+tag
# ll in Z | ll out Z
# -------------------- M_out [35,65] U [115,120]
# eu in Z | eu out Z
# M_in (70,110)
### ----- Initialize (samples):
plotter_ll=InitializePlotter(indir, addSig=False, addData=True,doRatio=False)
plotter_eu=InitializePlotter(indir, addSig=False, addData=True,doRatio=False, doElMu=True)
setcuts=SetCuts()
cuts_ll={'inclusive': setcuts.alphaCuts(Zmass='inclusive', zpt_cut=zpt_cut, met_cut=met_cut),
'in': setcuts.alphaCuts(isll=True, Zmass='in', zpt_cut=zpt_cut, met_cut=met_cut),
'out': setcuts.alphaCuts(isll=True, Zmass='out', zpt_cut=zpt_cut, met_cut=met_cut)}
cuts_eu={'inclusive': setcuts.alphaCuts(isll=False, Zmass='inclusive', zpt_cut=zpt_cut, met_cut=met_cut),
'in' : setcuts.alphaCuts(isll=False, Zmass='in', zpt_cut=zpt_cut, met_cut=met_cut),
'out': setcuts.alphaCuts(isll=False, Zmass='out', zpt_cut=zpt_cut, met_cut=met_cut)}
ROOT.gROOT.ProcessLine('.x ../src/tdrstyle.C')
### ----- Execute (plotting):
plotter_ll.Stack.drawStack('llnunu_l1_mass', cuts_ll['inclusive'], str(lumi*1000), 10, 0.0, 200.0, titlex = "M_{Z}^{ll}", units = "GeV",
output=tag+'_mll',outDir=outdir, separateSignal=True,
drawtex="", channel="")
plotter_eu.Stack.drawStack('elmununu_l1_mass', cuts_eu['inclusive'], str(lumi*1000), 10, 0.0, 200.0, titlex = "M_{Z}^{e#mu}", units = "GeV",
output=tag+'_melmu',outDir=outdir, separateSignal=True,
# ll in Z | ll out Z
# -------------------- M_out [35,65] U [115,120]
# eu in Z | eu out Z
# M_in (70,110)
### ----- Initialize (samples):
plotter_ll = InitializePlotter(addSig=False, addData=False, doRatio=False)
plotter_eu = InitializePlotter(addSig=False,
addData=False,
doRatio=False,
doElMu=True)
setcuts = SetCuts()
print "I am cuts_ll:"
cuts_ll = setcuts.alphaCuts(inclusive=True)
print "I am cuts_eu:"
cuts_eu = setcuts.alphaCuts(isll=False, inclusive=True)
exit(0)
### ----- Execute (plotting):
print product(lsChannel, lswhichregion)
for Channel, whichregion in product(lsChannel, lswhichregion):
cuts = mycuts.abcdCuts(Channel, whichregion)
print cuts
ROOT.gROOT.ProcessLine('.x tdrstyle.C')
ROOT.gStyle.SetPadBottomMargin(0.2)
ROOT.gStyle.SetPadLeftMargin(0.15)
# M_in (70,110)
### ----- Initialize (samples):
plotter_ll = InitializePlotter(indir,
addSig=False,
addData=True,
doRatio=False)
plotter_eu = InitializePlotter(indir,
addSig=False,
addData=True,
doRatio=False,
doElMu=True)
setcuts = SetCuts()
cuts_ll = {
'inclusive':
setcuts.alphaCuts(Zmass='inclusive', zpt_cut=zpt_cut, met_cut=met_cut),
'in':
setcuts.alphaCuts(isll=True, Zmass='in', zpt_cut=zpt_cut, met_cut=met_cut),
'out':
setcuts.alphaCuts(isll=True, Zmass='out', zpt_cut=zpt_cut, met_cut=met_cut)
}
cuts_eu = {
'inclusive':
setcuts.alphaCuts(isll=False,
Zmass='inclusive',
zpt_cut=zpt_cut,
met_cut=met_cut),
'in':
setcuts.alphaCuts(isll=False, Zmass='in', zpt_cut=zpt_cut,
met_cut=met_cut),
'out':
lswhichregion=[]
if whichregion=="": lswhichregion=['SR', 'VR']
else: lswhichregion.append(whichregion)
# ll in Z | ll out Z
# -------------------- M_out [35,65] U [115,120]
# eu in Z | eu out Z
# M_in (70,110)
### ----- Initialize (samples):
plotter_ll=InitializePlotter(addSig=False, addData=False,doRatio=False)
plotter_eu=InitializePlotter(addSig=False, addData=False,doRatio=False, doElMu=True)
setcuts=SetCuts()
print "I am cuts_ll:"
cuts_ll=setcuts.alphaCuts(inclusive=True)
print "I am cuts_eu:"
cuts_eu=setcuts.alphaCuts(isll=False, inclusive=True)
exit(0)
### ----- Execute (plotting):
print product(lsChannel, lswhichregion)
for Channel, whichregion in product(lsChannel, lswhichregion):
cuts=mycuts.abcdCuts(Channel, whichregion)
print cuts
ROOT.gROOT.ProcessLine('.x tdrstyle.C')
ROOT.gStyle.SetPadBottomMargin(0.2)
ROOT.gStyle.SetPadLeftMargin(0.15)