From c3fdd3fb605eadbebc8edf78ad5979a6a339e222 Mon Sep 17 00:00:00 2001
From: Remy Prechelt <prechelt@hawaii.edu>
Date: Wed, 3 Feb 2021 20:10:42 -1000
Subject: [PATCH] Add outputs to vertical_EAS example.
---
examples/vertical_EAS.cpp | 15 ++++++++++++---
1 file changed, 12 insertions(+), 3 deletions(-)
diff --git a/examples/vertical_EAS.cpp b/examples/vertical_EAS.cpp
index f012440e8..590a8e819 100644
--- a/examples/vertical_EAS.cpp
+++ b/examples/vertical_EAS.cpp
@@ -25,6 +25,8 @@
#include <corsika/framework/core/Cascade.hpp>
#include <corsika/framework/geometry/PhysicalGeometry.hpp>
+#include <corsika/output/OutputManager.hpp>
+
#include <corsika/media/Environment.hpp>
#include <corsika/media/FlatExponential.hpp>
#include <corsika/media/HomogeneousMedium.hpp>
@@ -188,6 +190,9 @@ int main(int argc, char** argv) {
ShowerAxis const showerAxis{injectionPos, (showerCore - injectionPos) * 1.5, env};
+ // create the output manager that we then register outputs with
+ OutputManager output("vertical_EAS_outputs");
+
// setup processes, decays and interactions
corsika::sibyll::Interaction sibyll;
@@ -232,8 +237,10 @@ int main(int argc, char** argv) {
LongitudinalProfile longprof{showerAxis};
Plane const obsPlane(showerCore, DirectionVector(rootCS, {0., 0., 1.}));
- ObservationPlane observationLevel(obsPlane, DirectionVector(rootCS, {1., 0., 0.}),
- "particles.dat");
+ ObservationPlane observationLevel(obsPlane, DirectionVector(rootCS, {1., 0., 0.}));
+
+ // register the observation plane with the output
+ output.add("obsplane", observationLevel);
corsika::urqmd::UrQMD urqmd;
InteractionCounter urqmdCounted{urqmd};
@@ -260,12 +267,14 @@ int main(int argc, char** argv) {
// define air shower object, run simulation
setup::Tracking tracking;
- Cascade EAS(env, tracking, sequence, stack);
+ Cascade EAS(env, tracking, sequence, output, stack);
// to fix the point of first interaction, uncomment the following two lines:
// EAS.forceInteraction();
+ output.startOfRun();
EAS.run();
+ output.endOfRun();
cut.showResults();
em_continuous.showResults();
--
GitLab