IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 9a85572f authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by Ralf Ulrich
Browse files

recovered testUrQMD from !205

parent afac3fd7
No related branches found
No related tags found
1 merge request!335Resolve "UrQMD cross-section regression"
......@@ -77,21 +77,19 @@ TEST_CASE("UrQMD") {
auto const& cs = *csPtr;
{ [[maybe_unused]] auto const& env_dummy = env; }
Code validProjectileCodes[] = {Code::PiPlus, Code::PiMinus, Code::Proton,
Code::Neutron, Code::KPlus, Code::KMinus,
Code::K0, Code::K0Bar, Code::K0Long};
Code validProjectileCodes[] = {Code::PiPlus, Code::PiMinus, Code::Proton,
Code::AntiProton, Code::AntiNeutron, Code::Neutron,
Code::KPlus, Code::KMinus, Code::K0,
Code::K0Bar, Code::K0Long};
for (auto code : validProjectileCodes) {
auto [stack, view] = setup::testing::setup_stack(
code, 0, 0, 100_GeV, (setup::Environment::BaseNodeType* const)nodePtr, cs);
CHECK(stack->getEntries() == 1);
CHECK(view->getEntries() == 0);
auto [stack, view] = setup::testing::setup_stack(code, 0, 0, 100_GeV, nodePtr, cs);
REQUIRE(stack->getEntries() == 1);
REQUIRE(view->getEntries() == 0);
// simple check whether the cross-section is non-vanishing
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Proton) / 1_mb > 0);
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Nitrogen) / 1_mb > 0);
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Oxygen) / 1_mb > 0);
CHECK(urqmd.getCrossSection(view->getProjectile(), Code::Argon) / 1_mb > 0);
// only nuclei with available tabluated data so far
REQUIRE(urqmd.getInteractionLength(stack->getNextParticle()) > 1_g / square(1_cm));
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment