IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 88539df0 authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by Ralf Ulrich
Browse files

increased CONEX output precision

(cherry picked from commit 7888ac08)
parent 545512cc
No related branches found
No related tags found
1 merge request!363Increased CONEX output precision
...@@ -260,14 +260,13 @@ namespace corsika { ...@@ -260,14 +260,13 @@ namespace corsika {
::conex::get_shower_hadron_(icuth, nX, Hadrons[0]); ::conex::get_shower_hadron_(icuth, nX, Hadrons[0]);
std::ofstream file{"conex_output.txt"}; std::ofstream file{"conex_output.txt"};
file << fmt::format("#{:>8} {:>13} {:>13} {:>13} {:>13} {:>13} {:>13} {:>13}\n", "X", file << fmt::format("#{:>10} {:>13} {:>13} {:>13} {:>13} {:>13} {:>13} {:>13}\n", "X",
"N", "dEdX", "Mu", "dMu", "Photon", "El", "Had"); "N", "dEdX", "Mu", "dMu", "Photon", "El", "Had");
for (int i = 0; i < nX; ++i) { for (int i = 0; i < nX; ++i) {
file << fmt::format( file << fmt::format(
" {:>8.2f} {:>13.3} {:>13.3} {:>13.3} {:>13.3} {:>13.3} {:>13.3} {:>13.3}\n", " {:>10.2f} {:.5e} {:.5e} {:.5e} {:.5e} {:.5e} {:.5e} {:.5e}\n",
X[i], N[i], dEdX[i], Mu[i], dMu[i], Photon[i], Electrons[i], Hadrons[i]); X[i], N[i], dEdX[i], Mu[i], dMu[i], Photon[i], Electrons[i], Hadrons[i]);
} }
file.close();
std::ofstream fitout{"conex_fit.txt"}; std::ofstream fitout{"conex_fit.txt"};
fitout << fitpars[1 - 1] << " # log10(eprima/eV)" << std::endl; fitout << fitpars[1 - 1] << " # log10(eprima/eV)" << std::endl;
...@@ -283,7 +282,6 @@ namespace corsika { ...@@ -283,7 +282,6 @@ namespace corsika {
fitout << fitpars[11 - 1] << " # phi" << std::endl; fitout << fitpars[11 - 1] << " # phi" << std::endl;
fitout << fitpars[12 - 1] << " # inelasticity 1st int." << std::endl; fitout << fitpars[12 - 1] << " # inelasticity 1st int." << std::endl;
fitout << fitpars[13 - 1] << " # ???" << std::endl; fitout << fitpars[13 - 1] << " # ???" << std::endl;
fitout.close();
} }
inline HEPEnergyType CONEXhybrid::getEnergyEM() const { return energy_em_; } inline HEPEnergyType CONEXhybrid::getEnergyEM() const { return energy_em_; }
......
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