38 int main(
int argc,
char **argv)
42 using namespace KM3NETDAQ;
54 JParser<> zap(
"Auxiliary program to store Monte Carlo true muon in format for subsequent fits.");
58 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
64 catch(
const exception& error) {
65 FATAL(error.what() << endl);
75 while (inputFile.hasNext()) {
77 STATUS(
"event: " << setw(10) << inputFile.getCounter() <<
'\r');
DEBUG(endl);
82 const Evt*
event = ps;
89 if(!take_electron) fermion = find_if(event->mc_trks.begin(),
event->mc_trks.end(),
is_neutrino);
90 else fermion = find_if(event->mc_trks.begin(),
event->mc_trks.end(),
is_electron);
92 if (fermion != event->mc_trks.end()) {
95 double time = converter.
putTime();
102 electron_dir *= shower_elongation;
103 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.
int main(int argc, char *argv[])
ROOT TTree parameter settings of various packages.
Synchronously read DAQ events and Monte Carlo events (and optionally other events).
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)...
Auxiliary class to convert DAQ hit time to/from Monte Carlo hit time.
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.
JDirection3D getDirection(const Vec &dir)
Get direction.
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.
JPosition3D getPosition(const Vec &pos)
Get position.
double putTime() const
Get Monte Carlo time minus DAQ/trigger time.
General purpose messaging.
Auxiliary include file for time conversion between DAQ/trigger hit and Monte Carlo hit...
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.
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.
const double getInverseSpeedOfLight()
Get inverse speed of light.
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.