26 int main(
int argc,
char **argv)
31 typedef JGrid<double> JGrid_t;
52 zap[
'f'] =
make_field(inputFile,
"oscillation probability table file");
54 zap[
'X'] =
make_field(abscissaType,
"abscissa type") =
58 zap[
'x'] =
make_field(abscissaGrid,
"abscissa binning");
59 zap[
'Y'] =
make_field(ordinateType,
"ordinate type") =
63 zap[
'y'] =
make_field(ordinateGrid,
"ordinate binning");
65 = JOscParameters(
false);
66 zap[
'#'] =
make_field(channel,
"oscillation channel")
67 = JOscChannel(JFlavour_t::MUON, JFlavour_t::MUON, JChargeParity_t::PARTICLE);
73 catch (
const exception& error) {
74 FATAL(error.what() << endl);
78 const JOscProbInterpolator<> interpolator(inputFile.c_str(),
parameters);
81 abscissaGrid.getSize(), abscissaGrid.getXmin(), abscissaGrid.getXmax(),
82 ordinateGrid.getSize(), ordinateGrid.getXmin(), ordinateGrid.getXmax());
84 const JOscillogram oscillogram(abscissaType, abscissaGrid,
85 ordinateType, ordinateGrid,
86 channel, &interpolator);
90 for (
unsigned int ix = 0; ix < abscissaGrid.getSize(); ++ix) {
91 for (
unsigned int iy = 0; iy < ordinateGrid.getSize(); ++iy) {
93 const double x = abscissaGrid.getX(ix);
94 const double y = ordinateGrid.getX(iy);
96 const double P = oscillogram.getP(ix, iy);
100 h0.SetBinContent(ix+1, iy+1, P);
Utility class to parse command line options.
int main(int argc, char *argv[])
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
Auxiliary data structure for floating point format specification.
I/O formatting auxiliaries.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
General purpose messaging.
then for APP in event gandalf start energy
Utility class to parse command line options.
#define DEBUG(A)
Message macros.