def __init__(self, useRoot=False ,basedir=None ,do_ddQCD=True):
if basedir is None:
self.basedir="/net/scratch_cms/institut_3a/materok/wprime/okt14/crab_projects/merged"
#self.basedir="/net/scratch_cms/institut_3a/materok/wprime/run1/crab_projects/merged"
else:
self.basedir=basedir
self.useRoot=useRoot
self.xs= ConfigObj("xsFile13TeV_25ns_wprime.cfg")
self.lumi=0
self.lumi+= 5882 #RunB
self.lumi+= 2646 #RunC
self.lumi+= 4353 #RunD
self.lumi+= 4049 #RunE
self.lumi+= 3160 #RunF
#self.lumi= 11048 #RunC + RunD + RunE
#self.lumi= 12881 #RunB + RunC + RunF
self.bghists=HistStorage(self.xs,self.lumi,path=self.basedir,useRoot=self.useRoot)
self.bghists.setDataDriven("dataDrivenQCD")
self.bghists.setDataDriven("dataDrivenTree")
self.bghists.setDataDriven("dataDrivenTreept")
self.bghists.setDataDriven("dataDrivenTreedecay_mode")
self.bghists.setDataDriven("dataDrivenTreeeta")
self.bghists.setDataDriven("dataDrivenTreelepton_n")
self.bghists.setDataDriven("dataDrivenTree15_pt_base")
self.bghists.setDataDriven("dataDrivenTree15_pt_baseFit")
self.bghists.setDataDriven("dataDrivenTree15_pt_decay_base")
self.bghists.setDataDriven("dataDrivenTree15_pt_decay_baseFit")
self.bghists.setDataDriven("dataDrivenTree15_pt_eta_baseFit")
self.bghists.setDataDriven("dataDrivenTree15_pt_eta_base")
self.bghists.setDataDriven("dataDrivenTree15_lepton_n_base")
self.bghists.setDataDriven("dataDrivenTree15_decay_mode_base")
#do_ddQCD=False
self.bglist=OrderedDict()
self.bglist["Diboson"]=[
"WW_TuneCUETP8M1_13TeV-pythia8",
"WZ_TuneCUETP8M1_13TeV-pythia8",
"ZZ_TuneCUETP8M1_13TeV-pythia8"
]
if not do_ddQCD:
self.bglist[r"$Z\rightarrow \nu\nu$"]=[
"ZJetsToNuNu_HT-100To200_13TeV-madgraph",
"ZJetsToNuNu_HT-200To400_13TeV-madgraph",
"ZJetsToNuNu_HT-400To600_13TeV-madgraph",
"ZJetsToNuNu_HT-600To800_13TeV-madgraph",
"ZJetsToNuNu_HT-800To1200_13TeV-madgraph",
"ZJetsToNuNu_HT-1200To2500_13TeV-madgraph",
"ZJetsToNuNu_HT-2500ToInf_13TeV-madgraph"
]
self.bglist[r"$\mathsf{Z/\gamma}$"]=[
"DYJetsToLL_M-5to50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8"
]
self.bglist[r'$\mathsf{t\bar{t}}$']=[
"TT_TuneCUETP8M1_13TeV-powheg-pythia8",
"ST_t-channel_antitop_4f_leptonDecays_13TeV-powheg-pythia8_TuneCUETP8M1",
"ST_t-channel_top_4f_leptonDecays_13TeV-powheg-pythia8_TuneCUETP8M1",
"ST_tW_antitop_5f_inclusiveDecays_13TeV-powheg-pythia8_TuneCUETP8M1",
"ST_tW_top_5f_inclusiveDecays_13TeV-powheg-pythia8_TuneCUETP8M1"
]
if do_ddQCD:
self.bglist["multijet"]=[
"dataDrivenQCD",
]
if not do_ddQCD:
self.bglist["QCD"]=[
"QCD_HT200to300_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"QCD_HT300to500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"QCD_HT500to700_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"QCD_HT700to1000_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"QCD_HT1000to1500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"QCD_HT1500to2000_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"QCD_HT2000toInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8"
]
#self.bglist[r'$\mathsf{W\rightarrow l \nu}$ $\mathsf{M_{W}<200 GeV}$']=[
#"WJetsToLNu_HT-100To200_13TeV_MG",
#"WJetsToLNu_HT-200To400_13TeV_MG",
#"WJetsToLNu_HT-400To600_13TeV_MG",
##"WJetsToLNu_HT-600To800_13TeV_MG",
#"WJetsToLNu_HT-600ToInf_13TeV_MG",
##"WJetsToLNu_HT-800To1200_13TeV_MG",
##"WJetsToLNu_HT-1200To2500_13TeV_MG",
##"WJetsToLNu_HT-2500ToInf_13TeV_MG",
##"WJetsToLNu_13TeV-FXFX_MC",
##"WJetsToLNu_13TeV_MG",
#"WJetsToLNu_13TeV_AM",
#]
#self.bglist[r'$\mathsf{W\rightarrow l \nu}$ MCatNLO']=[
###"WJetsToLNu_13TeV-FXFX_MC",
#"WJetsToLNu_13TeV_AM",
###"WJetsToLNu_13TeV-FXFX_MC",
##"WJetsToLNu_13TeVMLM_MG",
#]
self.bglist["WJets"]=[
"WJetsToLNu_HT-100To200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"WJetsToLNu_HT-200To400_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"WJetsToLNu_HT-400To600_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"WJetsToLNu_HT-600To800_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"WJetsToLNu_HT-800To1200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"WJetsToLNu_HT-1200To2500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8",
"WJetsToLNu_HT-2500ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8"
]
"""
self.bglist[r'$\mathsf{W\rightarrow \tau \nu}$']=[
"WToTauNu_M-100_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-200_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-500_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-1000_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-2000_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-3000_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-4000_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-5000_TuneCUETP8M1_13TeV-pythia8-tauola",
"WToTauNu_M-6000_TuneCUETP8M1_13TeV-pythia8-tauola"
]
"""
self.bghists.additionalWeight["dataDrivenQCD"]=0.938091502614
self.bghists.additionalWeight["TT_13TeV_AM"]=1./3.
self.bghists.additionalWeight["TT_13TeV_ext1_AM"]=1./3.
self.bghists.additionalWeight["TTJets_13TeV_AM"]=1./3.
self.bghists.additionalWeight["TT_13TeV_ext3_PH"]=1./2.
self.bghists.additionalWeight["TT_13TeV_ext4_PH"]=1./2.
self.bghists.additionalWeight["TT_Mtt-700to1000_13TeV_ext1_PH"]=1./2.
self.bghists.additionalWeight["TT_Mtt-700to1000_13TeV_PH"]=1./2.
self.bghists.additionalWeight["TT_Mtt-1000toInf_13TeV_ext1_PH"]=1./2.
self.bghists.additionalWeight["TT_Mtt-1000toInf_13TeV_ext2_PH"]=1./2.
self.bghists.additionalWeight["WWToLNuQQ_13TeV_ext1_PH"]=1./2.
self.bghists.additionalWeight["WWToLNuQQ_13TeV_PH"]=1./2.
#bghists.additionalWeight["WJetsToLNu_HT-100To200_13TeV_MG"]=0.5
#bghists.additionalWeight["TT_13TeV_PH"]=0.5
#bghists.additionalWeight["TT_13TeV_ext3_PH"]=0.5
#bghists.additionalWeight["ST_t-channel_4f_leptonDecays_13TeV_AM"]=0.5
#bghists.additionalWeight["DYJetsToLL_M-50_13TeV_ext1_AM"]=0.5
#bghists.additionalWeight["ST_t-channel_4f_leptonDecays_13TeV_ext1_AM"]=0.5
self.bgcolorList={}
self.bgcolorList[r'$\mathsf{W\rightarrow l \nu}$']=getRGBTColor(861)
self.bgcolorList[r'$\mathsf{W\rightarrow \tau \nu}$']=getRGBTColor(861)
self.bgcolorList["WJets"]="lightblue"
self.bgcolorList["W"]=getRGBTColor(861)
self.bgcolorList["W on shell"]="lightblue"
self.bgcolorList[r'$\mathsf{W\rightarrow l \nu}$ MCatNLO']="brown"
self.bgcolorList[r'$\mathsf{W\rightarrow l \nu}$ $\mathsf{M_{W}<200 GeV}$']="lightblue"
self.bgcolorList["W HT"]="orange"
self.bgcolorList["multijet"]=getRGBTColor(213)
self.bgcolorList["multijet MC"]=getRGBTColor(216)
self.bgcolorList["QCD from data"]=getRGBTColor(213)
self.bgcolorList["QCD jet"]=getRGBTColor(213)
self.bgcolorList["QCD MC"]="brown"
self.bgcolorList["QCD"]="brown"
self.bgcolorList["Z (inv)"]="darkblue"
self.bgcolorList["Top"]=getRGBTColor(616)
self.bgcolorList["$t\bar{t}$"]=getRGBTColor(616)
self.bgcolorList[r'$\mathsf{t\bar{t}}$']="orange"
self.bgcolorList["Diboson"]=getRGBTColor(413)
self.bgcolorList[r"$\mathsf{Z/\gamma}$"]=getRGBTColor(810)
self.bgcolorList["DY"]=getRGBTColor(810)
self.bgcolorList[r"$\mathsf{\gamma jets}$"]="yellow"
self.bgcolorList[r"$Z\rightarrow \nu\nu$"]="yellow"
# dummy signal
#setup signals
self.sghist=HistStorage(self.xs,self.lumi,path=self.basedir,useRoot=self.useRoot)
self.sgList=OrderedDict()
self.sgList["W' M=1TeV"]=["WprimeToTauNu_M-1000_TuneCUETP8M1_13TeV-pythia8-tauola"]
self.sgcolorlist={}
self.sgcolorlist["W' M=1TeV"]="darkred"
self.dat_hist=HistStorage(self.xs,self.lumi,path=self.basedir,isData=True,useRoot=self.useRoot)
#self.DataFiles=["TauRunC"]
#self.DataFiles=["TauRunD"]
#self.DataFiles=["TauRunE"]
#self.DataFiles=["Tau_Run2016E_PromptReco_v2","Tau_Run2016C_PromptReco_v2","Tau_Run2016D_PromptReco_v2"]
self.DataFiles=["Tau_Run2016B_PromptReco_v2","Tau_Run2016C_PromptReco_v2",
"Tau_Run2016D_PromptReco_v2","Tau_Run2016E_PromptReco_v2","Tau_Run2016F_PromptReco_v1"]
# output2015_7_30_14_0 -> data
# output2015_7_31_10_28 -> data
#
# output2015_7_31_11_21 -> data
#
# output2015_8_3_17_34 -> data
# output2015_8_4_17_55 -> MC
#
basedir="/net/scratch_cms/institut_3a/erdweg/public/13TeV_rpv/BACKGROUND/merged"
lumi = 40.8 + 16.3 + 25.028
# lumi = 10
lumisc = float(lumi)/float(1000)
xs= ConfigObj("/home/home1/institut_3a/keller/TAPAS/plotting/xs_Phys14.cfg")
bghists=HistStorage(xs,lumi,xstype=None,path=basedir)
bglist=OrderedDict()
bglist["DY"]=[
'ZToEE_NNPDF30_13TeV_M_50_120_PH-skimid1741',
'ZToEE_NNPDF30_13TeV_M_120_200_PH-skimid1736',
'ZToEE_NNPDF30_13TeV_M_200_400_PH-skimid1896',
'ZToEE_NNPDF30_13TeV_M_400_800_PH-skimid1730',
'ZToEE_NNPDF30_13TeV_M_800_1400_PH-skimid1740',
'ZToEE_NNPDF30_13TeV_M_1400_2300_PH-skimid1794',
'ZToEE_NNPDF30_13TeV_M_3500_4500_PH-skimid1770',
'ZToEE_NNPDF30_13TeV_M_6000_Inf_PH-skimid1726',
'ZToMuMu_NNPDF30_13TeV_M_50_120_PH-skimid1752',
'ZToMuMu_NNPDF30_13TeV_M_200_400_PH-skimid1727',
'ZToMuMu_NNPDF30_13TeV_M_400_800_PH-skimid1728',
if binning[-1] > plot_range[1]:
binning[-1] = plot_range[1]
res_file.Close()
if debug:
print(binning)
return binning
#
basedir="/disk1/erdweg/television/BACKGROUND/merged/"
# lumi = 40.8 + 16.3 + 25.028
lumi = 1620.957
# lumi = 981.610
lumisc = float(lumi)/float(1000)
xs= ConfigObj("/disk1/erdweg/plotting/xs_Phys14.cfg")
bghists=HistStorage(xs,lumi,xstype=None,path=basedir)
bglist=OrderedDict()
bglist["DY"]=[
'ZToEE_NNPDF30_13TeV_M_50_120_PH-skimid3044',
'ZToEE_NNPDF30_13TeV_M_120_200_PH-skimid3023',
'ZToEE_NNPDF30_13TeV_M_200_400_PH-skimid2867',
'ZToEE_NNPDF30_13TeV_M_400_800_PH-skimid2872',
'ZToEE_NNPDF30_13TeV_M_800_1400_PH-skimid3012',
'ZToEE_NNPDF30_13TeV_M_1400_2300_PH-skimid2854',
'ZToEE_NNPDF30_13TeV_M_2300_3500_PH-skimid2853',
'ZToEE_NNPDF30_13TeV_M_3500_4500_PH-skimid2857',
'ZToEE_NNPDF30_13TeV_M_4500_6000_PH-skimid2856',
'ZToEE_NNPDF30_13TeV_M_6000_Inf_PH-skimid2892',
'ZToMuMu_NNPDF30_13TeV_M_50_120_PH-skimid3048',
def main():
#get the list of bg and data files from the unusal config
sampleinfos, basedir = getSampleInfo("201705_mergedext.config")
lumi = sampleinfos.values()[-1]["lumi"]
useRoot = False
#we will use the same colors as the .c plotter for now
colors = [100, 90, 80, 70, 60, 95, 85, 75, 65, 55]
#reverse order and change to RGB
colors = [getRGBTColor(i) for i in colors[::-1]]
#make storage containers for the bgs and data
bghists = HistStorage(xs=sampleinfos,
xstype="BSM3G",
lumi=lumi,
path=basedir,
useRoot=useRoot)
dat_hist = HistStorage(path=basedir, isData=True, useRoot=useRoot)
#make some lists that we will put into the containers
colorList = {}
bglist = OrderedDict()
for sample in sampleinfos:
info = sampleinfos[sample]
sample = sample.replace(".root", "")
if info["type"] == "bg":
shortName = info["output"].replace(".root", "")
if "#" in shortName:
shortName = "$\mathsf{" + shortName.replace("#", "\\") + "}$"
if shortName not in bglist:
bglist[shortName] = [sample]
colorList[shortName] = colors.pop()
else:
bglist[shortName].append(sample)
#for data directly add the file
if info["type"] == "data":
dat_hist.addFile(sample)
for bg in bglist:
print("bglist['%s']=[" % (bg))
for ibg in bglist[bg]:
print("'ibg',")
print("]")
#now actually add the backgrounds
bghists.addFileList(bglist)
bghists.colorList = colorList
hists = [
"NDiMuonCombinations/MHT",
"NDiMuonCombinations/Meff",
]
histContainer = HistStorageContainer(bg=bghists, data=dat_hist)
#histContainer=HistStorageContainer(bg=bghists,sg=sghist,data=dat_hist)
binning = {
"_pt": 25,
"_mt": 25,
"_met": 25,
"MC_W_m_Gen": 10,
"_jet_ht": 10,
}
xranges = {
"_pt": [0, 1000],
"_jet_ht": [0, 1600],
"_mt": [150, 1600],
"boson_qt": [0, 600],
"_met": [0, 1300],
"_ET_MET": [0, 4],
"_DeltaPhi": [0.2, math.pi],
"MC_W_m_Gen": [0, 1300],
"relIso": [0, 0.05],
}
yranges = {
"_ET_MET": [1.01e-3, 1e8],
"_DeltaPhi": [1.01e-3, 1e8],
}
bghists.initStyle(style="bg")
#sghist.initStyle(style="sg")
for hist in hists:
##use different out dirs
dir_name = "test"
if not os.path.exists(dir_name):
os.mkdir(dir_name)
if "prof" in hist:
hist_style = sc.style_container(kind="Lines",
style='CMS',
useRoot=useRoot,
cmsPositon="upper left",
lumi=lumi,
cms=13)
else:
hist_style = sc.style_container(style='CMS',
useRoot=useRoot,
cmsPositon="upper left",
lumi=lumi,
cms=13)
#change the style if you want
#hist_style._error_bands_fcol=["grey","yellow"]
#hist_style._error_bands_ecol=["lightgrey","yellow"]
#hist_style._error_bands_alph=0.4
#hist_style.Set_error_line_width(0.0000000001)
hist_style.Set_poisson_error_off()
#hist_style.Set_minor_ticks()
#hist_style.Set_n_legend_columns(2)
#hist_style.Set_xerr()
#hist_style.Set_do_overflowbin()
#get the actual histograms
histContainer.getHist(hist)
#programmable binning
binf = getDictValue(hist, binning)
if binf is not None:
if isinstance(binf, list):
histContainer.rebin(vector=binf)
else:
histContainer.rebin(width=binf)
name = hist.replace("/", "")
#test = plotter(hist=histContainer.getBGList(),data_hist=histContainer.getData(),data=True,style=hist_style)
test = plotter(hist=histContainer.getBGList(), style=hist_style)
test.Add_data(histContainer.getData())
test.Add_plot('Ratio', pos=1, height=15)
test.Add_plot('DiffRatio_width_increase', pos=2, height=15)
mplt = test.create_plot()
test.SavePlot('%s/%s.png' % (dir_name, name))
test.SavePlot('%s/%s.pdf' % (dir_name, name))
def main():
#get the list of bg and data files from the unusal config
sampleinfos, basedir = getSampleInfo("201705_mergedext.config")
lumi = sampleinfos.values()[-1]["lumi"]
useRoot = False
#we will use the same colors as the .c plotter for now
colors = [100, 90, 80, 70, 60, 95, 85, 75, 65, 55]
#reverse order and change to RGB
colors = [getRGBTColor(i) for i in colors[::-1]]
#make storage containers for the bgs and data
bghists = HistStorage(xs=sampleinfos,
xstype="BSM3G",
lumi=lumi,
path=basedir,
useRoot=useRoot)
dat_hist = HistStorage(path=basedir, isData=True, useRoot=useRoot)
#make some lists that we will put into the containers
colorList = {}
bglist = OrderedDict()
#the idea is to have a readable dict:
bglist['DY+Jets'] = [
'DYJetsToLL_M-50_HT-70To100',
'DYJetsToLL_M-50_HT-600To800',
'DYJetsToLL_M-50_HT-400To600',
'DYJetsToLL_M-50_HT-800To1200',
'DYJetsToLL_M-50_HT-1200To2500',
'DYJetsToLL_M-50_HT-2500ToInf',
'DYJetsToLL_M-50_HT-200To400',
'DYJetsToLL_M-50_HT-0To70',
'DYJetsToLL_M-50_HT-100To200',
]
bglist['VV'] = [
'ZZ_TuneCUETP8M1_13TeV',
'WZ_TuneCUETP8M1_13TeV',
'WW_TuneCUETP8M1_13TeV',
]
bglist['SingleTop'] = [
'ST_tW_antitop_5f_inclusiveDecays_13TeV',
'ST_t-channel_top_4f_inclusiveDecays_13TeV',
'ST_tW_top_5f_inclusiveDecays_13TeV',
'ST_t-channel_antitop_4f_inclusiveDecays_13TeV',
]
bglist['$\mathsf{t\bar{t}}$'] = [
'TT_TuneCUETP8M2T4_13TeV',
]
bglist['W+Jets'] = [
'WJetsToLNu_HT-0To70',
]
##you can reverse the order
bglist = OrderedDict(reversed(list(bglist.items())))
#now actually add the backgrounds
bghists.addFileList(bglist)
bghists.colorList = colorList
hists = [
"NDiMuonCombinations/MHT",
"NDiMuonCombinations/Meff",
]
histContainer = HistStorageContainer(bg=bghists, data=dat_hist)
#histContainer=HistStorageContainer(bg=bghists,sg=sghist,data=dat_hist)
bghists.initStyle(style="bg")
for hist in hists:
#this is the magic line
histContainer.getHist(hist)
#now you can do what you want with the hists
allbg = sum(histContainer.getBGList())
allbg.Draw()
for bg in histContainer.getBGList():
bg.Draw("same")
data = histContainer.getData()
data.Draw("same p")
def main():
basedir="/disk1/erdweg/out/output2015_2_13_15_47/merged//"
lumi=1000
xs= ConfigObj("/disk1/erdweg/plotting/xs_Phys14.cfg")
bghists=HistStorage(xs,lumi,path=basedir)
# bghists.setDataDriven("dataDrivenQCD")
bglist=OrderedDict()
bglist["PU20bx25"]=['ZprimeToTauTau_M-5000_PU20bx25_PHYS14',
'ZprimeToTauTau_M-1000_PU20bx25_PHYS14',
]
bglist["AVE30BX50"]=['ZprimeToTauTau_M-1000_AVE30BX50_PHYS14',
'ZprimeToTauTau_M-5000_AVE30BX50_PHYS14',
]
bglist["PU40bx25"]=['ZprimeToTauTau_M-1000_PU40bx25_PHYS14',
'ZprimeToTauTau_M-5000_PU40bx25_PHYS14',
]
colorList={}
colorList["PU20bx25"]="lightblue"
colorList["AVE30BX50"]="lightgreen"
colorList["PU40bx25"]="darkmagenta"
bghists.addFileList(bglist)
dat_hist=HistStorage(xs,lumi,path=basedir,isData=True)
# dat_hist.addFile("Data")
hists=[#'Taus/h1_1_Tau_eta',
#'Taus/h1_1_Tau_phi',
#'Taus/h1_2_Tau_pt',
#'Taus/h1_Tau_pt_resolution_0_500',
'Eles/h1_Ele_pt_resolution_0_500',
'Eles/h1_Ele_pt_resolution_500_1000',
#'Eles/h1_2_Ele_pt',
'Muons/h1_Muon_pt_resolution_0_500',
'Muons/h1_Muon_pt_resolution_500_1000',
#'Muons/h1_2_Muon_pt',
#'emu/Stage_0/h1_0_emu_Mass',
#'mutau/Stage_0/h1_0_mutau_Mass',
#'mutau/Stage_6/h1_6_mutau_Mass',
#'Ctr/h1_Ctr_Vtx_unweighted',
#'Ctr/h1_Ctr_Vtx_weighted',
#'Ctr/h1_Ctr_Vtx_emu_unweighted',
#'Ctr/h1_Ctr_Vtx_emu_weighted',
#'Ctr/h1_Ctr_HT',
#'Ctr/h1_Ctr_pT_hat',
#'emu/Stage_0/h1_0_emu_Mass_resolution',
#'mutau/Stage_0/h1_0_mutau_Mass_resolution'
]
sghist=HistStorage(xs,lumi,path=basedir)
# sgName="$\mathsf{W' \, M=2.3\,TeV \cdot 0.02}$"
# sghist.additionalWeight={"WprimeToTauNu_M-2300_TuneZ2star_8TeV-pythia6-tauola_Summer12_DR53X-PU_S10_START53_V7A-v1SIM":0.02}
# sghist.addAllFiles(tag="WprimeToTauNu_M-2300",joinName=sgName)
# sghist.colorList={sgName :"darkred"}
histContainer=HistStorageContainer(sg=bghists,data=dat_hist,bg=sghist)
binning={
"Taus/h1_2_Tau_pt":10,
# "emu/Stage_0/h1_0_emu_Mass":range(200,300,20)+range(300,400,50)+range(400,1600,100)+range(1600,2000,200),
# "mutau/Stage_0/h1_0_mutau_Mass":range(200,300,20)+range(300,400,50)+range(400,1600,100)+range(1600,2000,200),
# "mutau/Stage_6/h1_6_mutau_Mass":range(200,300,20)+range(300,400,50)+range(400,1600,100)+range(1600,2000,200),
"_met_et":30,
}
xranges={
"Taus/h1_2_Tau_pt":[0,2000],
"emu/Stage_0/h1_0_emu_Mass":[0,2000],
"mutau/Stage_0/h1_0_mutau_Mass":[0,2000],
"mutau/Stage_6/h1_6_mutau_Mass":[0,2000],
"emu/Stage_0/h1_0_emu_Mass_resolution":[-2,2],
"mutau/Stage_0/h1_0_mutau_Mass_resolution":[-2,2],
"Taus/h1_Tau_pt_resolution_0_500":[-2,2],
"Eles/h1_Ele_pt_resolution_0_500":[-2,2],
"Eles/h1_Ele_pt_resolution_500_1000":[-2,2],
"Muons/h1_Muon_pt_resolution_0_500":[-2,2],
"Muons/h1_Muon_pt_resolution_500_1000":[-2,2],
"Ctr/h1_Ctr_Vtx_unweighted":[0,40],
"Ctr/h1_Ctr_Vtx_weighted":[0,40],
"Ctr/h1_Ctr_Vtx_emu_unweighted":[0,40],
"Ctr/h1_Ctr_Vtx_emu_weighted":[0,40],
}
bghists.initStyle(style="sg",colors=colorList)
# sghist.initStyle(style="sg")
for hist in hists:
print("Now plotting: " + hist)
histContainer.getHist(hist)
binf=getDictValue(hist,binning)
if binf is not None:
# dat_hist.rebin(width=binf)
bghists.rebin(width=binf)
bghists.colorList=colorList
# bghists.setStyle(bgcolors=colorList)
# dat_hist.getHistList()[0].SetTitle("data")
hist_style = sc.style_container(style = 'CMS', useRoot = False,cms=13,lumi=19700)
hist_style.Set_additional_text('plot dummy')
test = plotter(sig=bghists.getHistList(),style=hist_style)
# test.Add_data(dat_hist.getHistList()[0])
# test.Add_plot('Signi',pos=0, height=15)
#test.Add_plot('DiffRatio',pos=1, height=15)
# test.Add_plot('Signi',pos=2, height=15)
# test.Add_plot('Diff',pos=1, height=15)
# test.Add_plot('Ratio',pos=2, height=15)
#test.Add_error_hist([sys_hist_2,sys_hist], band_center = 'ref')
test._cms_text_x = 0.12
test._cms_text_y = 0.91
#
# mxrange=getDictValue(hist,xranges)
if hist in xranges.keys():
test.Set_axis(xmin=xranges[hist][0],xmax=xranges[hist][1],ymin=1e-3,ymax=1e3)
name=hist.replace("/","")
test.create_plot()
plot_gauss_fit(test.Get_axis1(), bghists.getHistList()[0], xmin = -0.5, xmax = 0.5)
# test.show_fig()
test.SavePlot('plots/%s.pdf'%(name))
return 42
#
# output2015_7_30_14_0 -> data
# output2015_7_31_10_28 -> data
#
# output2015_7_31_11_21 -> data
#
# output2015_8_3_17_34 -> data
# output2015_8_4_17_55 -> MC
# basedir="/disk1/erdweg/out/output2015_7_21_14_52/merged/"
basedir="/disk1/erdweg/out/output2015_8_19_18_16/merged/"
lumi = 40.028
lumisc = float(lumi)/float(1000)
xs= ConfigObj("/disk1/erdweg/plotting/xs_Phys14.cfg")
bghists=HistStorage(xs,lumi,xstype=None,path=basedir)
bglist=OrderedDict()
# bglist["DY"]=[
# 'DYJetsToLL_M-10to50_13TeV-FXFX_MC',
# 'DYJetsToLL_M-50_13TeV-FXFX_MC',
# 'DYJetsToLL_M-100to200_13TeV-FXFX_MC',
# 'DYJetsToLL_M-200to400_13TeV-FXFX_MC',
# 'DYJetsToLL_M-400to500_13TeV-FXFX_MC',
# 'DYJetsToLL_M-500to700_13TeV-FXFX_MC',
# 'DYJetsToLL_M-700to800_13TeV-FXFX_MC',
# 'DYJetsToLL_M-800to1000_13TeV-FXFXRaw_MC',
# 'DYJetsToLL_M-1000to1500_13TeV-FXFX_MC',
# 'DYJetsToLL_M-1500to2000_13TeV-FXFX_MC',
# 'DYJetsToLL_M-2000to3000_13TeV-FXFX_MC',
if binning[-1] > plot_range[1]:
binning[-1] = plot_range[1]
res_file.Close()
if debug:
print(binning)
return binning
#
basedir="/home/home1/institut_3a/mukherjee/PhysicsAnalysis/CMSSW_7_4_0/src/TAPAS/PxlAnalyzer/JobTry/BACKGROUND/merged/"
#lumi = 40.8 + 16.3 + 25.028
lumi = 1000
lumisc = float(lumi)/float(1)
xs= ConfigObj("/home/home1/institut_3a/mukherjee/TAPAS/DibosonPlots/xs_13TeV.cfg")
#####bghists=HistStorage(xs,lumi,xstype=None,path=basedir)
bghists=HistStorage(xs,lumi,path=basedir)
bglist=OrderedDict()
#bglist["DY"]=[
# 'ZToEE_NNPDF30_13TeV_M_50_120_PH-skimid1741',
# 'ZToEE_NNPDF30_13TeV_M_120_200_PH-skimid1736',
# 'ZToEE_NNPDF30_13TeV_M_200_400_PH-skimid1896',
# 'ZToEE_NNPDF30_13TeV_M_400_800_PH-skimid1730',
# 'ZToEE_NNPDF30_13TeV_M_800_1400_PH-skimid1740',
# 'ZToEE_NNPDF30_13TeV_M_1400_2300_PH-skimid1794',
# 'ZToEE_NNPDF30_13TeV_M_3500_4500_PH-skimid1770',
# 'ZToEE_NNPDF30_13TeV_M_6000_Inf_PH-skimid1726',
# 'ZToMuMu_NNPDF30_13TeV_M_50_120_PH-skimid1752',
# 'ZToMuMu_NNPDF30_13TeV_M_200_400_PH-skimid1727',
# 'ZToMuMu_NNPDF30_13TeV_M_400_800_PH-skimid1728',
