Jpp  pmt_effective_area_update
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Functions
JMechanics.cc File Reference

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.

Functions

int main (int argc, char **argv)
 

Detailed Description

Auxiliary application to fit mechanical constants.

Author
mdejong

Definition in file JMechanics.cc.

Function Documentation

int main ( int  argc,
char **  argv 
)

Definition at line 24 of file JMechanics.cc.

25 {
26  using namespace std;
27  using namespace JPP;
28 
29  string inputFile;
30  string outputFile;
31  int debug;
32 
33  try {
34 
35  JParser<> zap("Auxiliary application to fit mechanical constants.");
36 
37  zap['f'] = make_field(inputFile, "input file; see https://git.km3net.de/calibration/input_tables");
38  zap['o'] = make_field(outputFile) = "mechanics.root";
39  zap['d'] = make_field(debug) = 1;
40 
41  zap(argc, argv);
42  }
43  catch(const exception &error) {
44  FATAL(error.what() << endl);
45  }
46 
47 
48  const double error = 0.2;
49 
50  TGraphErrors g1;
51 
52  double xmax = 0.0;
53  double ymax = 0.0;
54 
55  {
56  ifstream in(inputFile.c_str());
57 
58  string buffer;
59 
60  getline(in, buffer);
61 
62  DEBUG("Header " << buffer << endl);
63 
64  const locale loc(in.getloc(), new JWhiteSpacesFacet(in.getloc(), ", \t\r\n"));
65 
66  in.imbue(loc);
67 
68  for (double x, y; in >> x >> y; ) {
69 
70  AddPoint(&g1, x, y, 0.0, error);
71 
72  if (x > xmax) {
73  xmax = x;
74  ymax = y;
75  }
76  }
77 
78  in.close();
79  }
80 
81  DEBUG("Number of points " << g1.GetN() << endl);
82 
83 
84  TF1 f1("f1", "[0]*(x + [1]*log(1.0 - [2]*x))");
85 
86  f1.SetParameter(0, ymax / xmax);
87  f1.SetParameter(1, 0.25 * xmax);
88  f1.SetParameter(2, 0.50 / xmax);
89 
90  g1.Fit(&f1);
91  g1.SetName("g1");
92 
93  cout << inputFile << endl;
94 
95  cout << "b [m] " << FIXED(7,3) << f1.GetParameter(1) << endl;
96  cout << "a [m^-1] " << FIXED(7,5) << f1.GetParameter(2) << endl;
97 
98  TFile out(outputFile.c_str(), "recreate");
99 
100  out << g1;
101 
102  out.Write();
103  out.Close();
104 }
Utility class to parse command line options.
Definition: JParser.hh:1500
void AddPoint(TGraph *g1, const Double_t x, const Double_t y)
Add point to TGraph.
Auxiliary data structure for floating point format specification.
Definition: JManip.hh:446
string outputFile
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1961
std::istream & getline(std::istream &in, JString &object)
Read string from input stream until end of line.
Definition: JString.hh:478
int debug
debug level
Definition: JSirene.cc:63
#define FATAL(A)
Definition: JMessage.hh:67
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
char * loc(char *orig)
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
Definition: JCanberra.sh:40
Double_t g1(const Double_t x)
Function.
Definition: JQuantiles.cc:25