#include "Params.hh"

Trk get_elong_trk( const Trk &trk, const double fraction, const double precision ){

  Trk trk_e = trk;
  trk_e.pos = trk.pos + trk.dir * geanz.getLength( trk.E, fraction, precision);
  trk_e.t   = trk.t   +           geanz.getLength( trk.E, fraction, precision) / c_light;

  return trk_e;

}

Trk propagate( Trk& trk, double length )
{
  trk.pos = trk.pos + trk.dir * length;
  trk.t   = trk.t + length / aa::c_light;
  return trk;
}