added TrackWriter to cascade_example for visualization

Documentation/Examples/CMakeLists.txt

@@ -28,6 +28,7 @@ target_link_libraries (cascade_example SuperStupidStack CORSIKAunits CORSIKAlogg
    ProcessSibyll
   CORSIKAcascade
   ProcessStackInspector
+  ProcessTrackWriter
   CORSIKAprocesses
   CORSIKAparticles
   CORSIKAgeometry

Documentation/Examples/cascade_example.cc

@@ -25,6 +25,8 @@
 #include <corsika/process/sibyll/Decay.h>
 #include <corsika/process/sibyll/Interaction.h>
 
+#include <corsika/process/track_writer/TrackWriter.h>
+
 #include <corsika/units/PhysicalUnits.h>
 
 #include <corsika/random/RNGManager.h>
@@ -219,14 +221,15 @@ int main() {
   corsika::process::sibyll::Interaction sibyll;
   corsika::process::sibyll::Decay decay;
   ProcessCut cut(8_GeV);
+  corsika::process::TrackWriter::TrackWriter trackWriter("tracks.dat");
 
   // assemble all processes into an ordered process list
-  const auto sequence = /*p0 <<*/ sibyll << decay << cut;
+  const auto sequence = /*p0 <<*/ sibyll << decay << cut << trackWriter;
 
   // setup particle stack, and add primary particle
   setup::Stack stack;
   stack.Clear();
-  const hep::EnergyType E0 = 1_TeV;
+  const hep::EnergyType E0 = 100_TeV;
   {
     auto particle = stack.NewParticle();
     particle.SetPID(Code::Proton);
@@ -235,7 +238,7 @@ int main() {
     particle.SetEnergy(E0);
     particle.SetMomentum(plab);
     particle.SetTime(0_ns);
-    Point p(rootCS, 0_m, 0_m, 10_km);
+    Point p(rootCS, 0_m, 0_m, 0_m);
     particle.SetPosition(p);
   }