DQMHistAnalysisSVDOccupancy.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 : DQMHistAnalysisSVDOccupancy.cc
// Description : module for DQM histogram analysis of SVD sensors occupancies
//-
 
 
#include <dqm/analysis/modules/DQMHistAnalysisSVDOccupancy.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(DQMHistAnalysisSVDOccupancy);
 
//-----------------------------------------------------------------
//                 Implementation
//-----------------------------------------------------------------
 
DQMHistAnalysisSVDOccupancyModule::DQMHistAnalysisSVDOccupancyModule()
  : DQMHistAnalysisSVDModule(false, true, true)
{
  //Parameter definition
  B2DEBUG(10, "DQMHistAnalysisSVDOccupancy: Constructor done.");
 
  setDescription("DQM Analysis Module that produces colored canvas for a straightforward interpretation of the SVD Data Quality.");
 
  addParam("occLevel_Error", m_occError, "Maximum Occupancy (%) allowed for safe operations (red)", double(3));
  addParam("occLevel_Warning", m_occWarning, "Occupancy (%) at WARNING level (orange)", double(2.5));
  addParam("occGroupIDsLevel_Error", m_groupIDsOccError, "Maximum group IDs Occupancy (%) allowed for safe operations (red)",
           double(2));
  addParam("occGroupIDsLevel_Warning", m_groupIDsOccWarning, "Group IDs Occupancy (%) at WARNING level (orange)", double(1.5));
  addParam("occLevel_Empty", m_occEmpty, "Maximum Occupancy (%) for which the sensor is considered empty", double(0));
  addParam("onlineOccLevel_Error", m_onlineOccError, "Maximum OnlineOccupancy (%) allowed for safe operations (red)", double(10));
  addParam("onlineOccLevel_Warning", m_onlineOccWarning, "OnlineOccupancy (%) at WARNING level (orange)", double(5));
  addParam("onlineOccLevel_Empty", m_onlineOccEmpty, "Maximum OnlineOccupancy (%) for which the sensor is considered empty",
           double(0));
  addParam("printCanvas", m_printCanvas, "if True prints pdf of the analysis canvas", bool(false));
  addParam("additionalPlots", m_additionalPlots, "Flag to produce additional plots",   bool(false));
  addParam("RPhiView", m_RPhiView, "Flag to produce RPhi view plots",   bool(true));
  addParam("RPhiViewId0", m_RPhiViewId0, "Flag to produce  RPhi view plots for Id0 group",   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:"));
}
 
DQMHistAnalysisSVDOccupancyModule::~DQMHistAnalysisSVDOccupancyModule() { }
 
void DQMHistAnalysisSVDOccupancyModule::initialize()
{
  B2DEBUG(10, "DQMHistAnalysisSVDOccupancy: initialized.");
 
  const VXD::GeoCache& geo = VXD::GeoCache::getInstance();
 
  //collect the list of all SVD Modules in the geometry here
  std::vector<VxdID> sensors = geo.getListOfSensors();
  for (const auto& aVxdID : sensors) {
    VXD::SensorInfoBase info = geo.getSensorInfo(aVxdID);
    // B2INFO("VXD " << aVxdID);
    if (info.getType() != VXD::SensorInfoBase::SVD) continue;
    m_SVDModules.push_back(aVxdID); // reorder, sort would be better
  }
  std::sort(m_SVDModules.begin(), m_SVDModules.end());  // back to natural order
 
  //occupancy chart chip
  m_cOccupancyChartChip = new TCanvas("SVDOccupancy/c_OccupancyChartChip");
 
  //strip occupancy per sensor
  if (m_additionalPlots)
    m_sensors = m_SVDModules.size();
  else
    m_sensors = 2;
 
  m_cStripOccupancyU = new TCanvas*[m_sensors];
  m_cStripOccupancyV = new TCanvas*[m_sensors];
  for (int i = 0; i < m_sensors; i++) {
    int tmp_layer = m_SVDModules[i].getLayerNumber();
    int tmp_ladder = m_SVDModules[i].getLadderNumber();
    int tmp_sensor = m_SVDModules[i].getSensorNumber();
    m_cStripOccupancyU[i] = new TCanvas(Form("SVDOccupancy/c_StripOccupancyU_%d_%d_%d", tmp_layer, tmp_ladder, tmp_sensor));
    m_cStripOccupancyV[i] = new TCanvas(Form("SVDOccupancy/c_StripOccupancyV_%d_%d_%d", tmp_layer, tmp_ladder, tmp_sensor));
  }
 
  gROOT->cd();
  m_cOccupancyU = new TCanvas("SVDAnalysis/c_SVDOccupancyU");
  m_cOccupancyV = new TCanvas("SVDAnalysis/c_SVDOccupancyV");
 
  m_cOnlineOccupancyU = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancyU");
  m_cOnlineOccupancyV = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancyV");
 
  m_cOccupancyUGroupId0 = new TCanvas("SVDAnalysis/c_SVDOccupancyUGroupId0");
  m_cOccupancyVGroupId0 = new TCanvas("SVDAnalysis/c_SVDOccupancyVGroupId0");
 
  if (m_RPhiView) {
    m_cOccupancyRPhiViewU = new TCanvas("SVDAnalysis/c_SVDOccupancyRPhiViewU", "", 800, 800);
    m_cOccupancyRPhiViewV = new TCanvas("SVDAnalysis/c_SVDOccupancyRPhiViewV", "", 800, 800);
 
    m_cOnlineOccupancyRPhiViewU = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancyRPhiViewU", "", 800, 800);
    m_cOnlineOccupancyRPhiViewV = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancyRPhiViewV", "", 800, 800);
 
    if (m_RPhiViewId0) {
      m_cOccupancyRPhiViewUGroupId0 = new TCanvas("SVDAnalysis/c_SVDOccupancyRPhiViewGroupId0U", "", 800, 800);
      m_cOccupancyRPhiViewVGroupId0 = new TCanvas("SVDAnalysis/c_SVDOccupancyRPhiViewGroupId0V", "", 800, 800);
    }
  }
 
  TString hName = getHistoNameFromCanvas(m_cOccupancyU->GetName(), "@view");
  m_hOccupancy =  new SVDSummaryPlots(hName.Data(), "Average OFFLINE Sensor Occupancy (%), @view/@side Side");
  m_hOccupancy->setStats(0);
 
  hName = getHistoNameFromCanvas(m_cOnlineOccupancyU->GetName(), "@view");
  m_hOnlineOccupancy = new SVDSummaryPlots(hName.Data(), "Average ONLINE Sensor Occupancy (%), @view/@side Side");
  m_hOnlineOccupancy->setStats(0);
 
  // change name by hand cos side index not at the end
  hName = getHistoNameFromCanvas("SVDAnalysis/c_SVDOccupancyGroupId0U", "@view");
  m_hOccupancyGroupId0 = new SVDSummaryPlots(hName.Data(),
                                             "Average OFFLINE Sensor Occupancy (%), @view/@side Side for cluster time group Id = 0, 1, 2 & 3");
  m_hOccupancyGroupId0->setStats(0);
 
  if (m_3Samples) {
    m_cOccupancyU3Samples = new TCanvas("SVDAnalysis/c_SVDOccupancy3SamplesU");
    m_cOccupancyV3Samples = new TCanvas("SVDAnalysis/c_SVDOccupancy3SamplesV");
 
    m_cOnlineOccupancyU3Samples = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancy3SamplesU");
    m_cOnlineOccupancyV3Samples = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancy3SamplesV");
 
    if (m_RPhiView) {
      m_cOccupancyRPhiViewU3Samples = new TCanvas("SVDAnalysis/c_SVDOccupancyRPhiView3SamplesU", "", 800, 800);
      m_cOccupancyRPhiViewV3Samples = new TCanvas("SVDAnalysis/c_SVDOccupancyRPhiView3SamplesV", "", 800, 800);
 
      m_cOnlineOccupancyRPhiViewU3Samples = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancyRPhiView3SamplesU", "", 800, 800);
      m_cOnlineOccupancyRPhiViewV3Samples = new TCanvas("SVDAnalysis/c_SVDOnlineOccupancyRPhiView3SamplesV", "", 800, 800);
    }
 
    hName = getHistoNameFromCanvas(m_cOccupancyU3Samples->GetName(), "@view");
    m_hOccupancy3Samples = new SVDSummaryPlots(hName.Data(), "Average OFFLINE Sensor Occupancy (%), @view/@side Side for 3 samples");
    m_hOnlineOccupancy->setStats(0);
 
    hName = getHistoNameFromCanvas(m_cOnlineOccupancyU3Samples->GetName(), "@view");
    m_hOnlineOccupancy3Samples = new SVDSummaryPlots(hName.Data(),
                                                     "Average ONLINE Sensor Occupancy (%), @view/@side Side for 3 samples");
    m_hOnlineOccupancy3Samples->setStats(0);
  }
 
 
  //register limits for EPICS
  registerEpicsPV(m_pvPrefix + "occupancyLimits", "occLimits");
  registerEpicsPV(m_pvPrefix + "occupancyOnlineLimits", "occOnlineLimits");
  registerEpicsPV(m_pvPrefix + "occupancyGoupIDsLimits", "occGroupIDsLimits");
}
 
void DQMHistAnalysisSVDOccupancyModule::beginRun()
{
  B2DEBUG(10, "DQMHistAnalysisSVDOccupancy: beginRun called.");
 
  m_cOccupancyChartChip->Clear();
  for (int i = 0; i < m_sensors; i++) {
    m_cStripOccupancyU[i]->Clear();
    m_cStripOccupancyV[i]->Clear();
  }
 
  // histo
  if (m_cOccupancyU)
    m_cOccupancyU->Clear();
  if (m_cOccupancyV)
    m_cOccupancyV->Clear();
 
  if (m_cOnlineOccupancyU)
    m_cOnlineOccupancyU->Clear();
  if (m_cOnlineOccupancyV)
    m_cOnlineOccupancyV->Clear();
 
  if (m_cOccupancyUGroupId0)
    m_cOccupancyUGroupId0->Clear();
  if (m_cOccupancyVGroupId0)
    m_cOccupancyVGroupId0->Clear();
 
  // RPhiView
  if (m_RPhiView) {
    if (m_cOccupancyRPhiViewU)
      m_cOccupancyRPhiViewU->Clear();
    if (m_cOccupancyRPhiViewV)
      m_cOccupancyRPhiViewV->Clear();
 
    if (m_cOnlineOccupancyRPhiViewU)
      m_cOnlineOccupancyRPhiViewU->Clear();
    if (m_cOnlineOccupancyRPhiViewV)
      m_cOnlineOccupancyRPhiViewV->Clear();
    if (m_RPhiViewId0) {
      if (m_cOccupancyRPhiViewUGroupId0)
        m_cOccupancyRPhiViewUGroupId0->Clear();
      if (m_cOccupancyRPhiViewVGroupId0)
        m_cOccupancyRPhiViewVGroupId0->Clear();
    }
  }
  // 3 samples
  if (m_3Samples) {
    m_cOccupancyU3Samples->Clear();
    m_cOccupancyV3Samples->Clear();
    m_cOnlineOccupancyU3Samples->Clear();
    m_cOnlineOccupancyV3Samples->Clear();
 
    if (m_RPhiView) {
      m_cOccupancyRPhiViewU3Samples->Clear();
      m_cOccupancyRPhiViewV3Samples->Clear();
      m_cOnlineOccupancyRPhiViewU3Samples->Clear();
      m_cOnlineOccupancyRPhiViewV3Samples->Clear();
    }
  }
 
  //Retrieve limits from EPICS
  double oocErrorLoOff = 0.;
  double oocErrorLoOn = 0.;
  double oocErrorLoIDs = 0.;
 
  double occWarningOff = 0.;
  double occWarningOn = 0.;
  double occWarningIDs = 0.;
 
  requestLimitsFromEpicsPVs("occLimits", oocErrorLoOff, occWarningOff, m_occWarning,  m_occError);
  requestLimitsFromEpicsPVs("occOnlineLimits", oocErrorLoOn, occWarningOn, m_onlineOccWarning,  m_onlineOccError);
  requestLimitsFromEpicsPVs("occGroupIDsLimits", oocErrorLoIDs, occWarningIDs, m_groupIDsOccWarning,  m_groupIDsOccError);
 
  B2DEBUG(10, " SVD occupancy thresholds taken from EPICS configuration file:");
  B2DEBUG(10, "  ONLINE OCCUPANCY: empty < " << m_onlineOccEmpty << " normal < " << m_onlineOccWarning << " warning < " <<
          m_onlineOccError <<
          " < error");
  B2DEBUG(10, "  OFFLINE OCCUPANCY: empty < " << m_occEmpty << " normal < " << m_occWarning << " warning < " << m_occError <<
          " < error with minimum statistics of " << m_occEmpty);
 
  B2DEBUG(10, "  Group IDs OCCUPANCY: empty < " << m_occEmpty << " normal < " << m_groupIDsOccWarning << " warning < " <<
          m_groupIDsOccError <<
          " < error with minimum statistics of " << m_occEmpty);
 
  // Create text panel
  //OFFLINE occupancy plots legend
  m_legProblem->Clear();
  m_legProblem->AddText("ERROR!");
  m_legProblem->AddText("at least one sensor with:");
  m_legProblem->AddText(Form("occupancy > %1.1f%%", m_occError));
 
  m_legWarning->Clear();
  m_legWarning->AddText("WARNING!");
  m_legWarning->AddText("at least one sensor with:");
  m_legWarning->AddText(Form("%1.1f%% < occupancy < %1.1f%%", m_occWarning, m_occError));
 
  m_legNormal->Clear();
  m_legNormal->AddText("OCCUPANCY WITHIN LIMITS");
  m_legNormal->AddText(Form("%1.1f%% < occupancy < %1.1f%%", m_occEmpty, m_occWarning));
 
  m_legLowStat->Clear();
  m_legLowStat->AddText("LOW STATISTICS");
  m_legLowStat->AddText("from at least one sensor");
 
  m_legEmpty->Clear();
  m_legEmpty->AddText("at least one sensor is emtpy");
 
  //ONLINE occupancy plots legend
  m_legOnlineProblem->Clear();
  m_legOnlineProblem->AddText("ERROR!");
  m_legOnlineProblem->AddText("at least one sensor with:");
  m_legOnlineProblem->AddText(Form("online occupancy > %1.1f%%", m_onlineOccError));
 
  m_legOnlineWarning->Clear();
  m_legOnlineWarning->AddText("WARNING!");
  m_legOnlineWarning->AddText("at least one sensor with:");
  m_legOnlineWarning->AddText(Form("%1.1f%% < online occupancy < %1.1f%%", m_onlineOccWarning, m_onlineOccError));
 
  m_legOnlineNormal->Clear();
  m_legOnlineNormal->AddText("OCCUPANCY WITHIN LIMITS");
  m_legOnlineNormal->AddText(Form("%1.1f%% < online occupancy < %1.1f%%", m_onlineOccEmpty, m_onlineOccWarning));
 
  //Group IDs occupancy plots legend
  m_legGroupIDsProblem->Clear();
  m_legGroupIDsProblem->AddText("ERROR!");
  m_legGroupIDsProblem->AddText("at least one sensor with:");
  m_legGroupIDsProblem->AddText(Form("group IDs occupancy > %1.1f%%", m_groupIDsOccError));
 
  m_legGroupIDsWarning->Clear();
  m_legGroupIDsWarning->AddText("WARNING!");
  m_legGroupIDsWarning->AddText("at least one sensor with:");
  m_legGroupIDsWarning->AddText(Form("%1.1f%% < group IDs occupancy < %1.1f%%", m_groupIDsOccWarning, m_groupIDsOccError));
 
  m_legGroupIDsNormal->Clear();
  m_legGroupIDsNormal->AddText("OCCUPANCY WITHIN LIMITS");
  m_legGroupIDsNormal->AddText(Form("%1.1f%% < group IDs occupancy < %1.1f%%", m_groupIDsOccEmpty, m_groupIDsOccWarning));
 
 
  m_occUstatus = good;
  m_occVstatus = good;
 
  m_occU3Samples = good;
  m_occV3Samples = good;
 
  m_occUGroupId0 = good;
  m_occVGroupId0 = good;
 
  m_onlineOccUstatus = good;
  m_onlineOccVstatus = good;
 
  m_onlineOccU3Samples = good;
  m_onlineOccV3Samples = good;
}
 
void DQMHistAnalysisSVDOccupancyModule::event()
{
  B2DEBUG(10, "DQMHistAnalysisSVDOccupancy: event called.");
 
  //find nEvents
  TH1* hnEvnts = findHist("SVDExpReco/SVDDQM_nEvents", true);
  if (hnEvnts == NULL) {
    B2INFO("no events, nothing to do here");
    return;
  } else {
    B2DEBUG(10, "SVDExpReco/SVDDQM_nEvents found");
  }
 
 
  TString tmp = hnEvnts->GetTitle();
  Int_t pos = tmp.Last('~');
  if (pos == -1) pos = 0;
 
  TString runID = tmp(pos, tmp.Length() - pos);
  B2INFO("DQMHistAnalysisSVDOccupancyModule::runID = " << runID);
  Float_t nEvents = hnEvnts->GetEntries();
 
  //occupancy chart
  auto hChart = findHist("SVDExpReco/SVDDQM_StripCountsChip");
 
  if (hChart != NULL) {
    m_hOccupancyChartChip.Clear();
    hChart->Copy(m_hOccupancyChartChip);
    m_hOccupancyChartChip.SetName("SVDOccupancyChart");
    m_hOccupancyChartChip.SetTitle(Form("SVD OFFLINE Occupancy per chip %s", runID.Data()));
    m_hOccupancyChartChip.Scale(1 / nEvents / 128);
    m_cOccupancyChartChip->cd();
    //    m_hOccupancyChartChip->SetStats(0);
    m_hOccupancyChartChip.Draw();
  }
  m_cOccupancyChartChip->Modified();
  m_cOccupancyChartChip->Update();
 
  if (m_printCanvas)
    m_cOccupancyChartChip->Print("c_OccupancyChartChip.pdf");
 
 
  //check MODULE OCCUPANCY online & offline
  //reset canvas color
  //update titles with exp and run number
  m_occUstatus = good;
  m_occVstatus = good;
 
  m_occU3Samples = good;
  m_occV3Samples = good;
 
  m_occUGroupId0 = good;
  m_occVGroupId0 = good;
 
  m_onlineOccUstatus = good;
  m_onlineOccVstatus = good;
 
  m_onlineOccU3Samples = good;
  m_onlineOccV3Samples = good;
 
  m_hOccupancy->reset();
  m_hOccupancy->setStats(0);
  m_hOccupancy->setRunID(runID);
 
  m_hOnlineOccupancy->reset();
  m_hOnlineOccupancy->setStats(0);
  m_hOnlineOccupancy->setRunID(runID);
 
  m_hOccupancyGroupId0->reset();
  m_hOccupancyGroupId0->setStats(0);
  m_hOccupancyGroupId0->setRunID(runID);
 
  if (m_3Samples) {
    m_hOccupancy3Samples->reset();
    m_hOccupancy3Samples->setStats(0);
    m_hOccupancy3Samples->setRunID(runID);
 
    m_hOnlineOccupancy3Samples->reset();
    m_hOnlineOccupancy3Samples->setStats(0);
    m_hOnlineOccupancy3Samples->setRunID(runID);
  }
 
  //set dedicate gStyle
  gStyle->SetOptStat(0);
  gStyle->SetPaintTextFormat("2.3f");
 
 
  for (unsigned int i = 0; i < m_SVDModules.size(); i++) {
    int tmp_layer = m_SVDModules[i].getLayerNumber();
    int tmp_ladder = m_SVDModules[i].getLadderNumber();
    int tmp_sensor = m_SVDModules[i].getSensorNumber();
 
    //look for U histogram - OFFLINE ZS
    TString tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_StripCountU", tmp_layer, tmp_ladder, tmp_sensor);
 
    auto hStripCountU = findHist(tmpname.Data());
    if (hStripCountU == NULL) {
      B2INFO("Occupancy U histogram not found");
      setOccStatus(-1, m_occUstatus);
    } else {
      Float_t occU = getOccupancy(hStripCountU->GetEntries(), tmp_layer, nEvents);
      m_hOccupancy->fill(m_SVDModules[i], 1, occU);
      setOccStatus(occU, m_occUstatus);
 
      //produce the occupancy plot
      if (m_additionalPlots) {
        m_hStripOccupancyU[i].Clear();
        hStripCountU->Copy(m_hStripOccupancyU[i]);
        m_hStripOccupancyU[i].Scale(1 / nEvents);
        m_hStripOccupancyU[i].SetName(Form("%d_%d_%d_OccupancyU", tmp_layer, tmp_ladder, tmp_sensor));
        m_hStripOccupancyU[i].SetTitle(Form("SVD Sensor %d_%d_%d U-Strip OFFLINE Occupancy vs Strip Number %s", tmp_layer, tmp_ladder,
                                            tmp_sensor, runID.Data()));
      } else {
        if (i == 0 || i == 1) {
          m_hStripOccupancyU[i].Clear();
          hStripCountU->Copy(m_hStripOccupancyU[i]);
          m_hStripOccupancyU[i].Scale(1 / nEvents);
          m_hStripOccupancyU[i].SetName(Form("%d_%d_%d_OccupancyU", tmp_layer, tmp_ladder, tmp_sensor));
          m_hStripOccupancyU[i].SetTitle(Form("SVD Sensor %d_%d_%d U-Strip OFFLINE Occupancy vs Strip Number %s", tmp_layer, tmp_ladder,
                                              tmp_sensor, runID.Data()));
        }
      }
    }
 
    if (m_3Samples) {
      //look for U histogram - OFFLINE ZS for 3 samples
      tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_Strip3CountU", tmp_layer, tmp_ladder, tmp_sensor);
 
      auto hStrip3CountU = findHist(tmpname.Data());
      if (hStrip3CountU == NULL) {
        B2INFO("Occupancy U histogram not found for 3 samples");
        setOccStatus(-1, m_occU3Samples);
 
      } else {
        Float_t occU = getOccupancy(hStrip3CountU->GetEntries(), tmp_layer, nEvents);
        m_hOccupancy3Samples->fill(m_SVDModules[i], 1, occU);
        setOccStatus(occU, m_occU3Samples);
      }
    }
 
    // groupId0 side U
    TString tmpnameGrpId0 = Form("SVDExpReco/SVDDQM_%d_%d_%d_StripCountSignalGroupIDsU", tmp_layer, tmp_ladder, tmp_sensor);
    auto hStripCountGroupIDsU = findHist(tmpnameGrpId0.Data());
    if (hStripCountGroupIDsU == NULL) {
      B2INFO("Occupancy U histogram for group Id0 not found");
      setOccStatus(-1, m_occUGroupId0, kGroupIDs);
 
    } else {
      Float_t occU = getOccupancy(hStripCountGroupIDsU->GetEntries(), tmp_layer, nEvents);
      m_hOccupancyGroupId0->fill(m_SVDModules[i], 1, occU);
      setOccStatus(occU, m_occUGroupId0, kGroupIDs);
    }
 
    //look for V histogram - OFFLINE ZS
    tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_StripCountV", tmp_layer, tmp_ladder, tmp_sensor);
 
    auto hStripCountV = findHist(tmpname.Data());
    if (hStripCountV == NULL) {
      B2INFO("Occupancy V histogram not found");
      setOccStatus(-1, m_occVstatus);
 
    } else {
      Float_t occV = getOccupancy(hStripCountV->GetEntries(), tmp_layer, nEvents, true);
      m_hOccupancy->fill(m_SVDModules[i], 0, occV);
      setOccStatus(occV, m_occVstatus);
 
      //produce the occupancy plot
      if (m_additionalPlots) {
        m_hStripOccupancyV[i].Clear();
        hStripCountV->Copy(m_hStripOccupancyV[i]);
        m_hStripOccupancyV[i].Scale(1 / nEvents);
        m_hStripOccupancyV[i].SetName(Form("%d_%d_%d_OccupancyV", tmp_layer, tmp_ladder, tmp_sensor));
        m_hStripOccupancyV[i].SetTitle(Form("SVD Sensor %d_%d_%d V-Strip OFFLINE Occupancy vs Strip Number %s", tmp_layer, tmp_ladder,
                                            tmp_sensor, runID.Data()));
      } else {
        if (i < (unsigned int)m_sensors) {
          m_hStripOccupancyV[i].Clear();
          hStripCountV->Copy(m_hStripOccupancyV[i]);
          m_hStripOccupancyV[i].Scale(1 / nEvents);
          m_hStripOccupancyV[i].SetName(Form("%d_%d_%d_OccupancyV", tmp_layer, tmp_ladder, tmp_sensor));
          m_hStripOccupancyV[i].SetTitle(Form("SVD Sensor %d_%d_%d V-Strip OFFLINE Occupancy vs Strip Number %s", tmp_layer, tmp_ladder,
                                              tmp_sensor, runID.Data()));
        }
      }
    }
 
    if (m_3Samples) {
      //look for V histogram - OFFLINE ZS for 3 samples
      tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_Strip3CountV", tmp_layer, tmp_ladder, tmp_sensor);
 
      auto hStrip3CountV = findHist(tmpname.Data());
      if (hStrip3CountV == NULL) {
        B2INFO("Occupancy V histogram not found");
        setOccStatus(-1, m_occV3Samples);
 
      } else {
        Float_t occV = getOccupancy(hStrip3CountV->GetEntries(), tmp_layer, nEvents, true);
        m_hOccupancy3Samples->fill(m_SVDModules[i], 0, occV);
        setOccStatus(occV, m_occV3Samples);
      }
    }
 
    // groupId0 side V
    tmpnameGrpId0 = Form("SVDExpReco/SVDDQM_%d_%d_%d_StripCountSignalGroupIDsV", tmp_layer, tmp_ladder, tmp_sensor);
 
    auto hStripCountGroupIDsV = findHist(tmpnameGrpId0.Data());
    if (hStripCountGroupIDsV == NULL) {
      B2INFO("Occupancy U histogram for group Id0 not found");
      setOccStatus(-1, m_occVGroupId0, kGroupIDs);
 
    } else {
      Float_t occV = getOccupancy(hStripCountGroupIDsV->GetEntries(), tmp_layer, nEvents, true);
      m_hOccupancyGroupId0->fill(m_SVDModules[i], 0, occV);
      setOccStatus(occV, m_occVGroupId0, kGroupIDs);
    }
 
    //look for V histogram - ONLINE ZS
    tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_OnlineZSStripCountV", tmp_layer, tmp_ladder, tmp_sensor);
 
    auto hOnlineStripCountV = findHist(tmpname.Data());
    if (hOnlineStripCountV == NULL) {
      B2INFO("OnlineOccupancy V histogram not found");
      setOccStatus(-1, m_onlineOccVstatus, kOnline);
 
    } else {
      Float_t onlineOccV = getOccupancy(hOnlineStripCountV->GetEntries(), tmp_layer, nEvents, true);
      m_hOnlineOccupancy->fill(m_SVDModules[i], 0, onlineOccV);
 
      for (int b = 1; b < hOnlineStripCountV->GetNbinsX() + 1; b++) {
        hOnlineStripCountV->SetBinContent(b, hOnlineStripCountV->GetBinContent(b) / nEvents * 100);
      }
      hOnlineStripCountV->GetYaxis()->SetTitle("ZS3 occupancy (%)");
      setOccStatus(onlineOccV, m_onlineOccVstatus, kOnline);
    }
 
    if (m_3Samples) {
      //look for V histogram - ONLINE ZS for 3 samples
      tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_OnlineZSStrip3CountV", tmp_layer, tmp_ladder, tmp_sensor);
 
      auto hOnlineStrip3CountV = findHist(tmpname.Data());
      if (hOnlineStrip3CountV == NULL) {
        B2INFO("OnlineOccupancy3 V histogram not found");
        setOccStatus(-1, m_onlineOccV3Samples, kOnline);
 
      } else {
        Float_t onlineOccV = getOccupancy(hOnlineStrip3CountV->GetEntries(), tmp_layer, nEvents, true);
        m_hOnlineOccupancy3Samples->fill(m_SVDModules[i], 0, onlineOccV);
 
        for (int b = 1; b < hOnlineStrip3CountV->GetNbinsX() + 1; b++) {
          hOnlineStrip3CountV->SetBinContent(b, hOnlineStrip3CountV->GetBinContent(b) / nEvents * 100);
        }
        hOnlineStrip3CountV->GetYaxis()->SetTitle("ZS3 occupancy (%)");
        setOccStatus(onlineOccV, m_onlineOccV3Samples, kOnline);
      }
    }
 
    //look for U histogram - ONLINE ZS
    tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_OnlineZSStripCountU", tmp_layer, tmp_ladder, tmp_sensor);
 
    auto hOnlineStripCountU = findHist(tmpname.Data());
    if (hOnlineStripCountU == NULL) {
      B2INFO("OnlineOccupancy U histogram not found");
      setOccStatus(-1, m_onlineOccUstatus, true);
 
    } else {
      Float_t onlineOccU = getOccupancy(hOnlineStripCountU->GetEntries(), tmp_layer, nEvents);
      m_hOnlineOccupancy->fill(m_SVDModules[i], 1, onlineOccU);
 
      for (int b = 1; b < hOnlineStripCountU->GetNbinsX() + 1; b++) {
        hOnlineStripCountU->SetBinContent(b, hOnlineStripCountU->GetBinContent(b) / nEvents * 100);
      }
      hOnlineStripCountU->GetYaxis()->SetTitle("ZS3 occupancy (%)");
      setOccStatus(onlineOccU, m_onlineOccUstatus, kOnline);
    }
 
    if (m_3Samples) {
      //look for U histogram - ONLINE ZS for 3 samples
      tmpname = Form("SVDExpReco/SVDDQM_%d_%d_%d_OnlineZSStrip3CountU", tmp_layer, tmp_ladder, tmp_sensor);
 
      auto hOnlineStrip3CountU = findHist(tmpname.Data());
      if (hOnlineStrip3CountU == NULL) {
        B2INFO("OnlineOccupancy3 U histogram not found");
        setOccStatus(-1, m_onlineOccU3Samples, kOnline);
 
      } else {
        Float_t onlineOccU = getOccupancy(hOnlineStrip3CountU->GetEntries(), tmp_layer, nEvents);
        m_hOnlineOccupancy3Samples->fill(m_SVDModules[i], 1, onlineOccU);
 
        for (int b = 1; b < hOnlineStrip3CountU->GetNbinsX() + 1; b++) {
          hOnlineStrip3CountU->SetBinContent(b, hOnlineStrip3CountU->GetBinContent(b) / nEvents * 100);
        }
        hOnlineStrip3CountU->GetYaxis()->SetTitle("ZS3 occupancy (%)");
        setOccStatus(onlineOccU, m_onlineOccU3Samples, kOnline);
      }
    }
 
    //update sensor occupancy canvas U and V
    if (m_additionalPlots) {
      m_cStripOccupancyU[i]->cd();
      m_hStripOccupancyU[i].Draw("histo");
      m_cStripOccupancyV[i]->cd();
      m_hStripOccupancyV[i].Draw("histo");
    } else {
      if (i == 0 || i == 1) {
        m_cStripOccupancyU[i]->cd();
        m_hStripOccupancyU[i].Draw("histo");
        m_cStripOccupancyV[i]->cd();
        m_hStripOccupancyV[i].Draw("histo");
      }
    }
  }
 
  //update summary offline occupancy U canvas
  updateCanvases(m_hOccupancy, m_cOccupancyU, m_cOccupancyRPhiViewU,  m_occUstatus, true);
 
  //update summary offline occupancy V canvas
  updateCanvases(m_hOccupancy, m_cOccupancyV, m_cOccupancyRPhiViewV,  m_occVstatus, false);
 
  //update summary offline occupancy U canvas for groupId0
  updateCanvases(m_hOccupancyGroupId0, m_cOccupancyUGroupId0, m_cOccupancyRPhiViewUGroupId0,  m_occUGroupId0, true, kGroupIDs);
 
  //update summary offline occupancy V canvas for groupId0
  updateCanvases(m_hOccupancyGroupId0, m_cOccupancyVGroupId0, m_cOccupancyRPhiViewVGroupId0,  m_occVGroupId0, false, kGroupIDs);
 
  //update summary online occupancy U canvas
  updateCanvases(m_hOnlineOccupancy, m_cOnlineOccupancyU, m_cOnlineOccupancyRPhiViewU,  m_onlineOccUstatus, true, kOnline);
 
  //update summary online occupancy V canvas
  updateCanvases(m_hOnlineOccupancy, m_cOnlineOccupancyV, m_cOnlineOccupancyRPhiViewV,  m_onlineOccVstatus, false, kOnline);
 
  if (m_3Samples) {
    //update summary offline occupancy U canvas for 3 samples
    updateCanvases(m_hOccupancy3Samples, m_cOccupancyU3Samples, m_cOccupancyRPhiViewU3Samples,  m_occU3Samples, true);
 
    //update summary offline occupancy V canvas for 3 samples
    updateCanvases(m_hOccupancy3Samples, m_cOccupancyV3Samples, m_cOccupancyRPhiViewV3Samples,  m_occV3Samples, false);
 
    //update summary online occupancy U canvas for 3 samples
    updateCanvases(m_hOnlineOccupancy3Samples, m_cOnlineOccupancyU3Samples, m_cOnlineOccupancyRPhiViewU3Samples,  m_onlineOccU3Samples,
                   true, kOnline);
 
    //update summary online occupancy V canvas for 3 samples
    updateCanvases(m_hOnlineOccupancy3Samples, m_cOnlineOccupancyV3Samples, m_cOnlineOccupancyRPhiViewV3Samples,  m_onlineOccV3Samples,
                   false, kOnline);
  }
 
  if (m_printCanvas) {
    m_cOccupancyU->Print("c_SVDOccupancyU.pdf");
    m_cOccupancyV->Print("c_SVDOccupancyV.pdf");
    m_cOnlineOccupancyU->Print("c_SVDOnlineOccupancyU.pdf");
    m_cOnlineOccupancyV->Print("c_SVDOnlineOccupancyV.pdf");
  }
}
 
void DQMHistAnalysisSVDOccupancyModule::endRun()
{
  B2DEBUG(10, "DQMHistAnalysisSVDOccupancy: endRun called");
}
 
 
void DQMHistAnalysisSVDOccupancyModule::terminate()
{
  B2DEBUG(10, "DQMHistAnalysisSVDOccupancy: terminate called");
 
  delete m_hOccupancy;
  delete m_hOnlineOccupancy;
  delete m_hOccupancyGroupId0;
  delete m_hOccupancy3Samples;
  delete m_hOnlineOccupancy3Samples;
 
  delete m_cOccupancyChartChip;
 
  delete m_cOccupancyU;
  delete m_cOccupancyV;
 
  delete m_cOnlineOccupancyU;
  delete m_cOnlineOccupancyV;
 
  delete m_cOccupancyU3Samples;
  delete m_cOccupancyV3Samples;
 
  delete m_cOnlineOccupancyU3Samples;
  delete m_cOnlineOccupancyV3Samples;
 
  delete m_cOccupancyUGroupId0;
  delete m_cOccupancyVGroupId0;
 
  if (m_RPhiView) {
    delete m_cOccupancyRPhiViewU;
    delete m_cOccupancyRPhiViewV;
 
    delete m_cOnlineOccupancyRPhiViewU;
    delete m_cOnlineOccupancyRPhiViewV;
 
    delete m_cOccupancyRPhiViewU3Samples;
    delete m_cOccupancyRPhiViewV3Samples;
 
    delete m_cOnlineOccupancyRPhiViewU3Samples;
    delete m_cOnlineOccupancyRPhiViewV3Samples;
    if (m_RPhiViewId0) {
      delete m_cOccupancyRPhiViewUGroupId0;
      delete m_cOccupancyRPhiViewVGroupId0;
    }
  }
 
  for (int module = 0; module < m_sensors; module++) {
    delete m_cStripOccupancyU[module];
    delete m_cStripOccupancyV[module];
  }
  delete m_cStripOccupancyU;
  delete m_cStripOccupancyV;
}
 
Float_t DQMHistAnalysisSVDOccupancyModule::getOccupancy(float entries, int tmp_layer, int nEvents, bool sideV)
{
  Int_t nStrips = 768;
  if (tmp_layer != 3 && sideV)
    nStrips = 512;
 
  return (entries / nStrips / nEvents * 100);
}
 
void DQMHistAnalysisSVDOccupancyModule::setOccStatus(float occupancy, svdStatus& occupancyStatus, int histoType)
{
  if (histoType == kOnline) {
    if (occupancy < 0)
      occupancyStatus = std::max(noStat, occupancyStatus);
    else if (occupancy <= m_onlineOccEmpty) {
      occupancyStatus = std::max(lowStat, occupancyStatus);
    } else if (occupancy < m_onlineOccWarning && occupancy >= m_onlineOccEmpty) {
      occupancyStatus = std::max(good, occupancyStatus);
    } else if (occupancy > m_onlineOccWarning && occupancy < m_onlineOccError) {
      occupancyStatus = std::max(warning, occupancyStatus);
    } else if (occupancy >= m_onlineOccError)  {
      occupancyStatus = std::max(error, occupancyStatus);
    }
  } else if (histoType == kOffline) {
    if (occupancy < 0)
      occupancyStatus = std::max(noStat, occupancyStatus);
    else if (occupancy <= m_occEmpty) {
      occupancyStatus = std::max(lowStat, occupancyStatus);
    } else if (occupancy < m_occWarning && occupancy >= m_occEmpty) {
      occupancyStatus = std::max(good, occupancyStatus);
    } else if (occupancy > m_occWarning && occupancy < m_occError) {
      occupancyStatus = std::max(warning, occupancyStatus);
    } else if (occupancy >= m_occError)  {
      occupancyStatus = std::max(error, occupancyStatus);
    }
  } else if (histoType == kGroupIDs) {
    if (occupancy < 0)
      occupancyStatus = std::max(noStat, occupancyStatus);
    else if (occupancy <= m_groupIDsOccEmpty) {
      occupancyStatus = std::max(lowStat, occupancyStatus);
    } else if (occupancy < m_groupIDsOccWarning && occupancy >= m_groupIDsOccEmpty) {
      occupancyStatus = std::max(good, occupancyStatus);
    } else if (occupancy > m_groupIDsOccWarning && occupancy < m_groupIDsOccError) {
      occupancyStatus = std::max(warning, occupancyStatus);
    } else if (occupancy >= m_groupIDsOccError)  {
      occupancyStatus = std::max(error, occupancyStatus);
    }
  }
}