fix EAS lenDec bug

NuSpaceSim · Jan 24, 2024 · e6fb9d3 · e6fb9d3
1 parent 0600453
commit e6fb9d3
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/nuspacesim/simulation/eas_optical/eas.py b/src/nuspacesim/simulation/eas_optical/eas.py
@@ -63,15 +63,15 @@ def altDec(self, beta, tauBeta, tauLorentz, u=None, *args, **kwargs):
 
         u = np.random.uniform(0, 1, len(beta)) if u is None else u
 
-        tDec = -tauLorentz * mean_Tau_life * np.log(u)
+        tDec = -tauLorentz * mean_Tau_life * np.log(u)  # seconds
 
-        lenDec = tDec * tauBeta * c.value
+        lenDec = 1000.0 * tDec * tauBeta * c.value  # km
 
         altDec = np.sqrt(
             R_earth.to(units.km).value ** 2
             + lenDec**2
             + 2.0 * R_earth.to(units.km).value * lenDec * np.sin(beta)
-        )
+        )  # km
 
         altDec -= R_earth.to(units.km).value