SiStripBadComponentInfo.cc

Go to the documentation of this file.

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

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include "CalibTracker/Records/interface/SiStripQualityRcd.h"
#include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
#include "DQM/SiStripMonitorClient/interface/SiStripUtility.h"
#include "DQM/SiStripMonitorClient/plugins/SiStripBadComponentInfo.h"

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

//
// -- Contructor
//
SiStripBadComponentInfo::SiStripBadComponentInfo(edm::ParameterSet const& pSet) : 
    bookedStatus_(false),
    nSubSystem_(6),
    qualityLabel_(pSet.getParameter<std::string>("StripQualityLabel"))
{ 
  // Create MessageSender
  LogDebug( "SiStripBadComponentInfo") << "SiStripBadComponentInfo::Deleting SiStripBadComponentInfo ";
}

SiStripBadComponentInfo::~SiStripBadComponentInfo()
{
  LogDebug("SiStripBadComponentInfo")
    << "SiStripBadComponentInfo::Deleting SiStripBadComponentInfo ";
}

//
// -- Read Condition
//
void
SiStripBadComponentInfo::checkBadComponents(edm::EventSetup const& eSetup)
{
  LogDebug("SiStripBadComponentInfo")
    << "SiStripBadComponentInfo:: Begining of Run";

  // Retrieve tracker topology from geometry
  eSetup.get<TrackerTopologyRcd>().get(tTopoHandle_);
  const TrackerTopology* const topo = tTopoHandle_.product();

 
  eSetup.get<SiStripQualityRcd>().get(qualityLabel_,siStripQuality_);
 
  auto const& BC = siStripQuality_->getBadComponentList();
  
  for (size_t i=0;i<BC.size();++i){
    int subdet=-999; int component=-999;

    //&&&&&&&&&&&&&&&&&
    // Single SubSyste
    //&&&&&&&&&&&&&&&&&
    int subDet = DetId(BC[i].detid).subdetId();
    if (subDet == StripSubdetector::TIB) {
      //&&&&&&&&&&&&&&&&&
      // TIB
      //&&&&&&&&&&&&&&&&&
      subdet = 3;
      component = topo->tibLayer(BC[i].detid);
    } else if (subDet == StripSubdetector::TID) {
      //&&&&&&&&&&&&&&&&&
      // TID
      //&&&&&&&&&&&&&&&&&
      if (topo->tidSide(BC[i].detid) == 2)
        subdet = 4;
      else
        subdet = 5;
      component = topo->tidWheel(BC[i].detid);
    } else if (subDet == StripSubdetector::TOB) {
      //&&&&&&&&&&&&&&&&&
      // TOB
      //&&&&&&&&&&&&&&&&&
      subdet = 6;
      component = topo->tobLayer(BC[i].detid);
    } else if (subDet == StripSubdetector::TEC) {
      //&&&&&&&&&&&&&&&&&
      // TEC
      //&&&&&&&&&&&&&&&&&
      if (topo->tecSide(BC[i].detid) == 2)
        subdet = 2;
      else
        subdet = 1;
      component = topo->tecWheel(BC[i].detid);
    }
    fillBadComponentMaps(subdet, component, BC[i]);
  }

  //&&&&&&&&&&&&&&&&&&
  // Single Strip Info
  //&&&&&&&&&&&&&&&&&&

  SiStripQuality::RegistryIterator rbegin =
    siStripQuality_->getRegistryVectorBegin();
  SiStripQuality::RegistryIterator rend =
    siStripQuality_->getRegistryVectorEnd();

  for (SiStripBadStrip::RegistryIterator rp = rbegin; rp != rend; ++rp) {
    uint32_t detid = rp->detid;
    int subdet = -999;
    int component = -999;
    DetId detectorId = DetId(detid);
    int subDet = detectorId.subdetId();
    if (subDet == StripSubdetector::TIB) {
      subdet = 3;
      component = topo->tibLayer(detid);
    } else if (subDet == StripSubdetector::TID) {
      if (topo->tidSide(detid) == 2)
        subdet = 5;
      else
        subdet = 4;
      component = topo->tidWheel(detid);
    } else if (subDet == StripSubdetector::TOB) {
      subdet = 6;
      component = topo->tobLayer(detid);
    } else if (subDet == StripSubdetector::TEC) {
      if (topo->tecSide(detid) == 2)
        subdet = 2;
      else
        subdet = 1;
      component = topo->tecWheel(detid);
    }

    SiStripQuality::Range sqrange =
      SiStripQuality::Range(siStripQuality_->getDataVectorBegin() + rp->ibegin,
                            siStripQuality_->getDataVectorBegin() + rp->iend);

    for (int it = 0; it < sqrange.second - sqrange.first; it++) {
      unsigned int range = siStripQuality_->decode(*(sqrange.first + it)).range;
      float val = (mapBadStrip.find(std::make_pair(subdet, component)) !=
                   mapBadStrip.end()) ?
                    mapBadStrip.at(std::make_pair(subdet, component)) :
                    0.;
      val += range;
      mapBadStrip[std::make_pair(subdet, component)] = val;
    }
  }
}

