DQMHistAnalysisSVDDose.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2021 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Ludovico Massaccesi                                      *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#ifdef _BELLE2_EPICS
#include "cadef.h"
#endif
 
#include <dqm/analysis/modules/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 : public DQMHistAnalysisModule {
  public:
    DQMHistAnalysisSVDDoseModule();
    virtual ~DQMHistAnalysisSVDDoseModule();
 
  private:
    typedef struct SensorGroup {
      TString nameSuffix; 
      TString titleSuffix; 
      const char* pvMiddle; 
      int nStrips; 
    } SensorGroup;
 
#ifdef _BELLE2_EPICS
 
    typedef struct MyPV {
      chid mychid; 
      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)
    {
      TString name = TString("occu_from_") + num->GetName();
      T* 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; 
    bool m_useEpics; 
    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<TCanvas*> m_c_occuHER; 
    std::vector<TCanvas*> m_c_occuLER1; 
    std::vector<TCanvas*> m_c_occuHER1; 
    std::vector<TH2F*> m_h_occuLER; 
    std::vector<TH2F*> m_h_occuHER; 
    std::vector<TH1F*> m_h_occuLER1; 
    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<TCanvas*> m_c_occuHERAll; 
    std::vector<TCanvas*> m_c_occuLER1All; 
    std::vector<TCanvas*> m_c_occuHER1All; 
    std::vector<TH2F*> m_h_occuLERAll; 
    std::vector<TH2F*> m_h_occuHERAll; 
    std::vector<TH1F*> m_h_occuLER1All; 
    std::vector<TH1F*> m_h_occuHER1All; 
#ifdef _BELLE2_EPICS
    std::vector<MyPV> m_myPVs; 
    double m_lastPVUpdate = -1.0; 
    chid m_timeSinceLastPVUpdateChan; 
    struct dbr_ctrl_enum m_stateCtrl; 
    chid m_stateChan; 
#endif
 
    static const std::vector<SensorGroup> c_sensorGroups;
  };
}