diff --git a/corsika/detail/modules/proposal/ContinuousProcess.inl b/corsika/detail/modules/proposal/ContinuousProcess.inl
index c22226ac5f80de5067151d00d832836d481d4f27..430d779c6062549b4e44d5142aca15c8dd85886b 100644
--- a/corsika/detail/modules/proposal/ContinuousProcess.inl
+++ b/corsika/detail/modules/proposal/ContinuousProcess.inl
@@ -128,11 +128,11 @@ namespace corsika::proposal {
     auto const energy = vP.getEnergy();
     auto const energy_lim = std::max(
         energy * 0.9, // either 10% relative loss max., or
-        get_energy_threshold(
-            code) // energy thresholds globally defined for individual particles
-            *
-            0.99 // need to go 1% below global e-cut to assure removal in ParticleCut. The
-                 // 1% does not matter since at cut-time the entire energy is removed.
+        get_kinetic_energy_threshold(code) +
+            get_mass(code) // energy thresholds globally defined for individual particles
+                * 0.99     // need to go 1% below global e-cut to assure removal in
+                       // ParticleCut. The 1% does not matter since at cut-time the entire
+                       // energy is removed.
     );
 
     // solving the track integral for giving energy lim
diff --git a/corsika/stack/VectorStack.hpp b/corsika/stack/VectorStack.hpp
index 8296c928d7b4bc36e3b672114657dc8d86c126ea..e6bc6a3a426373aa835108a8787551a295b3824d 100644
--- a/corsika/stack/VectorStack.hpp
+++ b/corsika/stack/VectorStack.hpp
@@ -127,10 +127,6 @@ namespace corsika {
      *
      * @{
      */
-    DirectionVector getDirection() const {
-      return this->getMomentum() / this->getEnergy();
-    }
-
     VelocityVector getVelocity() const {
       return this->getMomentum() / this->getEnergy() * constants::c;
     }