IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 4756a808 authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan:
Browse files

Merge branch 'conex_fit_file' into 'master'

Write CONEX fit parameters to file

See merge request AirShowerPhysics/corsika!267
parents 4d5021e2 1838600d
No related branches found
No related tags found
No related merge requests found
...@@ -164,6 +164,21 @@ void CONEXSourceCut::SolveCE() { ...@@ -164,6 +164,21 @@ void CONEXSourceCut::SolveCE() {
file << X[i] << " " << N[i] << " " << dEdX[i] << " " << Mu[i] << " " << dMu[i] << " " file << X[i] << " " << N[i] << " " << dEdX[i] << " " << Mu[i] << " " << dMu[i] << " "
<< Gamma[i] << " " << Electrons[i] << " " << Hadrons[i] << std::endl; << Gamma[i] << " " << Electrons[i] << " " << Hadrons[i] << std::endl;
} }
std::ofstream fitout{"conex_fit.txt"};
fitout << fitpars[1 - 1] << " # log10(eprima/eV)" << std::endl;
fitout << fitpars[2 - 1] << " # theta" << std::endl;
fitout << fitpars[3 - 1] << " # X1 (first interaction)" << std::endl;
fitout << fitpars[4 - 1] << " # Nmax" << std::endl;
fitout << fitpars[5 - 1] << " # X0" << std::endl;
fitout << fitpars[6 - 1] << " # P1" << std::endl;
fitout << fitpars[7 - 1] << " # P2" << std::endl;
fitout << fitpars[8 - 1] << " # P3" << std::endl;
fitout << fitpars[9 - 1] << " # chi^2 / sqrt(Nmax)" << std::endl;
fitout << fitpars[10 - 1] << " # Xmax" << std::endl;
fitout << fitpars[11 - 1] << " # phi" << std::endl;
fitout << fitpars[12 - 1] << " # inelasticity 1st int." << std::endl;
fitout << fitpars[13 - 1] << " # ???" << std::endl;
} }
CONEXSourceCut::CONEXSourceCut(geometry::Point center, CONEXSourceCut::CONEXSourceCut(geometry::Point center,
......
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