35 int main(
int argc,
char **argv)
50 zap[
'P'] =
make_field(oscProbTable,
"8D oscillation probability table file");
56 catch(
const exception& error) {
57 FATAL(error.what() << endl);
77 if (scanner->simul.size() > 0) {
78 STATUS(
"Scanning " << scanner->simul[0].program <<
" files..." << endl);
87 scanner->setFlux(multiFlux);
93 RIGHT(15) <<
"E [GeV]" <<
94 RIGHT(15) <<
"costh" <<
95 RIGHT(10) <<
"flux " <<
96 LEFT (30) <<
"[GeV^-1 * m^-2 * sr^-1 * s^-1]" << endl);
98 while (scanner->hasNext()) {
100 const Evt*
event = scanner->next();
103 const double E = neutrino.
E;
104 const double costh = -neutrino.
dir.
z / neutrino.
dir.
len();
105 const double weight = scanner->getWeight(*event);
107 STATUS(
LEFT (10) << scanner->getCounter() <<
FILL(10,
' ') <<
' ' <<
108 FIXED ( 5, 3) << right << E <<
FILL(10,
' ') <<
' ' <<
109 FIXED ( 5, 3) << right << costh <<
FILL(10,
' ') <<
' ' <<
Utility class to parse command line options.
int main(int argc, char *argv[])
Explicit template specialization for 6D interpolation table in:
Auxiliary class for interpolating an oscillation probability table in 8 dimensions in terms of: ...
const JHead & getHeader() const
Get header.
Data structure for 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.
I/O formatting auxiliaries.
#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...
General purpose messaging.
Auxiliary data structure for sequence of same character.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
Explicit template specialization for 2D interpolation table in:
Auxiliary base class for list of file names.
std::vector< filescanner_type >::iterator iterator
Utility class to parse command line options.
Auxiliary class for organising Monte Carlo file scanners associated with event weighters.
const Trk & get_neutrino(const Evt &evt)
Get incoming neutrino.
Auxiliaries for defining the range of iterations of objects.
The Trk class represents a Monte Carlo (MC) particle as well as a reconstructed track/shower.
Auxiliary data structure for floating point format specification.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
std::vector< double > weight
#define DEBUG(A)
Message macros.