From 59082738bf92abd53cf5895cdeefe31fa027742d Mon Sep 17 00:00:00 2001
From: Felix Riehn <friehn@lip.pt>
Date: Wed, 14 Feb 2024 18:57:01 +0100
Subject: [PATCH] use neutrino interaction
---
applications/c8_air_shower.cpp | 14 ++++++++++----
corsika/modules/Pythia8.hpp | 1 +
2 files changed, 11 insertions(+), 4 deletions(-)
diff --git a/applications/c8_air_shower.cpp b/applications/c8_air_shower.cpp
index ba60fd009..02840074b 100644
--- a/applications/c8_air_shower.cpp
+++ b/applications/c8_air_shower.cpp
@@ -147,7 +147,9 @@ int main(int argc, char** argv) {
->needs(opt_Z)
->check(CLI::Range(1, 58))
->group("Primary");
- app.add_option("-p,--pdg", "PDG code for primary.")
+ app.add_option("-p,--pdg",
+ "PDG code for primary (p=2212, gamma=22, e-=11, nu_e=12, mu-=13, "
+ "nu_mu=14, tau=15, nu_tau=16).")
->excludes(opt_A)
->excludes(opt_Z)
->group("Primary");
@@ -397,6 +399,9 @@ int main(int argc, char** argv) {
corsika::pythia8::Decay decayPythia;
+ // neutrino interactions
+ corsika::pythia8::NeutrinoInteraction neutrinoPrimaryPythia(E0);
+
// hadronic photon interactions in resonance region
corsika::sophia::InteractionModel sophia;
@@ -541,9 +546,9 @@ int main(int argc, char** argv) {
output.add("ZHS", zhs);
// assemble the final process sequence with radio
- auto sequence =
- make_sequence(stackInspect, hadronSequence, decayPythia, emCascade, emContinuous,
- coreas, zhs, longprof, observationLevel, thinning, cut);
+ auto sequence = make_sequence(stackInspect, neutrinoPrimaryPythia, hadronSequence,
+ decayPythia, emCascade, emContinuous, coreas, zhs,
+ longprof, observationLevel, thinning, cut);
/* === END: SETUP PROCESS LIST === */
@@ -554,6 +559,7 @@ int main(int argc, char** argv) {
Cascade EAS(env, tracking, sequence, output, stack);
// print our primary parameters all in one place
+ CORSIKA_LOG_INFO("Primary name: {}", beamCode);
if (app["--pdg"]->count() > 0) {
CORSIKA_LOG_INFO("Primary PDG ID: {}", app["--pdg"]->as<int>());
} else {
diff --git a/corsika/modules/Pythia8.hpp b/corsika/modules/Pythia8.hpp
index 9fdf773ad..ee5d18e74 100644
--- a/corsika/modules/Pythia8.hpp
+++ b/corsika/modules/Pythia8.hpp
@@ -10,3 +10,4 @@
#include <corsika/modules/pythia8/Decay.hpp>
#include <corsika/modules/pythia8/Interaction.hpp>
+#include <corsika/modules/pythia8/NeutrinoInteraction.hpp>
--
GitLab