diff --git a/corsika/detail/modules/energy_loss/BetheBlochPDG.inl b/corsika/detail/modules/energy_loss/BetheBlochPDG.inl index 1ac6d7098fb4cf41e84d4b6243f74e460e6daaa9..1511d5029789cafddf4e564eac602bc13cf149b8 100644 --- a/corsika/detail/modules/energy_loss/BetheBlochPDG.inl +++ b/corsika/detail/modules/energy_loss/BetheBlochPDG.inl @@ -184,9 +184,9 @@ namespace corsika { energy * 0.9, // either 10% relative loss max., or get_energy_threshold(vParticle.getPID()) // energy thresholds globally defined for // individual particles - * - 0.99999 // need to go 1% below global e-cut to assure removal in ParticleCut. The - // 1% does not matter since at cut-time the entire energy is removed. + * 0.99999 // need to go 1% below global e-cut to assure removal in + // ParticleCut. The 1% does not matter since at cut-time the entire + // energy is removed. ); auto const maxGrammage = (vParticle.getEnergy() - energy_lim) / dEdX; diff --git a/corsika/detail/modules/tracking/TrackingLeapFrogCurved.inl b/corsika/detail/modules/tracking/TrackingLeapFrogCurved.inl index fc66093a9780860ca1b15fa2257a96506708ac7b..cad673519f8de048c7df4bdaeb3a62e55be126b1 100644 --- a/corsika/detail/modules/tracking/TrackingLeapFrogCurved.inl +++ b/corsika/detail/modules/tracking/TrackingLeapFrogCurved.inl @@ -278,8 +278,8 @@ namespace corsika { std::vector<double> const deltaLs = solve_quadratic_real(denom, p, q); - CORSIKA_LOG_TRACE("deltaLs=[{}]", fmt::join(deltaLs, ", ")); - + CORSIKA_LOG_TRACE("deltaLs=[{}]", fmt::join(deltaLs, ", ")); + if (deltaLs.size() == 0) { return Intersections(std::numeric_limits<double>::infinity() * 1_s); }