IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 56d103b6 authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan:
Browse files

limit trajectory after end-point determination for later purposes

parent 0093a324
No related branches found
No related tags found
No related merge requests found
...@@ -113,6 +113,8 @@ namespace corsika::cascade { ...@@ -113,6 +113,8 @@ namespace corsika::cascade {
// std::visit(corsika::setup::ParticleUpdate<Particle>{particle}, step); // std::visit(corsika::setup::ParticleUpdate<Particle>{particle}, step);
particle.SetPosition(step.PositionFromArclength(min_distance)); particle.SetPosition(step.PositionFromArclength(min_distance));
// .... also update time, momentum, direction, ... // .... also update time, momentum, direction, ...
step.LimitEndTo(min_distance);
// apply all continuous processes on particle + track // apply all continuous processes on particle + track
corsika::process::EProcessReturn status = corsika::process::EProcessReturn status =
...@@ -129,7 +131,7 @@ namespace corsika::cascade { ...@@ -129,7 +131,7 @@ namespace corsika::cascade {
<< ((min_distance < distance_max) ? "yes" : "no") << std::endl; << ((min_distance < distance_max) ? "yes" : "no") << std::endl;
if (min_distance < distance_max) { // interaction to happen within geometric limit if (min_distance < distance_max) { // interaction to happen within geometric limit
// check weather decay or interaction limits this step // check whether decay or interaction limits this step
if (min_distance == distance_interact) { if (min_distance == distance_interact) {
std::cout << "collide" << std::endl; std::cout << "collide" << std::endl;
......
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