35 int main(
int argc,
char **argv)
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.
int main(int argc, char *argv[])
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
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.
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...
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.
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.