IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 9c746aba authored by Felix Riehn's avatar Felix Riehn
Browse files

Merge branch '272-boost-coordinate-system-in-process-sibyll-interaction' of...

Merge branch '272-boost-coordinate-system-in-process-sibyll-interaction' of gitlab.ikp.kit.edu:AirShowerPhysics/corsika into 272-boost-coordinate-system-in-process-sibyll-interaction
parents 405cbe35 5a04089c
No related branches found
No related tags found
No related merge requests found
...@@ -327,11 +327,12 @@ namespace corsika::process::sibyll { ...@@ -327,11 +327,12 @@ namespace corsika::process::sibyll {
Ecm_final += psib.GetEnergy(); Ecm_final += psib.GetEnergy();
} }
cout << "conservation (all GeV):" << endl cout << "conservation (all GeV):" << endl
<< "Ecm_initial(per nucleon)=" << Ecm / 1_GeV << " Ecm_final(per nucleon)=" << "Ecm_initial(per nucleon)=" << Ecm / 1_GeV
<< Ecm_final * 2. / (get_nwounded() + 1) / 1_GeV << endl << " Ecm_final(per nucleon)=" << Ecm_final * 2. / (get_nwounded() + 1) / 1_GeV << endl
<< "Elab_initial=" << Etot / 1_GeV << " Elab_final=" << Elab_final / 1_GeV << "Elab_initial=" << Etot / 1_GeV
<< " Elab_final=" << Elab_final / 1_GeV
<< " diff (%)=" << (Elab_final / Etot / get_nwounded() - 1) * 100 << " 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 << endl
<< "Plab_initial=" << (pProjectileLab / 1_GeV).GetComponents() << "Plab_initial=" << (pProjectileLab / 1_GeV).GetComponents()
<< ", Plab_final=" << (Plab_final / 1_GeV).GetComponents() << endl; << ", Plab_final=" << (Plab_final / 1_GeV).GetComponents() << endl;
......
...@@ -203,7 +203,7 @@ TEST_CASE("SibyllInterface", "[processes]") { ...@@ -203,7 +203,7 @@ TEST_CASE("SibyllInterface", "[processes]") {
[[maybe_unused]] const GrammageType length = model.GetInteractionLength(particle); [[maybe_unused]] const GrammageType length = model.GetInteractionLength(particle);
} }
SECTION("InteractionInterface - high energy") { SECTION("InteractionInterface - high energy") {
setup::Stack stack; setup::Stack stack;
const HEPEnergyType E0 = 60_EeV; const HEPEnergyType E0 = 60_EeV;
......
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