30int main(
int argc,
char **argv)
45 zap[
'P'] =
make_field(oscProbTable,
"oscillation probability table file");
51 catch(
const exception& error) {
52 FATAL(error.what() << endl);
67 scanner->setFlux(
flux);
69 if (scanner->simul.size() > 0) {
70 STATUS(
"Scanning " << scanner->simul[0].program <<
" files..." << endl);
74 RIGHT(15) <<
"E [GeV]" <<
' ' <<
75 RIGHT(15) <<
"costh" <<
' ' <<
76 RIGHT(40) <<
"flux [GeV^-1 * m^-2 * sr^-1 * s^-1]" <<
' ' <<
77 RIGHT(20) <<
"rate [Hz]" << endl);
79 while (scanner->hasNext()) {
81 const Evt*
event = scanner->next();
84 const double E = neutrino.
E;
85 const double costh = -neutrino.
dir.
z / neutrino.
dir.
len();
86 const double F =
flux.getFlux(*event);
87 const double weight = scanner->getWeight(*event);
89 STATUS(
LEFT (10) << scanner->getCounter() <<
' ' <<
90 FIXED (15, 3) << E <<
' ' <<
91 FIXED (15, 3) << costh <<
' ' <<
int main(int argc, char **argv)
General purpose messaging.
#define DEBUG(A)
Message macros.
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
I/O formatting auxiliaries.
ROOT TTree parameter settings of various packages.
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.
const Trk & get_neutrino(const Evt &evt)
Get incoming neutrino.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Auxiliary data structure for floating point format specification.
Implementation of atmospheric neutrino flux using official KM3NeT atmospheric flux function.
Helper class for oscillation probability calculators.
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Auxiliary class for organising Monte Carlo file scanners associated with event weighters.
std::vector< filescanner_type >::iterator iterator
Auxiliary base class for list of file names.
Auxiliary data structure for alignment of data.
Auxiliary data structure for floating point format specification.
The Trk class represents a Monte Carlo (MC) particle as well as a reconstructed track/shower.
double E
Energy [GeV] (either MC truth or reconstructed)
double len() const
Get length.