From b18f670445cb9bf9a9161aba40293e63fb9c04cf Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@tu-dortmund.de>
Date: Wed, 22 Apr 2020 17:25:53 +0200
Subject: [PATCH] allow K0Short interaction in UrQMD
---
Processes/UrQMD/UrQMD.cc | 9 +++------
1 file changed, 3 insertions(+), 6 deletions(-)
diff --git a/Processes/UrQMD/UrQMD.cc b/Processes/UrQMD/UrQMD.cc
index 09d171e26..a595d524e 100644
--- a/Processes/UrQMD/UrQMD.cc
+++ b/Processes/UrQMD/UrQMD.cc
@@ -292,10 +292,9 @@ corsika::process::EProcessReturn UrQMD::DoInteraction(SetupProjectile& vProjecti
rsys_.bdist = nucrad_(targetA) + nucrad_(1) + 2 * options_.CTParam[30 - 1];
rsys_.ebeam = (projectileEnergyLab - vProjectile.GetMass()) * (1 / 1_GeV);
- if (projectileCode == particles::Code::K0Long) {
+ if (projectileCode == particles::Code::K0Long ||
+ projectileCode == particles::Code::K0Short) {
projectileCode = fBooleanDist(fRNG) ? particles::Code::K0 : particles::Code::K0Bar;
- } else if (projectileCode == particles::Code::K0Short) {
- throw std::runtime_error("K0Short should not interact");
}
auto const [ityp, iso3] = ConvertToUrQMD(projectileCode);
@@ -438,9 +437,7 @@ std::pair<int, int> corsika::process::UrQMD::ConvertToUrQMD(
void UrQMD::readXSFile(std::string const& filename) {
std::ifstream file(filename, std::ios::in);
- if (!file.is_open()) {
- throw std::runtime_error(filename + " could not be opened.");
- }
+ if (!file.is_open()) { throw std::runtime_error(filename + " could not be opened."); }
std::string line;
--
GitLab