From fe278f6709f7fad993b41f82ba7b46f8dd1f8226 Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@tu-dortmund.de>
Date: Fri, 12 Jun 2020 17:19:59 +0200
Subject: [PATCH] use ShowerAxis in cascade_example
---
Documentation/Examples/cascade_example.cc | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc
index f5c09343b..7a2db1390 100644
--- a/Documentation/Examples/cascade_example.cc
+++ b/Documentation/Examples/cascade_example.cc
@@ -21,6 +21,7 @@
#include <corsika/environment/Environment.h>
#include <corsika/environment/HomogeneousMedium.h>
#include <corsika/environment/NuclearComposition.h>
+#include <corsika/environment/ShowerAxis.h>
#include <corsika/geometry/Sphere.h>
@@ -107,6 +108,8 @@ int main() {
rootCS, 0_m, 0_m,
height_atmosphere); // this is the CORSIKA 7 start of atmosphere/universe
+ ShowerAxis const showerAxis{injectionPos, Vector{rootCS, 0_m, 0_m, 120_km}, env};
+
{
auto elab2plab = [](HEPEnergyType Elab, HEPMassType m) {
return sqrt((Elab - m) * (Elab + m));
@@ -141,8 +144,7 @@ int main() {
process::particle_cut::ParticleCut cut(80_GeV);
process::track_writer::TrackWriter trackWriter("tracks.dat");
- process::energy_loss::EnergyLoss eLoss(
- injectionPos, geometry::Vector<dimensionless_d>(rootCS, {0, 0, -1}));
+ process::energy_loss::EnergyLoss eLoss{showerAxis};
// assemble all processes into an ordered process list
auto sequence = stackInspect << sibyll << sibyllNuc << decay << eLoss << cut
--
GitLab