39 using namespace KM3NETDAQ;
43 JMultipleFileScanner<Evt> inputFile;
52 JParser<> zap(
"Example program to analyse track fit results from AAnet formatted data.");
55 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
63 catch(
const exception& error) {
64 FATAL(error.what() << endl);
70 if (detectorFile !=
"") {
72 load(detectorFile, detector);
74 catch(
const JException& error) {
79 const JModuleRouter router(detector);
83 TH1D hx(
"hx", NULL, 100, -3.0, +2.3);
84 TH1D hd(
"hd", NULL, 100, 0.0, 10.0);
85 TH1D ht(
"ht", NULL, 100, -100.0, 100.0);
86 TH1D he(
"he", NULL, 100, -5.0, +5.0);
87 TH1D h1(
"h1", NULL, 100, -50.0, +50.0);
90 while (inputFile.hasNext()) {
92 STATUS(
"event: " << setw(10) << inputFile.getCounter() <<
'\r');
DEBUG(endl);
94 const Evt* evt = inputFile.next();
110 hx.Fill(log10(
getAngle(ta.getDirection(), tb.getDirection())));
111 hd.Fill((ta.getPosition() - tb.getPosition()).getLength());
112 ht.Fill( ta.getT() - tb.getT());
113 he.Fill(log10(tb.getE()/ta.getE()));
118 if (router.hasModule(i->dom_id)) {
120 const JHitL0 hit =
getHit(*i, router);
122 h1.Fill(hit.getT() - tb.getT(hit.getPosition()));
Utility class to parse command line options.
static const JGeane gWater(2.67e-1 *JTOOLS::DENSITY_SEA_WATER, 3.4e-4 *JTOOLS::DENSITY_SEA_WATER)
Function object for Energy loss of muon in sea water.
const Trk & get_muon(const Evt &evt)
Get first muon from the event tracklist.
double getAngle(const JFirst_t &first, const JSecond_t &second)
Get space angle between objects.
JTrack3E getTrack(const Trk &track)
Get track.
JLimit JLimit_t
Type definition of limit.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
void load(const JString &file_name, JDetector &detector)
Load detector from input file.
double getTimeSinceRTS(const int frame_index)
Get time in ns since last RTS for a given frame index.
End of muon fit applications.
bool has_muon(const Evt &evt)
Test whether given event has a muon.
Start of muon fit applications.
const JLimit & getLimit() const
Get limit.
JHitL0 getHit(const Hit &hit)
Get transformation.
#define DEBUG(A)
Message macros.