fix plotting variance (#263)

LPC-HH · Mar 1, 2025 · 7602989 · 7602989
1 parent a98aaf3
commit 7602989
Showing 1 changed file with 6 additions and 2 deletions.
diff --git 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