From 4fe1d9a018df60d78d4dedd439e102630a589641 Mon Sep 17 00:00:00 2001
From: Jean-Marco Alameddine <jean-marco.alameddine@tu-dortmund.de>
Date: Fri, 9 Dec 2022 17:16:33 +0100
Subject: [PATCH] Convert energy_production_threshold value from kinetic energy
to total energy before passing them to PROPOSAL
---
corsika/modules/proposal/ProposalProcessBase.hpp | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/corsika/modules/proposal/ProposalProcessBase.hpp b/corsika/modules/proposal/ProposalProcessBase.hpp
index 6780c3a6b..807d35819 100644
--- a/corsika/modules/proposal/ProposalProcessBase.hpp
+++ b/corsika/modules/proposal/ProposalProcessBase.hpp
@@ -46,9 +46,10 @@ namespace corsika::proposal {
corsika::units::si::HEPEnergyType
emCut) { //!< Stochastic losses smaller than the given cut
//!< will be handeled continuously.
- auto p_cut =
- std::make_shared<const PROPOSAL::EnergyCutSettings>(emCut / 1_MeV, 1, false);
- return PROPOSAL::GetStdCrossSections(T(), m, p_cut, true);
+ auto particle_def = T();
+ auto p_cut = std::make_shared<const PROPOSAL::EnergyCutSettings>(
+ particle_def.mass + emCut / 1_MeV, 1, false);
+ return PROPOSAL::GetStdCrossSections(particle_def, m, p_cut, true);
};
//!
--
GitLab