IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 37281802 authored by ralfulrich's avatar ralfulrich Committed by Felix Riehn
Browse files

added Felix suggestions

parent 186aefa5
No related branches found
No related tags found
No related merge requests found
...@@ -157,9 +157,11 @@ TEST_CASE("pythia process") { ...@@ -157,9 +157,11 @@ TEST_CASE("pythia process") {
SECTION("pythia interaction") { SECTION("pythia interaction") {
// feenableexcept(FE_INVALID); \todo how does this work nowadays //! feenableexcept(FE_INVALID); \todo how does this work nowadays
// this will be a p-p collision at sqrts=3.5TeV -> no problem for pythia
auto [stackPtr, secViewPtr] = setup::testing::setup_stack( auto [stackPtr, secViewPtr] = setup::testing::setup_stack(
Code::PiPlus, 0, 0, 100_GeV, (setup::Environment::BaseNodeType* const)nodePtr, Code::Proton, 0, 0, 7_TeV, (setup::Environment::BaseNodeType* const)nodePtr,
*csPtr); *csPtr);
auto& view = *secViewPtr; auto& view = *secViewPtr;
auto particle = stackPtr->first(); auto particle = stackPtr->first();
...@@ -167,6 +169,7 @@ TEST_CASE("pythia process") { ...@@ -167,6 +169,7 @@ TEST_CASE("pythia process") {
corsika::pythia8::Interaction model; corsika::pythia8::Interaction model;
model.doInteraction(view); model.doInteraction(view);
[[maybe_unused]] const GrammageType length = model.getInteractionLength(particle); [[maybe_unused]] const GrammageType length = model.getInteractionLength(particle);
CHECK(length / 1_kg * square(1_m) == Approx(82.2524)); CHECK(length / 1_kg * square(1_m) == Approx(43.04).margin(5e-1));
CHECK(view.getSize() == 38);
} }
} }
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