Jpp 20.0.0-195-g190c9e876
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 "TRootCanvas.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 195 of file JAAnet/JEvD.cc.

196{
197 using namespace std;
198 using namespace JPP;
199 using namespace KM3NETDAQ;
200
201 JSingleFileScanner<Evt> inputFile;
202 JLimit_t& numberOfEvents = inputFile.getLimit();
203 string pdfFile;
204 string outputFile;
205 JMuonGandalfParameters_t parameters;
206 int application;
207 JEventSelector event_selector;
208 JCanvas canvas;
209 bool batch;
210 struct : JStyle::JParameters {
211 double arrowSize = 0.003;
212 string arrowType = "|->";
213 double arrowScale = 250.0;
214 Width_t lineWidth = 2;
215 Style_t lineStyle = 1;
216 int nbinsX = 50;
217 int nbinsY = 250;
218 double T_ns = 0.0;
219 bool equalize = false;
220 } graphics;
221 string option;
222 int debug;
223
224
225 try {
226
227 JProperties properties = graphics.getProperties();
228
229 properties.insert(gmake_property(graphics.arrowSize));
230 properties.insert(gmake_property(graphics.arrowType));
231 properties.insert(gmake_property(graphics.arrowScale));
232 properties.insert(gmake_property(graphics.lineWidth));
233 properties.insert(gmake_property(graphics.lineStyle));
234 properties.insert(gmake_property(graphics.T_ns));
235 properties.insert(gmake_property(graphics.equalize));
236
237 parameters.numberOfPrefits = 1;
238
239 JParser<> zap("Program to display hit probabilities.");
240
241 zap['w'] = make_field(canvas, "size of canvas <nx>x<ny> [pixels]") = JCanvas(1200, 600);
242 zap['f'] = make_field(inputFile, "input file (output of JXXXMuonReconstruction.sh)");
243 zap['n'] = make_field(numberOfEvents) = JLimit::max();
244 zap['F'] = make_field(pdfFile);
245 zap['o'] = make_field(outputFile, "graphics output file name") = MAKE_STRING("display_" << WILDCARD << ".gif");
246 zap['@'] = make_field(parameters) = JPARSER::initialised();
247 zap['A'] = make_field(application) = JMUONENERGY, JMUONSTART, JMUONGANDALF, JMUONSIMPLEX, JMUONPREFIT;
248 zap['L'] = make_field(event_selector) = JPARSER::initialised();
249 zap['%'] = make_field(properties) = JPARSER::initialised();
250 zap['O'] = make_field(option, "draw option") = arrow_t, histogram_t;
251 zap['B'] = make_field(batch, "batch processing");
252 zap['d'] = make_field(debug) = 1;
253
254 zap(argc, argv);
255 }
256 catch(const exception& error) {
257 FATAL(error.what() << endl);
258 }
259
260 if (batch && outputFile == "") {
261 FATAL("Missing output file name " << outputFile << " in batch mode." << endl);
262 }
263
264 if (!batch && outputFile == "") {
265 outputFile = MAKE_STRING(WILDCARD << ".gif");
266 }
267
268 if (outputFile.find(WILDCARD) == string::npos) {
269 FATAL("Output file name " << outputFile << " has no wild card '" << WILDCARD << "'" << endl);
270 }
271
272 JHit::setSlewing(false);
273
274 JSummaryFileRouter summary(inputFile);
275
276 const JMuonPDF_t pdf(pdfFile, parameters.TTS_ns);
277
278 const JTimeRange T_ns(parameters.TMin_ns, parameters.TMax_ns);
279
280 typedef vector<JHitL0> JDataL0_t;
281 typedef vector<JHitW0> JDataW0_t;
282
283
284 // ROOT
285
286 gROOT->SetBatch(batch);
287
288 TApplication* tp = new TApplication("user", NULL, NULL);
289 TCanvas* cv = new TCanvas("display", "", canvas.x, canvas.y);
290
291 if (!batch) {
292 ((TRootCanvas *) cv->GetCanvasImp())->Connect("CloseWindow()", "TApplication", tp, "Terminate()");
293 }
294
295 unique_ptr<TStyle> gStyle(new JStyle("gplot", cv->GetWw(), cv->GetWh(), graphics));
296
297 gROOT->SetStyle("gplot");
298 gROOT->ForceStyle();
299
300 const size_t NUMBER_OF_PADS = 3;
301
302 cv->SetFillStyle(4000);
303 cv->SetFillColor(kWhite);
304
305 TPad* p1 = new TPad("p1", NULL, 0.0, 0.00, 1.0, 0.95);
306 TPad* p2 = new TPad("p2", NULL, 0.0, 0.95, 1.0, 1.00);
307
308 p1->Divide(NUMBER_OF_PADS, 1);
309
310 p1->Draw();
311 p2->Draw();
312
313 const double Dmax = 1000.0;
314 const double Rmin = 0.0;
315 const double Rmax = min(parameters.roadWidth_m, 0.4 * Dmax);
316 const double Tmin = min(parameters.TMin_ns, -10.0);
317 const double Tmax = max(parameters.TMax_ns, +100.0);
318 const double Amin = 0.002 * (Tmax - Tmin); // minimal arrow length [ns]
319 const double Amax = 0.8 * (Tmax - Tmin); // maximal arrow length [ns]
320 const double ymin = Tmin - (option == arrow_t ? 0.2 * Amax : 0.0);
321 const double ymax = Tmax + (option == arrow_t ? 0.5 * Amax : 0.0);
322
323 const string Xlabel[NUMBER_OF_PADS] = { "R [m]", "#phi [rad]", "z [m]" };
324 const double Xmin [NUMBER_OF_PADS] = { Rmin, -PI, -0.4 * Dmax };
325 const double Xmax [NUMBER_OF_PADS] = { Rmax, +PI, +0.4 * Dmax };
326
327 double Xs[NUMBER_OF_PADS];
328
329 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
330 Xs[i] = 0.003 * (Xmax[i] - Xmin[i]) * (0.5 * NUMBER_OF_PMTS); // x-offset arrow as function of PMT number
331 }
332
333 TH2D H2[NUMBER_OF_PADS];
334 TGraph G2[NUMBER_OF_PADS];
335
336 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
337
338 H2[i] = TH2D(MAKE_CSTRING("h" << i), NULL, graphics.nbinsX, Xmin[i] - Xs[i], Xmax[i] + Xs[i], graphics.nbinsY, ymin, ymax);
339
340 H2[i].GetXaxis()->SetTitle(Xlabel[i].c_str());
341 H2[i].GetYaxis()->SetTitle("#Deltat [ns]");
342
343 H2[i].GetXaxis()->CenterTitle(true);
344 H2[i].GetYaxis()->CenterTitle(true);
345
346 H2[i].SetStats(kFALSE);
347
348 G2[i].Set(2);
349
350 G2[i].SetPoint(0, H2[i].GetXaxis()->GetXmin(), 0.0);
351 G2[i].SetPoint(1, H2[i].GetXaxis()->GetXmax(), 0.0);
352
353 p1->cd(i+1);
354
355 H2[i].Draw("AXIS");
356 G2[i].Draw("SAME");
357 }
358
359
360 while (inputFile.hasNext()) {
361
362 cout << "event: " << setw(8) << inputFile.getCounter() << endl;
363
364 const Evt* evt = inputFile.next();
365
366 if (has_reconstructed_track<JPP_RECONSTRUCTION_TYPE>(*evt, rec_stages_range(application))) {
367
368 Trk fit = get_best_reconstructed_track<JPP_RECONSTRUCTION_TYPE>(*evt, rec_stages_range(application));
369
370 if (!event_selector(fit, *evt)) {
371 continue;
372 }
373
374 JDataL0_t dataL0;
375
376 for (const Hit& hit : evt->hits) {
377 dataL0.push_back(JHitL0(JDAQPMTIdentifier(hit.dom_id, hit.channel_id),
378 JAxis3D(getPosition(hit.pos), getDirection(hit.dir)),
379 JHit(hit.t, hit.tot)));
380 }
381
382 summary.update(JDAQChronometer(evt->det_id,
383 evt->run_id,
384 evt->frame_index,
385 JDAQUTCExtended(evt->t.GetSec(), evt->t.GetNanoSec() / 16)));
386
387 const time_converter converter = time_converter(*evt);
388
389 Trk muon; // Monte Carlo true muon
390
391 if (has_muon(*evt)) {
392
393 for (const auto& trk : evt->mc_trks) {
394 if (is_muon(trk)) {
395 if (trk.E > muon.E) {
396
397 muon = trk;
398 muon.t += converter.putTime();
399
400 setW(muon, JSTART_LENGTH_METRES, fabs(muon.len));
401 }
402 }
403 }
404 }
405
406
407 bool monte_carlo = false; // show Monte Carlo true muon
408
409 for (bool next = false; !next; ) {
410
411 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
412 H2[i].Reset();
413 }
414
415 Trk trk;
416
417 if (!monte_carlo)
418 trk = fit;
419 else
420 trk = muon;
421
422 JRotation3D R (getDirection(trk));
423 JLine1Z tz(getPosition (trk).rotate(R), trk.t);
424 JRange<double> Z_m;
425 /*
426 if (hasW(trk, JSTART_LENGTH_METRES)) {
427 Z_m = JRange<double>(0.0, getW(fit,JSTART_LENGTH_METRES)) + JRange<double>(parameters.ZMin_m, parameters.ZMax_m);
428 }
429 */
430 const JModel<JLine1Z> match(tz, parameters.roadWidth_m, T_ns, Z_m);
431
432 // hit selection based on fit result
433
434 JDataW0_t data;
435
436 for (JDataL0_t::const_iterator i = dataL0.begin(); i != dataL0.end(); ++i) {
437
438 const int type = 0;
439 const double QE = 1.0;
440 const double R_Hz = summary.getRate(i->getPMTIdentifier(), parameters.R_Hz);
441
442 JHitW0 hit(*i, type, QE, R_Hz);
443
444 hit.rotate(R);
445
446 if (match(hit)) {
447 data.push_back(hit);
448 }
449 }
450
451 // select first hit in PMT
452
453 sort(data.begin(), data.end(), JHitW0::compare);
454
455 JDataW0_t::iterator __end = unique(data.begin(), data.end(), equal_to<JDAQPMTIdentifier>());
456
457 double E_GeV = parameters.E_GeV;
458 /*
459 if (trk.E > 0.1) {
460 E_GeV = trk.E;
461 }
462 */
463
464 // move fit to geometrical center of hits
465
466 JRange<double> zs = JRange<double>::DEFAULT_RANGE();
467
468 for (JDataW0_t::iterator hit = data.begin(); hit != __end; ++hit) {
469
470 const double x = hit->getX() - tz.getX();
471 const double y = hit->getY() - tz.getY();
472 const double z = hit->getZ();
473 const double R = sqrt(x*x + y*y);
474
475 zs.include(z - R/getTanThetaC());
476 }
477
478 const double z0 = tz.getZ();
479 const double z1 = 0.5 * (zs.getLowerLimit() + zs.getUpperLimit());
480
481 tz.setZ(z1, getSpeedOfLight());
482
483
484 // graphics
485
486 ostringstream os;
487 vector<TArrow> arrow [NUMBER_OF_PADS];
488 vector<TMarker> marker[NUMBER_OF_PADS];
489
490 if (hasW(trk, JSTART_LENGTH_METRES) && getW(trk, JSTART_LENGTH_METRES) > 0.0) {
491
492 marker[2].push_back(TMarker(z0 - tz.getZ(), 0.0, kFullCircle));
493 marker[2].push_back(TMarker(z0 - tz.getZ() + getW(trk, JSTART_LENGTH_METRES), 0.0, kFullCircle));
494
495 static_cast<TAttMarker&>(marker[2][0]) = TAttMarker(kRed, kFullCircle, 0.7);
496 static_cast<TAttMarker&>(marker[2][1]) = TAttMarker(kRed, kFullCircle, 0.7);
497 }
498
499 DEBUG("trk: "
500 << FIXED(7,2) << tz.getX() << ' '
501 << FIXED(7,2) << tz.getY() << ' '
502 << FIXED(7,2) << tz.getZ() << ' '
503 << FIXED(12,2) << tz.getT() << endl);
504
505 double chi2 = 0;
506
507 for (JDataW0_t::const_iterator hit = data.begin(); hit != __end; ++hit) {
508
509 const double x = hit->getX() - tz.getX();
510 const double y = hit->getY() - tz.getY();
511 const double z = hit->getZ() - tz.getZ();
512 const double R = sqrt(x*x + y*y);
513
514 const double t1 = tz.getT() + (z + R * getTanThetaC()) * getInverseSpeedOfLight();
515
516 JDirection3D dir(hit->getDX(), hit->getDY(), hit->getDZ()); // PMT orientation
517
518 dir.rotate(JRotation3Z(-atan2(y,x))); // rotate PMT axis to x-z plane
519
520 const double theta = dir.getTheta();
521 const double phi = fabs(dir.getPhi()); // rotational symmetry of Cherenkov cone
522
523 //const double E = gWater.getE(E_GeV, z); // correct for energy loss
524 const double E = E_GeV;
525 const double dt = T_ns.constrain(hit->getT() - t1);
526
527 JMuonPDF_t::result_type H1 = pdf.calculate(E, R, theta, phi, dt);
528 JMuonPDF_t::result_type H0(hit->getR() * 1e-9, 0.0, T_ns);
529
530 if (H1.V >= parameters.VMax_npe) {
531 H1 *= parameters.VMax_npe / H1.V;
532 }
533
534 H1 += H0; // signal + background
535
536 chi2 += H1.getChi2() - H0.getChi2();
537
538 DEBUG("hit: "
539 << setw(8) << hit->getModuleID() << '.' << FILL(2,'0') << (int) hit->getPMTAddress() << FILL() << ' '
540 << SCIENTIFIC(8,2) << E << ' '
541 << FIXED(7,2) << R << ' '
542 << FIXED(7,4) << theta << ' '
543 << FIXED(7,4) << phi << ' '
544 << FIXED(7,3) << dt << ' '
545 << FIXED(7,3) << H1.getChi2() << ' '
546 << FIXED(7,3) << H0.getChi2() << endl);
547
548 const double derivative = H1.getDerivativeOfChi2() - H0.getDerivativeOfChi2();
549
550 double size = derivative * graphics.arrowScale; // size of arrow
551
552 if (fabs(size) < Amin) {
553 size = (size > 0.0 ? +Amin : -Amin);
554 } else if (fabs(size) > Amax) {
555 size = (size > 0.0 ? +Amax : -Amax);
556 }
557
558 const double X[NUMBER_OF_PADS] = { R, atan2(y,x), z - R/getTanThetaC() };
559
560 const double xs = (double) (NUMBER_OF_PMTS - 2 * hit->getPMTAddress()) / (double) NUMBER_OF_PMTS;
561
562 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
563
564 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());
565
566 a1.SetLineWidth(graphics.lineWidth);
567 a1.SetLineStyle(graphics.lineStyle);
568
569 arrow[i].push_back(a1);
570
571 H2[i].Fill(X[i], dt + graphics.T_ns);
572 }
573 }
574
575 if (graphics.equalize) {
576
577 double zmax = 0.0;
578
579 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
580 if (H2[i].GetMaximum() > zmax) {
581 zmax = H2[i].GetMaximum();
582 }
583 }
584
585 zmax *= 1.2;
586
587 for (size_t i = 0; i != NUMBER_OF_PADS; ++i) {
588 H2[i].SetMaximum(zmax);
589 }
590 }
591
592 os << FILL(6,'0') << evt->run_id << ":" << evt->frame_index << "/" << evt->trigger_counter << FILL();
593 os << " Q = " << FIXED(4,0) << trk.lik
594 << '/' << FIXED(4,0) << -chi2;
595 os << " E = " << SCIENTIFIC(7,1) << trk.E << " [GeV]";
596 os << " cos(#theta) = " << FIXED(6,3) << trk.dir.z;
597 os << " L = " << FIXED(6,2) << getW(trk, JSTART_LENGTH_METRES, 0.0) << " [m]";
598
599 if (monte_carlo)
600 os << " Monte Carlo";
601 else if (is_muon(muon))
602 os << " #Delta#alpha = " << FIXED(6,2) << getAngle(getDirection(muon), getDirection(trk)) << " [deg]";
603
604
605 // draw
606
607 TLatex title(0.05, 0.5, os.str().c_str());
608
609 title.SetTextAlign(12);
610 title.SetTextFont(42);
611 title.SetTextSize(0.6);
612
613 p2->cd();
614
615 title.Draw();
616
617 for (int i = 0; i != NUMBER_OF_PADS; ++i) {
618
619 p1->cd(i+1);
620
621 if (option == arrow_t) {
622
623 for (auto& a1 : arrow[i]) {
624 a1.Draw();
625 }
626
627 for (auto& m1 : marker[i]) {
628 m1.Draw();
629 }
630 }
631
632 if (option == histogram_t) {
633 H2[i].Draw("SAME");
634 }
635 }
636
637 cv->Update();
638
639
640 // action
641
642 if (batch) {
643
644 cv->SaveAs(replace(outputFile, WILDCARD, MAKE_STRING(inputFile.getCounter())).c_str());
645
646 next = true;
647
648 } else {
649
650 static int count = 0;
651
652 if (count++ == 0) {
653 cout << endl << "Type '?' for possible options." << endl;
654 }
655
656 for (bool user = true; user; ) {
657
658 cout << "\n> " << flush;
659
660 switch (JKeypress(true).get()) {
661
662 case '?':
663 cout << endl;
664 cout << "possible options: " << endl;
665 cout << 'q' << " -> " << "exit application" << endl;
666 cout << 'u' << " -> " << "update canvas" << endl;
667 cout << 's' << " -> " << "save graphics to file" << endl;
668 cout << 'M' << " -> " << "Monte Carlo true muon information" << endl;
669 cout << 'F' << " -> " << "fit information" << endl;
670 if (event_selector.is_valid()) {
671 cout << 'L' << " -> " << "reload event selector" << endl;
672 }
673 cout << 'r' << " -> " << "rewind input file" << endl;
674 cout << 'R' << " -> " << "switch to ROOT mode (quit ROOT to continue)" << endl;
675 cout << 'p' << " -> " << "print event information" << endl;
676 cout << ' ' << " -> " << "next event (as well as any other key)" << endl;
677 break;
678
679 case 'q':
680 cout << endl;
681 return 0;
682
683 case 'u':
684 cv->Update();
685 break;
686
687 case 's':
688 cv->SaveAs(replace(outputFile, WILDCARD, MAKE_STRING(inputFile.getCounter())).c_str());
689 break;
690
691 case 'M':
692 if (is_muon(muon))
693 monte_carlo = true;
694 else
695 ERROR(endl << "No Monte Carlo muon available." << endl);
696 user = false;
697 break;
698
699 case 'F':
700 monte_carlo = false;
701 user = false;
702 break;
703
704 case 'L':
705 if (event_selector.is_valid()) {
706 execute(MAKE_STRING("make -f " << getPath(argv[0]) << "/JMakeEventSelector libs"), 3);
707 event_selector.reload();
708 }
709 break;
710
711 case 'R':
712 tp->Run(kTRUE);
713 break;
714
715 case 'p':
716 cout << endl;
717 evt->print(cout);
718 if (debug >= debug_t) {
719 for (const auto& trk : evt->mc_trks) {
720 cout << "MC "; trk.print(cout); cout << endl;
721 }
722 for (const auto& trk : evt->trks) {
723 cout << "fit "; trk.print(cout); cout << endl;
724 }
725 for (const auto& hit : evt->hits) {
726 cout << "hit "; hit.print(cout); cout << endl;
727 }
728 }
729 break;
730
731 case 'r':
732 inputFile.rewind();
733
734 default:
735 next = true;
736 user = false;
737 break;
738 }
739 }
740 }
741 }
742 }
743 }
744 cout << endl;
745}
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:74
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:1139
JEEP::JProperties
Utility class to parse parameter values.
Definition JProperties.hh:503
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:25
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
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:27
JAANET::getDirection
JDirection3D getDirection(const Vec &dir)
Get direction.
Definition JAAnetToolkit.hh:221
JAANET::hasW
bool hasW(const Trk &trk, const int i)
Check availability of value.
Definition JAAnet/JEvD.cc:103
JAANET::has_muon
bool has_muon(const Evt &evt)
Test whether given event has a muon.
Definition JAAnetToolkit.hh:847
JAANET::setW
void setW(Trk &trk, const int i, const double value)
Set associated value.
Definition JAAnet/JEvD.cc:146
JAANET::getW
double getW(const Trk &track, const size_t index, const double value)
Get track information.
Definition JAAnetToolkit.hh:319
JAANET::getPosition
JPosition3D getPosition(const Vec &pos)
Get position.
Definition JAAnetToolkit.hh:185
JAANET::is_muon
bool is_muon(const Trk &track)
Test whether given track is a (anti-)muon.
Definition JAAnetToolkit.hh:369
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:48
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:652
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:611
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:85
JRECONSTRUCTION::JMuonGandalfParameters_t::TTS_ns
double TTS_ns
transition-time spread [ns]
Definition JMuonGandalfParameters_t.hh:84
JRECONSTRUCTION::JMuonGandalfParameters_t::TMin_ns
double TMin_ns
minimal time w.r.t. Cherenkov hypothesis [ns]
Definition JMuonGandalfParameters_t.hh:86
JRECONSTRUCTION::JMuonGandalfParameters_t::roadWidth_m
double roadWidth_m
road width [m]
Definition JMuonGandalfParameters_t.hh:81
JRECONSTRUCTION::JMuonGandalfParameters_t::TMax_ns
double TMax_ns
maximal time w.r.t. Cherenkov hypothesis [ns]
Definition JMuonGandalfParameters_t.hh:87
JRECONSTRUCTION::JMuonGandalfParameters_t::VMax_npe
double VMax_npe
maximum number of of photo-electrons
Definition JMuonGandalfParameters_t.hh:90
JRECONSTRUCTION::JMuonGandalfParameters_t::R_Hz
double R_Hz
default rate [Hz]
Definition JMuonGandalfParameters_t.hh:82
JRECONSTRUCTION::JMuonGandalfParameters_t::numberOfPrefits
size_t numberOfPrefits
number of prefits
Definition JMuonGandalfParameters_t.hh:83
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