IAP GITLAB

Skip to content
Snippets Groups Projects

Resolve "UrQMD cross-section regression"

Merged Maximilian Reininghaus requested to merge 403-urqmd-cross-section-regression into master
1 file
+ 9
11
Compare changes
  • Side-by-side
  • Inline
@@ -77,21 +77,19 @@ TEST_CASE("UrQMD") {
@@ -77,21 +77,19 @@ TEST_CASE("UrQMD") {
auto const& cs = *csPtr;
auto const& cs = *csPtr;
{ [[maybe_unused]] auto const& env_dummy = env; }
{ [[maybe_unused]] auto const& env_dummy = env; }
Code validProjectileCodes[] = {Code::PiPlus, Code::PiMinus, Code::Proton,
Code validProjectileCodes[] = {Code::PiPlus, Code::PiMinus, Code::Proton,
Code::Neutron, Code::KPlus, Code::KMinus,
Code::AntiProton, Code::AntiNeutron, Code::Neutron,
Code::K0, Code::K0Bar, Code::K0Long};
Code::KPlus, Code::KMinus, Code::K0,
 
Code::K0Bar, Code::K0Long};
for (auto code : validProjectileCodes) {
for (auto code : validProjectileCodes) {
auto [stack, view] = setup::testing::setup_stack(
auto [stack, view] = setup::testing::setup_stack(code, 0, 0, 100_GeV, nodePtr, cs);
code, 0, 0, 100_GeV, (setup::Environment::BaseNodeType* const)nodePtr, cs);
REQUIRE(stack->getEntries() == 1);
CHECK(stack->getEntries() == 1);
REQUIRE(view->getEntries() == 0);
CHECK(view->getEntries() == 0);
// simple check whether the cross-section is non-vanishing
// simple check whether the cross-section is non-vanishing
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Proton) / 1_mb > 0);
// only nuclei with available tabluated data so far
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Nitrogen) / 1_mb > 0);
REQUIRE(urqmd.getInteractionLength(stack->getNextParticle()) > 1_g / square(1_cm));
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Oxygen) / 1_mb > 0);
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Argon) / 1_mb > 0);
}
}
}
}
Loading