41int main(
int argc,
char **argv)
73 catch(
const exception& error) {
74 FATAL(error.what() << endl);
86 if (!oscProbTable.empty()) {
90 const JOscFlux oscFlux1(atmFlux, oscProb);
91 const JOscFlux oscFlux2(interpolator, oscProb);
110 log10Egrid.getSize(), log10Egrid.getXmin(), log10Egrid.getXmax(),
111 coszGrid .getSize(), coszGrid .getXmin(), coszGrid .getXmax()));
114 log10Egrid.getSize(), log10Egrid.getXmin(), log10Egrid.getXmax(),
115 coszGrid .getSize(), coszGrid .getXmin(), coszGrid .getXmax()));
117 DEBUG(
LEFT(15) <<
"E [GeV]" <<
RIGHT(10) <<
"cosz" <<
RIGHT(25) <<
"flux1" <<
RIGHT(52) <<
"flux2 [GeV^-1 * m^-2 * sr^-1 * s^-1]" << endl);
119 for (vector<int>::const_iterator type = nuTypes.cbegin(); type != nuTypes.cend(); ++type) {
121 for (
int i = 0; i < log10Egrid.getSize(); ++i) {
123 const double X = log10Egrid.getX(i);
125 for (
int j = 0;
j < coszGrid.getSize(); ++
j) {
127 const double cosz = coszGrid.getX(
j);
129 const double F0 = flux1(*type, X, cosz);
130 const double F1 = flux2(*type, X, cosz);
133 FIXED (10,3) << cosz <<
139 h0[*type]->SetBinContent(i+1,
j+1, F0);
140 h1[*type]->SetBinContent(i+1,
j+1, F1);
int main(int argc, char **argv)
Dynamic ROOT object management.
General purpose messaging.
#define DEBUG(A)
Message macros.
#define ASSERT(A,...)
Assert macro.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition of particle types.
I/O formatting auxiliaries.
ROOT TTree parameter settings of various packages.
Template specialisation for interpolator of azimuth-averaged Honda flux table.
Data structure for single set of oscillation parameters.
Template definition of a multi-dimensional oscillation probability interpolation table.
Utility class to parse command line options.
Auxiliary class to manage set of compatible ROOT objects (e.g. histograms) using unique keys.
void Write(TDirectory &out, const bool wm=false)
Write objects to file.
T pow(const T &x, const double y)
Power .
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Auxiliary data structure for floating point format specification.
Helper class for diffuse flux factor.
Specialisation of event-weight factor interface for atmospheric neutrino flux.
void configure(const pointer_type &p)
Configure event-weight factor.
Implementation of oscillated neutrino flux.
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Auxiliary base class for list of file names.
Auxiliary data structure for floating point format specification.