Particle references are const in getTime() method.

corsika/detail/framework/geometry/LeapFrogTrajectory.inl

@@ -55,7 +55,7 @@ namespace corsika {
   }
 
   template <typename Particle>
-  inline TimeType LeapFrogTrajectory::getTime(Particle& particle, double const u) const {
+  inline TimeType LeapFrogTrajectory::getTime(Particle const& particle, double const u) const {
     return particle.getTime() + getDuration(u);
   }
 

corsika/detail/framework/geometry/StraightTrajectory.inl

@@ -24,7 +24,7 @@ namespace corsika {
   }
 
   template <typename Particle>
-  inline TimeType StraightTrajectory::getTime(Particle& particle, double const u) const {
+  inline TimeType StraightTrajectory::getTime(Particle const& particle, double const u) const {
     return particle.getTime() + getDuration(u);//timeStep_ * u;
   }
 

corsika/framework/geometry/LeapFrogTrajectory.hpp

@@ -74,7 +74,7 @@ namespace corsika {
 
     ///! time at the start (u=0) or at the end (u=1) of the track of a particle
     template <typename Particle>
-    TimeType getTime(Particle& particle, double const u) const;
+    TimeType getTime(Particle const& particle, double const u) const;
 
     ///! total length along potentially bend trajectory
     LengthType getLength(double const u = 1) const;

corsika/framework/geometry/StraightTrajectory.hpp

@@ -86,7 +86,7 @@ namespace corsika {
 
     ///! time at the start (u=0) or at the end (u=1) of the track of a particle
     template <typename Particle>
-    TimeType getTime(Particle& particle, double const u) const;
+    TimeType getTime(Particle const& particle, double const u) const;
 
     ///! total length along potentially bend trajectory
     LengthType getLength(double const u = 1) const;

tests/framework/testGeometry.cpp

@@ -375,7 +375,7 @@ TEST_CASE("Geometry Trajectories") {
     StraightTrajectory base(line, t);
     CHECK(line.getPosition(t).getCoordinates() == base.getPosition(1.).getCoordinates());
     // test the getTime() method for straight trajectory
-    CHECK( base.getTime(particle1, 1) / 1_s == Approx(1));
+    CHECK(base.getTime(particle1, 1) / 1_s == Approx(1));
 
     CHECK((base.getDirection(0).getComponents(rootCS) -
            QuantityVector<dimensionless_d>{1, 0, 0})
@@ -442,9 +442,9 @@ TEST_CASE("Geometry Trajectories") {
     LeapFrogTrajectory base(pos, v0, B0, k, t);
 
     // test the getTime() method for trajectories
-    CHECK( (base.getTime(particle1,1) - t) / 1_s == 0 );
-    CHECK( base.getTime(particle1,0) / 1_s == Approx(0) );
-    CHECK( (base.getTime(particle1,0) + t) / 1_s == Approx(1e-12) );
+    CHECK((base.getTime(particle1,1) - t) / 1_s == Approx(0));
+    CHECK(base.getTime(particle1,0) / 1_s == Approx(0));
+    CHECK((base.getTime(particle1,0) + t) / 1_s == Approx(1e-12));
   }
 }