def monov_plot_nuisance_channels_combined():
#### Channels combined
# Nuisances
for year in [2017, 2018, "combined"]:
filler = {
"year": year,
}
outdir = './plots/combined_{year}/'.format(**filler)
diffnuis_file = 'diagnostics/diffnuisances_nominal_monov_{year}.root'.format(
**filler)
plot_nuis(diffnuis_file, outdir)
dataValidation("combined", "combinedW", category, ws_file, fitdiag_file,
outdir, lumi[year], year)
# Split by flavor
dataValidation("dimuon", "singlemuon", category, ws_file, fitdiag_file,
outdir, lumi[year], year)
dataValidation("dielectron", "singleelectron", category, ws_file,
fitdiag_file, outdir, lumi[year], year)
dataValidation("singleelectron", "gjets", category, ws_file, fitdiag_file,
outdir, lumi[year], year)
dataValidation("singlemuon", "gjets", category, ws_file, fitdiag_file,
outdir, lumi[year], year)
dataValidation("dielectron", "gjets", category, ws_file, fitdiag_file,
outdir, lumi[year], year)
dataValidation("dimuon", "gjets", category, ws_file, fitdiag_file, outdir,
lumi[year], year)
plot_nuis(diffnuis_file, outdir)
### Years fit together
filler = {"category": "monovtight"}
ws_file = "root/ws_monov_nominal_tight.root".format(**filler)
fitdiag_file = 'diagnostics/fitDiagnostics_nominal_{category}_combined.root'.format(
**filler)
diffnuis_file = 'diagnostics/diffnuisances_nominal_{category}_combined.root'.format(
**filler)
model_file = "root/combined_model_monov_nominal_tight.root".format(**filler)
for year in [2017, 2018]:
outdir = './plots/combined_{YEAR}/'.format(YEAR=year)
category = 'monovtight_{YEAR}'.format(YEAR=year)
for region in regions:
plotPreFitPostFit(region, category, ws_file, fitdiag_file, outdir,