32 int main(
int argc,
char **argv)
44 JParser<> zap(
"Example program to plot acoustic fit results.");
46 zap[
'f'] =
make_field(inputFile,
"input file (output of JKatoomba[.sh])");
47 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
53 catch(
const exception &error) {
54 FATAL(error.what() << endl);
58 TH1D
h1(
"chi2/NDF", NULL, 500, 0.0, 10.0);
69 while (inputFile.hasNext()) {
71 STATUS(
"event: " << setw(10) << inputFile.getCounter() <<
'\r');
DEBUG(endl);
73 const JEvt* evt = inputFile.next();
82 for (JEvt::const_iterator i = evt->begin(); i != evt->end(); ++i) {
85 const double tx = i->tx * 1.0e3;
86 const double ty = i->ty * 1.0e3;
87 const double vs = i->vs * 1.0e2;
88 const double ts = sqrt(tx*tx + ty*ty);
93 GO[id].put(t1, atan2(ty, tx));
110 out << JGraph(i->second,
MAKE_CSTRING(
"G[" <<
FILL(4,
'0') << i->first <<
"].orientation"));
114 out << JGraph(i->second,
MAKE_CSTRING(
"G[" <<
FILL(4,
'0') << i->first <<
"].amplitude"));
118 out << JGraph(i->second,
MAKE_CSTRING(
"G[" <<
FILL(4,
'0') << i->first <<
"].stretching"));
Utility class to parse command line options.
int npar
number of fit parameters
int main(int argc, char *argv[])
Data structure for graph data.
ROOT TTree parameter settings.
#define MAKE_CSTRING(A)
Make C-string.
then for HISTOGRAM in h0 h1
Dynamic ROOT object management.
Auxiliary data structure to build TGraph.
double UNIXTimeStop
stop time
Auxiliary class for defining the range of iterations of objects.
Auxiliary class to manage set of compatible ROOT objects (e.g. histograms) using unique keys...
I/O formatting auxiliaries.
double ndf
weighed number of degrees of freedom
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
General purpose messaging.
Auxiliary data structure for sequence of same character.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
double UNIXTimeStart
start time
General purpose class for object reading from a list of file names.
Utility class to parse command line options.
const JLimit & getLimit() const
Get limit.
#define DEBUG(A)
Message macros.
Double_t g1(const Double_t x)
Function.