27 int main(
int argc,
char **argv)
32 typedef JAbstractHistogram<Double_t> JHistogram_t;
46 zap[
'f'] =
make_field(inputFile,
"8D oscillation probability table file");
48 zap[
'x'] =
make_field(
X,
"histogram binning in log10(E)");
49 zap[
'y'] =
make_field(
Y,
"histogram binning in cos(theta)");
51 zap[
'#'] =
make_field(channel,
"Oscillation channel") = JOscChannel(JOscChannel::MUON, JOscChannel::MUON, JOscChannel::PARTICLE);
56 catch (
const exception& error) {
57 FATAL(error.what() << endl);
61 TH2D h0(
"h0",
"; log_{10}(E); cos(#theta)",
62 X.getNumberOfBins(),
X.getLowerLimit(),
X.getUpperLimit(),
63 Y.getNumberOfBins(),
Y.getLowerLimit(),
Y.getUpperLimit());
65 const JOscProbTable6D<> g6(inputFile.c_str());
68 const JOscProbInterpolator2D<> interpolator(&g2);
71 RIGHT(15) <<
"cos(theta)" <<
72 RIGHT(15) <<
"P" << endl);
74 for (Int_t ix = 1; ix <= h0.GetNbinsX(); ++ix) {
75 for (Int_t iy = 1; iy <= h0.GetNbinsY(); ++iy) {
77 const Double_t log10E = h0.GetXaxis()->GetBinCenter(ix);
78 const Double_t costh = h0.GetYaxis()->GetBinCenter(iy);
79 const Double_t
P = interpolator(channel, log10E, costh);
82 FIXED(15,3) << costh <<
85 h0.SetBinContent(ix, iy, 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
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
Auxiliary data structure for floating point format specification.
I/O formatting auxiliaries.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
T pow(const T &x, const double y)
Power .
then break fi done getCenter read X Y Z let X
Auxiliary methods to convert data members or return values of member methods of a set of objects to a...
General purpose messaging.
Utility class to parse command line options.
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG