Jpp 19.3.0-rc.3
the software that should make you happy
Loading...
Searching...
No Matches
Classes | Namespaces | Functions
JEvD.cc File Reference

Program to display hit probabilities. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <vector>
#include <algorithm>
#include <memory>
#include "TROOT.h"
#include "TApplication.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TH2D.h"
#include "TArrow.h"
#include "TLatex.h"
#include "TMarker.h"
#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/Evt.hh"
#include "km3net-dataformat/online/JDAQ.hh"
#include "km3net-dataformat/online/JDAQHeader.hh"
#include "km3net-dataformat/tools/time_converter.hh"
#include "JROOT/JStyle.hh"
#include "JROOT/JCanvas.hh"
#include "JROOT/JRootToolkit.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
#include "JTrigger/JHitL0.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JSingleFileScanner.hh"
#include "JSupport/JSummaryFileRouter.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JAAnet/JAAnetToolkit.hh"
#include "JPhysics/JPDF_t.hh"
#include "JFit/JLine1Z.hh"
#include "JFit/JModel.hh"
#include "JReconstruction/JHitW0.hh"
#include "JReconstruction/JMuonGandalfParameters_t.hh"
#include "JLang/JPredicate.hh"
#include "JLang/JComparator.hh"
#include "JMath/JMathToolkit.hh"
#include "JSystem/JKeypress.hh"
#include "JSystem/JProcess.hh"
#include "Jeep/JFunctionAdaptor.hh"
#include "Jeep/JProperties.hh"
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Classes

struct  JAANET::JEventSelector
 Event selector. More...
 

Namespaces

namespace  JAANET
 Extensions to Evt data format.
 

Functions

bool JAANET::hasW (const Trk &trk, const int i)
 Check availability of value.
 
double JAANET::getW (const Trk &trk, const int i)
 Get associated value.
 
double JAANET::getW (const Trk &trk, const int i, const double value)
 Get associated value.
 
void JAANET::setW (Trk &trk, const int i, const double value)
 Set associated value.
 
int main (int argc, char **argv)
 

Detailed Description

Program to display hit probabilities.

Author
mdejong

Definition in file JAAnet/JEvD.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 194 of file JAAnet/JEvD.cc.

195{
196 using namespace std;
197 using namespace JPP;
198 using namespace KM3NETDAQ;
199
200 JSingleFileScanner<Evt> inputFile;
201 JLimit_t& numberOfEvents = inputFile.getLimit();
202 string pdfFile;
203 string outputFile;
204 JMuonGandalfParameters_t parameters;
205 int application;
206 JEventSelector event_selector;
207 JCanvas canvas;
208 bool batch;
209 struct : JStyle::JParameters {
210 double arrowSize = 0.003;
211 string arrowType = "|->";
212 double arrowScale = 250.0;
213 Width_t lineWidth = 2;
214 Style_t lineStyle = 1;
215 int nbinsX = 50;
216 int nbinsY = 250;
217 double T_ns = 0.0;
218 } graphics;
219 string option;
220 int debug;
221
222
223 try {
224
225 JProperties properties = graphics.getProperties();
226
227 properties.insert(gmake_property(graphics.arrowSize));
228 properties.insert(gmake_property(graphics.arrowType));
229 properties.insert(gmake_property(graphics.arrowScale));
230 properties.insert(gmake_property(graphics.lineWidth));
231 properties.insert(gmake_property(graphics.lineStyle));
232 properties.insert(gmake_property(graphics.T_ns));
233
234 parameters.numberOfPrefits = 1;
235
236 JParser<> zap("Program to display hit probabilities.");
237
238 zap['w'] = make_field(canvas, "size of canvas <nx>x<ny> [pixels]") = JCanvas(1200, 600);
239 zap['f'] = make_field(inputFile, "input file (output of JXXXMuonReconstruction.sh)");
240 zap['n'] = make_field(numberOfEvents) = JLimit::max();
241 zap['P'] = make_field(pdfFile);
242 zap['o'] = make_field(outputFile, "graphics output file name") = MAKE_STRING("display_" << WILDCARD << ".gif");
243 zap['@'] = make_field(parameters) = JPARSER::initialised();
244 zap['A'] = make_field(application) = JMUONENERGY, JMUONSTART, JMUONGANDALF, JMUONSIMPLEX, JMUONPREFIT;
245 zap['L'] = make_field(event_selector) = JPARSER::initialised();
246 zap['%'] = make_field(properties) = JPARSER::initialised();
247 zap['O'] = make_field(option, "draw option") = arrow_t, histogram_t;
248 zap['B'] = make_field(batch, "batch processing");
249 zap['d'] = make_field(debug) = 1;
250
251 zap(argc, argv);
252 }
253 catch(const exception& error) {
254 FATAL(error.what() << endl);
255 }
256
257 if (batch && outputFile == "") {
258 FATAL("Missing output file name " << outputFile << " in batch mode." << endl);
259 }
260
261 if (!batch && outputFile == "") {
262 outputFile = MAKE_STRING(WILDCARD << ".gif");
263 }
264
265 if (outputFile.find(WILDCARD) == string::npos) {
266 FATAL("Output file name " << outputFile << " has no wild card '" << WILDCARD << "'" << endl);
267 }
268
269 JHit::setSlewing(false);
270
271 JSummaryFileRouter summary(inputFile);
272
273 const JMuonPDF_t pdf(pdfFile, parameters.TTS_ns);
274
275 const JTimeRange T_ns(parameters.TMin_ns, parameters.TMax_ns);
276
277 typedef vector<JHitL0> JDataL0_t;
278 typedef vector<JHitW0> JDataW0_t;
279
280
281 // ROOT
282
283 gROOT->SetBatch(batch);
284
285 TApplication* tp = new TApplication("user", NULL, NULL);
286 TCanvas* cv = new TCanvas("display", "", canvas.x, canvas.y);
287
288 unique_ptr<TStyle> gStyle(new JStyle("gplot", cv->GetWw(), cv->GetWh(), graphics));
289
290 gROOT->SetStyle("gplot");
291 gROOT->ForceStyle();
292
293 const size_t NUMBER_OF_PADS = 3;
294
295 cv->SetFillStyle(4000);
296 cv->SetFillColor(kWhite);
297
298 TPad* p1 = new TPad("p1", NULL, 0.0, 0.00, 1.0, 0.95);
299 TPad* p2 = new TPad("p2", NULL, 0.0, 0.95, 1.0, 1.00);
300
301 p1->Divide(NUMBER_OF_PADS, 1);
302
303 p1->Draw();
304 p2->Draw();
305
306 const double Dmax = 1000.0;
307 const double Rmin = 0.0;
308 const double Rmax = min(parameters.roadWidth_m, 0.4 * Dmax);
309 const double Tmin = min(parameters.TMin_ns, -10.0);
310 const double Tmax = max(parameters.TMax_ns, +100.0);
311 const double Amin = 0.002 * (Tmax - Tmin); // minimal arrow length [ns]
312 const double Amax = 0.8 * (Tmax - Tmin); // maximal arrow length [ns]
313 const double ymin = Tmin - (option == arrow_t ? 0.2 * Amax : 0.0);
314 const double ymax = Tmax + (option == arrow_t ? 0.5 * Amax : 0.0);
315
316 const string Xlabel[NUMBER_OF_PADS] = { "R [m]", "#phi [rad]", "z [m]" };
317 const double Xmin [NUMBER_OF_PADS] = { Rmin, -PI, -0.3 * Dmax };
318 const double Xmax [NUMBER_OF_PADS] = { Rmax, +PI, +0.3 * Dmax };
319
320 double Xs[NUMBER_OF_PADS];
321
322 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
323 Xs[i] = 0.003 * (Xmax[i] - Xmin[i]) * (0.5 * NUMBER_OF_PMTS); // x-offset arrow as function of PMT number
324 }
325
326 TH2D H2[NUMBER_OF_PADS];
327 TGraph G2[NUMBER_OF_PADS];
328
329 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
330
331 H2[i] = TH2D(MAKE_CSTRING("h" << i), NULL, graphics.nbinsX, Xmin[i] - Xs[i], Xmax[i] + Xs[i], graphics.nbinsY, ymin, ymax);
332
333 H2[i].GetXaxis()->SetTitle(Xlabel[i].c_str());
334 H2[i].GetYaxis()->SetTitle("#Deltat [ns]");
335
336 H2[i].GetXaxis()->CenterTitle(true);
337 H2[i].GetYaxis()->CenterTitle(true);
338
339 H2[i].SetStats(kFALSE);
340
341 G2[i].Set(2);
342
343 G2[i].SetPoint(0, H2[i].GetXaxis()->GetXmin(), 0.0);
344 G2[i].SetPoint(1, H2[i].GetXaxis()->GetXmax(), 0.0);
345
346 p1->cd(i+1);
347
348 H2[i].Draw("AXIS");
349 G2[i].Draw("SAME");
350 }
351
352
353 while (inputFile.hasNext()) {
354
355 cout << "event: " << setw(8) << inputFile.getCounter() << endl;
356
357 const Evt* evt = inputFile.next();
358
359 if (has_reconstructed_track<JPP_RECONSTRUCTION_TYPE>(*evt, rec_stages_range(application))) {
360
361 Trk fit = get_best_reconstructed_track<JPP_RECONSTRUCTION_TYPE>(*evt, rec_stages_range(application));
362
363 if (!event_selector(fit, *evt)) {
364 continue;
365 }
366
367 JDataL0_t dataL0;
368
369 for (const Hit& hit : evt->hits) {
370 dataL0.push_back(JHitL0(JDAQPMTIdentifier(hit.dom_id, hit.channel_id),
371 JAxis3D(getPosition(hit.pos), getDirection(hit.dir)),
372 JHit(hit.t, hit.tot)));
373 }
374
375 summary.update(JDAQChronometer(evt->det_id,
376 evt->run_id,
377 evt->frame_index,
378 JDAQUTCExtended(evt->t.GetSec(), evt->t.GetNanoSec() / 16)));
379
380 const time_converter converter = time_converter(*evt);
381
382 Trk muon; // Monte Carlo true muon
383
384 if (has_muon(*evt)) {
385
386 for (const auto& trk : evt->mc_trks) {
387 if (is_muon(trk)) {
388 if (trk.E > muon.E) {
389
390 muon = trk;
391 muon.t += converter.putTime();
392
393 setW(muon, JSTART_LENGTH_METRES, fabs(muon.len));
394 }
395 }
396 }
397 }
398
399
400 bool monte_carlo = false; // show Monte Carlo true muon
401
402 for (bool next = false; !next; ) {
403
404 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
405 H2[i].Reset();
406 }
407
408 Trk trk;
409
410 if (!monte_carlo)
411 trk = fit;
412 else
413 trk = muon;
414
415 JRotation3D R (getDirection(trk));
416 JLine1Z tz(getPosition (trk).rotate(R), trk.t);
417 JRange<double> Z_m;
418 /*
419 if (hasW(trk, JSTART_LENGTH_METRES)) {
420 Z_m = JRange<double>(0.0, getW(fit,JSTART_LENGTH_METRES)) + JRange<double>(parameters.ZMin_m, parameters.ZMax_m);
421 }
422 */
423 const JModel<JLine1Z> match(tz, parameters.roadWidth_m, T_ns, Z_m);
424
425 // hit selection based on fit result
426
427 JDataW0_t data;
428
429 for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
430
431 JHitW0 hit(*i, summary.getRate(i->getPMTIdentifier(), parameters.R_Hz));
432
433 hit.rotate(R);
434
435 if (match(hit)) {
436 data.push_back(hit);
437 }
438 }
439
440 // select first hit in PMT
441
442 sort(data.begin(), data.end(), JHitW0::compare);
443
444 JDataW0_t::iterator __end = unique(data.begin(), data.end(), equal_to<JDAQPMTIdentifier>());
445
446 double E_GeV = parameters.E_GeV;
447 /*
448 if (trk.E > 0.1) {
449 E_GeV = trk.E;
450 }
451 */
452
453 // move fit to geometrical center of hits
454
455 JRange<double> zs = JRange<double>::DEFAULT_RANGE();
456
457 for (JDataW0_t::iterator hit = data.begin(); hit != __end; ++hit) {
458
459 const double x = hit->getX() - tz.getX();
460 const double y = hit->getY() - tz.getY();
461 const double z = hit->getZ();
462 const double R = sqrt(x*x + y*y);
463
464 zs.include(z - R/getTanThetaC());
465 }
466
467 const double z0 = tz.getZ();
468 const double z1 = 0.5 * (zs.getLowerLimit() + zs.getUpperLimit());
469
470 tz.setZ(z1, getSpeedOfLight());
471
472
473 // graphics
474
475 ostringstream os;
476 vector<TArrow> arrow [NUMBER_OF_PADS];
477 vector<TMarker> marker[NUMBER_OF_PADS];
478
479 if (hasW(trk, JSTART_LENGTH_METRES) && getW(trk, JSTART_LENGTH_METRES) > 0.0) {
480
481 marker[2].push_back(TMarker(z0 - tz.getZ(), 0.0, kFullCircle));
482 marker[2].push_back(TMarker(z0 - tz.getZ() + getW(trk, JSTART_LENGTH_METRES), 0.0, kFullCircle));
483
484 static_cast<TAttMarker&>(marker[2][0]) = TAttMarker(kRed, kFullCircle, 0.7);
485 static_cast<TAttMarker&>(marker[2][1]) = TAttMarker(kRed, kFullCircle, 0.7);
486 }
487
488 DEBUG("trk: "
489 << FIXED(7,2) << tz.getX() << ' '
490 << FIXED(7,2) << tz.getY() << ' '
491 << FIXED(7,2) << tz.getZ() << ' '
492 << FIXED(12,2) << tz.getT() << endl);
493
494 double chi2 = 0;
495
496 for (JDataW0_t::const_iterator hit = data.begin(); hit != __end; ++hit) {
497
498 const double x = hit->getX() - tz.getX();
499 const double y = hit->getY() - tz.getY();
500 const double z = hit->getZ() - tz.getZ();
501 const double R = sqrt(x*x + y*y);
502
503 const double t1 = tz.getT() + (z + R * getTanThetaC()) * getInverseSpeedOfLight();
504
505 JDirection3D dir(hit->getDX(), hit->getDY(), hit->getDZ()); // PMT orientation
506
507 dir.rotate(JRotation3Z(-atan2(y,x))); // rotate PMT axis to x-z plane
508
509 const double theta = dir.getTheta();
510 const double phi = fabs(dir.getPhi()); // rotational symmetry of Cherenkov cone
511
512 //const double E = gWater.getE(E_GeV, z); // correct for energy loss
513 const double E = E_GeV;
514 const double dt = T_ns.constrain(hit->getT() - t1);
515
516 JMuonPDF_t::result_type H1 = pdf.calculate(E, R, theta, phi, dt);
517 JMuonPDF_t::result_type H0(hit->getR() * 1e-9, 0.0, T_ns);
518
519 if (H1.V >= parameters.VMax_npe) {
520 H1 *= parameters.VMax_npe / H1.V;
521 }
522
523 H1 += H0; // signal + background
524
525 chi2 += H1.getChi2() - H0.getChi2();
526
527 DEBUG("hit: "
528 << setw(8) << hit->getModuleID() << '.' << FILL(2,'0') << (int) hit->getPMTAddress() << FILL() << ' '
529 << SCIENTIFIC(8,2) << E << ' '
530 << FIXED(7,2) << R << ' '
531 << FIXED(7,4) << theta << ' '
532 << FIXED(7,4) << phi << ' '
533 << FIXED(7,3) << dt << ' '
534 << FIXED(7,3) << H1.getChi2() << ' '
535 << FIXED(7,3) << H0.getChi2() << endl);
536
537 const double derivative = H1.getDerivativeOfChi2() - H0.getDerivativeOfChi2();
538
539 double size = derivative * graphics.arrowScale; // size of arrow
540
541 if (fabs(size) < Amin) {
542 size = (size > 0.0 ? +Amin : -Amin);
543 } else if (fabs(size) > Amax) {
544 size = (size > 0.0 ? +Amax : -Amax);
545 }
546
547 const double X[NUMBER_OF_PADS] = { R, atan2(y,x), z - R/getTanThetaC() };
548
549 const double xs = (double) (NUMBER_OF_PMTS - 2 * hit->getPMTAddress()) / (double) NUMBER_OF_PMTS;
550
551 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
552
553 TArrow a1(X[i] + xs*Xs[i], dt + graphics.T_ns, X[i] + xs*Xs[i], dt + graphics.T_ns + size, graphics.arrowSize, graphics.arrowType.c_str());
554
555 a1.SetLineWidth(graphics.lineWidth);
556 a1.SetLineStyle(graphics.lineStyle);
557
558 arrow[i].push_back(a1);
559
560 H2[i].Fill(X[i], dt + graphics.T_ns);
561 }
562 }
563
564 os << FILL(6,'0') << evt->run_id << ":" << evt->frame_index << "/" << evt->trigger_counter << FILL();
565 os << " Q = " << FIXED(4,0) << trk.lik
566 << '/' << FIXED(4,0) << -chi2;
567 os << " E = " << SCIENTIFIC(7,1) << trk.E << " [GeV]";
568 os << " cos(#theta) = " << FIXED(6,3) << trk.dir.z;
569
570 if (monte_carlo)
571 os << " Monte Carlo";
572 else if (is_muon(muon))
573 os << " #Delta#alpha = " << FIXED(6,2) << getAngle(getDirection(muon), getDirection(trk)) << " [deg]";
574
575
576 // draw
577
578 TLatex title(0.05, 0.5, os.str().c_str());
579
580 title.SetTextAlign(12);
581 title.SetTextFont(42);
582 title.SetTextSize(0.6);
583
584 p2->cd();
585
586 title.Draw();
587
588 for (int i = 0; i != NUMBER_OF_PADS; ++i) {
589
590 p1->cd(i+1);
591
592 if (option == arrow_t) {
593
594 for (auto& a1 : arrow[i]) {
595 a1.Draw();
596 }
597
598 for (auto& m1 : marker[i]) {
599 m1.Draw();
600 }
601 }
602
603 if (option == histogram_t) {
604 H2[i].Draw("SAME");
605 }
606 }
607
608 cv->Update();
609
610
611 // action
612
613 if (batch) {
614
615 cv->SaveAs(replace(outputFile, WILDCARD, MAKE_STRING(inputFile.getCounter())).c_str());
616
617 next = true;
618
619 } else {
620
621 static int count = 0;
622
623 if (count++ == 0) {
624 cout << endl << "Type '?' for possible options." << endl;
625 }
626
627 for (bool user = true; user; ) {
628
629 cout << "\n> " << flush;
630
631 switch (JKeypress(true).get()) {
632
633 case '?':
634 cout << endl;
635 cout << "possible options: " << endl;
636 cout << 'q' << " -> " << "exit application" << endl;
637 cout << 'u' << " -> " << "update canvas" << endl;
638 cout << 's' << " -> " << "save graphics to file" << endl;
639 cout << 'M' << " -> " << "Monte Carlo true muon information" << endl;
640 cout << 'F' << " -> " << "fit information" << endl;
641 if (event_selector.is_valid()) {
642 cout << 'L' << " -> " << "reload event selector" << endl;
643 }
644 cout << 'r' << " -> " << "rewind input file" << endl;
645 cout << 'R' << " -> " << "switch to ROOT mode (quit ROOT to continue)" << endl;
646 cout << 'p' << " -> " << "print event information" << endl;
647 cout << ' ' << " -> " << "next event (as well as any other key)" << endl;
648 break;
649
650 case 'q':
651 cout << endl;
652 return 0;
653
654 case 'u':
655 cv->Update();
656 break;
657
658 case 's':
659 cv->SaveAs(replace(outputFile, WILDCARD, MAKE_STRING(inputFile.getCounter())).c_str());
660 break;
661
662 case 'M':
663 if (is_muon(muon))
664 monte_carlo = true;
665 else
666 ERROR(endl << "No Monte Carlo muon available." << endl);
667 user = false;
668 break;
669
670 case 'F':
671 monte_carlo = false;
672 user = false;
673 break;
674
675 case 'L':
676 if (event_selector.is_valid()) {
677 execute(MAKE_STRING("make -f " << getPath(argv[0]) << "/JMakeEventSelector libs"), 3);
678 event_selector.reload();
679 }
680 break;
681
682 case 'R':
683 tp->Run(kTRUE);
684 break;
685
686 case 'p':
687 cout << endl;
688 evt->print(cout);
689 if (debug >= debug_t) {
690 for (const auto& trk : evt->mc_trks) {
691 cout << "MC "; trk.print(cout); cout << endl;
692 }
693 for (const auto& trk : evt->trks) {
694 cout << "fit "; trk.print(cout); cout << endl;
695 }
696 for (const auto& hit : evt->hits) {
697 cout << "hit "; hit.print(cout); cout << endl;
698 }
699 }
700 break;
701
702 case 'r':
703 inputFile.rewind();
704
705 default:
706 next = true;
707 user = false;
708 break;
709 }
710 }
711 }
712 }
713 }
714 }
715 cout << endl;
716}
outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37
p1
TPaveText * p1
Definition JDrawModule3D.cc:36
DEBUG
#define DEBUG(A)
Message macros.
Definition JMessage.hh:62
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:72
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition JPrint.hh:72
MAKE_STRING
#define MAKE_STRING(A)
Make string.
Definition JPrint.hh:63
gmake_property
#define gmake_property(A)
macros to convert (template) parameter to JPropertiesElement object
Definition JProperties.hh:1293
JEEP::JProperties
Utility class to parse parameter values.
Definition JProperties.hh:501
JFIT::JLine1Z
Data structure for fit of straight line paralel to z-axis.
Definition JLine1Z.hh:29
JGEOMETRY3D::JAxis3D
Axis object.
Definition JAxis3D.hh:41
JGEOMETRY3D::JDirection3D
Data structure for direction in three dimensions.
Definition JDirection3D.hh:35
JGEOMETRY3D::JRotation3Z
Rotation around Z-axis.
Definition JRotation3D.hh:87
JPARSER::JParser
Utility class to parse command line options.
Definition JParser.hh:1698
JRECONSTRUCTION::JHitW0
Auxiliary class for a hit with background rate value.
Definition JHitW0.hh:23
JROOT::JCanvas
Data structure for size of TCanvas.
Definition JCanvas.hh:26
JROOT::JCanvas::y
int y
number of pixels in Y
Definition JCanvas.hh:99
JROOT::JCanvas::x
int x
number of pixels in X
Definition JCanvas.hh:98
JROOT::JStyle
Wrapper class around ROOT TStyle.
Definition JStyle.hh:24
JSUPPORT::JSingleFileScanner
Object reading from a list of files.
Definition JSingleFileScanner.hh:172
JSUPPORT::JSingleFileScanner::rewind
virtual void rewind() override
Rewind.
Definition JSingleFileScanner.hh:240
JSUPPORT::JSingleFileScanner::hasNext
virtual bool hasNext() override
Check availability of next element.
Definition JSingleFileScanner.hh:257
JSUPPORT::JSingleFileScanner::getCounter
counter_type getCounter() const
Get counter.
Definition JSingleFileScanner.hh:231
JSUPPORT::JSingleFileScanner::next
virtual const pointer_type & next() override
Get next element.
Definition JSingleFileScanner.hh:303
JSUPPORT::JSummaryFileRouter
File router for fast addressing of summary data.
Definition JSummaryFileRouter.hh:40
JSYSTEM::JKeypress
Enable unbuffered terminal input.
Definition JKeypress.hh:32
JTOOLS::JRange::DEFAULT_RANGE
static JRange< T, JComparator_t > DEFAULT_RANGE()
Default range.
Definition JRange.hh:555
JTOOLS::JRange::include
range_type & include(argument_type x)
Include given value to range.
Definition JRange.hh:397
JTOOLS::JRange::getLowerLimit
T getLowerLimit() const
Get lower limit.
Definition JRange.hh:202
JTOOLS::JRange::getUpperLimit
T getUpperLimit() const
Get upper limit.
Definition JRange.hh:213
JTRIGGER::JHitL0
Data structure for L0 hit.
Definition JHitL0.hh:31
JTRIGGER::JHit::setSlewing
static void setSlewing(const bool slewing)
Set slewing option.
Definition JTrigger/JHit.hh:121
KM3NETDAQ::JDAQUTCExtended
Data structure for UTC time.
Definition JDAQUTCExtended.hh:23
time_converter
Auxiliary class to convert DAQ hit time to/from Monte Carlo hit time.
Definition time_converter.hh:20
time_converter::putTime
double putTime() const
Get Monte Carlo time minus DAQ/trigger time.
Definition time_converter.hh:72
JMUONGANDALF
static const int JMUONGANDALF
Definition definitions/reconstruction.hh:16
JMUONPREFIT
static const int JMUONPREFIT
Definition definitions/reconstruction.hh:14
JMUONENERGY
static const int JMUONENERGY
Definition definitions/reconstruction.hh:17
JMUONSIMPLEX
static const int JMUONSIMPLEX
Definition definitions/reconstruction.hh:15
JMUONSTART
static const int JMUONSTART
Definition definitions/reconstruction.hh:18
JSTART_LENGTH_METRES
static const int JSTART_LENGTH_METRES
distance between projected positions on the track of optical modules for which the response does not ...
Definition fitparameters.hh:26
JAANET::getDirection
JDirection3D getDirection(const Vec &dir)
Get direction.
Definition JAAnetToolkit.hh:214
JAANET::hasW
bool hasW(const Trk &trk, const int i)
Check availability of value.
Definition JAAnet/JEvD.cc:102
JAANET::has_muon
bool has_muon(const Evt &evt)
Test whether given event has a muon.
Definition JAAnetToolkit.hh:765
JAANET::setW
void setW(Trk &trk, const int i, const double value)
Set associated value.
Definition JAAnet/JEvD.cc:145
JAANET::getW
double getW(const Trk &track, const size_t index, const double value)
Get track information.
Definition JAAnetToolkit.hh:288
JAANET::getPosition
JPosition3D getPosition(const Vec &pos)
Get position.
Definition JAAnetToolkit.hh:178
JAANET::is_muon
bool is_muon(const Trk &track)
Test whether given track is a (anti-)muon.
Definition JAAnetToolkit.hh:338
JEEP::debug_t
@ debug_t
debug
Definition JMessage.hh:29
JEEP::getPath
std::string getPath(const std::string &file_name)
Get path, i.e. part before last JEEP::PATHNAME_SEPARATOR if any.
Definition JeepToolkit.hh:148
JGEOMETRY3D::getAngle
double getAngle(const JQuaternion3D &first, const JQuaternion3D &second)
Get space angle between quanternions.
Definition JQuaternion3D.hh:1165
JLANG::replace
std::string replace(const std::string &input, const std::string &target, const std::string &replacement)
Replace tokens in string.
Definition JLangToolkit.hh:124
JMATH::PI
static const double PI
Mathematical constants.
Definition JMath/JConstants.hh:20
JPHYSICS::getInverseSpeedOfLight
const double getInverseSpeedOfLight()
Get inverse speed of light.
Definition JPhysics/JConstants.hh:123
JPHYSICS::getTanThetaC
double getTanThetaC()
Get average tangent of Cherenkov angle of water corresponding to group velocity.
Definition JPhysics/JConstants.hh:156
JPHYSICS::getSpeedOfLight
const double getSpeedOfLight()
Get speed of light.
Definition JPhysics/JConstants.hh:112
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:43
JRECONSTRUCTION::get_best_reconstructed_track
const JFit & get_best_reconstructed_track(const JEvt &evt, JTrackSelector_t selector, JQualitySorter_t comparator)
Get best reconstructed track.
Definition JReconstruction/JEvtToolkit.hh:502
JRECONSTRUCTION::has_reconstructed_track
bool has_reconstructed_track(const JEvt &evt, JTrackSelector_t selector)
Test whether given event has a track according selection.
Definition JReconstruction/JEvtToolkit.hh:461
KM3NETDAQ
KM3NeT DAQ data structures and auxiliaries.
Definition DataQueue.cc:39
KM3NETDAQ::WILDCARD
static const char WILDCARD
Definition JDAQTags.hh:56
Evt
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Definition Evt.hh:21
Evt::frame_index
int frame_index
from the raw data
Definition Evt.hh:29
Evt::run_id
int run_id
DAQ run identifier.
Definition Evt.hh:26
Evt::hits
std::vector< Hit > hits
list of hits
Definition Evt.hh:38
Evt::mc_trks
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
Definition Evt.hh:49
Evt::det_id
int det_id
detector identifier from DAQ
Definition Evt.hh:23
Evt::trigger_counter
ULong64_t trigger_counter
trigger counter
Definition Evt.hh:31
Evt::print
void print(std::ostream &out=std::cout) const
Print event.
Definition Evt.hh:74
Evt::trks
std::vector< Trk > trks
list of reconstructed tracks (can be several because of prefits,showers, etc).
Definition Evt.hh:39
Evt::t
TTimeStamp t
UTC time of the timeslice, or the event_time for MC. (default: 01 Jan 1970 00:00:00)
Definition Evt.hh:33
FILL
Auxiliary data structure for sequence of same character.
Definition JManip.hh:330
FIXED
Auxiliary data structure for floating point format specification.
Definition JManip.hh:448
Hit
Definition Hit.hh:10
Hit::dom_id
int dom_id
module identifier from the data (unique in the detector).
Definition Hit.hh:14
Hit::pos
Vec pos
hit position
Definition Hit.hh:25
Hit::print
void print(std::ostream &out=std::cout) const
Print hit.
Definition Hit.hh:60
Hit::dir
Vec dir
hit direction; i.e. direction of the PMT
Definition Hit.hh:26
Hit::channel_id
unsigned int channel_id
PMT channel id {0,1, .., 30} local to moduke.
Definition Hit.hh:15
Hit::tot
unsigned int tot
tot value as stored in raw data (int for pyroot)
Definition Hit.hh:17
Hit::t
double t
hit time (from tdc+calibration or MC truth)
Definition Hit.hh:23
JACOUSTICS::JHit
Acoustics hit.
Definition JAcoustics/JHit.hh:35
JACOUSTICS::JModel
Model for fit to acoustics data.
Definition JAcoustics/JModel.hh:480
JEEP::JFunctionAdaptorHelper::is_valid
bool is_valid() const
Check validity of function.
Definition JFunctionAdaptor.hh:56
JEEP::JFunctionAdaptorHelper::reload
void reload()
Reload function from shared library.
Definition JFunctionAdaptor.hh:180
JMuonPDF_t
Auxiliary data structure for muon PDF.
Definition JPDF_t.hh:135
JMuonPDF_t::result_type
JFunction1D_t::result_type result_type
Definition JPDF_t.hh:145
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JRECONSTRUCTION::JMuonGandalfParameters_t
Data structure for fit parameters.
Definition JMuonGandalfParameters_t.hh:23
JRECONSTRUCTION::JMuonGandalfParameters_t::E_GeV
double E_GeV
energy [GeV]
Definition JMuonGandalfParameters_t.hh:83
JRECONSTRUCTION::JMuonGandalfParameters_t::TTS_ns
double TTS_ns
transition-time spread [ns]
Definition JMuonGandalfParameters_t.hh:82
JRECONSTRUCTION::JMuonGandalfParameters_t::TMin_ns
double TMin_ns
minimal time w.r.t. Cherenkov hypothesis [ns]
Definition JMuonGandalfParameters_t.hh:84
JRECONSTRUCTION::JMuonGandalfParameters_t::roadWidth_m
double roadWidth_m
road width [m]
Definition JMuonGandalfParameters_t.hh:79
JRECONSTRUCTION::JMuonGandalfParameters_t::TMax_ns
double TMax_ns
maximal time w.r.t. Cherenkov hypothesis [ns]
Definition JMuonGandalfParameters_t.hh:85
JRECONSTRUCTION::JMuonGandalfParameters_t::VMax_npe
double VMax_npe
maximum number of of photo-electrons
Definition JMuonGandalfParameters_t.hh:88
JRECONSTRUCTION::JMuonGandalfParameters_t::R_Hz
double R_Hz
default rate [Hz]
Definition JMuonGandalfParameters_t.hh:80
JRECONSTRUCTION::JMuonGandalfParameters_t::numberOfPrefits
size_t numberOfPrefits
number of prefits
Definition JMuonGandalfParameters_t.hh:81
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
SCIENTIFIC
Auxiliary data structure for floating point format specification.
Definition JManip.hh:488
Trk
The Trk class represents a Monte Carlo (MC) particle as well as a reconstructed track/shower.
Definition Trk.hh:15
Trk::print
void print(std::ostream &out=std::cout) const
Print track.
Definition Trk.hh:182
Trk::dir
Vec dir
track direction
Definition Trk.hh:18
Trk::E
double E
Energy [GeV] (either MC truth or reconstructed)
Definition Trk.hh:20
Trk::t
double t
track time [ns] (when the particle is at pos )
Definition Trk.hh:19
Trk::len
double len
length, if applicable [m]
Definition Trk.hh:22
Trk::lik
double lik
likelihood or lambda value (for aafit, lambda)
Definition Trk.hh:23
Vec::z
double z
Definition Vec.hh:14
rec_stages_range
Range of reconstruction stages.
Definition tools/reconstruction.hh:26
Generated by doxygen 1.12.0