diff --git a/examples/cascade_examples/mars.cpp b/examples/cascade_examples/mars.cpp index 26029b80422cc9a28086ee7c265389245e83d004..d41b6f562301d95dd442796751ba74e5b96ea160 100644 --- a/examples/cascade_examples/mars.cpp +++ b/examples/cascade_examples/mars.cpp @@ -312,6 +312,15 @@ int main(int argc, char** argv) { HEPEnergyType const hadcut = 1_GeV; ParticleCut<SubWriter<decltype(dEdX)>> cut(emcut, emcut, hadcut, hadcut, true, dEdX); + // tell proposal that we are interested in all energy losses above the particle cut + set_energy_production_threshold(Code::Electron, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::Positron, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::Photon, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::MuMinus, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::MuPlus, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::TauMinus, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::TauPlus, std::min({emcut, hadcut})); + /* === START: SETUP PROCESS LIST === */ corsika::sibyll::Interaction sibyll{env}; InteractionCounter sibyllCounted(sibyll); diff --git a/examples/cascade_examples/water.cpp b/examples/cascade_examples/water.cpp index 9943b634e48466faf04a227953671aec49df6c07..add92b35a74ad495ae4e2d7e77920744832e8ebc 100644 --- a/examples/cascade_examples/water.cpp +++ b/examples/cascade_examples/water.cpp @@ -212,6 +212,15 @@ int main(int argc, char** argv) { HEPEnergyType const hadCut = eCut; ParticleCut<SubWriter<decltype(dEdX)>> cut(emCut, emCut, hadCut, hadCut, true, dEdX); + // tell proposal that we are interested in all energy losses above the particle cut + set_energy_production_threshold(Code::Electron, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::Positron, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::Photon, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::MuMinus, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::MuPlus, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::TauMinus, std::min({emcut, hadcut})); + set_energy_production_threshold(Code::TauPlus, std::min({emcut, hadcut})); + // hadronic interactions HEPEnergyType heHadronModelThreshold = std::pow(10, 1.9) * 1_GeV; corsika::sibyll::Interaction sibyll(env);