Do the actual processing.
286{
289
291 if (m_IsPhysicsRun == true) {
292
293
294 auto hist_comL1_GDLL1 =
getDelta(
"TRGGDL",
"hGDL_gdlL1TocomL1_all", 0,
true);
295 if (hist_comL1_GDLL1) {
296 double comL1_GDLL1_mean = 0.0;
297 hist_comL1_GDLL1->Draw();
298 comL1_GDLL1_mean = hist_comL1_GDLL1->GetMean();
299 B2DEBUG(1, "comL1_GDLL1_mean:" << comL1_GDLL1_mean);
300 setEpicsPV(
"comL1_GDLL1_mean", comL1_GDLL1_mean);
301 }
302
303
304 auto hist =
getDelta(
"EventT0DQMdir",
"m_histEventT0_TOP_hadron_L1_ECLTRG", 0,
true);
305 if (hist) {
306 double ECLTRG_timing_mean = 0.0;
307 hist->Draw();
308 ECLTRG_timing_mean = hist->GetMean();
309 B2DEBUG(1, "ECLTRG_timing_mean:" << ECLTRG_timing_mean);
310 setEpicsPV(
"ECLTRG_timing_mean", ECLTRG_timing_mean);
311 }
312
313
314
315 auto histCDCTRG =
getDelta(
"EventT0DQMdir",
"m_histEventT0_TOP_hadron_L1_CDCTRG", 0,
true);
316 if (histCDCTRG) {
317 double CDCTRG_timing_mean = 0.0;
318 histCDCTRG->Draw();
319 CDCTRG_timing_mean = histCDCTRG->GetMean();
320 B2DEBUG(1, "CDCTRG_timing_mean:" << CDCTRG_timing_mean);
321 setEpicsPV(
"CDCTRG_timing_mean", CDCTRG_timing_mean);
322 }
323
324
325 auto histTOPTRG =
getDelta(
"EventT0DQMdir",
"m_histEventT0_TOP_hadron_L1_TOPTRG", 0,
true);
326 if (histTOPTRG) {
327 double TOPTRG_timing_mean = 0.0;
328 histTOPTRG->Draw();
329 TOPTRG_timing_mean = histTOPTRG->GetMean();
330 B2DEBUG(1, "TOPTRG_timing_mean:" << TOPTRG_timing_mean);
331 setEpicsPV(
"TOPTRG_timing_mean", TOPTRG_timing_mean);
332 }
333
334
335 auto hist_ECLTRG_peak =
getDelta(
"TRGGRL",
"h_ECLL1", 0,
true);
336 if (hist_ECLTRG_peak) {
337
338 double ECLTRG_peak = 0.0;
339 hist_ECLTRG_peak->Draw();
340
341
342 ECLTRG_peak = hist_ECLTRG_peak->GetMean();
343
344
345
346
347
348
349
350
351
352 B2DEBUG(1, "ECLTRG_peak:" << ECLTRG_peak);
354
355 }
356
357
358 auto hist_ECLTRG_deadch =
getDelta(
"TRG",
"h_TCId", 0,
true);
359 if (hist_ECLTRG_deadch) {
360 hist_ECLTRG_deadch->Draw();
361 int numberOfBins = hist_ECLTRG_deadch->GetNbinsX();
362 double average = 0;
363 for (int i = 80; i < 512; i++) {
364 average += hist_ECLTRG_deadch->GetBinContent(i);
365 }
366 average = 1.0 * average / (512 - 80);
367 int ECLTRG_deadch = 0;
368
369
370
371 for (int i = 2; i <= numberOfBins - 1; i++) {
372 if (hist_ECLTRG_deadch->GetBinContent(i) <= 0.01 * average) { ECLTRG_deadch += 1; }
373 }
374 B2DEBUG(1, "ECLTRG_deadch:" << ECLTRG_deadch);
376 }
377
378
379 auto hist_TOPTRG_deadch =
getDelta(
"TRGGRL",
"h_slot_TOP", 0,
true);
380 if (hist_TOPTRG_deadch) {
381 hist_TOPTRG_deadch->Draw();
382 int numberOfBins = hist_TOPTRG_deadch->GetNbinsX();
383 int TOPTRG_deadch = 0;
384 for (int i = 2; i <= numberOfBins - 1; i++) {
385 if (hist_TOPTRG_deadch->GetBinContent(i) <= 0.01 * hist_TOPTRG_deadch->GetMaximum()) { TOPTRG_deadch += 1; }
386 }
387 B2DEBUG(1, "TOPTRG_deadch:" << TOPTRG_deadch);
389 }
390
391
392
393
394 auto hist_CDCTRG_2D_peak =
getDelta(
"TRGGRL",
"h_CDCL1", 0,
true);
395 if (hist_CDCTRG_2D_peak) {
396 double CDCTRG_2D_peak = 0.0;
397 hist_CDCTRG_2D_peak->Draw();
398 int bin_CDCTRG_2D_peak = hist_CDCTRG_2D_peak->GetMaximumBin();
399 CDCTRG_2D_peak = hist_CDCTRG_2D_peak->GetXaxis()->GetBinCenter(bin_CDCTRG_2D_peak);
400 B2DEBUG(1, "CDCTRG_2D_peak:" << CDCTRG_2D_peak);
402
403 }
404
405
406 auto hist_NN_peak =
getDelta(
"TRGGRL",
"h_CDCNNL1", 0,
true);
407 if (hist_NN_peak) {
408 double NN_peak = 0.0;
409 hist_NN_peak->Draw();
410 int bin_NN_peak = hist_NN_peak->GetMaximumBin();
411 NN_peak = hist_NN_peak->GetXaxis()->GetBinCenter(bin_NN_peak);
412 B2DEBUG(1, "NN_peak:" << NN_peak);
414
415 }
416
417
418 auto hist_CDCTRG_TSF_peak =
getDelta(
"TRGGRL",
"h_TSFL1", 0,
true);
419 if (hist_CDCTRG_TSF_peak) {
420 double CDCTRG_TSF_peak = 0.0;
421 hist_CDCTRG_TSF_peak->Draw();
422 int bin_CDCTRG_TSF_peak = hist_CDCTRG_TSF_peak->GetMaximumBin();
423 CDCTRG_TSF_peak = hist_CDCTRG_TSF_peak->GetXaxis()->GetBinCenter(bin_CDCTRG_TSF_peak);
424 B2DEBUG(1, "CDCTRG_TSF_peak:" << CDCTRG_TSF_peak);
425 setEpicsPV(
"CDCTRG_TSF_peak", CDCTRG_TSF_peak);
426
427 }
428
429
430 auto hist_CDCTRG_deadch =
getDelta(
"TRGCDCTNN",
"NeuroHWInTSID", 0,
true);
431 if (hist_CDCTRG_deadch) {
432 hist_CDCTRG_deadch->Draw();
433 int numberOfBins = hist_CDCTRG_deadch->GetNbinsX();
434
435 int CDCTRG_deadch = 0;
436 for (int i = 3; i <= numberOfBins; i++) {
437 if (hist_CDCTRG_deadch->GetBinContent(i) <= 0.01 * hist_CDCTRG_deadch->GetMaximum()) {CDCTRG_deadch += 1;}
438
439 }
440 B2DEBUG(1, "CDCTRG_deadch:" << CDCTRG_deadch);
442 }
443
444
445 auto hist_KLMTRG_peak =
getDelta(
"TRGGRL",
"h_KLML1", 0,
true);
446 if (hist_KLMTRG_peak) {
447 double KLMTRG_peak = 0.0;
448 hist_KLMTRG_peak->Draw();
449 int bin_KLMTRG_peak = hist_KLMTRG_peak->GetMaximumBin();
450 KLMTRG_peak = hist_KLMTRG_peak->GetXaxis()->GetBinCenter(bin_KLMTRG_peak);
451 B2DEBUG(1, "KLMTRG_peak:" << KLMTRG_peak);
453
454 }
455
456
457 auto hist_TOPTRG_peak =
getDelta(
"TRGGRL",
"h_TOPL1", 0,
true);
458 if (hist_TOPTRG_peak) {
459 double TOPTRG_peak = 0.0;
460 hist_TOPTRG_peak->Draw();
461 int bin_TOPTRG_peak = hist_TOPTRG_peak->GetMaximumBin();
462 TOPTRG_peak = hist_TOPTRG_peak->GetXaxis()->GetBinCenter(bin_TOPTRG_peak);
463 B2DEBUG(1, "TOPTRG_peak:" << TOPTRG_peak);
465
466 }
467
468
469 auto hist_ecltiming_cdctiming =
getDelta(
"TRGGDL",
"hGDL_ns_cdcToecl_all", 0,
true);
470 if (hist_ecltiming_cdctiming) {
471 double ecltiming_cdctiming = 0.0;
472 hist_ecltiming_cdctiming->Draw();
473 int bin_ecltiming_cdctiming = hist_ecltiming_cdctiming->GetMaximumBin();
474 ecltiming_cdctiming = hist_ecltiming_cdctiming->GetXaxis()->GetBinCenter(bin_ecltiming_cdctiming);
475 B2DEBUG(1, "ecltiming_cdctiming:" << ecltiming_cdctiming);
476 setEpicsPV(
"ecltiming_cdctiming", ecltiming_cdctiming);
477 }
478
479
480 auto hist_toptiming_ecltiming =
getDelta(
"TRGGDL",
"hGDL_ns_topToecl_all", 0,
true);
481 if (hist_toptiming_ecltiming) {
482 double toptiming_ecltiming = 0.0;
483 hist_toptiming_ecltiming->Draw();
484 int bin_toptiming_ecltiming = hist_toptiming_ecltiming->GetMaximumBin();
485 toptiming_ecltiming = hist_toptiming_ecltiming->GetXaxis()->GetBinCenter(bin_toptiming_ecltiming);
486 B2DEBUG(1, "toptiming_ecltiming:" << toptiming_ecltiming);
487 setEpicsPV(
"toptiming_ecltiming", toptiming_ecltiming);
488 }
489
490
491 auto hist_toptiming_cdctiming =
getDelta(
"TRGGDL",
"hGDL_ns_topTocdc_all", 0,
true);
492 if (hist_toptiming_cdctiming) {
493 double toptiming_cdctiming = 0.0;
494 hist_toptiming_cdctiming->Draw();
495 int bin_toptiming_cdctiming = hist_toptiming_cdctiming->GetMaximumBin();
496 toptiming_cdctiming = hist_toptiming_cdctiming->GetXaxis()->GetBinCenter(bin_toptiming_cdctiming);
497 B2DEBUG(1, "toptiming_cdctiming:" << toptiming_cdctiming);
498 setEpicsPV(
"toptiming_cdctiming", toptiming_cdctiming);
499 }
500
501
502 auto hist_gdll1_ecltiming =
getDelta(
"TRGGDL",
"hGDL_eclTogdlL1_all", 0,
true);
503 if (hist_gdll1_ecltiming) {
504 double gdll1_ecltiming = 0.0;
505 hist_gdll1_ecltiming->Draw();
506 int bin_gdll1_ecltiming = hist_gdll1_ecltiming->GetMaximumBin();
507 gdll1_ecltiming = hist_gdll1_ecltiming->GetXaxis()->GetBinCenter(bin_gdll1_ecltiming);
508 B2DEBUG(1, "gdll1_ecltiming:" << gdll1_ecltiming);
509 setEpicsPV(
"gdll1_ecltiming", gdll1_ecltiming);
510 }
511
512
513 auto hist_gdll1_cdctiming =
getDelta(
"TRGGDL",
"hGDL_cdcTogdlL1_all", 0,
true);
514 if (hist_gdll1_cdctiming) {
515 double gdll1_cdctiming = 0.0;
516 hist_gdll1_cdctiming->Draw();
517 int bin_gdll1_cdctiming = hist_gdll1_cdctiming->GetMaximumBin();
518 gdll1_cdctiming = hist_gdll1_cdctiming->GetXaxis()->GetBinCenter(bin_gdll1_cdctiming);
519 B2DEBUG(1, "gdll1_cdctiming:" << gdll1_cdctiming);
520 setEpicsPV(
"gdll1_cdctiming", gdll1_cdctiming);
521 }
522
523
524 auto hist_gdll1_toptiming =
getDelta(
"TRGGDL",
"hGDL_topTogdlL1_all", 0,
true);
525 if (hist_gdll1_toptiming) {
526 double gdll1_toptiming = 0.0;
527 hist_gdll1_toptiming->Draw();
528 int bin_gdll1_toptiming = hist_gdll1_toptiming->GetMaximumBin();
529 gdll1_toptiming = hist_gdll1_toptiming->GetXaxis()->GetBinCenter(bin_gdll1_toptiming);
530 B2DEBUG(1, "gdll1_toptiming:" << gdll1_toptiming);
531 setEpicsPV(
"gdll1_toptiming", gdll1_toptiming);
532 }
533
534
535 auto hist_barrel_klm_latency =
getDelta(
"TRGGDL",
"hGDL_itd_klm_hit_rise_all", 0,
true);
536 if (hist_barrel_klm_latency) {
537 double barrel_klm_latency = 0.0;
538 hist_barrel_klm_latency->Draw();
539 int bin_barrel_klm_latency = hist_barrel_klm_latency->GetMaximumBin();
540 barrel_klm_latency = hist_barrel_klm_latency->GetXaxis()->GetBinCenter(bin_barrel_klm_latency);
541 B2DEBUG(1, "barrel_klm_latency:" << barrel_klm_latency);
542 setEpicsPV(
"barrel_klm_latency", barrel_klm_latency);
543 }
544
545
546 auto hist_endcap_klm_latency =
getDelta(
"TRGGDL",
"hGDL_itd_eklm_hit_rise_all", 0,
true);
547 if (hist_endcap_klm_latency) {
548 double endcap_klm_latency = 0.0;
549 hist_endcap_klm_latency->Draw();
550 int bin_endcap_klm_latency = hist_endcap_klm_latency->GetMaximumBin();
551 endcap_klm_latency = hist_endcap_klm_latency->GetXaxis()->GetBinCenter(bin_endcap_klm_latency);
552 B2DEBUG(1, "endcap_klm_latency:" << endcap_klm_latency);
553 setEpicsPV(
"endcap_klm_latency", endcap_klm_latency);
554 }
555
556 bool m_IsPhysicsRun_ratio = (
getRunType() ==
"physics");
557 if (m_IsPhysicsRun_ratio == true) {
558
559
560 auto hist_software_trigger =
getDelta(
"softwaretrigger",
"skim", 0,
true);
561
562 if (hist_software_trigger) {
563 hist_software_trigger->Draw();
564
565 if (hist_software_trigger->GetBinContent(
566 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all")) != 0) {
567 double hadronb2_over_bhabha_all = 0.0;
568 hadronb2_over_bhabha_all = hist_software_trigger->GetBinContent(
569 hist_software_trigger->GetXaxis()->FindBin("accept_hadronb2")) / hist_software_trigger->GetBinContent(
570 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all"));
571
572 B2DEBUG(1, "hadronb2_over_bhabha_all:" << hadronb2_over_bhabha_all);
573 setEpicsPV(
"hadronb2_over_bhabha_all", hadronb2_over_bhabha_all);
574 }
575
576 if (hist_software_trigger->GetBinContent(
577 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all")) != 0) {
578 double mumu2trk_over_bhabha_all = 0.0;
579 mumu2trk_over_bhabha_all = hist_software_trigger->GetBinContent(
580 hist_software_trigger->GetXaxis()->FindBin("accept_mumu_2trk")) / hist_software_trigger->GetBinContent(
581 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all"));
582 B2DEBUG(1, "mumu2trk_over_bhabha_all:" << mumu2trk_over_bhabha_all);
583 setEpicsPV(
"mumu2trk_over_bhabha_all", mumu2trk_over_bhabha_all);
584 }
585
586 if (hist_software_trigger->GetBinContent(
587 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl")) != 0) {
588 double hadronb2_over_bhabha_trk_ecl = 0.0;
589 hadronb2_over_bhabha_trk_ecl = hist_software_trigger->GetBinContent(
590 hist_software_trigger->GetXaxis()->FindBin("accept_hadronb2")) / hist_software_trigger->GetBinContent(
591 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl"));
592
593 B2DEBUG(1, "hadronb2_over_bhabha_trk_ecl:" << hadronb2_over_bhabha_trk_ecl);
594 setEpicsPV(
"hadronb2_over_bhabha_trk_ecl", hadronb2_over_bhabha_trk_ecl);
595 }
596
597 if (hist_software_trigger->GetBinContent(
598 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl")) != 0) {
599 double mumu2trk_over_bhabha_trk_ecl = 0.0;
600 mumu2trk_over_bhabha_trk_ecl = hist_software_trigger->GetBinContent(
601 hist_software_trigger->GetXaxis()->FindBin("accept_mumu_2trk")) / hist_software_trigger->GetBinContent(
602 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl"));
603 B2DEBUG(1, "mumu2trk_over_bhabha_trk_ecl:" << mumu2trk_over_bhabha_trk_ecl);
604 setEpicsPV(
"mumu2trk_over_bhabha_trk_ecl", mumu2trk_over_bhabha_trk_ecl);
605 }
606
607 if (hist_software_trigger->GetBinContent(
608 hist_software_trigger->GetXaxis()->FindBin("accept_mumu_2trk")) != 0) {
609 double hadronb2_over_mumu2trk = 0.0;
610 hadronb2_over_mumu2trk = hist_software_trigger->GetBinContent(
611 hist_software_trigger->GetXaxis()->FindBin("accept_hadronb2")) / hist_software_trigger->GetBinContent(
612 hist_software_trigger->GetXaxis()->FindBin("accept_mumu_2trk"));
613 B2DEBUG(1, "hadronb2_over_mumu2trk:" << hadronb2_over_mumu2trk);
614 setEpicsPV(
"hadronb2_over_mumu2trk", hadronb2_over_mumu2trk);
615 }
616
617 if (hist_software_trigger->GetBinContent(
618 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all")) != 0) {
619 double mumu_tight_over_bhabha_all = 0.0;
620 mumu_tight_over_bhabha_all = hist_software_trigger->GetBinContent(
621 hist_software_trigger->GetXaxis()->FindBin("accept_mumutight")) / hist_software_trigger->GetBinContent(
622 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all"));
623 B2DEBUG(1, "mumu_tight_over_bhabha_all:" << mumu_tight_over_bhabha_all);
624 setEpicsPV(
"mumu_tight_over_bhabha_all", mumu_tight_over_bhabha_all);
625 }
626
627
628 if (hist_software_trigger->GetBinContent(
629 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all")) != 0) {
630 double gammagamma_over_bhabha_all = 0.0;
631 gammagamma_over_bhabha_all = hist_software_trigger->GetBinContent(
632 hist_software_trigger->GetXaxis()->FindBin("accept_gamma_gamma")) /
633 hist_software_trigger->GetBinContent(
634 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all"));
635 B2DEBUG(1, "gammagamma_over_bhabha_all:" << gammagamma_over_bhabha_all);
636 setEpicsPV(
"gammagamma_over_bhabha_all", gammagamma_over_bhabha_all);
637 }
638
639
640 if (hist_software_trigger->GetBinContent(
641 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all")) != 0) {
642 double tautau2trk_over_bhabha_all = 0.0;
643 tautau2trk_over_bhabha_all = hist_software_trigger->GetBinContent(
644 hist_software_trigger->GetXaxis()->FindBin("accept_tau_2trk")) / hist_software_trigger->GetBinContent(
645 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all"));
646 B2DEBUG(1, "tautau2trk_over_bhabha_all:" << tautau2trk_over_bhabha_all);
647 setEpicsPV(
"tautau2trk_over_bhabha_all", tautau2trk_over_bhabha_all);
648 }
649
650 if (hist_software_trigger->GetBinContent(
651 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all")) != 0) {
652 double hadron_over_bhabha_all = 0.0;
653 hadron_over_bhabha_all = hist_software_trigger->GetBinContent(
654 hist_software_trigger->GetXaxis()->FindBin("accept_hadron")) / hist_software_trigger->GetBinContent(
655 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_all"));
656 B2DEBUG(1, "hadron_over_bhabha_all:" << hadron_over_bhabha_all);
657 setEpicsPV(
"hadron_over_bhabha_all", hadron_over_bhabha_all);
658 }
659
660 if (hist_software_trigger->GetBinContent(
661 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl")) != 0) {
662 double mumu_tight_over_bhabha_trk_ecl = 0.0;
663 mumu_tight_over_bhabha_trk_ecl = hist_software_trigger->GetBinContent(
664 hist_software_trigger->GetXaxis()->FindBin("accept_mumutight")) / hist_software_trigger->GetBinContent(
665 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl"));
666 B2DEBUG(1, "mumu_tight_over_bhabha_trk_ecl:" << mumu_tight_over_bhabha_trk_ecl);
667 setEpicsPV(
"mumu_tight_over_bhabha_trk_ecl", mumu_tight_over_bhabha_trk_ecl);
668 }
669
670 if (hist_software_trigger->GetBinContent(
671 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl")) != 0) {
672 double gammagamma_over_bhabha_trk_ecl = 0.0;
673 gammagamma_over_bhabha_trk_ecl = hist_software_trigger->GetBinContent(
674 hist_software_trigger->GetXaxis()->FindBin("accept_gamma_gamma")) /
675 hist_software_trigger->GetBinContent(
676 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl"));
677 B2DEBUG(1, "gammagamma_over_bhabha_trk_ecl:" << gammagamma_over_bhabha_trk_ecl);
678 setEpicsPV(
"gammagamma_over_bhabha_trk_ecl", gammagamma_over_bhabha_trk_ecl);
679 }
680
681
682 if (hist_software_trigger->GetBinContent(
683 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl")) != 0) {
684 double tautau2trk_over_bhabha_trk_ecl = 0.0;
685 tautau2trk_over_bhabha_trk_ecl = hist_software_trigger->GetBinContent(
686 hist_software_trigger->GetXaxis()->FindBin("accept_tau_2trk")) / hist_software_trigger->GetBinContent(
687 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl"));
688 B2DEBUG(1, "tautau2trk_over_bhabha_trk_ecl:" << tautau2trk_over_bhabha_trk_ecl);
689 setEpicsPV(
"tautau2trk_over_bhabha_trk_ecl", tautau2trk_over_bhabha_trk_ecl);
690 }
691
692
693 if (hist_software_trigger->GetBinContent(
694 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl")) != 0) {
695 double hadron_over_bhabha_trk_ecl = 0.0;
696 hadron_over_bhabha_trk_ecl = hist_software_trigger->GetBinContent(
697 hist_software_trigger->GetXaxis()->FindBin("accept_hadron")) / hist_software_trigger->GetBinContent(
698 hist_software_trigger->GetXaxis()->FindBin("accept_bhabha_trk_ecl"));
699 B2DEBUG(1, "hadron_over_bhabha_trk_ecl:" << hadron_over_bhabha_trk_ecl);
700 setEpicsPV(
"hadron_over_bhabha_trk_ecl", hadron_over_bhabha_trk_ecl);
701 }
702
703 }
704
705 }
706
707 auto hist_nTSFHits_total =
getDelta(
"TRGGRL",
"h_wirecnt_sum", 0,
true);
708 if (hist_nTSFHits_total) {
709 double mean = hist_nTSFHits_total->GetMean();
710 B2DEBUG(1, "CDCTRG_meanTSFHits_total:" << mean);
712 }
713
714 auto hist_nTSFHits_clean =
getDelta(
"TRGGRL",
"h_wirecnt_sum_clean", 0,
true);
715 if (hist_nTSFHits_clean) {
716 double mean = hist_nTSFHits_clean->GetMean();
717 B2DEBUG(1, "CDCTRG_meanTSFHits_clean:" << mean);
719
720 double tail =
getThreshold(hist_nTSFHits_clean, 0.5, 0.02);
721 B2DEBUG(1, "CDCTRG_tailTSFHits_clean:" << tail);
723 }
724
725 auto hist_nTSFHits_injHER =
getDelta(
"TRGGRL",
"h_wirecnt_sum_injHER", 0,
true);
726 if (hist_nTSFHits_injHER) {
727 double mean = hist_nTSFHits_injHER->GetMean();
728 B2DEBUG(1, "CDCTRG_meanTSFHits_injHER:" << mean);
729 setEpicsPV(
"CDCTRG_meanTSFHits_injHER", mean);
730 }
731
732 auto hist_nTSFHits_injLER =
getDelta(
"TRGGRL",
"h_wirecnt_sum_injLER", 0,
true);
733 if (hist_nTSFHits_injLER) {
734 double mean = hist_nTSFHits_injLER->GetMean();
735 B2DEBUG(1, "CDCTRG_meanTSFHits_injLER:" << mean);
736 setEpicsPV(
"CDCTRG_meanTSFHits_injLER", mean);
737 }
738
739 auto hist_nTC_total =
getDelta(
"TRG",
"h_n_TChit_event_clkgrp", 0,
true);
740 if (hist_nTC_total) {
741 double mean = hist_nTC_total->GetMean();
742 B2DEBUG(1, "ECLTRG_meanTC_total:" << mean);
744 }
745
746 auto hist_nTC_clean =
getDelta(
"TRG",
"h_n_TChit_clean_clkgrp", 0,
true);
747 if (hist_nTC_clean) {
748 double mean = hist_nTC_clean->GetMean();
749 B2DEBUG(1, "ECLTRG_meanTC_clean:" << mean);
751
753 B2DEBUG(1, "ECLTRG_tailTC_clean:" << tail);
755 }
756
757 auto hist_nTC_injHER =
getDelta(
"TRG",
"h_n_TChit_injHER_clkgrp", 0,
true);
758 if (hist_nTC_injHER) {
759 double mean = hist_nTC_injHER->GetMean();
760 B2DEBUG(1, "ECLTRG_meanTC_injHER:" << mean);
762 }
763
764 auto hist_nTC_injLER =
getDelta(
"TRG",
"h_n_TChit_injLER_clkgrp", 0,
true);
765 if (hist_nTC_injLER) {
766 double mean = hist_nTC_injLER->GetMean();
767 B2DEBUG(1, "ECLTRG_meanTC_injLER:" << mean);
769 }
770
771 auto hist_nTC_FWD_total =
getDelta(
"TRG",
"h_n_TChit_FWD_event_clkgrp", 0,
true);
772 if (hist_nTC_FWD_total) {
773 double mean = hist_nTC_FWD_total->GetMean();
774 B2DEBUG(1, "ECLTRG_meanTC_FWD_total:" << mean);
776 }
777
778 auto hist_nTC_FWD_clean =
getDelta(
"TRG",
"h_n_TChit_FWD_clean_clkgrp", 0,
true);
779 if (hist_nTC_FWD_clean) {
780 double mean = hist_nTC_FWD_clean->GetMean();
781 B2DEBUG(1, "ECLTRG_meanTC_FWD_clean:" << mean);
783
785 B2DEBUG(1, "ECLTRG_tailTC_FWD_clean:" << tail);
787 }
788
789 auto hist_nTC_FWD_injHER =
getDelta(
"TRG",
"h_n_TChit_FWD_injHER_clkgrp", 0,
true);
790 if (hist_nTC_FWD_injHER) {
791 double mean = hist_nTC_FWD_injHER->GetMean();
792 B2DEBUG(1, "ECLTRG_meanTC_FWD_injHER:" << mean);
794 }
795
796 auto hist_nTC_FWD_injLER =
getDelta(
"TRG",
"h_n_TChit_FWD_injLER_clkgrp", 0,
true);
797 if (hist_nTC_FWD_injLER) {
798 double mean = hist_nTC_FWD_injLER->GetMean();
799 B2DEBUG(1, "ECLTRG_meanTC_FWD_injLER:" << mean);
801 }
802
803 auto hist_nTC_BRL_total =
getDelta(
"TRG",
"h_n_TChit_BRL_event_clkgrp", 0,
true);
804 if (hist_nTC_BRL_total) {
805 double mean = hist_nTC_BRL_total->GetMean();
806 B2DEBUG(1, "ECLTRG_meanTC_BRL_total:" << mean);
808 }
809
810 auto hist_nTC_BRL_clean =
getDelta(
"TRG",
"h_n_TChit_BRL_clean_clkgrp", 0,
true);
811 if (hist_nTC_BRL_clean) {
812 double mean = hist_nTC_BRL_clean->GetMean();
813 B2DEBUG(1, "ECLTRG_meanTC_BRL_clean:" << mean);
815
817 B2DEBUG(1, "ECLTRG_tailTC_BRL_clean:" << tail);
819 }
820
821 auto hist_nTC_BRL_injHER =
getDelta(
"TRG",
"h_n_TChit_BRL_injHER_clkgrp", 0,
true);
822 if (hist_nTC_BRL_injHER) {
823 double mean = hist_nTC_BRL_injHER->GetMean();
824 B2DEBUG(1, "ECLTRG_meanTC_BRL_injHER:" << mean);
826 }
827
828 auto hist_nTC_BRL_injLER =
getDelta(
"TRG",
"h_n_TChit_BRL_injLER_clkgrp", 0,
true);
829 if (hist_nTC_BRL_injLER) {
830 double mean = hist_nTC_BRL_injLER->GetMean();
831 B2DEBUG(1, "ECLTRG_meanTC_BRL_injLER:" << mean);
833 }
834
835 auto hist_nTC_BWD_total =
getDelta(
"TRG",
"h_n_TChit_BWD_event_clkgrp", 0,
true);
836 if (hist_nTC_BWD_total) {
837 double mean = hist_nTC_BWD_total->GetMean();
838 B2DEBUG(1, "ECLTRG_meanTC_BWD_total:" << mean);
840 }
841
842 auto hist_nTC_BWD_clean =
getDelta(
"TRG",
"h_n_TChit_BWD_clean_clkgrp", 0,
true);
843 if (hist_nTC_BWD_clean) {
844 double mean = hist_nTC_BWD_clean->GetMean();
845 B2DEBUG(1, "ECLTRG_meanTC_BWD_clean:" << mean);
847
849 B2DEBUG(1, "ECLTRG_tailTC_BWD_clean:" << tail);
851 }
852
853 auto hist_nTC_BWD_injHER =
getDelta(
"TRG",
"h_n_TChit_BWD_injHER_clkgrp", 0,
true);
854 if (hist_nTC_BWD_injHER) {
855 double mean = hist_nTC_BWD_injHER->GetMean();
856 B2DEBUG(1, "ECLTRG_meanTC_BWD_injHER:" << mean);
858 }
859
860 auto hist_nTC_BWD_injLER =
getDelta(
"TRG",
"h_n_TChit_BWD_injLER_clkgrp", 0,
true);
861 if (hist_nTC_BWD_injLER) {
862 double mean = hist_nTC_BWD_injLER->GetMean();
863 B2DEBUG(1, "ECLTRG_meanTC_BWD_injLER:" << mean);
865 }
866
867
868
870
871
873 }
874
875}
static void UpdateCanvas(const std::string &name, bool updated=true)
Mark canvas as updated (or not)
static const std::string & getRunType(void)
Get the list of the reference histograms.
TH1 * getDelta(const std::string &fullname, int n=0, bool onlyIfUpdated=true)
Get Delta histogram.
void setEpicsPV(const std::string &keyname, double value)
Write value to a EPICS PV.
static double getThreshold(const TH1 *hist, const double lower_bound, const double widthFraction=0.01)
A utility function to find the location of the tail of histograms.
TCanvas * m_canvas
Final Canvas.