35 int main(
int argc,
char **argv)
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.
int main(int argc, char *argv[])
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...
I/O formatting auxiliaries.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Implementation of event-weight factor for multiple particle types.
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.
Data structure for single set of oscillation parameters.
std::vector< filescanner_type >::iterator iterator
Utility class to parse command line options.
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.
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.
#define DEBUG(A)
Message macros.