From 0a72c1adaa187c5f57c34ddbda5bd215070ba26c Mon Sep 17 00:00:00 2001 From: Javier Duarte Date: Sat, 1 Mar 2025 12:58:35 -0800 Subject: [PATCH 1/4] try 20 gev bins --- src/HH4b/postprocessing/PostProcess.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/HH4b/postprocessing/PostProcess.py b/src/HH4b/postprocessing/PostProcess.py index dc133ec9..48a9f535 100644 --- a/src/HH4b/postprocessing/PostProcess.py +++ b/src/HH4b/postprocessing/PostProcess.py @@ -1383,7 +1383,7 @@ def postprocess_run3(args): # different for glopart-v2 elif args.txbb == "glopart-v2": fom_window_by_mass["H2PNetMass"] = [110, 155] # use wider range for FoM scan - blind_window_by_mass["H2PNetMass"] = [110, 140] # only blind 3 bins + blind_window_by_mass["H2PNetMass"] = [100, 140] # only blind 2 bins # different for pnet-v12 elif args.txbb == "pnet-v12": fom_window_by_mass["H2PNetMass"] = [120, 150] @@ -1700,7 +1700,7 @@ def postprocess_run3(args): parser.add_argument( "--mass-bins", type=int, - default=10, + default=20, choices=[10, 15, 20], help="width of mass bins", ) From bf9feff96da82442962cd785d26762d1e29e1d65 Mon Sep 17 00:00:00 2001 From: Javier Duarte Date: Sat, 1 Mar 2025 13:01:38 -0800 Subject: [PATCH 2/4] try 20 gev bins --- src/HH4b/combine/run_blinded_hh4b.sh | 2 +- src/HH4b/postprocessing/PlotFits.py | 8 ++++---- 2 files changed, 5 insertions(+), 5 deletions(-) diff --git a/src/HH4b/combine/run_blinded_hh4b.sh b/src/HH4b/combine/run_blinded_hh4b.sh index 683c0cc6..27c7906d 100755 --- a/src/HH4b/combine/run_blinded_hh4b.sh +++ b/src/HH4b/combine/run_blinded_hh4b.sh @@ -157,7 +157,7 @@ fi # freeze qcd params in blinded bins setparamsblinded="" freezeparamsblinded="" -for bin in {5..7} +for bin in {2..3} do setparamsblinded+="${CMS_PARAMS_LABEL}_tf_dataResidual_Bin${bin}=0," freezeparamsblinded+="${CMS_PARAMS_LABEL}_tf_dataResidual_Bin${bin}," diff --git a/src/HH4b/postprocessing/PlotFits.py b/src/HH4b/postprocessing/PlotFits.py index 5544a5b6..f64ca5db 100644 --- a/src/HH4b/postprocessing/PlotFits.py +++ b/src/HH4b/postprocessing/PlotFits.py @@ -66,17 +66,17 @@ def plot_fits(args): fit_shape_var = ShapeVar( "H2Msd", r"Jet 2 $m_\mathrm{SD}$ (GeV)", - [16, 60, 220], + [8, 60, 220], reg=False, - blind_window=[110, 140], + blind_window=[100, 140], ) else: fit_shape_var = ShapeVar( "H2PNetMass", r"Jet 2 $m_\mathrm{reg}$ (GeV)", - [16, 60, 220], + [8, 60, 220], reg=True, - blind_window=[110, 140], + blind_window=[100, 140], ) shape_vars = [fit_shape_var] From 50b92b45b459d2431a232f892a327f1b72c99ce5 Mon Sep 17 00:00:00 2001 From: Javier Duarte Date: Sat, 1 Mar 2025 15:44:14 -0800 Subject: [PATCH 3/4] 20 GeV --- src/HH4b/plotting.py | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/HH4b/plotting.py b/src/HH4b/plotting.py index a29cadba..7b880d67 100644 --- a/src/HH4b/plotting.py +++ b/src/HH4b/plotting.py @@ -813,9 +813,8 @@ def get_variances(bg_hist): # if math.isinf(yhist[5]): # blind! - yhist[5] = 0 - yhist[6] = 0 - yhist[7] = 0 + yhist[2] = 0 + yhist[3] = 0 hep.histplot( yhist, From 590df10394b5d75cc217d4015af409f351c08890 Mon Sep 17 00:00:00 2001 From: Javier Duarte Date: Sat, 1 Mar 2025 12:55:40 -0800 Subject: [PATCH 4/4] fix plotting variance --- src/HH4b/plotting.py | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/src/HH4b/plotting.py b/src/HH4b/plotting.py index 7b880d67..179a5260 100644 --- a/src/HH4b/plotting.py +++ b/src/HH4b/plotting.py @@ -587,8 +587,12 @@ def get_variances(bg_hist): else "Stat. MC Uncertainty" ) - bg_tot = sum([hists[sample, :] for sample in bg_keys]) - bg_err_tot_mcstat = np.sqrt(bg_tot.variances()) + # this version has an issue: + # bg_tot no longer weighted, returns None for variances + # bg_tot = sum([hists[sample, :] for sample in bg_keys]) + # bg_err_tot_mcstat = np.sqrt(bg_tot.variances()) + # compute summed variance manually + bg_err_tot_mcstat = np.sqrt(sum([hists[sample, :].variances() for sample in bg_keys])) # print("mcstat ",bg_err_tot_mcstat) plot_shaded = False