From eef1529de5d9c9cc82b145fd4ad626e17878ce9e Mon Sep 17 00:00:00 2001 From: ralfulrich <ralf.ulrich@kit.edu> Date: Sun, 6 Jun 2021 22:54:12 +0200 Subject: [PATCH] sibyl tests --- corsika/detail/modules/sibyll/NuclearInteraction.inl | 9 ++++----- tests/modules/testSibyll.cpp | 2 +- 2 files changed, 5 insertions(+), 6 deletions(-) diff --git a/corsika/detail/modules/sibyll/NuclearInteraction.inl b/corsika/detail/modules/sibyll/NuclearInteraction.inl index ccdd56b2d..a75207f1c 100644 --- a/corsika/detail/modules/sibyll/NuclearInteraction.inl +++ b/corsika/detail/modules/sibyll/NuclearInteraction.inl @@ -499,18 +499,17 @@ namespace corsika::sibyll { // in pf in common fragments, neglected fragm_(kATarget, kAProj, nIntProj, impactPar, nFragments, AFragments); - // this should not occur but well :) + // this should not occur but well :) (LCOV_EXCL_START) if (nFragments > (int)getMaxNFragments()) throw std::runtime_error("Number of nuclear fragments in NUCLIB exceeded!"); + // (LCOV_EXCL_STOP) CORSIKA_LOG_DEBUG("number of fragments: " + std::to_string(nFragments)); - for (int j = 0; j < nFragments; ++j) - CORSIKA_LOG_DEBUG("fragment {}: A={} px={} py={} pz={}", j, AFragments[j], - fragments_.ppp[j][0], fragments_.ppp[j][1], fragments_.ppp[j][2]); - CORSIKA_LOG_DEBUG("adding nuclear fragments to particle stack.."); // put nuclear fragments on corsika stack for (int j = 0; j < nFragments; ++j) { + CORSIKA_LOG_DEBUG("fragment {}: A={} px={} py={} pz={}", j, AFragments[j], + fragments_.ppp[j][0], fragments_.ppp[j][1], fragments_.ppp[j][2]); Code specCode; const auto nuclA = AFragments[j]; // get Z from stability line diff --git a/tests/modules/testSibyll.cpp b/tests/modules/testSibyll.cpp index 5e23c07e2..58e65000e 100644 --- a/tests/modules/testSibyll.cpp +++ b/tests/modules/testSibyll.cpp @@ -283,7 +283,7 @@ TEST_CASE("SibyllInterface", "[processes]") { SECTION("NuclearInteractionInterface") { auto [stack, viewPtr] = - setup::testing::setup_stack(Code::Nucleus, 4, 2, 500_GeV, + setup::testing::setup_stack(Code::Nucleus, 8, 4, 900_GeV, (setup::Environment::BaseNodeType* const)nodePtr, cs); setup::StackView& view = *viewPtr; auto particle = stack->first(); -- GitLab