50 zap[
'P'] =
make_field(oscProbTable,
"oscillation probability table file");
56 catch(
const exception& error) {
57 FATAL(error.what() << endl);
74 if (scanner->simul.size() > 0) {
75 STATUS(
"Scanning " << scanner->simul[0].program <<
" files..." << endl);
84 scanner->setFlux(multiFlux);
90 RIGHT(15) <<
"E [GeV]" <<
91 RIGHT(15) <<
"costh" <<
92 RIGHT(10) <<
"flux " <<
93 LEFT (30) <<
"[GeV^-1 * m^-2 * sr^-1 * s^-1]" << endl);
95 while (scanner->hasNext()) {
97 const Evt*
event = scanner->next();
100 const double E = neutrino.
E;
101 const double costh = -neutrino.
dir.
z / neutrino.
dir.
len();
102 const double weight = scanner->getWeight(*event);
104 STATUS(
LEFT (10) << scanner->getCounter() <<
FILL(10,
' ') <<
' ' <<
105 FIXED ( 5, 3) << right << E <<
FILL(10,
' ') <<
' ' <<
106 FIXED ( 5, 3) << right << costh <<
FILL(10,
' ') <<
' ' <<
Utility class to parse command line options.
const JHead & getHeader() const
Get header.
Data structure for single set of oscillation parameters.
*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)...
Auxiliary data structure for floating point format specification.
JEvtWeightFactorFunction< JFunction_t, JFlux > make_fluxFunction(const JFunction_t &flux)
Auxiliary method for creating an interface to a flux function.
double E
Energy [GeV] (either MC truth or reconstructed)
double len() const
Get length.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Implementation of event-weight factor for multiple particle types.
Implementation of atmospheric neutrino flux using official KM3NeT atmospheric flux function...
Auxiliary data structure for sequence of same character.
Auxiliary base class for list of file names.
std::vector< filescanner_type >::iterator iterator
Auxiliary class for organising Monte Carlo file scanners associated with event weighters.
const Trk & get_neutrino(const Evt &evt)
Get incoming neutrino.
The Trk class represents a Monte Carlo (MC) particle as well as a reconstructed track/shower.
Auxiliary data structure for floating point format specification.
Template definition of a multi-dimensional oscillation probability interpolation table.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
#define DEBUG(A)
Message macros.