PXDDAQDQMModule.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.                  *
 **************************************************************************/
 
#include <pxd/modules/pxdDQM/PXDDAQDQMModule.h>
#include <pxd/dataobjects/PXDDAQStatus.h>
#include <rawdata/dataobjects/RawSVD.h>
#include <mdst/dataobjects/EventLevelTriggerTimeInfo.h>
 
#include <TDirectory.h>
#include <boost/format.hpp>
#include <string>
 
using namespace std;
using namespace Belle2;
using namespace Belle2::PXD;
using namespace Belle2::PXD::PXDError;
using namespace Belle2::VXD;
using boost::format;
 
//-----------------------------------------------------------------
//                 Register the Module
//-----------------------------------------------------------------
REG_MODULE(PXDDAQDQM);
 
//-----------------------------------------------------------------
//                 Implementation
//-----------------------------------------------------------------
 
PXDDAQDQMModule::PXDDAQDQMModule() : HistoModule(), m_vxdGeometry(VXD::GeoCache::getInstance())
{
  //Set module properties
  setDescription("Monitor DAQ errors");
  setPropertyFlags(c_ParallelProcessingCertified);
  addParam("histogramDirectoryName", m_histogramDirectoryName, "Name of the directory where histograms will be placed",
           std::string("pxdDAQ"));
}
 
