allow non-strict monotony of efficiency vector

Author: lubynets

PWGHF/D2H/Macros/compute_fraction_cutvar.py

@@ -251,9 +251,11 @@ def main(config):
                 all_vectors_monotonous = MinimisationStatus.MonotonyViolation
                 print("\0\33[33mWARNING! main(): the raw yield uncertainties vector is not monotonous. Check the input for stability.\0\33[0m")
                 print(f"raw yield uncertainties vector elements = {unc_rawy}\n")
-            if not (np.all(effp[1:] > effp[:-1]) or np.all(effp[1:] < effp[:-1])):
+            if not (np.all(effp[1:] >= effp[:-1]) or np.all(effp[1:] <= effp[:-1])):
+                print(f"prompt efficiency vector elements = {effp}\n")
                 sys.exit("\33[31mFatal error: the prompt efficiency vector is not monotonous. Check the input. Exit.\33[0m")
-            if not (np.all(effnp[1:] > effnp[:-1]) or np.all(effnp[1:] < effnp[:-1])):
+            if not (np.all(effnp[1:] >= effnp[:-1]) or np.all(effnp[1:] <= effnp[:-1])):
+                print(f"nonprompt efficiency vector elements = {effnp}\n")
                 sys.exit("\33[31mFatal error: the nonprompt efficiency vector is not monotonous. Check the input. Exit.\33[0m")
 
         minimiser = CutVarMinimiser(rawy, effp, effnp, unc_rawy, unc_effp, unc_effnp)