36 int main(
int argc,
char **argv)
40 using namespace KM3NETDAQ;
52 JParser<> zap(
"Auxiliary program to store Monte Carlo true muon in format for subsequent fits.");
56 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
62 catch(
const exception& error) {
63 FATAL(error.what() << endl);
73 while (inputFile.hasNext()) {
75 STATUS(
"event: " << setw(10) << inputFile.getCounter() <<
'\r');
DEBUG(endl);
80 const Evt*
event = ps;
87 if(!take_electron) fermion = find_if(event->mc_trks.begin(),
event->mc_trks.end(),
is_neutrino);
88 else fermion = find_if(event->mc_trks.begin(),
event->mc_trks.end(),
is_electron);
90 if (fermion != event->mc_trks.end()) {
93 double time = converter.
putTime();
100 electron_dir *= shower_elongation;
101 position.
add(electron_dir);
bool is_electron(const Trk &track)
Test whether given track is a (anti-)electron.
Utility class to parse command line options.
ROOT TTree parameter settings.
Synchronously read DAQ events and Monte Carlo events (and optionally other events).
double putTime() const
Get Monte Carlo minus DAQ/trigger hit time.
Recording of objects on file according a format that follows from the file name extension.
Auxiliary class to synchronously read DAQ events and Monte Carlo events (and optionally other events)...
JTime & add(const JTime &value)
Addition operator.
Head getHeader(const JMultipleFileScanner_t &file_list)
Get Monte Carlo header.
bool is_neutrino(const Trk &track)
Test whether given track is a neutrino.
Auxiliary class for defining the range of iterations of objects.
double getE() const
Get energy.
Data structure for vector in three dimensions.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
static const JGeanz geanz(1.85, 0.62, 0.54)
Function object for longitudinal EM-shower profile.
General purpose messaging.
JFit getFit(const JHistory &history, const JTrack3D &track, const double Q, const int NDF, const double energy=0.0, const int status=0)
Get fit.
Data structure for set of track fit results.
double getMaximum(const double E) const
Get depth of shower maximum.
General purpose class for object reading from a list of file names.
Utility class to parse command line options.
JDirection3D getDirection(const Vec &v)
Get direction.
Auxiliary class to convert DAQ/trigger hit time to/from Monte Carlo hit time.
const JLimit & getLimit() const
Get limit.
Longitudinal emission profile EM-shower.
General purpose class for multiple pointers.
JVector3D & add(const JVector3D &vector)
Add vector.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
#define DEBUG(A)
Message macros.
JPosition3D getPosition(const Vec &v)
Get position.
int main(int argc, char *argv[])