35 int main(
int argc,
char **argv)
39 using namespace KM3NETDAQ;
50 JParser<> zap(
"Auxiliary program to store Monte Carlo true muon in format for subsequent fits.");
52 zap[
'f'] =
make_field(inputFile,
"output of JTriggerEfficiency[.sh]");
54 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
59 catch(
const exception& error) {
60 FATAL(error.what() << endl);
71 while (inputFile.hasNext()) {
73 STATUS(
"event: " << setw(10) << inputFile.getCounter() <<
'\r');
DEBUG(endl);
78 const Evt*
event = ps;
85 if (muon != event->mc_trks.end()) {
Utility class to parse command line options.
int main(int argc, char *argv[])
ROOT TTree parameter settings of various packages.
JTrack3E getTrack(const Trk &track)
Get track.
Synchronously read DAQ events and Monte Carlo events (and optionally other events).
bool has_neutrino(const Evt &evt)
Test whether given event has an incoming neutrino.
bool is_muon(const Trk &track)
Test whether given track is a (anti-)muon.
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.
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
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...
Vec pos
postion of the track at time t [m]
JFit getFit(const int id, const JMODEL::JString &string)
Get fit parameters of string.
General purpose class for object reading from a list of file names.
Utility class to parse command line options.
const JLimit & getLimit() const
Get limit.
const Trk & get_neutrino(const Evt &evt)
Get incoming neutrino.
General purpose class for multiple pointers.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
#define DEBUG(A)
Message macros.