From 5770f80c5ece4e99e9775754a96446f47fbdd694 Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@kit.edu>
Date: Fri, 21 Dec 2018 19:18:34 +0100
Subject: [PATCH] added TrackWriter to cascade_example for visualization
---
Documentation/Examples/CMakeLists.txt | 1 +
Documentation/Examples/cascade_example.cc | 9 ++++++---
2 files changed, 7 insertions(+), 3 deletions(-)
diff --git a/Documentation/Examples/CMakeLists.txt b/Documentation/Examples/CMakeLists.txt
index 5a5bd6553..cd7e52513 100644
--- a/Documentation/Examples/CMakeLists.txt
+++ b/Documentation/Examples/CMakeLists.txt
@@ -28,6 +28,7 @@ target_link_libraries (cascade_example SuperStupidStack CORSIKAunits CORSIKAlogg
ProcessSibyll
CORSIKAcascade
ProcessStackInspector
+ ProcessTrackWriter
CORSIKAprocesses
CORSIKAparticles
CORSIKAgeometry
diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc
index 232e7d002..0e6ce843b 100644
--- a/Documentation/Examples/cascade_example.cc
+++ b/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);
}
--
GitLab