diff --git a/src/nuspacesim/simulation/eas_optical/shower_properties.py b/src/nuspacesim/simulation/eas_optical/shower_properties.py
index 9cf9bed..818a84b 100644
--- a/src/nuspacesim/simulation/eas_optical/shower_properties.py
+++ b/src/nuspacesim/simulation/eas_optical/shower_properties.py
@@ -136,7 +136,7 @@ def gaisser_hillas_particle_count(
     Xmax = Xm + (70.0 * np.log10(Eshow * 1.0e-8))
     λ = p1 + p2 * Xmask + p3 * Xmask * Xmask
     λ[λ > 100.0] = 100.0
-    λ[λ < 1.0] = 1.0
+    λ[λ < 1.0e-5] = 1.0e-5
 
     # Parametric Form Parameters
     x = (Xmask - X0) / λ