SiStripCertificationInfo.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
#include "DQM/SiStripMonitorClient/interface/SiStripUtility.h"
#include "SiStripCertificationInfo.h"

#include "DataFormats/Histograms/interface/DQMToken.h"

//Run Info
#include "CondFormats/RunInfo/interface/RunInfo.h"

#include <iostream>
#include <iomanip>
#include <cstdio>
#include <string>
#include <sstream>
#include <cmath>

SiStripCertificationInfo::SiStripCertificationInfo(edm::ParameterSet const&) {
  consumes<DQMToken, edm::InRun>(edm::InputTag("siStripOfflineAnalyser", "DQMGenerationSiStripAnalyserRun"));
  consumes<DQMToken, edm::InLumi>(edm::InputTag("siStripOfflineAnalyser", "DQMGenerationSiStripAnalyserLumi"));
  detCablingToken_ = esConsumes<edm::Transition::BeginRun>();
  tTopoToken_ = esConsumes<edm::Transition::EndRun>();
  runInfoToken_ = esConsumes<edm::Transition::BeginRun>();
}

void SiStripCertificationInfo::beginRun(edm::Run const& run, edm::EventSetup const& eSetup) {
  edm::LogInfo("SiStripCertificationInfo") << "SiStripCertificationInfo:: Begining of Run";
  detCabling_ = &eSetup.getData(detCablingToken_);

  constexpr int siStripFedIdMin{FEDNumbering::MINSiStripFEDID};
  constexpr int siStripFedIdMax{FEDNumbering::MAXSiStripFEDID};

  if (eSetup.tryToGet<RunInfoRcd>()) {
    if (auto sumFED = eSetup.getHandle(runInfoToken_)) {
      for (auto const fedID : sumFED->m_fed_in) {
        if (fedID >= siStripFedIdMin && fedID <= siStripFedIdMax)
          ++nFEDConnected_;
      }
      LogDebug("SiStripDcsInfo") << " SiStripDcsInfo :: Connected FEDs " << nFEDConnected_;
    }
  }

  auto& dqm_store = *edm::Service<DQMStore>{};
  bookSiStripCertificationMEs(dqm_store);
  fillDummySiStripCertification(dqm_store);
}
//
// -- Book MEs for SiStrip Sertification fractions
//
void SiStripCertificationInfo::bookSiStripCertificationMEs(DQMStore& dqm_store) {
  if (sistripCertificationBooked_)
    return;

  dqm_store.cd();
  std::string strip_dir = "";
  SiStripUtility::getTopFolderPath(dqm_store, "SiStrip", strip_dir);
  if (!strip_dir.empty())
    dqm_store.setCurrentFolder(strip_dir + "/EventInfo");
  else
    dqm_store.setCurrentFolder("SiStrip/EventInfo");

  SiStripCertification = dqm_store.bookFloat("CertificationSummary");

  std::string hname = "CertificationReportMap";
  std::string htitle = "SiStrip Certification for Good Detector Fraction";
  SiStripCertificationSummaryMap = dqm_store.book2D(hname, htitle, 6, 0.5, 6.5, 9, 0.5, 9.5);
  SiStripCertificationSummaryMap->setAxisTitle("Sub Detector Type", 1);
  SiStripCertificationSummaryMap->setAxisTitle("Layer/Disc Number", 2);
  int ibin = 0;
  for (auto const& pr : SubDetMEsMap) {
    ++ibin;
    auto const& det = pr.first;
    SiStripCertificationSummaryMap->setBinLabel(ibin, det);
  }

  SubDetMEs local_mes;
  std::string tag;
  dqm_store.cd();
  if (!strip_dir.empty())
    dqm_store.setCurrentFolder(strip_dir + "/EventInfo/CertificationContents");
  else
    dqm_store.setCurrentFolder("SiStrip/EventInfo/CertificationContents");
  tag = "TIB";

  local_mes.folder_name = "TIB";
  local_mes.subdet_tag = "TIB";
  local_mes.n_layer = 4;
  local_mes.det_fractionME = dqm_store.bookFloat("SiStrip_" + tag);
  SubDetMEsMap.emplace(tag, local_mes);

  tag = "TOB";
  local_mes.folder_name = "TOB";
  local_mes.subdet_tag = "TOB";
  local_mes.n_layer = 6;
  local_mes.det_fractionME = dqm_store.bookFloat("SiStrip_" + tag);
  SubDetMEsMap.emplace(tag, local_mes);

  tag = "TECF";
  local_mes.folder_name = "TEC/PLUS";
  local_mes.subdet_tag = "TEC+";
  local_mes.n_layer = 9;
  local_mes.det_fractionME = dqm_store.bookFloat("SiStrip_" + tag);
  SubDetMEsMap.emplace(tag, local_mes);

  tag = "TECB";
  local_mes.folder_name = "TEC/MINUS";
  local_mes.subdet_tag = "TEC-";
  local_mes.n_layer = 9;
  local_mes.det_fractionME = dqm_store.bookFloat("SiStrip_" + tag);
  SubDetMEsMap.emplace(tag, local_mes);

  tag = "TIDF";
  local_mes.folder_name = "TID/PLUS";
  local_mes.subdet_tag = "TID+";
  local_mes.n_layer = 3;
  local_mes.det_fractionME = dqm_store.bookFloat("SiStrip_" + tag);
  SubDetMEsMap.emplace(tag, local_mes);

  tag = "TIDB";
  local_mes.folder_name = "TID/MINUS";
  local_mes.subdet_tag = "TID-";
  local_mes.n_layer = 3;
  local_mes.det_fractionME = dqm_store.bookFloat("SiStrip_" + tag);
  SubDetMEsMap.emplace(tag, local_mes);

  dqm_store.cd();
  if (!strip_dir.empty())
    dqm_store.setCurrentFolder(strip_dir + "/EventInfo");

  sistripCertificationBooked_ = true;
  dqm_store.cd();
}

