diff --git a/benchmarks/imaging_shower_ML/scripts/epcut_scan.py b/benchmarks/imaging_shower_ML/scripts/epcut_scan.py
index 611fe55efd3f456572e4b9abadd8529c594ffb79..6250aebd2da59e3c78f5607a634e65e11df229d7 100644
--- a/benchmarks/imaging_shower_ML/scripts/epcut_scan.py
+++ b/benchmarks/imaging_shower_ML/scripts/epcut_scan.py
@@ -255,8 +255,7 @@ if __name__ == '__main__':
 
     # rejection power
     rej_pow = npi/nres
-    # NOTE: assumed half of them are pions, so need a factor of sqrt(2)
-    rej_err = ((rej_pow - npi/nres_hi)*np.sqrt(2), (npi/nres_lo - rej_pow)*np.sqrt(2))
+    rej_err = ((rej_pow - npi/nres_hi), (npi/nres_lo - rej_pow))
 
     fig, ax1 = plt.subplots(figsize=(8, 8))
     ax2 = ax1.twinx()