diff --git a/flamedisx/non_asymptotic_inference.py b/flamedisx/non_asymptotic_inference.py
index ea916968..e3ace9bb 100644
--- a/flamedisx/non_asymptotic_inference.py
+++ b/flamedisx/non_asymptotic_inference.py
@@ -787,8 +787,7 @@ def get_median_disco_asymptotic(self, sigma_level=3):
 
             disco_sig_curves = np.stack(disco_sig_curves, axis=0)
             median_disco_sigs = [np.median(disco_sigs) for disco_sigs in disco_sig_curves]
-
-            median_crossing_point = self.critical_disco_value(median_disco_sigs, mus, 3)
+            median_crossing_point = self.critical_disco_value(median_disco_sigs, mus, sigma_level)
             medians[signal_source] = median_crossing_point
 
         return medians