compare plots

andreypz · Aug 26, 2024 · cad663a · cad663a
1 parent fc8185d
commit cad663a
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/pocket_coffea/utils/plot_utils.py b/pocket_coffea/utils/plot_utils.py
@@ -15,6 +15,7 @@
 import matplotlib.pyplot as plt
 from matplotlib.pyplot import cm
 import mplhep as hep
+from mplhep.error_estimation import poisson_interval
 from cycler import cycler
 
 from omegaconf import OmegaConf
@@ -696,9 +697,9 @@ def get_datamc_ratio(self, cat):
                 hden = hden * (1./den_integral)        
 
         ratio = hnum.values()/hden.values()
-        # Total uncertainy of num x den :
-        #ratio_variance = np.power(ratio,2)*( hnum.variances()*np.power(hnum.values(), -2) + hden.variances()*np.power(hden.values(), -2))
-        # Only the uncertainty of hnum (DATA) propagated:
+        # Total uncertainy propagation of num / den :
+        # ratio_variance = np.power(ratio,2)*( hnum.variances()*np.power(hnum.values(), -2) + hden.variances()*np.power(hden.values(), -2))
+        # Only the uncertainty of num (DATA) propagated:
         ratio_variance = hnum.variances()*np.power(hden.values(), -2)
 
         print(dir(hep))