DQMHistAnalysisSVDDose.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#ifdef _BELLE2_EPICS
#include "cadef.h"
#endif
 
#include <dqm/core/DQMHistAnalysis.h>
#include <TCanvas.h>
#include <TString.h>
#include <TPaveText.h>
#include <TH2.h>
#include <TH1.h>
#include <string>
#include <vector>
 
namespace Belle2 {
  class DQMHistAnalysisSVDDoseModule final : public DQMHistAnalysisModule {
  public:
    DQMHistAnalysisSVDDoseModule();
 
  private:
    typedef struct SensorGroup {
      TString nameSuffix; 
      TString titleSuffix; 
      const char* pvMiddle; 
      int nStrips; 
    } SensorGroup;
 
#ifdef _BELLE2_EPICS
    typedef struct MyPV {
      chid mychid = nullptr; 
      double lastNHits = 0.0; 
      double lastNEvts = 0.0; 
    } MyPV;
#endif

    void initialize() override final;

    void beginRun() override final;

    void event() override final;

    void endRun() override final;
 
    void updateCanvases(); 

    template<typename T>
    inline T* findHistT(TString name) { return dynamic_cast<T*>(findHist(name.Data())); }

    template<typename T>
    static T* divide(T* num, T* den, float scale = 1.0f, T* res = nullptr)
    {
      if (!res) {
        TString name = TString("occu_from_") + num->GetName();
        res = (T*)num->Clone(name);
      }
      for (int i = 0; i < num->GetNcells(); i++) {
        float n = num->GetBinContent(i), d = den->GetBinContent(i), e = num->GetBinError(i);
        res->SetBinContent(i, d ? scale * n / d : 0.0f);
        res->SetBinError(i, d ? scale * e / d : 0.0f);
      }
      return res;
    }

    static void carryOverflowOver(TH1F* h);
 
    // Steerable data members (parameters)
    std::string m_pvPrefix; 
    double m_epicsUpdateSeconds; 
    std::string m_pvSuffix; 
    std::string m_deltaTPVSuffix; 
    std::string m_statePVSuffix; 
 
    // Data members for outputs
    MonitoringObject* m_monObj = nullptr; 
    TPaveText* m_legend = nullptr; 
    // Canvases & output histos for Poisson trigger (TTYP_POIS) events
    std::vector<TCanvas*> m_c_instOccu; 
    std::vector<TCanvas*> m_c_occuLER; 
    std::vector<TH2F*> m_h_occuLER; 
    std::vector<TCanvas*> m_c_occuHER; 
    std::vector<TH2F*> m_h_occuHER; 
    std::vector<TCanvas*> m_c_occuLER1; 
    std::vector<TH1F*> m_h_occuLER1; 
    std::vector<TCanvas*> m_c_occuHER1; 
    std::vector<TH1F*> m_h_occuHER1; 
    // Canvases & output histos for all events
    std::vector<TCanvas*> m_c_instOccuAll; 
    std::vector<TCanvas*> m_c_occuLERAll; 
    std::vector<TH2F*> m_h_occuLERAll; 
    std::vector<TCanvas*> m_c_occuHERAll; 
    std::vector<TH2F*> m_h_occuHERAll; 
    std::vector<TCanvas*> m_c_occuLER1All; 
    std::vector<TH1F*> m_h_occuLER1All; 
    std::vector<TCanvas*> m_c_occuHER1All; 
    std::vector<TH1F*> m_h_occuHER1All; 
 
#ifdef _BELLE2_EPICS
    std::vector<MyPV> m_myPVs; 
    double m_lastPVUpdate = -1.0; 
    chid m_timeSinceLastPVUpdateChan = nullptr; 
    struct dbr_ctrl_enum m_stateCtrl; 
    chid m_stateChan = nullptr; 
#endif

    static const std::vector<SensorGroup> c_sensorGroups;
  };

}