IAP GITLAB

Skip to content
Snippets Groups Projects
Commit abcf2348 authored by ralfulrich's avatar ralfulrich
Browse files

output

parent 528d16d2
No related branches found
No related tags found
1 merge request!369Resolve "Trajectories crossing observation plane"
Pipeline #4970 passed
...@@ -118,10 +118,12 @@ TEST_CASE("Pythia8Interface", "modules") { ...@@ -118,10 +118,12 @@ TEST_CASE("Pythia8Interface", "modules") {
corsika::pythia8::Decay decay(particleList); corsika::pythia8::Decay decay(particleList);
CORSIKA_LOG_INFO("stack: {} {}", stack.asString(), particle.asString());
[[maybe_unused]] const TimeType time = decay.getLifetime(particle); [[maybe_unused]] const TimeType time = decay.getLifetime(particle);
double const gamma = particle.getEnergy() / get_mass(Code::PiPlus); double const gamma = particle.getEnergy() / get_mass(Code::PiPlus);
CHECK(time == get_lifetime(Code::PiPlus) * gamma); CHECK(time == get_lifetime(Code::PiPlus) * gamma);
decay.doDecay(*secViewPtr); decay.doDecay(*secViewPtr);
CORSIKA_LOG_INFO("piplus->{}", stack.asString());
CHECK(stack.getEntries() == 3); // piplus, muplu, numu CHECK(stack.getEntries() == 3); // piplus, muplu, numu
auto const pSum = sumMomentum(view, cs); auto const pSum = sumMomentum(view, cs);
CHECK((pSum - plab).getNorm() / 1_GeV == Approx(0).margin(1e-4)); CHECK((pSum - plab).getNorm() / 1_GeV == Approx(0).margin(1e-4));
......
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