DQMHistAnalysisTOP.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
 
#include <dqm/core/DQMHistAnalysis.h>
 
#include <TCanvas.h>
#include <TH1.h>
#include <TH2.h>
#include <TLine.h>
#include <TPaveText.h>
#include <THStack.h>
#include <TLegend.h>
#include <vector>
#include <string>
 
namespace Belle2 {
  class DQMHistAnalysisTOPModule final : public DQMHistAnalysisModule {
 
  public:
 
    DQMHistAnalysisTOPModule();
 
    ~DQMHistAnalysisTOPModule();
 
    void initialize() override final;
 
    void beginRun() override final;
 
    void event() override final;
 
    void endRun() override final;
 
    void terminate() override final;
 
 
  private:
 
    enum EAlarmStates {
      c_Gray = 0,    // undefined
      c_Green = 1,   // good
      c_Yellow = 2,  // warning
      c_Red = 3      // alarm
    };
 
    void updateWindowVsSlotCanvas();
 
    void updateEventMonitorCanvas();
 
    const TH1F* makeDeadAndHotFractionsPlot();
 
 
    void makePhotonYieldsAndBGRatePlots(const TH1F* activeFraction);
 
    void makeJunkFractionPlot();
 
    void setZAxisRange(const std::string& name, double scale);
 
    void makeBGSubtractedTimimgPlot(const std::string& name);
 
    void setMiraBelleVariables(const std::string& variableName, const TH1* histogram);
 
    int getAlarmState(double value, const std::vector<double>& alarmLevels, bool bigRed = true) const;
 
    int getAlarmColor(unsigned alarmState) const
    {
      if (alarmState < m_alarmColors.size()) return m_alarmColors[alarmState];
      return m_alarmColors[c_Gray];
    }
 
    void setAlarmLines(const std::vector<double>& alarmLevels, double xmin, double xmax, std::vector<TLine*>& alarmLines,
                       bool bigRed = true);
 
    void setAlarmLines();
 
    double getMean(const TH2* h);
 
    void setEpicsVariables();
 
    void updateLimits();
 
    void setIncludedBoardstacks(const std::vector<std::string>& excludedBoardstacks);
 
    // module parameters
 
    std::vector<int> m_asicWindowsBand = {215, 235}; 
    std::vector<double> m_asicWindowsAlarmLevels = {0.002, 0.02}; 
    std::vector<double> m_eventMonitorAlarmLevels = {1e-4, 2e-3}; 
    std::vector<double> m_junkHitsAlarmLevels = {0.05, 0.25}; 
    std::vector<double> m_deadChannelsAlarmLevels = {0.1, 0.35}; 
    std::vector<double> m_backgroundAlarmLevels = {5.0, 10.0}; 
    std::vector<double> m_photonYieldsAlarmLevels = {15.0, 25.0}; 
    std::vector<std::string> m_excludedBoardstacks; 
    std::string  m_pvPrefix;  
    // other
 
    std::vector<int> m_alarmColors = {kGray, kGreen, kYellow, kRed}; 
    std::vector<bool> m_includedBoardstacks; 
    std::map<std::string, int> m_bsmap;  
    bool m_IsNullRun = false; 
    TH1D* m_windowFractions = nullptr; 
    double m_totalWindowFraction = 0;  
    TH1D* m_photonYields = nullptr; 
    TCanvas* m_c_photonYields = nullptr; 
    TH1D* m_backgroundRates = nullptr; 
    TCanvas* m_c_backgroundRates = nullptr; 
    TH1F* m_hotFraction = nullptr; 
    TH1F* m_deadFraction = nullptr; 
    TH1F* m_activeFraction = nullptr; 
    THStack* m_stack = nullptr;  
    TLegend* m_legend = nullptr; 
    TCanvas* m_c_deadAndHot = nullptr; 
    TH1F* m_junkFraction = nullptr; 
    TCanvas* m_c_junkFraction = nullptr; 
    std::vector<TLine*> m_asicWindowsBandLines; 
    std::vector<TLine*> m_verticalLines; 
    std::vector<TLine*> m_junkHitsAlarmLines; 
    std::vector<TLine*> m_deadChannelsAlarmLines; 
    std::vector<TLine*> m_backgroundAlarmLines; 
    std::vector<TLine*> m_photonYieldsAlarmLines; 
    TPaveText* m_text1 = nullptr; 
    TPaveText* m_text2 = nullptr; 
    TPaveText* m_text3 = nullptr; 
    MonitoringObject* m_monObj = nullptr; 
  };
} // end namespace Belle2