kraus
--- a/src/Algorithms/ParallelTTracker.cpp
+++ b/src/Algorithms/ParallelTTracker.cpp
+        if (pathLength_m > stepSizes_m.getZStop()) {
+            ++ stepSizes_m;

-            dtAllTracks_m.pop();
-            localTrackSteps_m.pop();
-            zStop_m.pop();
+            if (stepSizes_m.reachedEnd()) return;

+            dtCurrentTrack_m = stepSizes_m.getdT();
            changeDT();
        }

        double speed = euclidean_norm(itsBunch_m->RefPartP_m * Physics::c / Util::getGamma(itsBunch_m->RefPartP_m));
        if (std::abs(pathLength_m - zstart_m) <=  0.5 * itsBunch_m->getdT() * speed) {
-            double tau = (pathLength_m - zstart_m) / speed;
+            double tau = (zstart_m - pathLength_m) / speed;