New getTime() method for start and end of track

corsika/detail/framework/geometry/LeapFrogTrajectory.inl

@@ -54,6 +54,11 @@ namespace corsika {
     return step / 2 * (correction + 1);
   }
 
+  template <typename Particle>
+  inline TimeType LeapFrogTrajectory::getTime(Particle& particle, double const u) const {
+    return particle.getTime() + getDuration(u);
+  }
+
   inline LengthType LeapFrogTrajectory::getLength(double const u) const {
     return getDuration(u) * initialVelocity_.getNorm();
   }

corsika/detail/framework/geometry/StraightTrajectory.inl

@@ -23,6 +23,11 @@ namespace corsika {
     return u * timeStep_;
   }
 
+  template <typename Particle>
+  inline TimeType StraightTrajectory::getTime(Particle& particle, double const u) const {
+    return particle.getTime() + getDuration(u);//timeStep_ * u;
+  }
+
   inline LengthType StraightTrajectory::getLength(double const u) const {
     if (timeLength_ == 0_s) return 0_m;
     if (timeStep_ == std::numeric_limits<TimeType::value_type>::infinity() * 1_s)

corsika/framework/geometry/LeapFrogTrajectory.hpp

@@ -72,6 +72,10 @@ namespace corsika {
     ///! duration along potentially bend trajectory
     TimeType getDuration(double const u = 1) const;
 
+    ///! 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;
+
     ///! total length along potentially bend trajectory
     LengthType getLength(double const u = 1) const;
 

corsika/framework/geometry/StraightTrajectory.hpp

@@ -84,6 +84,10 @@ namespace corsika {
     ///! duration along potentially bend trajectory
     TimeType getDuration(double const u = 1) const;
 
+    ///! 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;
+
     ///! total length along potentially bend trajectory
     LengthType getLength(double const u = 1) const;