nJs = ["NJ"] #"NJ68","NJ9i"]
#nJs = ["NJ68","NJ9i"]
basename = os.path.basename(pattern)
print basename
if basename == 'LT':
for nB in nBs:
for nJ in nJs:
patternnew = pattern + "*" + nB + "*" + nJ
basenamenew = os.path.basename(patternnew)
mask = basenamenew.replace("*", "X_")
print "Plots for", patternnew
yds = yp.YieldStore("lepYields")
yds.addFromFiles(patternnew, ("lep", "sele"))
yds.showStats()
canvs = []
# Samples
samps = ["EWK"] #,"TT","TTdiLep","TTsemiLep","WJets"]
for samp in samps:
# RCS MB
yp.colorDict[samp + "_Rcs_MB"] = yp.kBlue
hRcsMB = yp.makeSampHisto(yds, samp, "Rcs_MB",
samp + "_Rcs_MB")
#hRcsMB.SetTitle("R_{CS} (MB)")
'''
if len(sys.argv) > 1:
pattern = sys.argv[1]
print '# pattern is', pattern
else:
print "No pattern given!"
exit(0)
'''
pattern = ""
#BinMask LTX_HTX_NBX_NJX for canvas names
basename = os.path.basename(pattern)
mask = basename.replace("*", "X_")
# Define storage
yds1 = yp.YieldStore("DilepCorrs")
#path1 = "Yields/systs/btag/hadFlavour/fixXsec/allSF_noPU/meth1/merged/"
#path1 = "Yields/wData/jecv7_fixSR/lumi2p25fb/unblind/allSF_noPU/meth1A/merged/"
path1 = "Yields_DilepCorrYesNo/grid-dilep/merged/"
yds1.addFromFiles(path1, ("lep", "sele"))
yds2 = yp.YieldStore("StdWeights")
#path2 = "Yields/systs/btag/hadFlavour/fixXsec/allSF_noPU/meth2/merged/"
#path2 = "Yields/wData/jecv7/lumi2p25fb/unblind/allSF_noPU/meth1A/merged/"
path2 = "Yields_DilepCorrYesNo/grid/merged/"
yds2.addFromFiles(path2, ("lep", "sele"))
sysCols = [2, 4, 7, 8, 3, 9, 6] + range(10, 50) #[1,2,3] + range(4,10)
# canvs and hists
hists = []
storeDict = True
pckname = "pickles/sigSysts_fix_" + mask + ".pckz"
if storeDict == True and os.path.exists(pckname):
print "#Loading saved yields from pickle:", pckname
import cPickle as pickle
import gzip
yds = pickle.load(gzip.open(pckname, "rb"))
else:
print "#Reading yields from files!"
# Define storage
yds = yp.YieldStore("Sele")
#paths = []
## Add files
#scalePath = "Yields/signal/systs/scale/T1tttt/normScale2/allSF_noPU/meth1A/merged/"; paths.append(scalePath)
#isrPath = "Yields/signal/systs/ISR/T1tttt/allSF_noPU_v2/meth1A/merged/"; paths.append(isrPath)
#puPath = "Yields/signal/systs/pileup/T1tttt/allSF_noPU_fix/meth1A/merged/"; paths.append(puPath)
#btagPath = "Yields/signal/systs/btag/T1tttt/allSF_noPU_fixLepSF/meth1A/merged/"; paths.append(btagPath)
#jecPath = "Yields/signal/systs/JEC/allSF_noSF/merged/"; paths.append(jecPath)
## for central values
#normPath = "Yields/signal/fixSR/lumi2p3fb/jPt3TeV/merged/"; paths.append(normPath)
# Add all signal systematics
paths = glob('{}/signal_*/merged/'.format(pattern))
# Add central value
if len(sys.argv) > 1:
pattern = sys.argv[1]
print '# pattern is', pattern
else:
print "No pattern given!"
exit(0)
#BinMask LTX_HTX_NBX_NJX for canvas names
basename = os.path.basename(pattern)
mask = basename.replace("*","X_")
## Create Yield Storage
#ydsAntiEle = yp.YieldStore("EleAnti")
#ydsAntiEle.addFromFiles(pattern,("ele","anti"))
ydsSele = yp.YieldStore("Sele")
ydsSele.addFromFiles(pattern,("lep","sele"))
ydsSeleEle = yp.YieldStore("EleSele")
ydsSeleEle.addFromFiles(pattern,("ele","sele"))
ydsSeleMu = yp.YieldStore("MuSele")
ydsSeleMu.addFromFiles(pattern,("mu","sele"))
ydsAntiMu = yp.YieldStore("MuAnti")
ydsAntiMu.addFromFiles(pattern,("mu","anti"))
# Category
cats = ["CR_MB"]#,"CR_SB"]
for cat in cats:
'''
if len(sys.argv) > 1:
pattern = sys.argv[1]
print '# pattern is', pattern
else:
print "No pattern given!"
exit(0)
'''
pattern = ""
#BinMask LTX_HTX_NBX_NJX for canvas names
basename = os.path.basename(pattern)
mask = basename.replace("*","X_")
# Define storage
yds1 = yp.YieldStore("Method1")
#path1 = "Yields/systs/btag/hadFlavour/fixXsec/allSF_noPU/meth1/merged/"
#path1 = "Yields/wData/jecv7_fixSR/lumi2p25fb/unblind/allSF_noPU/meth1A/merged/"
path1 = "Yields/signal/fixSR/lumi2p25fb/allSF_noPU/mainNJ/merged/"
yds1.addFromFiles(path1,("lep","sele"))
yds2 = yp.YieldStore("Method2")
#path2 = "Yields/systs/btag/hadFlavour/fixXsec/allSF_noPU/meth2/merged/"
#path2 = "Yields/wData/jecv7/lumi2p25fb/unblind/allSF_noPU/meth1A/merged/"
path2 = "Yields/signal/fixSR/lumi2p25fb/allSF_noPU_cutMET3TeV/mergeNJ6p/"
yds2.addFromFiles(path2,("lep","sele"))
sysCols = [2,4,7,8,3,9,6] + range(10,50)#[1,2,3] + range(4,10)
# canvs and hists
hists = []
pckname = "pickles/bkgSysts"+mask+".pck"
if storeDict == True and os.path.exists(pckname):
print "#Loading saved yields from pickle!"
import cPickle as pickle
ydsSyst = pickle.load( open( pckname, "rb" ) )
ydsSyst.showStats()
else:
print "#Reading yields from files!"
# Define storage
ydsSyst = yp.YieldStore("Sele")
paths = []
# Add files
tptPath = "Yields/systs/topPt/MC/allSF_noPU/meth1A/merged/"; paths.append(tptPath)
puPath = "Yields/systs/PU/MC/allSF/meth1A/merged/"; paths.append(puPath)
wxsecPath = "Yields/systs/wXsec/MC/allSF_noPU/meth1A/merged/"; paths.append(wxsecPath)
ttvxsecPath = "Yields/systs/TTVxsec/MC/allSF_noPU/meth1A/merged/"; paths.append(ttvxsecPath)
wpolPath = "Yields/systs/Wpol/MC/allSF_noPU/meth1A/merged/"; paths.append(wpolPath)
dlConstPath = "Yields/systs/DLConst/merged/"; paths.append(dlConstPath)
dlSlopePath = "Yields/systs/DLSlope/merged/"; paths.append(dlSlopePath)
jerPath = "Yields/systs/JER/merged/"; paths.append(jerPath)
jerNoPath = "Yields/systs/JER_YesNo/merged/"; paths.append(jerNoPath)
btagPath = "Yields/systs/btag/hadFlavour/fixXsec/allSF_noPU/meth1A/merged/"; paths.append(btagPath)
jecPath = "Yields/systs/JEC/MC/allSF_noPU/meth1A/merged/"; paths.append(jecPath)
pckname = "pickles/bkgSysts" + mask + ".pck"
pcksig = "pickles/sigSysts_fix_" + mask + ".pckz"
if storeDict == True and os.path.exists(pckname):
print "#Loading saved yields from pickle!"
import cPickle as pickle
ydsSyst = pickle.load(open(pckname, "rb"))
#ydsSyst.showStats()
else:
print "#Reading yields from files!"
# Define storage
ydsSyst = yp.YieldStore("Sele")
paths = []
# Add files
tptPath = "Yields/systs/topPt/MC/allSF_noPU/meth1A/merged/"
paths.append(tptPath)
puPath = "Yields/systs/PU/MC/allSF/meth1A/merged/"
paths.append(puPath)
wxsecPath = "Yields/systs/wXsec/MC/allSF_noPU/meth1A/merged/"
paths.append(wxsecPath)
ttvxsecPath = "Yields/systs/TTVxsec/MC/allSF_noPU/meth1A/merged/"
paths.append(ttvxsecPath)
wpolPath = "Yields/systs/Wpol/MC/allSF_noPU/meth1A/merged/"
paths.append(wpolPath)
dlConstPath = "Yields/systs/DLConst/merged/"
paths.append(dlConstPath)
'''
if len(sys.argv) > 1:
pattern = sys.argv[1]
print '# pattern is', pattern
else:
print "No pattern given!"
exit(0)
'''
pattern = ""
#BinMask LTX_HTX_NBX_NJX for canvas names
basename = os.path.basename(pattern)
mask = basename.replace("*", "X_")
# Define storage
ydsMC = yp.YieldStore("MC")
pathMC = "Yields/wData/jecv7_fixSR/lumi2p25fb/unblind/allSF_noPU/meth1A/merged/"
ydsMC.addFromFiles(pathMC, ("lep", "sele"))
## Store dict in pickle file
storeDict = True
pckname = "pickles/allSigCentral.pckz"
if storeDict == True and os.path.exists(pckname):
print "#Loading saved yields from pickle:", pckname
import cPickle as pickle
import gzip
ydsSig = pickle.load(gzip.open(pckname, "rb"))
else:
ydsSig = yp.YieldStore("Signal")