Make more explicit the percentiles are for Normal dist

scikit-hep · Jan 20, 2021 · 26ada72 · 26ada72
1 parent b8207b9
commit 26ada72
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/pyhf/infer/calculators.py b/src/pyhf/infer/calculators.py
@@ -771,10 +771,11 @@ def expected_pvalues(self, sig_plus_bkg_distribution, bkg_only_distribution):
         # c.f. Issue #815, PR #817
         import numpy as np
 
-        _percentiles = [2.27501319, 15.86552539, 50.0, 84.13447461, 97.72498681]
+        # percentiles for -2, -1, 0, 1, 2 standard deviations of the Normal distribution
+        normal_percentiles = [2.27501319, 15.86552539, 50.0, 84.13447461, 97.72498681]
         pvalues_exp_band = np.percentile(
             tb.tolist(pvalues),
-            _percentiles,
+            normal_percentiles,
             axis=0,
         ).T.tolist()
         return pvalues_exp_band