IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 138fefe9 authored by André Schmidt's avatar André Schmidt Committed by ralfulrich
Browse files

magMaxLength is not a decay

parent 98d72bf5
No related branches found
No related tags found
No related merge requests found
...@@ -230,8 +230,7 @@ namespace corsika::cascade { ...@@ -230,8 +230,7 @@ namespace corsika::cascade {
auto magneticfield = corsika::geometry::Vector(fEnvironment.GetCoordinateSystem(), auto magneticfield = corsika::geometry::Vector(fEnvironment.GetCoordinateSystem(),
0_uT, 50_uT, 0_uT); 0_uT, 50_uT, 0_uT);
geometry::Vector<SpeedType::dimension_type> const velocityVerticalMag = geometry::Vector<SpeedType::dimension_type> const velocityVerticalMag =
velocity - magneticfield * velocity.dot(magneticfield) / velocity - velocity.parallelProjectionOnto(magneticfield);
(magneticfield.GetSquaredNorm());
LengthType const gyroradius = LengthType const gyroradius =
vParticle.GetEnergy() * velocityVerticalMag.GetNorm() * 1_V / vParticle.GetEnergy() * velocityVerticalMag.GetNorm() * 1_V /
(corsika::units::constants::cSquared * abs(chargeNumber) * (corsika::units::constants::cSquared * abs(chargeNumber) *
...@@ -308,7 +307,7 @@ namespace corsika::cascade { ...@@ -308,7 +307,7 @@ namespace corsika::cascade {
TStackView secondaries(vParticle); TStackView secondaries(vParticle);
if (min_distance != distance_max) { if (min_distance != distance_max && min_distance != magMaxLength) {
/* /*
Create SecondaryView object on Stack. The data container Create SecondaryView object on Stack. The data container
remains untouched and identical, and 'projectil' is identical remains untouched and identical, and 'projectil' is identical
......
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