IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 7552cc35 authored by Maximilian Sackel's avatar Maximilian Sackel
Browse files

remove scaling factor from em component maximal step length

parent 370d68ad
No related branches found
No related tags found
1 merge request!245Include proposal process rebase
Pipeline #2201 canceled
...@@ -78,7 +78,7 @@ namespace corsika::process::proposal { ...@@ -78,7 +78,7 @@ namespace corsika::process::proposal {
auto energy_loss = TotalEnergyLoss(vP, dX); auto energy_loss = TotalEnergyLoss(vP, dX);
if (vP.GetChargeNumber() != 0) Scatter(vP, energy_loss, dX); if (vP.GetChargeNumber() != 0) Scatter(vP, energy_loss, dX);
vP.SetEnergy(vP.GetEnergy() - energy_loss); 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()); vP.SetMomentum(vP.GetMomentum() * vP.GetEnergy() / vP.GetMomentum().GetNorm());
return process::EProcessReturn::eOk; return process::EProcessReturn::eOk;
} }
...@@ -94,8 +94,7 @@ namespace corsika::process::proposal { ...@@ -94,8 +94,7 @@ namespace corsika::process::proposal {
auto grammage = get<DISPLACEMENT>(c->second)->SolveTrackIntegral( auto grammage = get<DISPLACEMENT>(c->second)->SolveTrackIntegral(
vP.GetEnergy() / 1_MeV, energy_lim / 1_MeV) * vP.GetEnergy() / 1_MeV, energy_lim / 1_MeV) *
1_g / square(1_cm); 1_g / square(1_cm);
return vP.GetNode()->GetModelProperties().ArclengthFromGrammage(vT, grammage) * return vP.GetNode()->GetModelProperties().ArclengthFromGrammage(vT, grammage);
1.0001;
} }
} // namespace corsika::process::proposal } // namespace corsika::process::proposal
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