IAP GITLAB

Skip to content
Snippets Groups Projects
Commit c4b8af2e authored by ralfulrich's avatar ralfulrich
Browse files

fixed wrong logic in BetheBlochPDG

parent f3270a8c
No related branches found
No related tags found
No related merge requests found
...@@ -141,14 +141,16 @@ namespace corsika { ...@@ -141,14 +141,16 @@ namespace corsika {
inline ProcessReturn BetheBlochPDG::doContinuous(setup::Stack::particle_type& p, inline ProcessReturn BetheBlochPDG::doContinuous(setup::Stack::particle_type& p,
setup::Trajectory const& t, setup::Trajectory const& t,
bool const limitStep) { bool const) {
// if this step was limiting the CORSIKA stepping, the particle is lost // if this step was limiting the CORSIKA stepping, the particle is lost
/* see Issue https://gitlab.ikp.kit.edu/AirShowerPhysics/corsika/-/issues/389
if (limitStep) { if (limitStep) {
fillProfile(t, p.getEnergy()); fillProfile(t, p.getEnergy());
p.setEnergy(p.getMass()); p.setEnergy(p.getMass());
return ProcessReturn::ParticleAbsorbed; return ProcessReturn::ParticleAbsorbed;
} }
*/
if (p.getChargeNumber() == 0) return ProcessReturn::Ok; if (p.getChargeNumber() == 0) return ProcessReturn::Ok;
......
...@@ -101,8 +101,10 @@ int main(int argc, char** argv) { ...@@ -101,8 +101,10 @@ int main(int argc, char** argv) {
CORSIKA_LOG_INFO("vertical_EAS"); CORSIKA_LOG_INFO("vertical_EAS");
if (argc < 4) { if (argc < 4) {
std::cerr << "usage: vertical_EAS <A> <Z> <energy/GeV> [seed]" << std::endl; std::cerr << "usage: vertical_EAS <A> <Z> <energy/GeV> [seed] \n"
std::cerr << " if no seed is given, a random seed is chosen" << std::endl; " if A=0, Z is interpreted as PDG code \n"
" if no seed is given, a random seed is chosen \n"
<< std::endl;
return 1; return 1;
} }
feenableexcept(FE_INVALID); feenableexcept(FE_INVALID);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment