From 48b6c2d92af3b3cd85fd06c440d13aab5eb9e32e Mon Sep 17 00:00:00 2001
From: Andre Schmidt <schmidt-a@iklx288.ikp.kit.edu>
Date: Fri, 16 Oct 2020 18:23:34 +0200
Subject: [PATCH] sth
---
Processes/ObservationPlane/ObservationPlane.cc | 11 ++++++++---
1 file changed, 8 insertions(+), 3 deletions(-)
diff --git a/Processes/ObservationPlane/ObservationPlane.cc b/Processes/ObservationPlane/ObservationPlane.cc
index f1fb72638..f624c8e3a 100644
--- a/Processes/ObservationPlane/ObservationPlane.cc
+++ b/Processes/ObservationPlane/ObservationPlane.cc
@@ -70,13 +70,18 @@ LengthType ObservationPlane::MaxStepLength(setup::Stack::ParticleType const& vPa
auto const* currentLogicalVolumeNode = vParticle.GetNode();
auto magneticfield = currentLogicalVolumeNode->GetModelProperties().GetMagneticField(vParticle.GetPosition());
auto direction = trajectory.GetV0().normalized();
-
- std::cout << " " << plane_.GetNormal().dot(trajectory.GetV0().cross(magneticfield)) << std::endl;
- if (chargeNumber != 0 && abs(plane_.GetNormal().dot(trajectory.GetV0().cross(magneticfield))) * 1_s / 1_m / 1_T > 1e-5) {
+ if (chargeNumber != 0 && abs(plane_.GetNormal().dot(trajectory.GetV0().cross(magneticfield))) * 1_s / 1_m / 1_T > 1e-6) {
auto const* currentLogicalVolumeNode = vParticle.GetNode();
auto magneticfield = currentLogicalVolumeNode->GetModelProperties().GetMagneticField(vParticle.GetPosition());
auto k = chargeNumber * corsika::units::constants::c * 1_eV / (vParticle.GetMomentum().norm() * 1_V);
+
+ if (direction.dot(plane_.GetNormal()) * direction.dot(plane_.GetNormal()) -
+ (plane_.GetNormal().dot(trajectory.GetR0() - plane_.GetCenter()) *
+ plane_.GetNormal().dot(direction.cross(magneticfield)) * 2 * k) < 0)
+ {
+ return std::numeric_limits<double>::infinity() * 1_m;
+ }
LengthType MaxStepLength1 =
( sqrt(direction.dot(plane_.GetNormal()) * direction.dot(plane_.GetNormal()) -
--
GitLab