From 7552cc35b29ce939789743948e2c7fb0d1105a03 Mon Sep 17 00:00:00 2001
From: Maximilian Sackel <maximilian.sackel@tu-dortmund.de>
Date: Wed, 30 Sep 2020 12:48:11 +0000
Subject: [PATCH] remove scaling factor from em component maximal step length
---
Processes/Proposal/ContinuousProcess.cc | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/Processes/Proposal/ContinuousProcess.cc b/Processes/Proposal/ContinuousProcess.cc
index 582a59882..49427f03e 100644
--- a/Processes/Proposal/ContinuousProcess.cc
+++ b/Processes/Proposal/ContinuousProcess.cc
@@ -78,7 +78,7 @@ namespace corsika::process::proposal {
auto energy_loss = TotalEnergyLoss(vP, dX);
if (vP.GetChargeNumber() != 0) Scatter(vP, energy_loss, dX);
vP.SetEnergy(vP.GetEnergy() - energy_loss);
- if (vP.GetEnergy() < cut.GetECut()) return process::EProcessReturn::eParticleAbsorbed;
+ if (vP.GetEnergy() <= cut.GetECut()) return process::EProcessReturn::eParticleAbsorbed;
vP.SetMomentum(vP.GetMomentum() * vP.GetEnergy() / vP.GetMomentum().GetNorm());
return process::EProcessReturn::eOk;
}
@@ -94,8 +94,7 @@ namespace corsika::process::proposal {
auto grammage = get<DISPLACEMENT>(c->second)->SolveTrackIntegral(
vP.GetEnergy() / 1_MeV, energy_lim / 1_MeV) *
1_g / square(1_cm);
- return vP.GetNode()->GetModelProperties().ArclengthFromGrammage(vT, grammage) *
- 1.0001;
+ return vP.GetNode()->GetModelProperties().ArclengthFromGrammage(vT, grammage);
}
} // namespace corsika::process::proposal
--
GitLab