From 479e610ae3eae6937a87805fa7f45206dfe0a9ec Mon Sep 17 00:00:00 2001
From: ralfulrich <ralf.ulrich@kit.edu>
Date: Fri, 15 Oct 2021 13:35:59 +0200
Subject: [PATCH] fixed unit test values
---
tests/modules/testSibyll.cpp | 22 ++++++++++++----------
1 file changed, 12 insertions(+), 10 deletions(-)
diff --git a/tests/modules/testSibyll.cpp b/tests/modules/testSibyll.cpp
index c3b10af49..73e812b5c 100644
--- a/tests/modules/testSibyll.cpp
+++ b/tests/modules/testSibyll.cpp
@@ -263,23 +263,25 @@ TEST_CASE("SibyllInterface", "modules") {
SECTION("NuclearInteractionInterface") {
- HEPMomentumType const P0 = 2500_GeV; // per nucleon
- MomentumVector plab = MomentumVector(cs, {P0, 0_eV, 0_eV}); // per nucleon
+ HEPMomentumType const P0 = 50_TeV;
+ MomentumVector plab = MomentumVector(cs, {P0, 0_eV, 0_eV});
corsika::sibyll::InteractionModel hmodel;
NuclearInteractionModel model(hmodel, *env);
- HEPEnergyType const ElabNuc =
- sqrt(static_pow<2>(P0 * 8) + static_pow<2>(get_nucleus_mass(8, 4))) / 8;
- HEPEnergyType const sqrtSnn = sqrt((ElabNuc + constants::nucleonMass + P0) *
- (ElabNuc + constants::nucleonMass - P0));
- model.doInteraction(view, getCOMboost(ElabNuc, plab, cs), Code::Iron, Code::Oxygen,
+ Code const pid = Code::Oxygen;
+ size_t const A = get_nucleus_A(pid);
+ HEPEnergyType const Elab = sqrt(static_pow<2>(P0) + static_pow<2>(get_mass(pid)));
+ HEPEnergyType const sqrtSnn = sqrt((Elab / A + constants::nucleonMass + P0 / A) *
+ (Elab / A + constants::nucleonMass - P0 / A));
+ model.doInteraction(view, getCOMboost(Elab, plab, cs), Code::Oxygen, Code::Oxygen,
sqrtSnn);
- CrossSectionType const cx = model.getCrossSection(Code::Iron, Code::Oxygen, sqrtSnn);
+ CrossSectionType const cx =
+ model.getCrossSection(Code::Oxygen, Code::Oxygen, sqrtSnn);
// Felix, are those changes OK? Below are the checks before refactory-2020
// CHECK(length / 1_g * 1_cm * 1_cm == Approx(44.2).margin(.1));
// CHECK(view.getSize() == 11);
- CHECK(cx / 1_mb == Approx(1900).margin(100)); // this is not physics validation
+ CHECK(cx / 1_mb == Approx(1300).margin(300)); // this is not physics validation
// CHECK(view.getSize() == 20); // also sibyll not stable wrt. to compiler changes
- CHECK(view.getSize() == Approx(300).margin(90)); // this is not physics validation
+ CHECK(view.getSize() == Approx(100).margin(90)); // this is not physics validation
}
}
--
GitLab