RecoGeneralQa.h

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] */
 
#ifndef ALGO_QA_RECOGENERALQA_H
#define ALGO_QA_RECOGENERALQA_H 1
 
#include "CbmDefs.h"
#include "algo/base/HistogramSender.h"
#include "algo/qa/QaData.h"
 
#include <StorableTimeslice.hpp>
 
namespace cbm::algo::qa
{
  struct GeneralQaData {
    double total_time_ms       = 0.0;  // Total Processing time in ms
    double unpack_time_ms      = 0.0;  // Total unpacking time in ms
    double sts_reco_time_ms    = 0.0;  // Total STS reco time in ms
    double tof_reco_time_ms    = 0.0;  // Total TOF reco time in ms
    double trd_reco_time_ms    = 0.0;  // Total TRD reco time in ms
    double ca_time_ms          = 0.0;  // Total CA tracking time in ms
    double digimultrig_time_ms = 0.0;  // Total time used for the digi multiplicity trigger in ms
    double hitmultrig_time_ms  = 0.0;  // Total time used for the hit multiplicity trigger in ms
    double v0trig_time_ms      = 0.0;  // Total time used for the V0 trigger in ms
    double evtbuild_time_ms    = 0.0;  // Total time used for the event building in ms
 
    std::map<std::string, size_t> digis_nb = {};  // total size of timeslice contents for each detector
 
    void Reset()
    {
      total_time_ms       = 0.0;  // Total Processing time in ms
      unpack_time_ms      = 0.0;  // Total unpacking time in ms
      sts_reco_time_ms    = 0.0;  // Total STS reco time in ms
      tof_reco_time_ms    = 0.0;  // Total TOF reco time in ms
      trd_reco_time_ms    = 0.0;  // Total TRD reco time in ms
      ca_time_ms          = 0.0;  // Total CA tracking time in ms
      digimultrig_time_ms = 0.0;  // Total time used for the digi multiplicity trigger in ms
      hitmultrig_time_ms  = 0.0;  // Total time used for the hit multiplicity trigger in ms
      v0trig_time_ms      = 0.0;  // Total time used for the V0 trigger in ms
      evtbuild_time_ms    = 0.0;  // Total time used for the event building in ms
 
      digis_nb.clear();  // total size of timeslice contents for each detector
    }
  };

  class RecoGeneralQa {
   public:
    RecoGeneralQa(const uint64_t& runStartTimeNs, std::shared_ptr<HistogramSender> pSender);

    RecoGeneralQa() = delete;

    RecoGeneralQa(const RecoGeneralQa&) = delete;

    RecoGeneralQa(RecoGeneralQa&&) = delete;

    RecoGeneralQa& operator=(const RecoGeneralQa&) = delete;

    RecoGeneralQa& operator=(RecoGeneralQa&&) = delete;

    void operator()(const fles::Timeslice& ts, const GeneralQaData& qaData);
 
   private:  // methods
   private:  // members
    uint64_t fRunStartTimeNs;
    std::shared_ptr<HistogramSender> fpSender = nullptr;  
    qa::Data fQaData{"Reco"};                             
    bool fInitNotDone = true;
 
    // Constants
    static const int32_t kNbTsPerBinCount = 100;   // 100 TS duration per bin for raw counts
    static const int32_t kNbTsPerBinFrac  = 1000;  // 1000 TS duration per bin for fractions
    static const int32_t kNbTsPerBinTime  = 8;     // 8 TS duration per bin for processing time evo
    static const int32_t kNbMsPerBinTime  = 20;    // 10 ms per bin for processing time total
 
    // ---- Histograms
    qa::H1D* fphTimeslicesCountEvo     = nullptr;  
    qa::H1D* fphTimeslicesFractionEvo  = nullptr;  
    qa::H2D* fphTotalProcessingTimeEvo = nullptr;  
 
    qa::Prof1D* fphUnpProcTimeFractEvo = nullptr;  
    qa::Prof1D* fphStsProcTimeFractEvo = nullptr;  
    qa::Prof1D* fphTofProcTimeFractEvo = nullptr;  
    qa::Prof1D* fphTrdProcTimeFractEvo = nullptr;  
    qa::Prof1D* fphCaProcTimeFractEvo  = nullptr;  
 
    qa::Prof1D* fphDigiMulTrigProcTimeFractEvo = nullptr;  
    qa::Prof1D* fphHitMulTrigProcTimeFractEvo  = nullptr;  
    qa::Prof1D* fphV0TrigProcTimeFractEvo      = nullptr;  
    qa::Prof1D* fphEvtBuildProcTimeFractEvo    = nullptr;  
 
    // FIXME: replace with standard CBM enumerator? but maybe then too many plots from deprecated detectors
    std::vector<std::string> fDetsList = {"BMON", "MVD", "STS", "RICH", "MUCH", "TRD", "TRD2D", "TOF", "FSD"};
    std::map<std::string, qa::H1D*> fphDigisNbPerDet = {};
  };
}  // namespace cbm::algo::qa
 
#endif /* ALGO_QA_RECOGENERALQA_H */