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:
  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/MonitorElement.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( 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);
      }
    }
  }
}

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);