From 5a04089c6ec0be7a7f88db67af6ab5063f3e3035 Mon Sep 17 00:00:00 2001
From: Felix Riehn <felix@matilda>
Date: Sat, 1 Aug 2020 14:02:00 +0100
Subject: [PATCH] changed interaction print outs
---
Processes/Sibyll/Interaction.cc | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/Processes/Sibyll/Interaction.cc b/Processes/Sibyll/Interaction.cc
index 385e35ac2..a3ddcd492 100644
--- a/Processes/Sibyll/Interaction.cc
+++ b/Processes/Sibyll/Interaction.cc
@@ -330,11 +330,12 @@ namespace corsika::process::sibyll {
Ecm_final += psib.GetEnergy();
}
cout << "conservation (all GeV):" << endl
- << "Ecm_initial=" << Ecm / 1_GeV << " Ecm_final=" << Ecm_final / 1_GeV
- << endl
- << "Elab_initial=" << eProjectileLab / 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 / eProjectileLab / get_nwounded() - 1) * 100
+ << " diff (%)=" << (Elab_final / Etot / get_nwounded() - 1) * 100
+ << " E in nucleons=" << constants::nucleonMass * get_nwounded() / 1_GeV
<< endl
<< "Plab_initial=" << (pProjectileLab / 1_GeV).GetComponents()
<< ", Plab_final=" << (Plab_final / 1_GeV).GetComponents() << endl;
--
GitLab