void SiStripCertificationInfo::analyze(edm::Event const& event, edm::EventSetup const& eSetup) {}

void SiStripCertificationInfo::endLuminosityBlock(edm::LuminosityBlock const& lumiSeg, edm::EventSetup const& iSetup) {
  edm::LogInfo("SiStripDaqInfo") << "SiStripDaqInfo::endLuminosityBlock";

  if (nFEDConnected_ > 0) {
    auto& dqm_store = *edm::Service<DQMStore>{};
    fillSiStripCertificationMEsAtLumi(dqm_store);
  }
}

void SiStripCertificationInfo::endRun(edm::Run const& run, edm::EventSetup const& eSetup) {
  edm::LogInfo("SiStripCertificationInfo") << "SiStripCertificationInfo:: End Run";

  if (nFEDConnected_ > 0) {
    auto& dqm_store = *edm::Service<DQMStore>{};
    fillSiStripCertificationMEs(dqm_store, eSetup);
  }
}

//
// --Fill SiStrip Certification
//
void SiStripCertificationInfo::fillSiStripCertificationMEs(DQMStore& dqm_store, edm::EventSetup const& eSetup) {
  if (!sistripCertificationBooked_) {
    edm::LogError("SiStripCertificationInfo")
        << " SiStripCertificationInfo::fillSiStripCertificationMEs : MEs missing ";
    return;
  }

  //Retrieve tracker topology from geometry
  const auto tTopo = &eSetup.getData(tTopoToken_);

  resetSiStripCertificationMEs(dqm_store);
  std::string mdir = "MechanicalView";
  dqm_store.cd();
  if (!SiStripUtility::goToDir(dqm_store, mdir))
    return;
  std::string mechanical_dir = dqm_store.pwd();
  uint16_t nDetTot = 0;
  uint16_t nFaultyTot = 0;
  uint16_t nSToNTot = 0;
  float sToNTot = 0.0;
  SiStripFolderOrganizer folder_organizer;
  int xbin = 0;
  for (auto const& [name, subDetME] : SubDetMEsMap) {
    ++xbin;
    MonitorElement* me = subDetME.det_fractionME;
    if (!me)
      continue;
    std::string tag = subDetME.subdet_tag;
    std::string bad_module_folder = mechanical_dir + "/" + subDetME.folder_name + "/" + "BadModuleList";
    std::vector<MonitorElement*> faulty_detMEs = dqm_store.getContents(bad_module_folder);

    uint16_t ndet_subdet = 0;
    uint16_t nfaulty_subdet = 0;
    int nlayer = subDetME.n_layer;
    int ybin = 0;
    for (int ilayer = 0; ilayer < nlayer; ilayer++) {
      uint16_t ndet_layer = detCabling_->connectedNumber(tag, ilayer + 1);
      ndet_subdet += ndet_layer;
      ybin++;
      uint16_t nfaulty_layer = 0;
      for (auto me : faulty_detMEs) {
        if (me->kind() != MonitorElement::Kind::INT)
          continue;
        if (me->getIntValue() == 0)
          continue;
        uint32_t detId = atoi(me->getName().c_str());
        std::pair<std::string, int32_t> det_layer_pair = folder_organizer.GetSubDetAndLayer(detId, tTopo, false);
        if (abs(det_layer_pair.second) == ilayer + 1)
          nfaulty_layer++;
      }

      nfaulty_subdet += nfaulty_layer;
      float fraction_layer = -1.0;
      if (ndet_layer > 0)
        fraction_layer = 1 - ((nfaulty_layer * 1.0) / ndet_layer);
      if (SiStripCertificationSummaryMap)
        SiStripCertificationSummaryMap->Fill(xbin, ilayer + 1, fraction_layer);
    }
    if (ybin <= SiStripCertificationSummaryMap->getNbinsY()) {
      for (int k = ybin + 1; k <= SiStripCertificationSummaryMap->getNbinsY(); k++)
        SiStripCertificationSummaryMap->Fill(xbin, k, -1.0);
    }
    float fraction_subdet = -1.0;
    if (ndet_subdet > 0)
      fraction_subdet = 1 - ((nfaulty_subdet * 1.0) / ndet_subdet);
    // Check S/N status flag and use the minimum between the two
    std::string full_path = mechanical_dir.substr(0, mechanical_dir.find_last_of('/')) +
                            "/EventInfo/reportSummaryContents/SiStrip_SToNFlag_" + name;
    MonitorElement* me_ston = dqm_store.get(full_path);
    me->Reset();
    if (me_ston && me_ston->kind() == MonitorElement::Kind::REAL) {
      float ston_flg = me_ston->getFloatValue();
      sToNTot += ston_flg;
      nSToNTot++;
      me->Fill(fminf(fraction_subdet, ston_flg));
    } else
      me->Fill(fraction_subdet);
    nDetTot += ndet_subdet;
    nFaultyTot += nfaulty_subdet;
  }
  float fraction_global = -1.0;
  if (nDetTot > 0)
    fraction_global = 1.0 - ((nFaultyTot * 1.0) / nDetTot);
  float ston_frac_global = 1.0;
  if (nSToNTot > 0)
    ston_frac_global = sToNTot / nSToNTot;
  SiStripCertification->Fill(fminf(fraction_global, ston_frac_global));
}
//
// --Fill SiStrip Certification
//
void SiStripCertificationInfo::resetSiStripCertificationMEs(DQMStore& dqm_store) {
  if (!sistripCertificationBooked_) {
    bookSiStripCertificationMEs(dqm_store);
  }
  assert(sistripCertificationBooked_);

  SiStripCertification->Reset();
  for (auto const& pr : SubDetMEsMap) {
    pr.second.det_fractionME->Reset();
  }
  SiStripCertificationSummaryMap->Reset();
}
//
// -- Fill Dummy SiStrip Certification
//
void SiStripCertificationInfo::fillDummySiStripCertification(DQMStore& dqm_store) {
  resetSiStripCertificationMEs(dqm_store);
  if (sistripCertificationBooked_) {
    SiStripCertification->Fill(-1.0);
    for (auto const& pr : SubDetMEsMap) {
      pr.second.det_fractionME->Reset();
      pr.second.det_fractionME->Fill(-1.0);
    }

    for (int xbin = 1; xbin < SiStripCertificationSummaryMap->getNbinsX() + 1; xbin++) {
      for (int ybin = 1; ybin < SiStripCertificationSummaryMap->getNbinsY() + 1; ybin++) {
        SiStripCertificationSummaryMap->Fill(xbin, ybin, -1.0);
      }
    }
  }
}
//
// --Fill SiStrip Certification
//
void SiStripCertificationInfo::fillSiStripCertificationMEsAtLumi(DQMStore& dqm_store) {
  if (!sistripCertificationBooked_) {
    edm::LogError("SiStripCertificationInfo")
        << " SiStripCertificationInfo::fillSiStripCertificationMEsAtLumi : MEs missing ";
    return;
  }
  resetSiStripCertificationMEs(dqm_store);
  dqm_store.cd();
  std::string strip_dir = "";
  SiStripUtility::getTopFolderPath(dqm_store, "SiStrip", strip_dir);
  if (strip_dir.empty())
    strip_dir = "SiStrip";

  std::string full_path;
  float dcs_flag = 1.0;
  float dqm_flag = 1.0;
  for (auto const& [type, subDetME] : SubDetMEsMap) {
    full_path = strip_dir + "/EventInfo/DCSContents/SiStrip_" + type;
    MonitorElement* me_dcs = dqm_store.get(full_path);
    if (me_dcs && me_dcs->kind() == MonitorElement::Kind::REAL)
      dcs_flag = me_dcs->getFloatValue();
    full_path = strip_dir + "/EventInfo/reportSummaryContents/SiStrip_" + type;
    MonitorElement* me_dqm = dqm_store.get(full_path);
    if (me_dqm && me_dqm->kind() == MonitorElement::Kind::REAL)
      dqm_flag = me_dqm->getFloatValue();
    subDetME.det_fractionME->Reset();
    subDetME.det_fractionME->Fill(fminf(dqm_flag, dcs_flag));
  }
  dcs_flag = 1.0;
  dqm_flag = 1.0;
  full_path = strip_dir + "/EventInfo/reportSummary";
  MonitorElement* me_dqm = dqm_store.get(full_path);
  if (me_dqm && me_dqm->kind() == MonitorElement::Kind::REAL)
    dqm_flag = me_dqm->getFloatValue();
  full_path = strip_dir + "/EventInfo/DCSSummary";
  MonitorElement* me_dcs = dqm_store.get(full_path);
  if (me_dcs && me_dcs->kind() == MonitorElement::Kind::REAL)
    dcs_flag = me_dcs->getFloatValue();
  SiStripCertification->Reset();
  SiStripCertification->Fill(fminf(dqm_flag, dcs_flag));
}
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(SiStripCertificationInfo);