DQMHistAnalysisSVDUnpacker.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 : DQMHistAnalysisSVDUnpacker.cc
// Description : module for DQM histogram analysis of SVD unpacker fromat error
//-
 
 
#include <dqm/analysis/modules/DQMHistAnalysisSVDUnpacker.h>
#include <vxd/geometry/GeoCache.h>
 
#include <TROOT.h>
#include <TStyle.h>
#include <TString.h>
#include <TAxis.h>
 
#include <TMath.h>
#include <iostream>
 
using namespace std;
using namespace Belle2;
 
//-----------------------------------------------------------------
//                 Register the Module
//-----------------------------------------------------------------
REG_MODULE(DQMHistAnalysisSVDUnpacker);
 
//-----------------------------------------------------------------
//                 Implementation
//-----------------------------------------------------------------
 
DQMHistAnalysisSVDUnpackerModule::DQMHistAnalysisSVDUnpackerModule()
  : DQMHistAnalysisSVDModule()
{
  //Parameter definition
  B2DEBUG(10, "DQMHistAnalysisSVDUnpacker: Constructor done.");
 
  setDescription("DQM Analysis Module that produces colored canvas for a straightforward interpretation of the SVD Data Quality.");
 
  addParam("unpackerErrorLevel", m_unpackError, "Maximum bin_content/ # events allowed before throwing ERROR", double(0.00001));
  addParam("printCanvas", m_printCanvas, "if True prints pdf of the analysis canvas", bool(false));
  addParam("samples3", m_3Samples, "if True 3 samples histograms analysis is performed", bool(false));
  addParam("PVPrefix", m_pvPrefix, "PV Prefix", std::string("SVD:"));
}
 
DQMHistAnalysisSVDUnpackerModule::~DQMHistAnalysisSVDUnpackerModule() { }
 
void DQMHistAnalysisSVDUnpackerModule::initialize()
{
  B2DEBUG(10, "DQMHistAnalysisSVDUnpacker: initialized.");
 
  m_legProblem->Clear();
  m_legProblem->AddText("ERROR!!");
 
  gROOT->cd();
  m_cUnpacker = new TCanvas("SVDAnalysis/c_SVDDataFormat");
  m_cUnpacker->SetGrid(1);
 
  //register limits for EPICS
  registerEpicsPV(m_pvPrefix + "UnpackError", "UnpackError");
}
 
 
void DQMHistAnalysisSVDUnpackerModule::beginRun()
{
  B2DEBUG(10, "DQMHistAnalysisSVDUnpacker: beginRun called.");
  m_cUnpacker->Clear();
 
  //Retrieve limits from EPICS
 
  double unpackWarnLo = 0.;
  double unpackWarnUp = 0.;
  double unpackErrorLo = 0.;
  requestLimitsFromEpicsPVs("UnpackError", unpackErrorLo, unpackWarnLo, unpackWarnUp,  m_unpackError);
  B2DEBUG(10, " SVD unpack error threshold taken from EPICS configuration file:");
  B2DEBUG(10, "  DATA UNPACK: error > " << m_unpackError);
}
 
void DQMHistAnalysisSVDUnpackerModule::event()
{
  B2DEBUG(10, "DQMHistAnalysisSVDUnpacker: event called.");
 
  //check data errors
  auto hErr = findHist("SVDUnpacker/DQMErrorEventsHisto");
 
  Float_t ratio = 0;
 
  if (hErr != NULL) {
    Float_t events = hErr->GetBinContent(1) + hErr->GetBinContent(2);
    Float_t errors = hErr->GetBinContent(2);
 
    if (events > 0) {
      ratio = errors / events;
      setEpicsPV("UnpackError", ratio);
    }
  } else {
    B2INFO("Histogram SVDUnpacker/DQMErrorEventsHisto from SVDUnpackerDQM not found!");
    return;
  }
 
  TString tmp = hErr->GetTitle();
  Int_t pos = tmp.Last('~');
  Int_t pos1 = tmp.Last('#');
  if (pos == -1)  pos = 0;
  if (pos1 == -1) pos1 = 0;
 
  TString runID = tmp(pos, pos1 - pos);
  B2INFO("DQMHistAnalysisSVDUnpackerModule::runID = " << runID);
 
  //check DATA FORMAT
  auto h = findHist("SVDUnpacker/DQMUnpackerHisto");
 
  if (h != NULL) {
    bool hasError = false;
 
    h->SetTitle(Form("SVD Data Format Monitor %s", runID.Data()));
 
    if (ratio > m_unpackError)
      hasError = true;
 
    if (!hasError) {
      m_cUnpacker->cd();
      h->Draw("colztext");
      h->SetStats(0);
      setStatusOfCanvas(good, m_cUnpacker, false);
    } else {
      m_cUnpacker->cd();
      h->Draw("colztext");
      h->SetStats(0);
      setStatusOfCanvas(error, m_cUnpacker, true);
    }
  } else {
    B2INFO("Histogram SVDUnpacker/DQMUnpackerHisto from SVDUnpackerDQM not found!");
    m_cUnpacker->cd();
    colorizeCanvas(m_cUnpacker, c_StatusDefault);
    setStatusOfCanvas(noStat, m_cUnpacker);
  }
 
  if (m_printCanvas)
    m_cUnpacker->Print("c_SVDDataFormat.pdf");
}
 
void DQMHistAnalysisSVDUnpackerModule::endRun()
{
  B2DEBUG(10, "DQMHistAnalysisSVDUnpacker: endRun called");
}
 
 
void DQMHistAnalysisSVDUnpackerModule::terminate()
{
  B2DEBUG(10, "DQMHistAnalysisSVDUnpacker: terminate called");
 
  delete m_cUnpacker;
}