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