25int main(
int argc,
char **argv)
51 zap[
'f'] =
make_field(inputFile,
"oscillation probability table file");
53 zap[
'X'] =
make_field(abscissaType,
"abscissa type") =
57 zap[
'x'] =
make_field(abscissaGrid,
"abscissa binning");
58 zap[
'Y'] =
make_field(ordinateType,
"ordinate type") =
62 zap[
'y'] =
make_field(ordinateGrid,
"ordinate binning");
63 zap[
'@'] =
make_field(parameters,
"oscillation parameters")
65 zap[
'#'] =
make_field(channel,
"oscillation channel")
66 =
JOscChannel(JFlavour_t::MUON, JFlavour_t::MUON, JChargeParity_t::PARTICLE);
72 catch (
const exception& error) {
73 FATAL(error.what() << endl);
80 abscissaGrid.getSize(), abscissaGrid.getXmin(), abscissaGrid.getXmax(),
81 ordinateGrid.getSize(), ordinateGrid.getXmin(), ordinateGrid.getXmax());
83 const JOscillogram oscillogram(abscissaType, abscissaGrid,
84 ordinateType, ordinateGrid,
85 channel, interpolator);
89 for (
int ix = 0; ix < abscissaGrid.getSize(); ++ix) {
90 for (
int iy = 0; iy < ordinateGrid.getSize(); ++iy) {
92 const double x = abscissaGrid.getX(ix);
93 const double y = ordinateGrid.getX(iy);
95 const double P = oscillogram.
getP(ix, iy);
99 h0.SetBinContent(ix+1, iy+1, P);