DQMHistAnalysisDeltaEpicsMonObjExample.cc

/**************************************************************************
 * 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.                  *
 **************************************************************************/
//+
// File : DQMHistAnalysisDeltaEpicsMonObjExample.cc
// Description : Test Module for Delta Histogram Access
//-
 
 
#include <dqm/analysis/modules/DQMHistAnalysisDeltaEpicsMonObjExample.h>
#include <TROOT.h>
 
using namespace Belle2;
 
//-----------------------------------------------------------------
//                 Register the Module
//-----------------------------------------------------------------
REG_MODULE(DQMHistAnalysisDeltaEpicsMonObjExample);
 
//-----------------------------------------------------------------
//                 Implementation
//-----------------------------------------------------------------
 
DQMHistAnalysisDeltaEpicsMonObjExampleModule::DQMHistAnalysisDeltaEpicsMonObjExampleModule()
  : DQMHistAnalysisModule()
{
  // This module CAN NOT be run in parallel!
  setDescription("Example module for using delta histogramming and sending values to EPICS");
 
  // Parameter definition
  addParam("histogramDirectoryName", m_histogramDirectoryName, "Name of Histogram dir", std::string("test"));
  addParam("histogramName", m_histogramName, "Name of Histogram", std::string("testHist"));
  addParam("PVPrefix", m_pvPrefix, "PV Prefix", std::string("DQM:TEST:"));
  B2DEBUG(1, "DQMHistAnalysisDeltaEpicsMonObjExample: Constructor done.");
 
}
 
void DQMHistAnalysisDeltaEpicsMonObjExampleModule::initialize()
{
  B2DEBUG(1, "DQMHistAnalysisDeltaEpicsMonObjExample: initialized.");
 
  gROOT->cd(); // this seems to be important before creating new Canvas, or strange things happen
  m_canvas = new TCanvas((m_histogramDirectoryName + "/c_Test").data());
 
  m_monObj = getMonitoringObject("test");
  m_monObj->addCanvas(m_canvas);
 
  // delta parameters, two examples:
  addDeltaPar(m_histogramDirectoryName, m_histogramName, HistDelta::c_Events,
              100000, 1); // update each 100k events (from daq histogram)
  // addDeltaPar(m_histogramDirectoryName, m_histogramName, HistDelta::c_Entries, 10000, 1); // update each 10000 entries
 
  registerEpicsPV(m_pvPrefix + "mean", "mean");
  registerEpicsPV(m_pvPrefix + "width", "width");
  registerEpicsPV(m_pvPrefix + "alarm", "alarm");
 
  m_meanLowerAlarmLine = new TLine();
  m_meanLowerAlarmLine->SetLineColor(kRed + 3);
  m_meanLowerAlarmLine->SetLineWidth(3);
  m_meanLowerAlarmLine->SetVertical(true);
  m_meanLowerWarnLine = new TLine();
  m_meanLowerWarnLine->SetLineColor(kOrange - 3);
  m_meanLowerWarnLine->SetLineWidth(3);
  m_meanLowerWarnLine->SetVertical(true);
  m_meanUpperWarnLine = new TLine();
  m_meanUpperWarnLine->SetLineColor(kOrange - 3);
  m_meanUpperWarnLine->SetLineWidth(3);
  m_meanUpperWarnLine->SetVertical(true);
  m_meanUpperAlarmLine = new TLine();
  m_meanUpperAlarmLine->SetLineColor(kRed + 3);
  m_meanUpperAlarmLine->SetLineWidth(3);
  m_meanUpperAlarmLine->SetVertical(true);
}
 
void DQMHistAnalysisDeltaEpicsMonObjExampleModule::beginRun()
{
  B2DEBUG(1, "DQMHistAnalysisDeltaEpicsMonObjExample: beginRun called.");
 
  double unused = NAN;
  m_meanLowerAlarm = -1.0; // this is some way to set default values, which are used
  m_meanLowerWarn = -0.9; // when none are set in epics PV directly.
  m_meanUpperWarn = 0.9;
  m_meanUpperAlarm = 1.0;
  requestLimitsFromEpicsPVs("mean", m_meanLowerAlarm, m_meanLowerWarn, m_meanUpperWarn, m_meanUpperAlarm);
 
  if (!std::isnan(m_meanLowerAlarm)) {
    m_meanLowerAlarmLine->SetX1(m_meanLowerAlarm);
    m_meanLowerAlarmLine->SetX2(m_meanLowerAlarm);
  }
  if (!std::isnan(m_meanLowerWarn)) {
    m_meanLowerWarnLine->SetX1(m_meanLowerWarn);
    m_meanLowerWarnLine->SetX2(m_meanLowerWarn);
  }
  if (!std::isnan(m_meanUpperWarn)) {
    m_meanUpperWarnLine->SetX1(m_meanUpperWarn);
    m_meanUpperWarnLine->SetX2(m_meanUpperWarn);
  }
  if (!std::isnan(m_meanUpperAlarm)) {
    m_meanUpperAlarmLine->SetX1(m_meanUpperAlarm);
    m_meanUpperAlarmLine->SetX2(m_meanUpperAlarm);
  }
 
  requestLimitsFromEpicsPVs("width", unused, unused, m_widthUpperWarn, m_widthUpperAlarm);
}
 