void
SiStripBadComponentInfo::endRun(edm::Run const& run,
                                edm::EventSetup const& eSetup)
{
  LogDebug("SiStripBadComponentInfo") << "SiStripBadComponentInfo:: End of Run";
  checkBadComponents(eSetup);
}

void
SiStripBadComponentInfo::dqmEndJob(DQMStore::IBooker& ibooker,
                                   DQMStore::IGetter& igetter)
{
  LogDebug("SiStripBadComponentInfo") << "SiStripBadComponentInfo::dqmEndRun";
  bookBadComponentHistos(ibooker, igetter);
  createSummary(badAPVME_, mapBadAPV);
  createSummary(badFiberME_, mapBadFiber);
  createSummary(badStripME_, mapBadStrip);
}
//
// -- Book MEs for SiStrip Dcs Fraction
//
void
SiStripBadComponentInfo::bookBadComponentHistos(DQMStore::IBooker& ibooker,
                                                DQMStore::IGetter& igetter)
{
  if (!bookedStatus_) {
    std::string strip_dir = "";
    ibooker.cd();
    //    SiStripUtility::getTopFolderPath(dqmStore_, "SiStrip", strip_dir);
    if (igetter.dirExists("SiStrip")) {
      ibooker.cd("SiStrip");
      strip_dir = ibooker.pwd();
    }
    edm::LogInfo("SiStripBadComponentInfo")
      << "SiStripBadComponentInfo::bookBadComponentHistos ==> " << strip_dir
      << " " << ibooker.pwd() << std::endl;
    if (!strip_dir.empty())
      ibooker.setCurrentFolder(strip_dir + "/EventInfo");
    else
      ibooker.setCurrentFolder("SiStrip/EventInfo");

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

    std::string hname, htitle;
    hname = "BadAPVMap";
    htitle = "SiStrip Bad APVs";
    badAPVME_ = ibooker.book2D(
      hname, htitle, nSubSystem_, 0.5, nSubSystem_ + 0.5, 10, 0.5, 10.5);
    badAPVME_->setAxisTitle("Sub Detector Type", 1);
    badAPVME_->setAxisTitle("Layer/Disc Number", 2);

    hname = "BadFiberMap";
    htitle = "SiStrip Bad Fibers";
    badFiberME_ = ibooker.book2D(
      hname, htitle, nSubSystem_, 0.5, nSubSystem_ + 0.5, 10, 0.5, 10.5);
    badFiberME_->setAxisTitle("Sub Detector Type", 1);
    badFiberME_->setAxisTitle("Layer/Disc Number", 2);

    hname = "BadStripMap";
    htitle = "SiStrip Bad Strips";
    badStripME_ = ibooker.book2D(
      hname, htitle, nSubSystem_, 0.5, nSubSystem_ + 0.5, 10, 0.5, 10.5);
    badStripME_->setAxisTitle("Sub Detector Type", 1);
    badStripME_->setAxisTitle("Layer/Disc Number", 2);

    std::vector<std::string> names;
    names.push_back("TECB");
    names.push_back("TECF");
    names.push_back("TIB");
    names.push_back("TIDB");
    names.push_back("TIDF");
    names.push_back("TOB");

    for (unsigned int i = 0; i < names.size(); i++) {
      badAPVME_->setBinLabel(i + 1, names[i]);
      badFiberME_->setBinLabel(i + 1, names[i]);
      badStripME_->setBinLabel(i + 1, names[i]);
    }

    bookedStatus_ = true;
    ibooker.cd();
  }
}

void SiStripBadComponentInfo::fillBadComponentMaps(int xbin,int component,SiStripQuality::BadComponent const& BC){

  auto index = std::make_pair(xbin, component);

  if (BC.BadApvs) {
    int ntot = std::bitset<16>(BC.BadApvs & 0x3f).count();
    float val =
      (mapBadAPV.find(index) != mapBadAPV.end()) ? mapBadAPV.at(index) : 0.;
    val += ntot;
    mapBadAPV[index] = val;
  }
  if (BC.BadFibers) {
    int ntot = std::bitset<16>(BC.BadFibers & 0x7).count();
    float val = (mapBadFiber.find(index) != mapBadFiber.end()) ?
                  mapBadFiber.at(index) :
                  0.;
    val += ntot;
    mapBadFiber[index] = val;
  }
}
void
SiStripBadComponentInfo::createSummary(
  MonitorElement* me,
  const std::map<std::pair<int, int>, float>& map)
{
  for (int i = 1; i < nSubSystem_ + 1; i++) {
    float sum = 0.0;
    for (int k = 1; k < me->getNbinsY(); k++) {
      auto index = std::make_pair(i, k);
      if (map.find(index) != map.end()) {
        me->setBinContent(i, k, map.at(index)); // fill the layer/wheel bins
        sum += map.at(index);
      }
    }
    me->setBinContent(
      i, me->getNbinsY(), sum); // fill the summary bin (last one)
  }
}
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(SiStripBadComponentInfo);