From a152a56add46a1e64cd77b3d081900a7b1dd506b Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@kit.edu>
Date: Wed, 17 Apr 2019 17:26:41 -0300
Subject: [PATCH] set At also for meson projectile to avoid NaN
---
Processes/UrQMD/UrQMD.cc | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/Processes/UrQMD/UrQMD.cc b/Processes/UrQMD/UrQMD.cc
index 4caedb88c..adee6aafc 100644
--- a/Processes/UrQMD/UrQMD.cc
+++ b/Processes/UrQMD/UrQMD.cc
@@ -92,9 +92,8 @@ corsika::process::EProcessReturn UrQMD::DoInteraction(SetupProjectile& projectil
cascinit_(sys_.Zp, sys_.Ap, id);
} else {
inputs_.prspflg = 1;
- int const Ap =
- 1; // what value to use here for non-baryons??? see CONEX UrQMD interface
- rsys_.bdist = nucrad_(Ap) + nucrad_(Atarget) + 2 * options_.CTParam[30 - 1];
+ sys_.Ap = 1; // even for non-baryons this has to be set, see vanilla UrQMD.f
+ rsys_.bdist = nucrad_(sys_.Ap) + nucrad_(Atarget) + 2 * options_.CTParam[30 - 1];
auto const [ityp, iso3] = ConvertToUrQMD(projectileCode);
// todo: conversion of K_long/short into strong eigenstates;
--
GitLab