Jpp 19.3.0-rc.2
the software that should make you happy
Loading...
Searching...
No Matches
Functions
JMonitorMultiplicity.cc File Reference

Auxiliary program to analyse coincidence multiplicity;. More...

#include <iostream>
#include <string>
#include "TFile.h"
#include "TH1D.h"
#include "TH2D.h"
#include "TMath.h"
#include "TROOT.h"
#include "TVectorD.h"
#include "JAAnet/JHead.hh"
#include "JCalibrate/JCalibrateK40.hh"
#include "km3net-dataformat/online/JDAQ.hh"
#include "JDAQ/JDAQEvaluator.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JDetectorAddressMap.hh"
#include "JDetector/JDetectorSupportkit.hh"
#include "JDetector/JModuleAddressMap.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
#include "JMath/JMathToolkit.hh"
#include "Jeep/JTimer.hh"
#include "JLang/JObjectMultiplexer.hh"
#include "JLang/JString.hh"
#include "JROOT/JManager.hh"
#include "JGizmo/JGizmoToolkit.hh"
#include "JROOT/JROOTClassSelector.hh"
#include "JSupport/JMeta.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JTreeScanner.hh"
#include "JSupport/JTriggerParametersSupportkit.hh"
#include "JTrigger/JHit.hh"
#include "JTrigger/JHitR0.hh"
#include "JTrigger/JMatchL0.hh"
#include "JTrigger/JSummaryRouter.hh"
#include "JTrigger/JSuperFrame1D.hh"
#include "JTrigger/JSuperFrame2D.hh"
#include "JTools/JRange.hh"

Go to the source code of this file.

Functions

int main (int argc, char **argv)
 

Detailed Description

Auxiliary program to analyse coincidence multiplicity;.

Author
mlincetto, mkarel

Definition in file JMonitorMultiplicity.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 52 of file JMonitorMultiplicity.cc.

