SiStripPopConBadComponentsDQM.cc

Go to the documentation of this file.

#include "DQMOffline/CalibTracker/plugins/SiStripDQMPopConSourceHandler.h"
#include "CondFormats/SiStripObjects/interface/SiStripBadStrip.h"

class SiStripPopConBadComponentsHandlerFromDQM : public SiStripDQMPopConSourceHandler<SiStripBadStrip> {
public:
  typedef dqm::legacy::MonitorElement MonitorElement;
  typedef dqm::legacy::DQMStore DQMStore;

  explicit SiStripPopConBadComponentsHandlerFromDQM(const edm::ParameterSet& iConfig);
  ~SiStripPopConBadComponentsHandlerFromDQM() override;
  // interface methods: implemented in template
  void initES(const edm::EventSetup&) override;
  void dqmEndJob(DQMStore::IBooker& booker, DQMStore::IGetter& getter) override;
  SiStripBadStrip* getObj() const override;

protected:
  std::string getMetaDataString() const override;

private:
  edm::FileInPath fp_;
  SiStripBadStrip m_obj;
  const TrackerTopology* trackerTopo_;
};

#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"

SiStripPopConBadComponentsHandlerFromDQM::SiStripPopConBadComponentsHandlerFromDQM(const edm::ParameterSet& iConfig)
    : SiStripDQMPopConSourceHandler<SiStripBadStrip>(iConfig),
      fp_{iConfig.getUntrackedParameter<edm::FileInPath>(
          "file", edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"))} {
  edm::LogInfo("SiStripBadComponentsDQMService") << "[SiStripBadComponentsDQMService::SiStripBadComponentsDQMService]";
}

SiStripPopConBadComponentsHandlerFromDQM::~SiStripPopConBadComponentsHandlerFromDQM() {
  edm::LogInfo("SiStripBadComponentsDQMService") << "[SiStripBadComponentsDQMService::~SiStripBadComponentsDQMService]";
}

void SiStripPopConBadComponentsHandlerFromDQM::initES(const edm::EventSetup& setup) {
  edm::ESHandle<TrackerTopology> tTopo;
  setup.get<TrackerTopologyRcd>().get(tTopo);
  trackerTopo_ = tTopo.product();
}

std::string SiStripPopConBadComponentsHandlerFromDQM::getMetaDataString() const {
  std::stringstream ss;
  ss << SiStripDQMPopConSourceHandler<SiStripBadStrip>::getMetaDataString();
  getObj()->printSummary(ss, trackerTopo_);
  return ss.str();
}

namespace {
  void getModuleFolderList(SiStripPopConBadComponentsHandlerFromDQM::DQMStore::IGetter& getter,
                           const std::string& pwd,
                           std::vector<std::string>& mfolders) {
    if (std::string::npos != pwd.find("module_")) {
      //    std::string mId = pwd.substr(pwd.find("module_")+7, 9);
      mfolders.push_back(pwd);
    } else {
      for (const auto& subdir : getter.getSubdirs()) {
        getter.cd(subdir);
        getModuleFolderList(getter, subdir, mfolders);
        getter.cd();
        getter.cd(pwd);
      }
    }
  }
}  // namespace

void SiStripPopConBadComponentsHandlerFromDQM::dqmEndJob(DQMStore::IBooker&, DQMStore::IGetter& getter) {
  //*LOOP OVER THE LIST OF SUMMARY OBJECTS TO INSERT IN DB*//

  m_obj = SiStripBadStrip();

  SiStripDetInfoFileReader reader(fp_.fullPath());

  getter.cd();

  const std::string mechanicalview_dir = "MechanicalView";
  if (!getter.dirExists(mechanicalview_dir))
    return;

  const std::vector<std::string> subdet_folder = {"TIB", "TOB", "TEC/side_1", "TEC/side_2", "TID/side_1", "TID/side_2"};

  int nDetsTotal = 0;
  int nDetsWithErrorTotal = 0;
  for (const auto& im : subdet_folder) {
    const std::string dname = mechanicalview_dir + "/" + im;
    getter.cd();
    if (!getter.dirExists(dname))
      continue;
    getter.cd(dname);

    std::vector<std::string> module_folders;
    getModuleFolderList(getter, dname, module_folders);
    int nDets = module_folders.size();

    int nDetsWithError = 0;
    const std::string bad_module_folder = dname + "/" + "BadModuleList";
    getter.cd();
    if (getter.dirExists(bad_module_folder)) {
      for (const MonitorElement* me : getter.getContents(bad_module_folder)) {
        nDetsWithError++;
        std::cout << me->getName() << " " << me->getIntValue() << std::endl;
        uint32_t detId = std::stoul(me->getName());
        short flag = me->getIntValue();

        std::vector<unsigned int> theSiStripVector;

        unsigned short firstBadStrip = 0, NconsecutiveBadStrips = 0;
        unsigned int theBadStripRange;

        // for(std::vector<uint32_t>::const_iterator is=BadApvList_.begin(); is!=BadApvList_.end(); ++is){

        //   firstBadStrip=(*is)*128;
        NconsecutiveBadStrips = reader.getNumberOfApvsAndStripLength(detId).first * 128;

        theBadStripRange = m_obj.encode(firstBadStrip, NconsecutiveBadStrips, flag);

        LogDebug("SiStripBadComponentsDQMService") << "detid " << detId << " \t"
                                                   << ", flag " << flag << std::endl;

        theSiStripVector.push_back(theBadStripRange);
        // }

        SiStripBadStrip::Range range(theSiStripVector.begin(), theSiStripVector.end());
        if (!m_obj.put(detId, range)) {
          edm::LogError("SiStripBadFiberBuilder")
              << "[SiStripBadFiberBuilder::analyze] detid already exists" << std::endl;
        }
      }
    }
    nDetsTotal += nDets;
    nDetsWithErrorTotal += nDetsWithError;
  }
  getter.cd();
}

SiStripBadStrip* SiStripPopConBadComponentsHandlerFromDQM::getObj() const { return new SiStripBadStrip(m_obj); }

#include "FWCore/Framework/interface/MakerMacros.h"
#include "DQMOffline/CalibTracker/plugins/SiStripPopConDQMEDHarvester.h"
using SiStripPopConBadComponentsDQM = SiStripPopConDQMEDHarvester<SiStripPopConBadComponentsHandlerFromDQM>;
DEFINE_FWK_MODULE(SiStripPopConBadComponentsDQM);