25int main(
int argc,
char **argv)
36 JParser<> zap(
"Auxiliary application to fit mechanical constants.");
38 zap[
'f'] =
make_field(inputFile,
"input file; see https://git.km3net.de/calibration/input_tables");
44 catch(
const exception &error) {
45 FATAL(error.what() << endl);
49 const double error = 0.2;
57 ifstream in(inputFile.c_str());
63 DEBUG(
"Header " << buffer << endl);
69 for (
double x, y; in >> x >> y; ) {
71 AddPoint(&
g1, x, y, 0.0, error);
82 DEBUG(
"Number of points " <<
g1.GetN() << endl);
85 TF1 f1(
"f1",
"[0]*(x + [1]*log(1.0 - [2]*x))");
87 f1.SetParameter(0, ymax / xmax);
88 f1.SetParameter(1, 0.25 * xmax);
89 f1.SetParameter(2, 0.50 / xmax);
94 cout << inputFile << endl;
96 cout <<
"T_x = " <<
FIXED(5,2) << f1.GetParameter(0) <<
" [s^2/m^2] * v_x^2" << endl;
98 cout <<
"b [m] " <<
FIXED(7,3) << f1.GetParameter(1) << endl;
99 cout <<
"a [m^-1] " <<
FIXED(7,5) << f1.GetParameter(2) << endl;