variables = {"tau_nc_ratio": {"Title": "Neutral-Charged asymmetry", "Log": False, "VariableWidth": False}}
legendPositions = {"tau_nc_ratio": [0.3, 0.2, 0.7, 0.4]}
plots = []
for var, options in variables.items():
for iso in isos:
for sample in samples:
histoIso = histosIso[sample]
histoInvertIso = histosInvertIso[sample]
plot = ComparisonPlot()
plot.name = var + "_" + sample + "_" + iso[0] + "_" + iso[1]
plot.plotDir = plotDir
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
#
cfgIso = deepcopy(configIso)
cfgIso.unitScaling = True
cfgIso.binWidthScaling = options["VariableWidth"]
cfgIso.xTitle = options["Title"]
plot.addHisto(histoIso[0], histoIso[1].format(VAR=var, ISO=iso[0]), cfgIso)
#
cfgInvertIso = deepcopy(configInvertIso)
cfgInvertIso.unitScaling = True
cfgInvertIso.binWidthScaling = options["VariableWidth"]
cfgInvertIso.xTitle = options["Title"]
plot.addHisto(histoInvertIso[0], histoInvertIso[1].format(VAR=var, ISO=iso[1]), cfgInvertIso)
plot.plot()
plots.append(plot)
plots = []
## Compare uncertainties with several numbers of tries
for name,sys in systematics.items():
for sample,histo in histos.items():
plot = ComparisonPlot()
plot.name = "statUncertainties_mvis_"+sample+'_'+name
plot.logy = False
plot.legendPosition = legendPosition
#
config1 = copy(configRawStat)
config1.xTitle = "m_{vis} [GeV]"
config1.yTitle = 'Relative uncertainty'
config1.legend = "Stat. unc."
histo1 = errorHisto(histo[0],name+'/'+histo[1])
plot.addHisto(histo1, config1)
#
config2 = copy(configFactorStat4)
config2.xTitle = "m_{vis} [GeV]"
config2.legend = "Factor unc. 10"
config2.yTitle = 'Relative uncertainty'
histo2 = fakeFactorErrorHisto(histo[0],histo[1], sys[10:20])
plot.addHisto(histo2, config2)
#
config3 = copy(configFactorStat3)
config3.xTitle = "m_{vis} [GeV]"
config3.legend = "Factor unc. 50"
config3.yTitle = 'Relative uncertainty'
histo3 = fakeFactorErrorHisto(histo[0],histo[1], sys[0:50])
plot.addHisto(histo3, config3)
#
plots = []
## Compare uncertainties with several numbers of tries
for name, sys in systematics.items():
for sample, histo in histos.items():
plot = ComparisonPlot()
plot.name = "statUncertainties_mvis_" + sample + '_' + name
plot.logy = False
plot.legendPosition = legendPosition
#
config1 = copy(configRawStat)
config1.xTitle = "m_{vis} [GeV]"
config1.yTitle = 'Relative uncertainty'
config1.legend = "Stat. unc."
histo1 = errorHisto(histo[0], name + '/' + histo[1])
plot.addHisto(histo1, config1)
#
config2 = copy(configFactorStat4)
config2.xTitle = "m_{vis} [GeV]"
config2.legend = "Factor unc. 10"
config2.yTitle = 'Relative uncertainty'
histo2 = fakeFactorErrorHisto(histo[0], histo[1], sys[10:20])
plot.addHisto(histo2, config2)
#
config3 = copy(configFactorStat3)
config3.xTitle = "m_{vis} [GeV]"
config3.legend = "Factor unc. 50"
config3.yTitle = 'Relative uncertainty'
histo3 = fakeFactorErrorHisto(histo[0], histo[1], sys[0:50])
plot.addHisto(histo3, config3)
#
"tau_decayMode":[0.2, 0.7, 0.55, 0.9],
"tau_pdgId":[0.2, 0.7, 0.55, 0.9],
}
plots = []
for var,options in variables.items():
for nameRef,histoRef in histosRef.items():
for name,histo in histos.items():
plot = ComparisonPlot()
plot.name = var+"_"+nameRef+"_"+name
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
configRef = copy(configZMuMu)
configRef.unitScaling = True
configRef.binWidthScaling = options["VariableWidth"]
configRef.xTitle = options["Title"]
plot.addHisto(histoRef[0],histoRef[1].format(VAR=var), configRef)
config = copy(configMuTau)
config.unitScaling = True
config.binWidthScaling = options["VariableWidth"]
config.xTitle = options["Title"]
config.legend = name+" "+config.legend
plot.addHisto(histo[0],histo[1].format(VAR=var), config)
plot.plot()
plots.append(plot)
plots = []
for var, options in variables.items():
for selection in selectionsOSSSInvertIso:
plot = ComparisonPlot()
plot.name = var + "_" + selection + '_MTbins'
plot.plotDir = 'plots/WJets/'
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
for bin in mtbins:
config = copy(HistoConfigs.configMTbins[bin])
config.unitScaling = True
config.binWidthScaling = options["VariableWidth"]
config.xTitle = options["Title"]
#config.legend = 'MT bin {}'.format(bin)
plot.addHisto(inputFileName,
histoNames.format(VAR=var, SEL=selection, BIN=bin),
config)
plot.plot()
plots.append(plot)
## Compare distributions in MT bins for OS/SS regions
plots = []
for var, options in variables.items():
for bin in mtbins:
## Invert iso
plot = ComparisonPlot()
plot.name = var + "_MTBin_" + str(bin) + '_InvertIso_OS_SS'
plot.plotDir = 'plots/WJets/'
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
for selection in selectionsOSSSInvertIso:
## Variables and selections
variables = {
"tau_nc_ratio":{"Title":"Neutral-Charged asymmetry", "Log":False, "VariableWidth":False},
}
legendPositions = {
"tau_nc_ratio":[0.3, 0.2, 0.7, 0.4],
}
plots = []
for var,options in variables.items():
plot = ComparisonPlot()
plot.name = var+"_backgrounds"
plot.plotDir = plotDir
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
for sample,histo in histos.items():
config = copy(configBackgrounds[sample])
config.unitScaling = True
config.binWidthScaling = options["VariableWidth"]
config.xTitle = options["Title"]
plot.addHisto(histo[0],histo[1].format(VAR=var), config)
plot.plot()
plots.append(plot)
## Compare OS/SS distributions in different MT bins
plots = []
for var,options in variables.items():
for selection in selectionsOSSSInvertIso:
plot = ComparisonPlot()
plot.name = var+"_"+selection+'_MTbins'
plot.plotDir = 'plots/WJets/'
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
for bin in mtbins:
config = copy(HistoConfigs.configMTbins[bin])
config.unitScaling = True
config.binWidthScaling = options["VariableWidth"]
config.xTitle = options["Title"]
#config.legend = 'MT bin {}'.format(bin)
plot.addHisto(inputFileName,histoNames.format(VAR=var,SEL=selection,BIN=bin), config)
plot.plot()
plots.append(plot)
## Compare distributions in MT bins for OS/SS regions
plots = []
for var,options in variables.items():
for bin in mtbins:
## Invert iso
plot = ComparisonPlot()
plot.name = var+"_MTBin_"+str(bin)+'_InvertIso_OS_SS'
plot.plotDir = 'plots/WJets/'
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
for fakefactor,sample in fakefactors.items():
for name,histo in histos.items():
filename = histo['File'].format(DIR=histoDir,SAMPLE=sample,VERSION=version)
histoTrue = histo['HistoTrue'].format(VAR=var)
histoEst = histo['HistoEst'].format(FACTOR=fakefactor,VAR=var)
#
plot = ComparisonPlot()
plot.plotDir = 'results/'
plot.name = var+"_"+name+'_'+fakefactor+"_"+sample
plot.xRange = [0.,5.9999]
plot.yRatioRange = [0.5,4.5]
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
configTrue = copy(HistoConfigs.configTrue)
#configTrue.unitScaling = True
configTrue.binWidthScaling = options["VariableWidth"]
configTrue.xTitle = options["Title"]
configTrue.legend = 'True'
plot.addHisto(filename,histoTrue, configTrue)
configEst = copy(HistoConfigs.configEst)
#configEst.unitScaling = True
configEst.binWidthScaling = options["VariableWidth"]
configEst.xTitle = options["Title"]
configEst.legend = 'Est.'
plot.addHisto(filename,histoEst, configEst)
plot.plot()
plot.plot_ratio()
plots.append(plot)
}
plots = []
for var, options in variables.items():
for iso in isos:
for sample in samples:
histoIso = histosIso[sample]
histoInvertIso = histosInvertIso[sample]
plot = ComparisonPlot()
plot.name = var + "_" + sample + "_" + iso[0] + "_" + iso[1]
plot.plotDir = plotDir
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
#
cfgIso = deepcopy(configIso)
cfgIso.unitScaling = True
cfgIso.binWidthScaling = options["VariableWidth"]
cfgIso.xTitle = options["Title"]
plot.addHisto(histoIso[0], histoIso[1].format(VAR=var, ISO=iso[0]),
cfgIso)
#
cfgInvertIso = deepcopy(configInvertIso)
cfgInvertIso.unitScaling = True
cfgInvertIso.binWidthScaling = options["VariableWidth"]
cfgInvertIso.xTitle = options["Title"]
plot.addHisto(histoInvertIso[0],
histoInvertIso[1].format(VAR=var,
ISO=iso[1]), cfgInvertIso)
plot.plot()
plots.append(plot)
for fakefactor, sample in fakefactors.items():
for name, histo in histos.items():
filename = histo['File'].format(DIR=histoDir,
SAMPLE=sample,
VERSION=version)
histoTrue = histo['HistoTrue'].format(VAR=var)
histoEst = histo['HistoEst'].format(FACTOR=fakefactor, VAR=var)
#
plot = ComparisonPlot()
plot.plotDir = 'results/'
plot.name = var + "_" + name + '_' + fakefactor + "_" + sample
plot.xRange = [0., 5.9999]
plot.yRatioRange = [0.5, 4.5]
plot.logy = options["Log"]
plot.legendPosition = legendPositions[var]
configTrue = copy(HistoConfigs.configTrue)
#configTrue.unitScaling = True
configTrue.binWidthScaling = options["VariableWidth"]
configTrue.xTitle = options["Title"]
configTrue.legend = 'True'
plot.addHisto(filename, histoTrue, configTrue)
configEst = copy(HistoConfigs.configEst)
#configEst.unitScaling = True
configEst.binWidthScaling = options["VariableWidth"]
configEst.xTitle = options["Title"]
configEst.legend = 'Est.'
plot.addHisto(filename, histoEst, configEst)
plot.plot()
plot.plot_ratio()
plots.append(plot)
plots = []
## Compare uncertainties with several numbers of tries
for name, sys in systematics.items():
for sample, histo in histos1.items():
plot = ComparisonPlot()
plot.name = "statUncertainties_mvis_" + sample + "_" + name
plot.logy = False
plot.legendPosition = legendPosition
#
config1 = copy(configRawStat)
config1.xTitle = "m_{vis} [GeV]"
config1.yTitle = "Relative uncertainty"
config1.legend = "Stat. unc. in anti-iso"
histo1 = errorHisto([h[0] for h in histo], [name + "/" + h[1] for h in histo])
plot.addHisto(histo1, config1)
#
# config2 = copy(configFactorStat4)
# config2.xTitle = "m_{vis} [GeV]"
# config2.legend = "Factor unc. 10"
# config2.yTitle = 'Relative uncertainty'
# histo2 = fakeFactorErrorHisto([h[0] for h in histo],[h[1] for h in histo] sys[10:20])
# plot.addHisto(histo2, config2)
#
config3 = copy(configFactorStat3)
config3.xTitle = "m_{vis} [GeV]"
config3.legend = "Factor unc., 50 toys"
config3.yTitle = "Relative uncertainty"
histo3 = fakeFactorErrorHisto([h[0] for h in histo], [h[1] for h in histo], sys[0:50])
plot.addHisto(histo3, config3)
#
