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