From 1838600d2c38286db8d6eefc1ce2f529c4a81c11 Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@tu-dortmund.de>
Date: Tue, 14 Jul 2020 22:04:43 +0200
Subject: [PATCH] write fitparams to file
---
Processes/CONEXSourceCut/CONEXSourceCut.cc | 15 +++++++++++++++
1 file changed, 15 insertions(+)
diff --git a/Processes/CONEXSourceCut/CONEXSourceCut.cc b/Processes/CONEXSourceCut/CONEXSourceCut.cc
index 8039d8ed8..94059e886 100644
--- a/Processes/CONEXSourceCut/CONEXSourceCut.cc
+++ b/Processes/CONEXSourceCut/CONEXSourceCut.cc
@@ -164,6 +164,21 @@ void CONEXSourceCut::SolveCE() {
file << X[i] << " " << N[i] << " " << dEdX[i] << " " << Mu[i] << " " << dMu[i] << " "
<< 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,
--
GitLab