From cf90d1f54e14539259d038d700e12fda25c4412f Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@tu-dortmund.de>
Date: Sun, 26 Apr 2020 02:26:03 +0200
Subject: [PATCH] UrQMD impact parameter range from C7 for O,N,Ar
---
Processes/UrQMD/UrQMD.cc | 13 ++++++++++++-
1 file changed, 12 insertions(+), 1 deletion(-)
diff --git a/Processes/UrQMD/UrQMD.cc b/Processes/UrQMD/UrQMD.cc
index a595d524e..927054df1 100644
--- a/Processes/UrQMD/UrQMD.cc
+++ b/Processes/UrQMD/UrQMD.cc
@@ -289,9 +289,20 @@ corsika::process::EProcessReturn UrQMD::DoInteraction(SetupProjectile& vProjecti
} else {
inputs_.prspflg = 1;
sys_.Ap = 1; // even for non-baryons this has to be set, see vanilla UrQMD.f
- rsys_.bdist = nucrad_(targetA) + nucrad_(1) + 2 * options_.CTParam[30 - 1];
rsys_.ebeam = (projectileEnergyLab - vProjectile.GetMass()) * (1 / 1_GeV);
+ switch (targetCode) {
+ case particles::Code::Nitrogen:
+ case particles::Code::Oxygen:
+ rsys_.bdist = 6;
+ break;
+ case particles::Code::Argon:
+ rsys_.bdist = 7;
+ break;
+ default:
+ rsys_.bdist = nucrad_(targetA) + nucrad_(1) + 2 * options_.CTParam[30 - 1];
+ }
+
if (projectileCode == particles::Code::K0Long ||
projectileCode == particles::Code::K0Short) {
projectileCode = fBooleanDist(fRNG) ? particles::Code::K0 : particles::Code::K0Bar;
--
GitLab