diff --git a/Processes/Sibyll/Interaction.cc b/Processes/Sibyll/Interaction.cc index bb4e93d80cd8b95dbcdf1becdb2495daff9c8593..d738dc81caa233416e74b1ca56d5e6e3d5662473 100644 --- a/Processes/Sibyll/Interaction.cc +++ b/Processes/Sibyll/Interaction.cc @@ -327,12 +327,11 @@ namespace corsika::process::sibyll { Ecm_final += psib.GetEnergy(); } cout << "conservation (all GeV):" << endl - << "Ecm_initial(per nucleon)=" << Ecm / 1_GeV - << " Ecm_final(per nucleon)=" << Ecm_final * 2. / (get_nwounded() + 1) / 1_GeV << endl - << "Elab_initial=" << Etot / 1_GeV - << " Elab_final=" << Elab_final / 1_GeV + << "Ecm_initial(per nucleon)=" << Ecm / 1_GeV << " Ecm_final(per nucleon)=" + << Ecm_final * 2. / (get_nwounded() + 1) / 1_GeV << endl + << "Elab_initial=" << Etot / 1_GeV << " Elab_final=" << Elab_final / 1_GeV << " diff (%)=" << (Elab_final / Etot / get_nwounded() - 1) * 100 - << " E in nucleons=" << constants::nucleonMass * get_nwounded() / 1_GeV + << " E in nucleons=" << constants::nucleonMass * get_nwounded() / 1_GeV << endl << "Plab_initial=" << (pProjectileLab / 1_GeV).GetComponents() << ", Plab_final=" << (Plab_final / 1_GeV).GetComponents() << endl; diff --git a/Processes/Sibyll/testSibyll.cc b/Processes/Sibyll/testSibyll.cc index 9aad0251e4324bec303d01c0e143353f628c58dd..18c0b94c55ca529232f3d84e52793b8f4892031a 100644 --- a/Processes/Sibyll/testSibyll.cc +++ b/Processes/Sibyll/testSibyll.cc @@ -203,7 +203,7 @@ TEST_CASE("SibyllInterface", "[processes]") { [[maybe_unused]] const GrammageType length = model.GetInteractionLength(particle); } - SECTION("InteractionInterface - high energy") { + SECTION("InteractionInterface - high energy") { setup::Stack stack; const HEPEnergyType E0 = 60_EeV;