diff --git a/Processes/Sibyll/Interaction.cc b/Processes/Sibyll/Interaction.cc
index 00ac5976813da70061aa9feca2f42d7df84651b8..bb4e93d80cd8b95dbcdf1becdb2495daff9c8593 100644
--- a/Processes/Sibyll/Interaction.cc
+++ b/Processes/Sibyll/Interaction.cc
@@ -327,11 +327,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;