Auxiliary application to fit mechanical constants.
More...
#include <string>
#include <iostream>
#include <iomanip>
#include <fstream>
#include "TROOT.h"
#include "TFile.h"
#include "TGraphErrors.h"
#include "TF1.h"
#include "JLang/JWhiteSpacesFacet.hh"
#include "JROOT/JRootToolkit.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
Go to the source code of this file.
|
int | main (int argc, char **argv) |
|
Auxiliary application to fit mechanical constants.
- Author
- mdejong
Definition in file JMechanics.cc.
int main |
( |
int |
argc, |
|
|
char ** |
argv |
|
) |
| |
Definition at line 24 of file JMechanics.cc.
35 JParser<> zap(
"Auxiliary application to fit mechanical constants.");
37 zap[
'f'] =
make_field(inputFile,
"input file; see https://git.km3net.de/calibration/input_tables");
43 catch(
const exception &error) {
44 FATAL(error.what() << endl);
48 const double error = 0.2;
56 ifstream
in(inputFile.c_str());
62 DEBUG(
"Header " << buffer << endl);
64 const locale
loc(
in.getloc(),
new JWhiteSpacesFacet(
in.getloc(),
", \t\r\n"));
68 for (
double x, y;
in >> x >> y; ) {
81 DEBUG(
"Number of points " <<
g1.GetN() << endl);
84 TF1 f1(
"f1",
"[0]*(x + [1]*log(1.0 - [2]*x))");
86 f1.SetParameter(0, ymax / xmax);
87 f1.SetParameter(1, 0.25 * xmax);
88 f1.SetParameter(2, 0.50 / xmax);
93 cout << inputFile << endl;
95 cout <<
"b [m] " <<
FIXED(7,3) << f1.GetParameter(1) << endl;
96 cout <<
"a [m^-1] " <<
FIXED(7,5) << f1.GetParameter(2) << endl;
Utility class to parse command line options.
void AddPoint(TGraph *g1, const Double_t x, const Double_t y)
Add point to TGraph.
Auxiliary data structure for floating point format specification.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
std::istream & getline(std::istream &in, JString &object)
Read string from input stream until end of line.
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Double_t g1(const Double_t x)
Function.