IAP GITLAB

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

changed interaction print outs

parent 68ef9aa4
No related branches found
No related tags found
2 merge requests!234WIP: Initial example of python as script language from C++,!204Resolve "boost & coordinate system in process::sibyll::Interaction"
Pipeline #1984 passed with warnings
...@@ -330,11 +330,12 @@ namespace corsika::process::sibyll { ...@@ -330,11 +330,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=" << Ecm / 1_GeV << " Ecm_final=" << Ecm_final / 1_GeV << "Ecm_initial(per nucleon)=" << Ecm / 1_GeV
<< endl << " Ecm_final(per nucleon)=" << Ecm_final * 2. / (get_nwounded() + 1) / 1_GeV << endl
<< "Elab_initial=" << eProjectileLab / 1_GeV << "Elab_initial=" << Etot / 1_GeV
<< " Elab_final=" << Elab_final / 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 << 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;
......
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