demote non-monotonous efficiencies to warning

Author: lubynets

PWGHF/D2H/Macros/compute_fraction_cutvar.py

@@ -279,25 +279,21 @@ def main(config):
             if not (np.all(rawy[1:] > rawy[:-1]) or np.all(rawy[1:] < rawy[:-1])):
                 all_vectors_monotonous = MinimisationStatus.MonotonyViolation
                 print(
-                    "\0\33[33mWARNING! main(): the raw yield vector is not monotonous. Check the input for stability.\0\33[0m"
+                    "\0\33[33m\nWARNING! main(): the raw yield vector is not monotonous. Check the input for stability.\n\0\33[0m"
                 )
                 print(f"raw yield vector elements = {rawy}\n")
             if not (np.all(unc_rawy[1:] > unc_rawy[:-1]) or np.all(unc_rawy[1:] < unc_rawy[:-1])):
                 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"
+                    "\0\33[33m\nWARNING! main(): the raw yield uncertainties vector is not monotonous. Check the input for stability.\n\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("\0\33[33m\nWARNING! the prompt efficiency vector is not monotonous. Check the input.\0\33[0m")
                 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("\0\33[33m\nWARNING! the nonprompt efficiency vector is not monotonous. Check the input.\0\33[0m")
                 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)
         status = minimiser.minimise_system(cfg["minimisation"]["correlated"])