From 6beecfbc4262426e7cfecd8b12fddca98c7cb1e5 Mon Sep 17 00:00:00 2001
From: Felix Riehn <felix@matilda>
Date: Thu, 19 Jan 2023 12:45:41 +0100
Subject: [PATCH] added sophia to water example
---
examples/water.cpp | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/examples/water.cpp b/examples/water.cpp
index 73accbb53..d4b04737a 100644
--- a/examples/water.cpp
+++ b/examples/water.cpp
@@ -32,6 +32,7 @@
#include <corsika/modules/Pythia8.hpp>
#include <corsika/modules/Random.hpp>
#include <corsika/modules/Sibyll.hpp>
+#include <corsika/modules/Sophia.hpp>
#include <corsika/modules/UrQMD.hpp>
#include <corsika/modules/tracking/TrackingStraight.hpp>
@@ -58,6 +59,7 @@ void registerRandomStreams(int seed) {
RNGManager<>::getInstance().registerRandomStream("cascade");
RNGManager<>::getInstance().registerRandomStream("qgsjet");
RNGManager<>::getInstance().registerRandomStream("sibyll");
+ RNGManager<>::getInstance().registerRandomStream("sophia");
RNGManager<>::getInstance().registerRandomStream("epos");
RNGManager<>::getInstance().registerRandomStream("pythia");
RNGManager<>::getInstance().registerRandomStream("urqmd");
@@ -241,9 +243,11 @@ int main(int argc, char** argv) {
}};
auto decaySequence = make_sequence(decayPythia, decaySibyll);
+ corsika::sophia::InteractionModel sophia;
+
// EM process
- corsika::proposal::Interaction emCascade(env, sibyll.getHadronInteractionModel(),
- heHadronModelThreshold);
+ corsika::proposal::Interaction emCascade(
+ env, sophia, sibyll.getHadronInteractionModel(), heHadronModelThreshold);
corsika::proposal::ContinuousProcess<SubWriter<decltype(dEdX)>> emContinuous(env, dEdX);
// total physics list
--
GitLab