void PXDDAQDQMModule::defineHisto()
{
  TDirectory* oldDir = gDirectory;
  if (m_histogramDirectoryName != "") {
    oldDir->mkdir(m_histogramDirectoryName.c_str());// do not rely on return value, might be ZERO
    oldDir->cd(m_histogramDirectoryName.c_str());
  }
 
  hDAQErrorEvent = new TH1D("PXDDAQError", "PXDDAQError/Event;;Count", ONSEN_USED_TYPE_ERR, 0, ONSEN_USED_TYPE_ERR);
  hDAQErrorDHC = new TH2D("PXDDAQDHCError", "PXDDAQError/DHC;DHC ID;", 16, 0, 16, ONSEN_USED_TYPE_ERR, 0, ONSEN_USED_TYPE_ERR);
  hDAQErrorDHE = new TH2D("PXDDAQDHEError", "PXDDAQError/DHE;DHE ID;", 64, 0, 64, ONSEN_USED_TYPE_ERR, 0, ONSEN_USED_TYPE_ERR);
  hDAQUseableModule = new TH1D("PXDDAQUseableModule", "PXDDAQUseableModule/DHE;DHE ID;", 64, 0, 64);
  hDAQNotUseableModule = new TH1D("PXDDAQNotUseableModule", "PXDDAQNotUseableModule/DHE;DHE ID;", 64, 0, 64);
  hDAQDHPDataMissing = new TH1D("PXDDAQDHPDataMissing", "PXDDAQDHPDataMissing/DHE*DHP;DHE ID;", 64 * 4, 0, 64);
  hDAQEndErrorDHC = new TH2D("PXDDAQDHCEndError", "PXDDAQEndError/DHC;DHC ID;", 16, 0, 16, 32, 0, 32);
  hDAQEndErrorDHE = new TH2D("PXDDAQDHEEndError", "PXDDAQEndError/DHE;DHE ID;", 64, 0, 64, 4 * 2 * 8, 0, 4 * 2 * 8);
 
  // histograms might get unreadable, but, if necessary, you can zoom in anyways.
  // we could use full alphanumeric histograms, but then, the labels would change (in the worst case) depending on observed errors
  // and ... the histogram would contain NO labels if there is NO error ... confusing.
  // ... an we would have to use alphanumeric X axis (DHE ID, DHC ID), too)
  for (int i = 0; i < ONSEN_USED_TYPE_ERR; i++) {
    const char* label = getPXDBitErrorName(i).c_str();
    hDAQErrorEvent->GetXaxis()->SetBinLabel(i + 1, label);
    hDAQErrorDHE->GetYaxis()->SetBinLabel(i + 1, label);
    hDAQErrorDHC->GetYaxis()->SetBinLabel(i + 1, label);
  }
 
  hDAQErrorEvent->LabelsOption("v"); // rotate the labels.
 
  std::vector<VxdID> sensors = m_vxdGeometry.getListOfSensors();
  for (VxdID& avxdid : sensors) {
    VXD::SensorInfoBase info = m_vxdGeometry.getSensorInfo(avxdid);
    if (info.getType() != VXD::SensorInfoBase::PXD) continue;
    //Only interested in PXD sensors
 
    TString buff = (std::string)avxdid;
    TString buffus = buff;
    buffus.ReplaceAll(".", "_");
 
//     string s = str(format("DHE %d:%d:%d (DHH ID %02Xh)") % num1 % num2 % num3 % i);
//     string s2 = str(format("_%d.%d.%d") % num1 % num2 % num3);
 
    hDAQDHETriggerGate[avxdid] = new TH1D("PXDDAQDHETriggerGate_" + buffus,
                                          "TriggerGate DHE " + buff + "; Trigger Gate; Counts", 192, 0, 192);
    hDAQDHEReduction[avxdid] = new TH1D("PXDDAQDHEDataReduction_" + buffus, "Data Reduction DHE " + buff + "; Raw/Red; Counts", 200, 0,
                                        40);// If max changed, check overflow copy below
    hDAQCM[avxdid] = new TH2D("PXDDAQCM_" + buffus, "Common Mode on DHE " + buff + "; Gate+Chip*192; Common Mode", 192 * 4, 0, 192 * 4,
                              64, 0, 64);
    hDAQCM2[avxdid] = new TH1D("PXDDAQCM2_" + buffus, "Common Mode on DHE " + buff + "; Common Mode", 64, 0, 64);
  }
  for (int i = 0; i < 16; i++) {
    hDAQDHCReduction[i] = new TH1D(("PXDDAQDHCDataReduction_" + str(format("%d") % i)).c_str(),
                                   ("Data Reduction DHC " + str(format(" %d") % i) + "; Raw/Red; Counts").c_str(), 200, 0,
                                   40); // If max changed, check overflow copy below
  }
//   hDAQErrorEvent->LabelsDeflate("X");
//   hDAQErrorEvent->LabelsOption("v");
//   hDAQErrorEvent->SetStats(0);
  hEODBAfterInjLER  = new TH1I("PXDEODBInjLER", "PXDEODBInjLER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hEODBAfterInjHER  = new TH1I("PXDEODBInjHER", "PXDEODBInjHER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hCM63AfterInjLER  = new TH1I("PXDCM63InjLER", "PXDCM63InjLER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hCM63AfterInjHER  = new TH1I("PXDCM63InjHER", "PXDCM63InjHER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hTruncAfterInjLER  = new TH1I("PXDTruncInjLER", "PXDTruncInjLER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hTruncAfterInjHER  = new TH1I("PXDTruncInjHER", "PXDTruncInjHER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hMissAfterInjLER  = new TH1I("PXDMissInjLER", "PXDMissInjLER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hMissAfterInjHER  = new TH1I("PXDMissInjHER", "PXDMissInjHER/Time;Time in #mus;Events/Time (5 #mus bins)", 4000, 0, 20000);
  hEODBTrgDiff  = new TH1I("PXDEODBTrgDiff", "PXDEODBTrgDiff/DiffTime;DiffTime in #mus;Events/Time (1 #mus bins)", 2000, 0, 2000);
  hCM63TrgDiff  = new TH1I("PXDCM63TrgDiff", "PXDCM63TrgDiff/DiffTime;DiffTime in #mus;Events/Time (1 #mus bins)", 2000, 0, 2000);
  hTruncTrgDiff  = new TH1I("PXDTruncTrgDiff", "PXDTruncTrgDiff/DiffTime;DiffTime in #mus;Events/Time (1 #mus bins)", 2000, 0, 2000);
  hMissTrgDiff  = new TH1I("PXDMissTrgDiff", "PXDMissTrgDiff/DiffTime;DiffTime in #mus;Events/Time (1 #mus bins)", 2000, 0, 2000);
 
  hDAQStat  = new TH1D("PXDDAQStat", "PXDDAQStat", 20, 0, 20);
  auto xa = hDAQStat->GetXaxis();
  if (xa) {
    // underflow: number of events -> for normalize
    xa->SetBinLabel(0 + 1, "EODB/HLT rej"); // event of doom or rejected
    xa->SetBinLabel(1 + 1, "Trunc 8%");
    xa->SetBinLabel(2 + 1, "HER Trunc");
    xa->SetBinLabel(3 + 1, "LER Trunc");
    xa->SetBinLabel(4 + 1, "CM63");
    xa->SetBinLabel(5 + 1, "HER CM63");
    xa->SetBinLabel(6 + 1, "LER CM63");
    xa->SetBinLabel(7 + 1, "HER CM63>1ms");
    xa->SetBinLabel(8 + 1, "LER CM63>1ms");
    xa->SetBinLabel(9 + 1, "HER Trunc>1ms");
    xa->SetBinLabel(10 + 1, "LER Trunc>1ms");
    xa->SetBinLabel(11 + 1, "MissFrame");
    xa->SetBinLabel(12 + 1, "Timeout");
    xa->SetBinLabel(13 + 1, "Link Down");
    xa->SetBinLabel(14 + 1, "Mismatch");
    xa->SetBinLabel(15 + 1, "HER MissFrame");
    xa->SetBinLabel(16 + 1, "LER MissFrame");
    xa->SetBinLabel(17 + 1, "HER MissFrame>1ms");
    xa->SetBinLabel(18 + 1, "LER MissFrame>1ms");
    xa->SetBinLabel(19 + 1, "unused");
  }
  // cd back to root directory
  oldDir->cd();
}
 
void PXDDAQDQMModule::initialize()
{
  REG_HISTOGRAM
  m_storeDAQEvtStats.isRequired();
  m_rawSVD.isOptional(); 
  m_EventLevelTriggerTimeInfo.isRequired();
 
}
 
void PXDDAQDQMModule::beginRun()
{
  hDAQErrorEvent->Reset();
  hDAQErrorDHC->Reset();
  hDAQErrorDHE->Reset();
  hDAQUseableModule->Reset();
  hDAQNotUseableModule->Reset();
  hDAQDHPDataMissing->Reset();
  hDAQEndErrorDHC->Reset();
  hDAQEndErrorDHE->Reset();
  for (auto& it : hDAQDHETriggerGate) if (it.second) it.second->Reset();
  for (auto& it : hDAQDHCReduction) if (it.second) it.second->Reset();
  for (auto& it : hDAQDHEReduction) if (it.second) it.second->Reset();
  for (auto& it : hDAQCM) if (it.second) it.second->Reset();
  for (auto& it : hDAQCM2) if (it.second) it.second->Reset();
  if (hCM63AfterInjLER) hCM63AfterInjLER->Reset();
  if (hCM63AfterInjHER) hCM63AfterInjHER->Reset();
  if (hTruncAfterInjLER) hTruncAfterInjLER->Reset();
  if (hTruncAfterInjHER) hTruncAfterInjHER->Reset();
  if (hMissAfterInjLER) hMissAfterInjLER->Reset();
  if (hMissAfterInjHER) hMissAfterInjHER->Reset();
  hDAQStat->Reset();
}
 
void PXDDAQDQMModule::event()
{
  hDAQErrorEvent->Fill(-1);// Event counter
  hDAQStat->Fill(-1); // to normalize to the number of events
  hDAQDHPDataMissing->Fill(-1); // to normalize to the number of events
  hDAQErrorDHC->Fill(-1, -1); // to normalize to the number of events
  hDAQErrorDHE->Fill(-1, -1); // to normalize to the number of events
  for (auto& it : hDAQCM2) if (it.second) it.second->Fill(-1); // to normalize to the number of events
  for (auto& it : hDAQCM) if (it.second) it.second->Fill(-1, -1); // to normalize to the number of events
 
  bool eodbFlag = m_rawSVD.getEntries() == 0;
 
  bool truncFlag = false; // flag events which are DHE truncated
  bool nolinkFlag = false; // flag events which are DHE truncated
  bool missingFlag = false; // flag events where frame is missing
  bool timeoutFlag = false; // flag events where frame timeout
  bool mismatchFlag = false; // flag events where trig mismatched
  bool cm63Flag = false; // flag event which are CM63 truncated
 
  B2DEBUG(20, "Iterate PXD DAQ Status");
  auto evt = *m_storeDAQEvtStats;
  auto evt_emask = evt.getErrorMask();
  for (int i = 0; i < ONSEN_MAX_TYPE_ERR; i++) {
    if (evt_emask[i]) hDAQErrorEvent->Fill(getPXDBitErrorName(i).c_str(), 1);
  }
  B2DEBUG(20, "Iterate PXD Packets, Err " << evt_emask);
  for (auto& pkt : evt) {
    B2DEBUG(20, "Iterate PXD DHC in Pkt " << pkt.getPktIndex());
    for (auto& dhc : pkt) {
      hDAQErrorDHC->Fill(dhc.getDHCID(), -1);// normalize
      auto dhc_emask = dhc.getErrorMask();
      for (int i = 0; i < ONSEN_MAX_TYPE_ERR; i++) {
        if (dhc_emask[i]) hDAQErrorDHC->Fill(dhc.getDHCID(), i);
      }
      unsigned int cmask = dhc.getEndErrorInfo();
      for (int i = 0; i < 32; i++) {
        unsigned int mask = (1u << i);
        if ((cmask & mask) == mask) hDAQEndErrorDHC->Fill(dhc.getDHCID(), i);
      }
      if (hDAQDHCReduction[dhc.getDHCID()]) {
        float red = dhc.getRedCnt() ? float(dhc.getRawCnt()) / dhc.getRedCnt() : 0.;
        B2DEBUG(98, "==DHC " << dhc.getDHCID() << "(Raw)" << dhc.getRawCnt() << " / (Red)" << dhc.getRedCnt() << " = " << red);
        if (red >= 40.) red = 39.999999999; // Bad, bad workaround. but we want to see the overflows
        hDAQDHCReduction[dhc.getDHCID()]->Fill(red);
      }
      B2DEBUG(20, "Iterate PXD DHE in DHC " << dhc.getDHCID() << " , Err " << dhc_emask);
      for (auto& dhe : dhc) {
        hDAQErrorDHE->Fill(dhe.getDHEID(), -1);// normalize
        auto dhe_emask = dhe.getErrorMask();
        B2DEBUG(20, "DHE " << dhe.getDHEID() << " , Err " << dhe_emask);
        for (int i = 0; i < ONSEN_MAX_TYPE_ERR; i++) {
          if (dhe_emask[i]) hDAQErrorDHE->Fill(dhe.getDHEID(), i);
        }
        if (dhe.isUsable()) {
          hDAQUseableModule->Fill(dhe.getDHEID());
        } else {
          hDAQNotUseableModule->Fill(dhe.getDHEID());
        }
        for (int i = 0; i < 4; i++) {
          if ((dhe.getDHPFoundMask() & (1 << i)) == 0) hDAQDHPDataMissing->Fill(dhe.getDHEID() + i * 0.25);
        }
        for (auto& dhp : dhe) {
          truncFlag |= dhp.getTruncated(); // new firmware workaround flag
        }
        unsigned int emask = dhe.getEndErrorInfo();
        // TODO differentiate between link-lost and truncation
        for (int i = 0; i < 4 * 2; i++) {
          auto sm = (emask >> i * 4) & 0xF;
          if (sm >= 8) sm = 7; // clip unknown to 7, as value >6 undefined for now
          if (sm > 0) hDAQEndErrorDHE->Fill(dhe.getDHEID(), i * 8 + sm); // we dont want to fill noerror=0
          missingFlag |= sm == 0x1; // missing
          timeoutFlag |= sm == 0x2; // timeout
          nolinkFlag |= sm == 0x3; // link down
          // 4 is DHP masked
          mismatchFlag |= sm == 0x5; // start/end mismatch
          truncFlag |= sm == 0x6; // trunc because of size
        }
 
        if (hDAQDHETriggerGate[dhe.getSensorID()]) hDAQDHETriggerGate[dhe.getSensorID()]->Fill(dhe.getTriggerGate());
        if (hDAQDHEReduction[dhe.getSensorID()]) {
          float red = dhe.getRedCnt() ? float(dhe.getRawCnt()) / dhe.getRedCnt() : 0.;
          B2DEBUG(98, "==DHE " << dhe.getSensorID() << "(Raw)" << dhe.getRawCnt() << " / (Red)" << dhe.getRedCnt() << " = " << red);
          if (red >= 40.) red = 39.999999999; // Bad, bad workaround. but we want to see the overflows
          hDAQDHEReduction[dhe.getSensorID()]->Fill(red);
        }
        for (auto cm = dhe.cm_begin(); cm < dhe.cm_end(); ++cm) {
          // uint8_t, uint16_t, uint8_t ; tuple of Chip ID (2 bit), Row (10 bit), Common Mode (6 bit)
          if (hDAQCM[dhe.getSensorID()]) hDAQCM[dhe.getSensorID()]->Fill(std::get<0>(*cm) * 192 + std::get<1>(*cm) / 4, std::get<2>(*cm));
          if (hDAQCM2[dhe.getSensorID()]) hDAQCM2[dhe.getSensorID()]->Fill(std::get<2>(*cm));
          //if (hDAQCM2[dhe.getSensorID()] && std::get<1>(*cm) < 768 - 2) {
          // Deactivated again as the clean-up effect was insignificant for beam data
          // ignore the always-on rows (second to last: 766 and last 767)
          // row in CM array is already in readout-direction and not V-ID
          // hDAQCM2[dhe.getSensorID()]->Fill(std::get<2>(*cm));
          //}
          cm63Flag |= 63 == std::get<2>(*cm);
        }
      }
    }
  }
 
 
 
  // And check if the stored data is valid
  if (m_EventLevelTriggerTimeInfo.isValid() and m_EventLevelTriggerTimeInfo->isValid()) {
 
    double lasttrig = m_EventLevelTriggerTimeInfo->getTimeSincePrevTrigger() / 127.; //  127MHz clock ticks to us, inexact rounding
    if (eodbFlag && hEODBTrgDiff) hEODBTrgDiff->Fill(lasttrig);
    if (cm63Flag && hCM63TrgDiff) hCM63TrgDiff->Fill(lasttrig);
    if (truncFlag && hTruncTrgDiff) hTruncTrgDiff->Fill(lasttrig);
    if (missingFlag && hMissTrgDiff) hMissTrgDiff->Fill(lasttrig);
 
    // check time overflow, too long ago
    if (m_EventLevelTriggerTimeInfo->hasInjection()) {
      // get last injection time
      double difference = m_EventLevelTriggerTimeInfo->getTimeSinceLastInjection() / 127.; //  127MHz clock ticks to us, inexact rounding
      if (m_EventLevelTriggerTimeInfo->isHER()) {
        if (eodbFlag) {
          if (hEODBAfterInjHER) hEODBAfterInjHER->Fill(difference);
        }
        if (cm63Flag) {
          hDAQStat->Fill(5); // sum CM63 after HER
          if (difference > 1000) hDAQStat->Fill(7); // sum CM63 after HER, but outside injections, 1ms
          if (hCM63AfterInjHER) hCM63AfterInjHER->Fill(difference);
        }
        if (truncFlag) {
          hDAQStat->Fill(2); // sum truncs after HER
          if (difference > 1000) hDAQStat->Fill(9); // sum truncs after HER, but outside injections, 1ms
          if (hTruncAfterInjHER) hTruncAfterInjHER->Fill(difference);
        }
        if (missingFlag) {
          hDAQStat->Fill(15); // sum missframe after HER
          if (difference > 1000) hDAQStat->Fill(17); // sum missframe after HER, but outside injections, 1ms
          if (hMissAfterInjHER) hMissAfterInjHER->Fill(difference);
        }
      } else {
        if (eodbFlag) {
          if (hEODBAfterInjLER) hEODBAfterInjLER->Fill(difference);
        }
        if (cm63Flag) {
          hDAQStat->Fill(6); // sum CM63 after LER
          if (difference > 1000) hDAQStat->Fill(8); // sum CM63 after LER, but outside injections, 1ms
          if (hCM63AfterInjLER) hCM63AfterInjLER->Fill(difference);
        }
        if (truncFlag) {
          hDAQStat->Fill(3); // sum truncs after LER
          if (difference > 1000) hDAQStat->Fill(10); // sum truncs after LER, but outside injections, 1ms
          if (hTruncAfterInjLER) hTruncAfterInjLER->Fill(difference);
        }
        if (missingFlag) {
          hDAQStat->Fill(16); // sum missframe after LER
          if (difference > 1000) hDAQStat->Fill(18); // sum missframe after LER, but outside injections, 1ms
          if (hMissAfterInjLER) hMissAfterInjLER->Fill(difference);
        }
      }
    }
  }
 
// make some nice statistics
  if (truncFlag) hDAQStat->Fill(1);
  if (cm63Flag) hDAQStat->Fill(4);
  if (missingFlag) hDAQStat->Fill(11);
  if (timeoutFlag) hDAQStat->Fill(12);
  if (nolinkFlag) hDAQStat->Fill(13);
  if (mismatchFlag) hDAQStat->Fill(14);
 
// Check Event-of-doom-busted or otherwise HLT rejected events
  if (eodbFlag) hDAQStat->Fill(0);
}