void DQMHistAnalysisDeltaEpicsMonObjExampleModule::endRun()
{
  B2DEBUG(1, "DQMHistAnalysisDeltaEpicsMonObjExample: endRun called.");
  // MiraBelle export code should run at end of Run
  // but it still "remembers" the state from last event call.
  doHistAnalysis(true);
}
 
void DQMHistAnalysisDeltaEpicsMonObjExampleModule::event()
{
  B2DEBUG(1, "DQMHistAnalysisDeltaEpicsMonObjExample: event called.");
  doHistAnalysis(false);
}
 
void DQMHistAnalysisDeltaEpicsMonObjExampleModule::doHistAnalysis(bool forMiraBelle)
{
  double data_mean = 0.0;
  double data_width = 0.0;
  EStatus status = c_StatusDefault;
 
  m_canvas->Clear();
  m_canvas->cd(0);
 
  // more handy to have it as a full name, but find/get functions
  // can work with one or two parameters
  // std::string fullname = m_histogramDirectoryName + "/" + m_histogramName;
 
  // In the following, enable one
 
  // get most recent delta, but only if updated since last call, nullptr otherwise
  auto hist =  getDelta(m_histogramDirectoryName, m_histogramName);
  // get basic histogram (run integrated up) but only if updated since last call
  // auto hist = findHist(m_histogramDirectoryName, m_histogramName, true);// only if updated
 
  // the following cases do not make sense, unless you want to achieve something special.
  // get most recent delta even if not updated
  // auto hist =  getDelta(m_histogramDirectoryName, m_histogramName,false);// even if no update
  // get basic histogram (run integrated up) even if not updated
  // auto hist = findHist(m_histogramDirectoryName, m_histogramName, false);// even if no update
 
  if (hist != nullptr) {
    data_mean = hist->GetMean();
    data_width = hist->GetRMS();
 
    // your conditions for warn and alarm
    // here we would have to check against NaN, too ... just too lazy to write the code in this example
    // maybe other condition, e.g. too low to judge -> 0
    status = makeStatus(true,
                        data_mean > m_meanUpperWarn || data_mean < m_meanLowerWarn || data_width > m_widthUpperWarn,
                        data_mean > m_meanUpperAlarm || data_mean < m_meanLowerAlarm || data_width > m_widthUpperAlarm);
 
    hist->Draw("hist");
 
    // line position need to be adjusted to histogram content
    double y1 = hist->GetMaximum();
    double y2 = hist->GetMinimum();
    auto ym = y1 + (y1 - y2) * 0.05;
    m_meanLowerAlarmLine->SetY1(ym);
    m_meanLowerWarnLine->SetY1(ym);
    m_meanUpperWarnLine->SetY1(ym);
    m_meanUpperAlarmLine->SetY1(ym);
    auto yn = y2 + (y1 - y2) * 0.05;
    m_meanLowerAlarmLine->SetY2(yn);
    m_meanLowerWarnLine->SetY2(yn);
    m_meanUpperWarnLine->SetY2(yn);
    m_meanUpperAlarmLine->SetY2(yn);
 
    if (!std::isnan(m_meanLowerAlarm)) m_meanLowerAlarmLine->Draw();
    if (!std::isnan(m_meanLowerWarn)) m_meanLowerWarnLine->Draw();
    if (!std::isnan(m_meanUpperWarn)) m_meanUpperWarnLine->Draw();
    if (!std::isnan(m_meanUpperAlarm)) m_meanUpperAlarmLine->Draw();
 
    colorizeCanvas(m_canvas, status);
  }
 
  // Tag canvas as updated ONLY if things have changed.
  UpdateCanvas(m_canvas->GetName(), hist != nullptr);
  // Remark: if you do not tag, you may have an empty Canvas (as it is cleared above)
 
  // this if left over from jsroot, may not be needed anymore (to check)
  m_canvas->Update();
 
  // actually, we would prefer to only update the epics variable
  // if the main histogram or delta~ has changed.
 
  if (hist) {
    if (forMiraBelle) {
      // currently, monObj are only used at end of run, thus this is not necessary
      // but, if in the future we want to use finer granularity, this code needs change.
      m_monObj->setVariable("mean", data_mean);
      m_monObj->setVariable("width", data_width);
    } else {
      // we do not want to fill epics again at end of run, as this would be identical to the last event called before.
      setEpicsPV("mean", data_mean);
      setEpicsPV("width", data_width);
      setEpicsPV("alarm", int(status));
      // until now, changes are buffered only locally, not waiting for them being written out to performance reason
      // framework will take care after all analysis modules are done
    }
  }
}
 
void DQMHistAnalysisDeltaEpicsMonObjExampleModule::terminate()
{
  B2DEBUG(1, "DQMHistAnalysisDeltaEpicsMonObjExample: terminate called");
}