diff --git a/Processes/UrQMD/UrQMD.cc b/Processes/UrQMD/UrQMD.cc index 0c9f632657c41b7f60722ff0f4ae0f7e1d499f60..196fcdf10c746148b7f7d520e4f18069ab1f2caf 100644 --- a/Processes/UrQMD/UrQMD.cc +++ b/Processes/UrQMD/UrQMD.cc @@ -293,20 +293,9 @@ 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;