53{
54 using namespace std;
55 using namespace KM3NETDAQ;
56 using namespace JPP;
57
58 //-----------------------------------------------------------
59 // parameter interface
60 //-----------------------------------------------------------
61
62 JMultipleFileScanner<JDAQTimesliceTypes_t> inputFile;
63 JLimit_t& numberOfEvents = inputFile.getLimit();
64 string outputFile;
65 string detectorFile;
66 Double_t TMax_ns;
67 JRange<double> rateVeto_Hz;
68 JRange<double> totVeto_ns;
69 JROOTClassSelector selector;
70 int debug;
71 int filterLevel;
72 int muteChannel;
73 bool twoDim;
74 bool monitorOccupancy;
75 bool notJoin;
76
77 try {
78
79 JParser<> zap("Auxiliary program to estimate PMT and hit multiplicities.");
80
81 zap['f'] = make_field(inputFile);
82 zap['o'] = make_field(outputFile) = "monitormultiplicity.root";
83 zap['a'] = make_field(detectorFile);
84 zap['n'] = make_field(numberOfEvents) = JLimit::max();
85 zap['T'] = make_field(TMax_ns) = 10.0; // [ns]
86 zap['V'] = make_field(rateVeto_Hz) = JRange<double>(0, 1e5); // [Hz]
87 zap['t'] = make_field(totVeto_ns) = JRange<double>(0, 1e4); // min and max value of ToT of hits to be monitored [ns]
88 zap['C'] = make_field(selector) = getROOTClassSelection<JDAQTimesliceTypes_t>();
89 zap['d'] = make_field(debug) = 1;
90 zap['M'] = make_field(muteChannel) = -1;
91 zap['F'] = make_field(filterLevel, "0 = raw data; 1 = filtered data; 2+ = only clean frames") = 1;
92 zap['D'] = make_field(twoDim, "2D mode for background subtraction");
93 zap['O'] = make_field(monitorOccupancy, "produces 2D PMT vs multiplicity plots");
94 zap['j'] = make_field(notJoin, "do not join consecutive hits on same PMT (diagnostic purpose)");
95
96 zap(argc, argv);
97 }
98 catch(const exception &error) {
99 FATAL(error.what() << endl);
100 }
101
102
103 //-----------------------------------------------------------
104 // check the input parameters
105 //-----------------------------------------------------------
106
107 JTriggerParameters parameters;
108 bool hasTriggerParameters = true;
109
110 try {
111 parameters = getTriggerParameters(inputFile);
112 NOTICE("Get trigger parameters from input." << endl);
113 }
114 catch(const JException& error) {
115 ERROR("No trigger parameters from input." << endl);
116 hasTriggerParameters = false;
117
118 }
119
120 // keep track of prescale factor for proper MC normalization.
121 int prescale = 1;
122
123 if (hasTriggerParameters) {
124
125 // for L1 (most common for sea analysis) check that time window is not larger than in trigger
126
127 if (parameters.writeL1.prescale > 0 && (selector == "JDAQTimeslice" || selector == "JDAQTimesliceL1") &&
128 parameters.TMaxLocal_ns < TMax_ns) {
129 FATAL("TMax_ns (-T) " << TMax_ns << " ns is larger than in the trigger " << parameters.TMaxLocal_ns << " ns." << endl);
130 }
131
132 if (selector == "JDAQTimeslice") {
133 // Jpp v8 legacy compatibility
134 prescale = (parameters.writeL1.prescale != 0) ? parameters.writeL1.prescale : parameters.writeL0.prescale;
135 }
136
137 if (selector == "JDAQTimesliceL0") {
138 prescale = parameters.writeL0.prescale;
139 }
140
141 if (selector == "JDAQTimesliceL1") {
142 prescale = parameters.writeL1.prescale;
143 }
144
145 if (selector == "JDAQTimesliceSN") {
146 prescale = parameters.writeSN.prescale;
147 }
148
149 }
150
151 if (prescale == 0) {
152 FATAL("[R] According to trigger parameters, the " << selector << " stream should be empty.");
153 } else {
154 NOTICE("[R] Prescale factor is " << prescale << endl);
155 }
156
157 if (!totVeto_ns.is_valid()) {
158 FATAL("Invalid time over threshold " << totVeto_ns << endl);
159 }
160
161 //-----------------------------------------------------------
162 // load the detector file
163 //-----------------------------------------------------------
164
165 JDetector detector;
166
167 try {
168 load(detectorFile, detector);
169 }
170 catch(const JException& error) {
171 FATAL(error);
172 }
173
174 if (detector.empty())
175 FATAL("Empty detector." << endl);
176
177 const JModuleRouter router(detector);
178 JSummaryRouter summaryRouter;
179
180 int detSize = getNumberOfModules(detector);
181
182 int detID = detector.getID();
183
184 JDetectorAddressMap& detectorAddressMap = getDetectorAddressMap(detID);
185
186
187 //-----------------------------------------------------------
188 // summary of input parameters
189 //-----------------------------------------------------------
190
191 NOTICE("[R] Frame time [ms] " << getFrameTime() * 1e-6 << endl);
192 NOTICE("[R] Tmax [ns] " << TMax_ns << endl);
193 NOTICE("[R] Rate range [Hz] " << rateVeto_Hz << endl);
194 NOTICE("[R] ToT range [ns] " << totVeto_ns << endl);
195 NOTICE("[R] Timeslice stream " << selector << endl);
196 NOTICE("[R] Detector file has " << detSize << " DOMs." << endl);
197
198 //-----------------------------------------------------------
199 // initialise histograms
200 //-----------------------------------------------------------
201
202 // nominal livetime + effective livetime for each DOM
203 TH1D* h_livetime = new TH1D("LiveTime", "Livetime", 1 + detSize, 0.0, 1.0 + detSize);
204
205 int ibin = 1;
206 h_livetime->GetXaxis()->SetBinLabel(ibin++, "Nominal");
207 for (JDetector::iterator module = detector.begin(); module != detector.end(); ++module) {
208 h_livetime->GetXaxis()->SetBinLabel(ibin++, TString(getLabel(*module)));
209 }
210
211 h_livetime->GetYaxis()->SetTitle("Livetime (s)");
212
213
214 // multiplicity histograms
215 const int nx = 1 + NUMBER_OF_PMTS;
216 const double xmin = -0.5;
217 const double xmax = nx - 0.5;
218
219 typedef JManager<TString, TH1D> JManager_t;
220 typedef JManager<TString, TH2D> JManager2D_t;
221 JManager_t H(new TH1D("%_H", NULL, nx, xmin, xmax)); H->Sumw2();
222 JManager_t P(new TH1D("%_P", NULL, nx, xmin, xmax)); P->Sumw2();
223
224 JManager2D_t HT(new TH2D("%_HT", NULL, nx, xmin, xmax, 100, -TMax_ns, +TMax_ns));
225
226 TH2D* CO_proto = new TH2D("%_CO", NULL, NUMBER_OF_PMTS, 0.5, 0.5 + NUMBER_OF_PMTS, NUMBER_OF_PMTS, -0.5, -0.5 + NUMBER_OF_PMTS);
227
228 if (monitorOccupancy) {
229 JModuleAddressMap memo = getDetectorAddressMap<JKM3NeT_t>().getDefaultModuleAddressMap();
230 setAxisLabels(CO_proto->GetYaxis(), memo);
231 }
232
233 JManager2D_t CO(CO_proto);
234
235 // benchmark histogram for time distribution of high-multiplicities
236 const int sn_mul_min = 6;
237 TH1D* time_distr = new TH1D("T_distr", "Coincidence (M > 6) distribution in timeslice", 100, 0, 1e8);
238
239 //-----------------------------------------------------------
240 // output file and auxiliary data structures
241 //-----------------------------------------------------------
242
243 // output file
244 TFile out(outputFile.c_str(), "recreate");
245
246 // auxiliary counters
247
248 int count_HRV = 0; // high-rate-veto counter
249 int count_FAF = 0; // fifo-almost-full counter
250 int count_URV = 0; // user-rate-veto counter
251
252 int treeMismatchCount = 0;
253
254 //-----------------------------------------------------------
255 // file pre-processing
256 //-----------------------------------------------------------
257
258 JTreeScanner<JDAQSummaryslice, JDAQEvaluator> summaryFile(inputFile);
259
260 JObjectMultiplexer<JDAQTimesliceTypes_t, JDAQTimeslice> in(inputFile, selector);
261
262 vector<JHitR0> filteredData;
263
264 typedef JSuperFrame2D<JHitR0> JSuperFrame2D_t;
265 typedef JSuperFrame1D<JHitR0> JSuperFrame1D_t;
266 const JMatchL0<JHitR0> match(TMax_ns);
267
268 counter_type counter = 0;
269
270 // process Monte Carlo aanet header
271 bool hasMCHeader = true ;
272 bool isMupage = false;
273 double liveTimeMC = 0 ;
274 double liveTimeMCErr = 0 ;
275
276 JHead mc_header;
277
278 try {
279 mc_header = (JHead)getHeader(inputFile);
280 }
281 catch (const JException& error) {
282 hasMCHeader = false;
283 }
284
285 NOTICE("[R] File source is " << (hasMCHeader ? "MC" : "DAQ") << "." << endl);
286 if (hasMCHeader) {
287 liveTimeMC = mc_header.livetime.numberOfSeconds;
288 liveTimeMCErr = mc_header.livetime.errorOfSeconds;
289 if (liveTimeMC > 0) {
290 isMupage = true;
291 NOTICE("[R] mupage live time is (" << liveTimeMC << " +/- " << liveTimeMCErr << ") s." << endl);
292 } else {
293 NOTICE("[R] aanet header found but mupage live time is zero. Likely not a mupage file.");
294 }
295 }
296
297 // retrieving and exporting run number to tree metadata;
298
299 int runNumber = 0;
300
301 if (summaryFile.hasNext()) {
302 runNumber = summaryFile.getEntry(0)->getRunNumber();
303 NOTICE("[R] Processing run #" << runNumber << endl);
304 summaryFile.rewind();
305 }
306
307 //-----------------------------------------------------------
308 // file processing: loop over timeslices
309 //-----------------------------------------------------------
310
311 for ( ; in.hasNext() && counter != inputFile.getLimit(); ++counter) {
312
313 STATUS("Entry: " << setw(10) << counter << '\r');
314
315 const JDAQTimeslice* timeslice = in.next();
316 const Long64_t index = summaryFile.find(*timeslice);
317 JDAQSummaryslice* summary = (index != -1 ? summaryFile.getEntry(index) : NULL);
318
319 summaryRouter.update(summary);
320
321 // check correspondence between summary data and timeslice data
322 if (summary != NULL) {
323 if (timeslice->getFrameIndex() != summary->getFrameIndex()) {
324 ++treeMismatchCount;
325 ERROR("[R>T] Frame indices do not match [counter=" << counter << ", timeslice=" << timeslice->getFrameIndex() << ", summaryslice=" << summary->getFrameIndex() << "]" << endl);
326 if (treeMismatchCount > 1) {
327 FATAL("ROOT TTrees not aligned. Run #" << runNumber << endl);
328 }
329 continue;
330 }
331 }
332
333 //-----------------------------------------------------------
334 // timeslice processing: loop over superframes
335 //-----------------------------------------------------------
336 for (JDAQTimeslice::const_iterator super_frame = timeslice->begin(); super_frame != timeslice->end(); ++super_frame) {
337
338 // empty superframes are not skipped because SN timeslices could have empty superframes even if the module was active;
339 // if DAQ error or FIFO full -> frame counts as dead time;
340 // if frame is missing but no errors -> frame counts as live time (SN stream has many empty frames);
341
342 int moduleID = super_frame->getModuleID();
343
344 if (!router.hasModule(moduleID)) {
345 // module is not mapped in detector file, nothing can be done
346 continue;
347 }
348
349 JModuleAddressMap moduleAddressMap = detectorAddressMap.get(moduleID);
350
351 //----------------------------------------------------------
352 // retrieval of summary data from superframe or summaryframe
353 //----------------------------------------------------------
354
355 vector<double> rate_Hz(NUMBER_OF_PMTS, 0.0);
356 JDAQFrameStatus status;
357
358 if (summary == NULL) {
359 WARNING("[R>T>F] Missing summary for timeslice." << endl);
360 // PMT rates estimated from timeslice data
361 for (JDAQSuperFrame::const_iterator i = super_frame->begin(); i != super_frame->end(); ++i) {
362 rate_Hz[i->getPMT()] += 1e9 / getFrameTime();
363 }
364
365 // frame status retrieved from superframe
366 status = super_frame->getDAQFrameStatus();
367
368 } else {
369
370 // summary frame lookup
371
372 JDAQSummaryFrame summary_frame;
373
374 if (!summaryRouter.hasSummaryFrame(super_frame->getModuleIdentifier())) {
375 summaryRouter.print(cout);
376 FATAL("[R>D>F] Summary frame not found, module = " << super_frame->getModuleID() << endl);
377 } else {
378 summary_frame = summaryRouter.getSummaryFrame(super_frame->getModuleID());
379 }
380
381 // PMT rates recovered from summaryslice data
382 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
383 rate_Hz[i] = summary_frame.getRate(i);
384 }
385
386 // frame status retrieved from summaryframe
387 status = summary_frame.getDAQFrameStatus();
388
389 }
390
391 //---------------------------
392 // processing of summary data
393 //---------------------------
394
395 // DAQ error (e.g. missing UDP trailer) -> discard frame
396 if (!status.testDAQStatus()) {
397 WARNING("[R>D>F] DAQ status not okay in timeslice #" << timeslice->getFrameIndex() << ", DOM=" << super_frame->getModuleID() << endl);
398 continue;
399 }
400
401 vector<bool> veto(NUMBER_OF_PMTS, false);
402 bool moduleLiveTime = 0;
403
404 // default status of module is alive if at least one PMT is active
405 for (vector<double>::const_iterator r = rate_Hz.begin(); !moduleLiveTime && r != rate_Hz.end(); ++r) {
406 if ((*r) > 0) {
407 moduleLiveTime = 1;
408 }
409 }
410
411 int frame_HRV_count = status.countHighRateVeto();
412 int frame_FAF_count = status.countFIFOStatus();
413 int frame_URV_count = 0; // user rate veto
414
415 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
416 if ((veto[i] = (veto[i] || !rateVeto_Hz(rate_Hz[i])))) {
417 frame_URV_count++;
418 }
419 }
420
421 count_HRV += frame_HRV_count;
422 count_FAF += frame_FAF_count;
423 count_URV += frame_URV_count;
424
425 int frame_VETO_count = frame_HRV_count + frame_FAF_count + frame_URV_count;
426
427
428 if (filterLevel >= 2 && frame_VETO_count > 0) {
429 moduleLiveTime = 0;
430 fill(veto.begin(), veto.end(), true);
431 } else if (filterLevel >= 1 && frame_HRV_count + frame_FAF_count > 0) {
432 for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
433 veto[i] = (veto[i] || (status.testFIFOStatus(i) || status.testHighRateVeto(i)));
434 }
435 }
436
437 // recover DOM in case the dead channel corresponds to the muted channel
438
439 if (muteChannel != -1) {
440
441 // veto conditions are not exclusive
442 int mute_VETO_count =
443 status.testHighRateVeto(muteChannel) +
444 status.testFIFOStatus(muteChannel ) +
445 !rateVeto_Hz(rate_Hz[muteChannel] );
446
447 if (mute_VETO_count == frame_VETO_count) {
448 moduleLiveTime = 1;
449 fill(veto.begin(), veto.end(), false);
450 }
451
452 veto[muteChannel] = true;
453
454 }
455
456 // get address and module for current superframe DOM;
457 const JModuleAddress& address = router.getAddress(super_frame->getModuleID());
458 const JModule& module = detector.getModule(address);
459 const TString moduleLabel = getLabel(module);
460
461 // acknowledge DOM is live in this timeslice for analysis purpose;
462 if (moduleLiveTime) {
463 h_livetime->Fill(moduleLabel, getFrameTime() * 1e-9 * moduleLiveTime);
464 }
465
466 //--------------------------------------------------------------------
467 // superframe processing: hit joining, calibration, filter, processing
468 //--------------------------------------------------------------------
469
470 JSuperFrame2D_t& buffer = JSuperFrame2D_t::demultiplex(*super_frame, router.getModule(super_frame->getModuleID()));
471
472 if (!notJoin) {
473 for (JSuperFrame2D_t::iterator i = buffer.begin(); i != buffer.end(); ++i) {
474 i->join(match);
475 }
476 }
477
478 JSuperFrame1D_t& data = JSuperFrame1D_t::multiplex(buffer);
479
480 filteredData.clear();
481
482 for (vector<JHitR0>::const_iterator h = data.begin(); h != data.end(); ++h) {
483
484 const int pmt = h->getPMT();
485
486 // veto and calibration;
487 if (!veto[pmt] &&
488 rateVeto_Hz(rate_Hz[pmt]) &&
489 totVeto_ns(h->getToT()) ) {
490 filteredData.push_back(*h);
491 }
492 }
493
494 // processing signal;
495
496 if (filteredData.size() > 1) {
497
498 TH1D* h_hit = H[moduleLabel];
499 TH1D* h_pmt = P[moduleLabel];
500
501 TH2D* h2_hit = HT[moduleLabel];
502 TH2D* h2_co = CO[moduleLabel];
503
504 sort(filteredData.begin(), filteredData.end());
505
506 //-----------------------------------------------------------
507 // superframe processing: loop on calibrated hits
508 //-----------------------------------------------------------
509 for (vector<JHitR0>::const_iterator p = filteredData.begin(); p != filteredData.end(); ) {
510
511 vector<JHitR0>::const_iterator q = p;
512
513 // coincidence building using a fixed time window
514 while (++q != filteredData.end() && q->getT() - p->getT() <= TMax_ns ) {}
515
516 // multiplicity estimation
517 int hit_multiplicity = distance(p,q);
518 set<int> pmthit;
519 for (vector<JHitR0>::const_iterator __p = p ; __p != q; ++__p ) { pmthit.insert(__p->getPMT()) ; }
520 int pmt_multiplicity = pmthit.size() ;
521
522 // histogram filling
523 h_pmt->Fill(pmt_multiplicity);
524 h_hit->Fill(hit_multiplicity);
525
526 if (twoDim && hit_multiplicity > 1) {
527 const double W = factorial(hit_multiplicity, 2);
528 for (vector<JHitR0>::const_iterator __p = p; __p != q; ++__p) {
529 for (vector<JHitR0>::const_iterator __q = __p; ++__q != q; ) {
530 double dt = JCombinatorics::getSign(__p->getPMT(), __q->getPMT()) * (__q->getT() - __p->getT());
531 h2_hit->Fill(hit_multiplicity, dt, 1.0/W);
532 }
533 }
534 }
535
536 if (monitorOccupancy) {
537 for (set<int>::const_iterator r = pmthit.begin(); r != pmthit.end(); ++r) {
538 JPMTPhysicalAddress pmtAddress = moduleAddressMap.getAddressTranslator(*r);
539 TString pmtLabel(pmtAddress.toString());
540 h2_co->Fill(pmt_multiplicity, pmtLabel, 1.0/pmt_multiplicity);
541 }
542 }
543
544 // benchmarking
545 if (hit_multiplicity >= sn_mul_min) {
546 time_distr->Fill(p->getT());
547 }
548
549 p = q;
550
551 }
552
553 }
554
555 }
556
557 }
558
559 NOTICE("[R] " << counter << " timeslices processed." << endl);
560
561 if (isMupage) {
562 h_livetime->Scale(liveTimeMC / (getFrameTime() * 1e-9 * counter * prescale));
563 }
564
565 // nominal live time;
566 double nominalLiveTime = (isMupage ? liveTimeMC : counter * getFrameTime() * 1e-9) / prescale;
567 h_livetime->Fill("Nominal", nominalLiveTime);
568
569 out.cd();
570
571 int nLiveDOMs = H.size();
572
573 h_livetime->Write();
574
575 // writing hit multiplicities is redundant given new options
576 // H.Write(out);
577
578 P.Write(out);
579
580 if (twoDim) {
581 HT.Write(out);
582 }
583
584 if (monitorOccupancy) {
585 CO.Write(out);
586 }
587
588 time_distr->Write();
589
590 double rate_HRV = count_HRV / (1.0 * counter);
591 double rate_FAF = count_FAF / (1.0 * counter);
592 double rate_URV = count_URV / (1.0 * counter);
593 double biolumIndex = (rate_HRV + rate_FAF) / (nLiveDOMs * NUMBER_OF_PMTS);
594 NOTICE("[R] Average HRV count per timeslice: " << rate_HRV << endl);
595 NOTICE("[R] Average FIFO-almost-full count per timeslice: " << rate_FAF << endl);
596 NOTICE("[R] Average user rate veto count per timeslice: " << rate_URV << endl);
597 NOTICE("[R] " << H.size() << " DOMs were active in the run." << endl);
598 NOTICE("[R] Bioluminescence proxy index: " << biolumIndex << endl);
599
600 putObject(out, JMeta(argc,argv));
601
602 out.Close();
603
604 NOTICE(endl);
605
606 return (counter ? 0 : 1);
607
608}
outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37
STATUS
#define STATUS(A)
Definition JMessage.hh:63
NOTICE
#define NOTICE(A)
Definition JMessage.hh:64
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:72
WARNING
#define WARNING(A)
Definition JMessage.hh:65
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
distance
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
Definition PhysicsEvent.hh:434
JAANET::JHead
Monte Carlo run header.
Definition JHead.hh:1236
JAANET::JHead::livetime
JAANET::livetime livetime
Definition JHead.hh:1604
JDETECTOR::JDetectorAddressMap
Lookup table for PMT addresses in detector.
Definition JDetectorAddressMap.hh:26
JDETECTOR::JDetectorAddressMap::get
const JModuleAddressMap & get(const int id) const
Get module address map.
Definition JDetectorAddressMap.hh:69
JDETECTOR::JDetectorAddressMap::getDefaultModuleAddressMap
virtual const JModuleAddressMap & getDefaultModuleAddressMap() const =0
Get default module address map.
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JModuleAddressMap
Lookup table for PMT addresses in optical module.
Definition JModuleAddressMap.hh:88
JDETECTOR::JModuleAddressMap::getAddressTranslator
const JPMTAddressTranslator & getAddressTranslator(const int tdc) const
Get PMT address translator.
Definition JModuleAddressMap.hh:257
JDETECTOR::JModuleAddress
Address of module in detector data structure.
Definition JModuleAddress.hh:32
JDETECTOR::JModuleRouter
Router for direct addressing of module data in detector data structure.
Definition JModuleRouter.hh:38
JDETECTOR::JModule
Data structure for a composite optical module.
Definition JModule.hh:75
JDETECTOR::JPMTPhysicalAddress
Data structure for PMT physical address.
Definition JPMTPhysicalAddress.hh:31
JDETECTOR::JPMTPhysicalAddress::toString
std::string toString() const
Convert PMT physical address to string.
Definition JPMTPhysicalAddress.hh:141
JLANG::JException
General exception.
Definition JException.hh:24
JLANG::JObjectMultiplexer
Auxiliary class for multiplexing object iterators.
Definition JObjectMultiplexer.hh:35
JPARSER::JParser
Utility class to parse command line options.
Definition JParser.hh:1698
JROOT::JManager
Auxiliary class to manage set of compatible ROOT objects (e.g. histograms) using unique keys.
Definition JManager.hh:47
JSUPPORT::JMultipleFileScanner
General purpose class for object reading from a list of file names.
Definition JMultipleFileScanner.hh:306
JSUPPORT::JSummaryRouter
Router for fast addressing of summary data in KM3NETDAQ::JDAQSummaryslice data structure as a functio...
Definition JSupport/JSummaryRouter.hh:41
JSUPPORT::JSummaryRouter::hasSummaryFrame
bool hasSummaryFrame(const JDAQModuleIdentifier &module) const
Has summary frame.
Definition JSupport/JSummaryRouter.hh:204
JSUPPORT::JSummaryRouter::getSummaryFrame
const JDAQSummaryFrame & getSummaryFrame() const
Get default summary frame.
Definition JSupport/JSummaryRouter.hh:72
JSUPPORT::JSummaryRouter::update
void update(const JDAQSummaryslice *ps)
Update router.
Definition JSupport/JSummaryRouter.hh:141
JSUPPORT::JTreeScanner
Template definition for direct access of elements in ROOT TChain.
Definition JTreeScanner.hh:698
JTOOLS::JCombinatorics::getSign
static int getSign(const int first, const int second)
Sign of pair of indices.
Definition JCombinatorics.hh:215
JTOOLS::JRange
Range of values.
Definition JRange.hh:42
JTOOLS::JRange::is_valid
bool is_valid() const
Check validity of range.
Definition JRange.hh:311
JTRIGGER::JMatchL0
L0 match criterion.
Definition JMatchL0.hh:29
JTRIGGER::JSuperFrame1D
1-dimensional frame with time calibrated data from one optical module.
Definition JSuperFrame1D.hh:40
JTRIGGER::JSuperFrame2D
2-dimensional frame with time calibrated data from one optical module.
Definition JSuperFrame2D.hh:46
KM3NETDAQ::JDAQChronometer::getFrameIndex
int getFrameIndex() const
Get frame index.
Definition JDAQChronometer.hh:132
KM3NETDAQ::JDAQFrameStatus::testFIFOStatus
bool testFIFOStatus() const
Test FIFO status.
Definition JDAQFrameStatus.hh:245
KM3NETDAQ::JDAQFrameStatus::countHighRateVeto
int countHighRateVeto() const
Count high-rate veto status.
Definition JDAQFrameStatus.hh:224
KM3NETDAQ::JDAQFrameStatus::getDAQFrameStatus
const JDAQFrameStatus & getDAQFrameStatus() const
Get DAQ frame status.
Definition JDAQFrameStatus.hh:80
KM3NETDAQ::JDAQFrameStatus::countFIFOStatus
int countFIFOStatus() const
Count FIFO status.
Definition JDAQFrameStatus.hh:268
KM3NETDAQ::JDAQFrameStatus::testHighRateVeto
bool testHighRateVeto() const
Test high-rate veto status.
Definition JDAQFrameStatus.hh:201
KM3NETDAQ::JDAQFrameStatus::testDAQStatus
bool testDAQStatus() const
Test DAQ status of packets.
Definition JDAQFrameStatus.hh:135
KM3NETDAQ::JDAQHit
Hit data structure.
Definition JDAQHit.hh:35
KM3NETDAQ::JDAQSummaryFrame
Data storage class for rate measurements of all PMTs in one module.
Definition JDAQSummaryFrame.hh:335
KM3NETDAQ::JDAQSummaryFrame::getRate
double getRate(const int tdc, const double factor=1.0) const
Get count rate.
Definition JDAQSummaryFrame.hh:468
GAUSS_LEGENDRE::xmax
const double xmax
Definition JQuadrature.cc:24
GAUSS_LEGENDRE::xmin
const double xmin
Definition JQuadrature.cc:23
JDETECTOR::getLabel
std::string getLabel(const JLocation &location)
Get module label for monitoring and other applications.
Definition JLocation.hh:247
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition JDetectorToolkit.hh:508
JDETECTOR::getDetectorAddressMap
JDetectorAddressMap & getDetectorAddressMap()
Get detector address map.
Definition JDetectorSupportkit.hh:542
JGIZMO::setAxisLabels
void setAxisLabels(TAxis *axis, const JModuleAddressMap &memo)
Set axis with PMT address labels.
Definition JGizmoToolkit.hh:1525
JMATH::factorial
long long int factorial(const long long int n)
Determine factorial.
Definition JMathToolkit.hh:43
JPHYSICS::H
static const double H
Planck constant [eV s].
Definition JPhysics/JConstants.hh:29
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:43
JROOT::putObject
bool putObject(TDirectory &dir, const TObject &object)
Write object to ROOT directory.
Definition JRootFileWriter.hh:41
JROOT::getROOTClassSelection
std::set< JROOTClassSelector > getROOTClassSelection(const bool option=false)
Get ROOT class selection.
Definition JROOTClassSelector.hh:318
JROOT::counter_type
Long64_t counter_type
Type definition for counter.
Definition JROOT/JCounter.hh:24
JSUPPORT::getTriggerParameters
JTriggerParameters getTriggerParameters(const JMultipleFileScanner_t &file_list)
Get trigger parameters.
Definition JTriggerParametersSupportkit.hh:148
JSUPPORT::getHeader
Head getHeader(const JMultipleFileScanner_t &file_list)
Get Monte Carlo header.
Definition JMonteCarloFileSupportkit.hh:547
JTRIGGER::getNumberOfModules
static const JModuleCounter getNumberOfModules
Function object to count unique modules.
Definition JModuleCounter.hh:67
KM3NETDAQ
KM3NeT DAQ data structures and auxiliaries.
Definition DataQueue.cc:39
KM3NETDAQ::getFrameTime
double getFrameTime()
Get frame time duration.
Definition JDAQClock.hh:162
KM3NETDAQ::NUMBER_OF_PMTS
static const int NUMBER_OF_PMTS
Total number of PMTs in module.
Definition JDAQ.hh:26
JAANET::detector
Detector file.
Definition JHead.hh:227
JAANET::livetime::numberOfSeconds
double numberOfSeconds
Live time [s].
Definition JHead.hh:896
JAANET::livetime::errorOfSeconds
double errorOfSeconds
Uncertainty on live time [s].
Definition JHead.hh:897
JROOT::JROOTClassSelector
Auxiliary class to select ROOT class based on class name.
Definition JROOTClassSelector.hh:38
JSUPPORT::JLimit
Auxiliary class for defining the range of iterations of objects.
Definition JLimit.hh:45
JSUPPORT::JLimit::max
static counter_type max()
Get maximum counter value.
Definition JLimit.hh:128
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72
Generated by doxygen 1.12.0