RecoGeneralQa.cxx

Go to the documentation of this file.

/* Copyright (C) 2024 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: P.-A. Loizeau [committer] */
 
#include "algo/qa/RecoGeneralQa.h"
 
#include <cmath>
 
using cbm::algo::qa::Data;
using cbm::algo::qa::H1D;
 
namespace cbm::algo::qa
{
  RecoGeneralQa::RecoGeneralQa(const uint64_t& runStartTimeNs, std::shared_ptr<HistogramSender> pSender)
    : fRunStartTimeNs(runStartTimeNs)
    , fpSender(pSender)
  {
  }
 
  // ---   Execution   --------------------------------------------------------
  void RecoGeneralQa::operator()(const fles::Timeslice& ts, const GeneralQaData& qaData)
  {
    if (!fpSender.get()) {
      return;
    }
 
    if (fInitNotDone) {
      double_t dBegAxisX  = 0.0;
      double_t dSizeTsSec = 0.128;  
      if (ts.num_components() != 0 && ts.num_microslices(0) > 1) {
        dSizeTsSec = ((ts.descriptor(0, 1).idx - ts.descriptor(0, 0).idx) * ts.num_core_microslices()) * 1.0e-9;
      }
      auto cName = "processed_ts";
      auto cTitl = "Statistics of TS processed online";
      auto canv  = CanvasConfig(cName, cTitl, 2, 1);
      {
        auto pad              = PadConfig();
        auto name             = "timeslices_count_evo";
        auto titl             = "Number of processed TS vs time in run; time in run [s]; Nb TS []";
        int32_t nbBins        = std::ceil(7200.0 / (dSizeTsSec * kNbTsPerBinCount));  // a bit more than 2h range
        double_t dEndAxisX    = nbBins * dSizeTsSec * kNbTsPerBinCount;
        fphTimeslicesCountEvo = fQaData.MakeObj<H1D>(name, titl, nbBins, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphTimeslicesCountEvo, "hist");
        canv.AddPadConfig(pad);
      }
      {
        auto pad                 = PadConfig();
        auto name                = "timeslices_fraction_evo";
        auto titl                = "Fraction of TS processed vs time in run; time in run [s]; Processed TS []";
        int32_t nbBins           = std::ceil(7200.0 / (dSizeTsSec * kNbTsPerBinFrac));  // a bit more than 2h range
        double_t dEndAxisX       = nbBins * dSizeTsSec * kNbTsPerBinFrac;
        fphTimeslicesFractionEvo = fQaData.MakeObj<H1D>(name, titl, nbBins, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphTimeslicesFractionEvo, "hist");
        canv.AddPadConfig(pad);
      }
      fQaData.AddCanvasConfig(canv);
 
      auto cNameB = "processing_time";
      auto cTitlB = "Statistics of TS processing time";
      auto canvB  = CanvasConfig(cNameB, cTitlB, 1, 1);
      {
        auto pad           = PadConfig();
        auto name          = "processing_time_evo";
        auto titl          = "Total processing time per TS vs time in run; time in run [s]; Proc. time [ms]";
        int32_t nbBinsX    = std::ceil(1800.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        int32_t nbBinsY    = std::ceil(60000.0 / kNbMsPerBinTime);  // 1 minute max range
        double_t dEndAxisY = nbBinsY * kNbMsPerBinTime;
        fphTotalProcessingTimeEvo =
          fQaData.MakeObj<H2D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX, nbBinsY, 0.0, dEndAxisY);
        pad.RegisterHistogram(fphTotalProcessingTimeEvo, "colz");
        canvB.AddPadConfig(pad);
      }
      fQaData.AddCanvasConfig(canvB);
 
      auto cNameC = "processing_time_fractions";
      auto cTitlC = "Details of TS processing time";
      auto canvC  = CanvasConfig(cNameC, cTitlC, 3, 3);
      {
        auto pad  = PadConfig();
        auto name = "unp_proc_time_frac_evo";
        auto titl = "Fraction of Total proc. time for unpacking per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX        = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX     = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphUnpProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphUnpProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad        = PadConfig();
        auto name       = "sts_proc_time_frac_evo";
        auto titl       = "Fraction of Total proc. time for STS reco per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX     = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphStsProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphStsProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad        = PadConfig();
        auto name       = "tof_proc_time_frac_evo";
        auto titl       = "Fraction of Total proc. time for TOF reco per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX     = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphTofProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphTofProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad        = PadConfig();
        auto name       = "trd_proc_time_frac_evo";
        auto titl       = "Fraction of Total proc. time for TRD reco per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX     = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphTrdProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphTrdProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad              = PadConfig();
        auto name             = "ca_proc_time_frac_evo";
        auto titl             = "Fraction of Total proc. time for CA per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX       = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX    = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphCaProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphCaProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad  = PadConfig();
        auto name = "digmul_proc_time_frac_evo";
        auto titl =
          "Fraction of Total proc. time for Digi Mult. trigger per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX    = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphDigiMulTrigProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphDigiMulTrigProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad  = PadConfig();
        auto name = "hitmul_proc_time_frac_evo";
        auto titl =
          "Fraction of Total proc. time for Hit Mult. trigger per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX               = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX            = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphHitMulTrigProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphHitMulTrigProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad  = PadConfig();
        auto name = "v0_proc_time_frac_evo";
        auto titl = "Fraction of Total proc. time for V0 trigger per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX           = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX        = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphV0TrigProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphV0TrigProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      {
        auto pad  = PadConfig();
        auto name = "evtbuild_proc_time_frac_evo";
        auto titl = "Fraction of Total proc. time for event building per TS vs time in run; time in run [s]; Fract. []";
        int32_t nbBinsX             = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX          = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphEvtBuildProcTimeFractEvo = fQaData.MakeObj<Prof1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphEvtBuildProcTimeFractEvo, "hist");
        canvC.AddPadConfig(pad);
      }
      fQaData.AddCanvasConfig(canvC);
 
 
      auto canvDigisNb = CanvasConfig("digis_nb_evo", "Digis Nb per detector", 3, 3);
      for (auto system : fDetsList) {
        auto pad                 = PadConfig();
        auto name                = "digis_nb_evo_" + system;
        auto titl                = "Digis Nb for " + system + " vs time in run; time in run [s]; Digis []";
        int32_t nbBinsX          = std::ceil(3600.0 / (dSizeTsSec * kNbTsPerBinTime));  // a bit more than 1h range
        double_t dEndAxisX       = nbBinsX * dSizeTsSec * kNbTsPerBinTime;
        fphDigisNbPerDet[system] = fQaData.MakeObj<H1D>(name, titl, nbBinsX, dBegAxisX, dEndAxisX);
        pad.RegisterHistogram(fphDigisNbPerDet[system], "hist");
        canvDigisNb.AddPadConfig(pad);
      }
      fQaData.AddCanvasConfig(canvDigisNb);
 
      fQaData.Init(fpSender);
      fInitNotDone = false;
    }
 
    fphTimeslicesCountEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9);
    fphTimeslicesFractionEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9, 1.0 / kNbTsPerBinFrac);
 
    fphTotalProcessingTimeEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9, qaData.total_time_ms);
    fphUnpProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                 qaData.unpack_time_ms / qaData.total_time_ms);
    fphStsProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                 qaData.sts_reco_time_ms / qaData.total_time_ms);
    fphTofProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                 qaData.tof_reco_time_ms / qaData.total_time_ms);
    fphTrdProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                 qaData.trd_reco_time_ms / qaData.total_time_ms);
    fphCaProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9, qaData.ca_time_ms / qaData.total_time_ms);
 
    fphDigiMulTrigProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                         qaData.digimultrig_time_ms / qaData.total_time_ms);
    fphHitMulTrigProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                        qaData.hitmultrig_time_ms / qaData.total_time_ms);
    fphV0TrigProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                    qaData.v0trig_time_ms / qaData.total_time_ms);
    fphEvtBuildProcTimeFractEvo->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9,
                                      qaData.evtbuild_time_ms / qaData.total_time_ms);
 
    for (auto system : fDetsList) {
      auto it = qaData.digis_nb.find(system);
      if (qaData.digis_nb.end() != it) {
        fphDigisNbPerDet[system]->Fill((ts.start_time() - fRunStartTimeNs) * 1e-9, it->second);
      }
    }
 
    fQaData.Send(fpSender);  // Send and reset!
  }
 
}  // namespace cbm::algo::qa