IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 39e8aed6 authored by Nikos Karastathis's avatar Nikos Karastathis :ocean: Committed by Maximilian Reininghaus
Browse files

fix hybrid_MC.cpp

parent 3182d3a0
No related branches found
No related tags found
1 merge request!426Resolve "Cascade: Problems with Multiple Scattering (in combination with tracking)"
...@@ -23,6 +23,7 @@ ...@@ -23,6 +23,7 @@
#include <corsika/framework/core/Logging.hpp> #include <corsika/framework/core/Logging.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp> #include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/core/Cascade.hpp> #include <corsika/framework/core/Cascade.hpp>
#include <corsika/framework/core/Step.hpp>
#include <corsika/framework/core/EnergyMomentumOperations.hpp> #include <corsika/framework/core/EnergyMomentumOperations.hpp>
#include <corsika/framework/random/RNGManager.hpp> #include <corsika/framework/random/RNGManager.hpp>
...@@ -117,14 +118,14 @@ public: ...@@ -117,14 +118,14 @@ public:
* @param particle * @param particle
* @return ProcessReturn * @return ProcessReturn
*/ */
template <typename TParticle, typename TTrack> template <typename TParticle>
ProcessReturn doContinuous(TParticle const& particle, TTrack const&, bool const) { ProcessReturn doContinuous(Step<TParticle> const& step, bool const) {
auto const delta = particle.getPosition() - plane_.getCenter(); auto const delta = step.getParticlePre().getPosition() - plane_.getCenter();
auto const n = plane_.getNormal(); auto const n = plane_.getNormal();
auto const proj = n.dot(delta); auto const proj = n.dot(delta);
if (proj < -1_m) { if (proj < -1_m) {
CORSIKA_LOG_INFO("particle {} failes: proj={}, delta={}, p={}", particle.asString(), CORSIKA_LOG_INFO("particle {} failes: proj={}, delta={}, p={}", step.getParticlePre().asString(),
proj, delta, particle.getPosition()); proj, delta, step.getPositionPost());
throw std::runtime_error("particle below obs level"); throw std::runtime_error("particle below obs level");
} }
return ProcessReturn::Ok; return ProcessReturn::Ok;
......
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