53 zap[
'P'] =
make_field(oscProbTable,
"oscillation probability table file");
59 catch(
const exception& error) {
60 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.
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Implementation of event-weight factor interface.
*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.
double E
Energy [GeV] (either MC truth or reconstructed)
double len() const
Get length.
JEvtWeightFactorFunction< JAtmosphericNeutrinoFlux, JFlux > make_atmosphericNeutrinoFluxFunction(const std::string &oscProbTableFile, const JOscParametersInterface< double > &oscParameters)
Auxiliary method for creating an interface to an atmospheric neutrino flux function using an oscillat...
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Implementation of event-weight factor for multiple particle types.
Auxiliary data structure for sequence of same character.
Auxiliary base class for list of file names.
Data structure for single set of oscillation parameters.
std::vector< filescanner_type >::iterator iterator
const JHead & getHeader() const
Get header.
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.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
#define DEBUG(A)
Message macros.