make it compile and work

Documentation/Examples/cascade_example.cc

@@ -12,6 +12,7 @@
 #include <corsika/cascade/Cascade.h>
 #include <corsika/process/ProcessSequence.h>
 #include <corsika/process/stack_inspector/StackInspector.h>
+#include <corsika/process/tracking_line/TrackingLine.h>
 
 #include <corsika/setup/SetupStack.h>
 #include <corsika/setup/SetupTrajectory.h>
@@ -40,7 +41,7 @@ public:
   ProcessSplit() {}
 
   template <typename Particle>
-  double MinStepLength(Particle& p) const {
+  double MinStepLength(Particle& p, setup::Trajectory&) const {
 
     // beam particles for sibyll : 1, 2, 3 for p, pi, k
     // read from cross section code table
@@ -91,7 +92,7 @@ public:
     return next_step;
   }
 
-  template <typename Particle, typename Trajectory, typename Stack>
+  template <typename Particle, typename Stack>
   EProcessReturn DoContinuous(Particle&, Trajectory&, Stack&) const {
     // corsika::utls::ignore(p);
     return EProcessReturn::eOk;
@@ -237,12 +238,13 @@ double s_rndm_(int&) {
 
 int main() {
 
-  stack_inspector::StackInspector<setup::Stack, setup::Trajectory> p0(true);
+  tracking_line::TrackingLine<setup::Stack> tracking;
+  stack_inspector::StackInspector<setup::Stack> p0(true);
   ProcessSplit p1;
   const auto sequence = p0 + p1;
   setup::Stack stack;
 
-  corsika::cascade::Cascade EAS(sequence, stack);
+  corsika::cascade::Cascade EAS(tracking, sequence, stack);
 
   stack.Clear();
   auto particle = stack.NewParticle();

Framework/Cascade/testCascade.cc

@@ -13,6 +13,7 @@
 
 #include <corsika/process/ProcessSequence.h>
 #include <corsika/process/stack_inspector/StackInspector.h>
+#include <corsika/process/tracking_line/TrackingLine.h>
 
 #include <corsika/stack/super_stupid/SuperStupidStack.h>
 
@@ -72,26 +73,10 @@ public:
 private:
 };
 
-template <typename Stack>
-class NewtonTracking { // naja.. not yet
-  typedef typename Stack::ParticleType Particle;
-
-public:
-  void Init() {}
-  corsika::setup::Trajectory GetTrack(Particle& p) {
-    corsika::geometry::Vector<SpeedType::dimension_type> v = p.GetDirection();
-    corsika::geometry::Line traj(p.GetPosition(), v);
-    {
-      CoordinateSystem& rootCS = RootCoordinateSystem::GetInstance().GetRootCS();
-      cout << v.GetComponents(rootCS) << endl;
-    }
-    return corsika::geometry::Trajectory<corsika::geometry::Line>(traj, 100_ns);
-  }
-};
 
 TEST_CASE("Cascade", "[Cascade]") {
 
-  NewtonTracking<setup::Stack> tracking;
+  tracking_line::TrackingLine<setup::Stack> tracking;
 
   stack_inspector::StackInspector<setup::Stack> p0(true);
   ProcessSplit p1;
@@ -108,7 +93,7 @@ TEST_CASE("Cascade", "[Cascade]") {
   particle.SetEnergy(E0);
   particle.SetPosition(Point(rootCS, {0_m, 0_m, 10_km}));
   particle.SetMomentum(
-      corsika::stack::super_stupid::MomentumVector(rootCS, {0_GeV, 0_GeV, -1_MeV}));
+		       corsika::stack::super_stupid::MomentumVector(rootCS, {0*newton*second, 0*newton*second, -1*newton*second}));
   EAS.Init();
   EAS.Run();
 

Processes/StackInspector/StackInspector.h

@@ -16,9 +16,6 @@
 
 #include <corsika/setup/SetupTrajectory.h>
 
-// namespace corsika::setup {
-// class BaseTrajectory;
-//}
 
 namespace corsika::process {