SiStripPopConPedestalsDQM.cc

Go to the documentation of this file.

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

class SiStripPopConPedestalsHandlerFromDQM : public SiStripDQMPopConSourceHandler<SiStripPedestals> {
public:
  typedef dqm::legacy::MonitorElement MonitorElement;
  typedef dqm::legacy::DQMStore DQMStore;

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

private:
  edm::FileInPath fp_;
  std::string MEDir_;
  SiStripPedestals m_obj;
};

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

SiStripPopConPedestalsHandlerFromDQM::SiStripPopConPedestalsHandlerFromDQM(const edm::ParameterSet& iConfig,
                                                                           edm::ConsumesCollector&&)
    : SiStripDQMPopConSourceHandler<SiStripPedestals>(iConfig),
      fp_{iConfig.getUntrackedParameter<edm::FileInPath>("file",
                                                         edm::FileInPath(SiStripDetInfoFileReader::kDefaultFile))},
      MEDir_{iConfig.getUntrackedParameter<std::string>("ME_DIR", "DQMData")} {
  edm::LogInfo("SiStripPedestalsDQMService") << "[SiStripPedestalsDQMService::SiStripPedestalsDQMService]";
}

SiStripPopConPedestalsHandlerFromDQM::~SiStripPopConPedestalsHandlerFromDQM() {
  edm::LogInfo("SiStripPedestalsDQMService") << "[SiStripPedestalsDQMService::~SiStripPedestalsDQMService]";
}

void SiStripPopConPedestalsHandlerFromDQM::dqmEndJob(DQMStore::IBooker&, DQMStore::IGetter& getter) {
  std::cout << "SiStripPedestalsDQMService::readPedestals" << std::endl;

  m_obj = SiStripPedestals();

  const auto detInfo = SiStripDetInfoFileReader::read(fp_.fullPath());

  // getter.cd(iConfig_.getUntrackedParameter<std::string>("ME_DIR"));
  getter.cd();

  uint32_t stripsPerApv = 128;

  // Get the full list of monitoring elements
  // const std::vector<MonitorElement*>& MEs = getter.getAllContents(iConfig_.getUntrackedParameter<std::string>("ME_DIR","DQMData"));

  // Take a copy of the vector
  std::vector<MonitorElement*> MEs = getter.getAllContents(MEDir_);
  // Remove all but the MEs we are using
  MEs.erase(std::remove_if(MEs.begin(),
                           MEs.end(),
                           [](const MonitorElement* ME) -> bool {
                             return std::string::npos == ME->getName().find("PedsPerStrip__det__");
                           }),
            MEs.end());

  // The histograms are one per DetId, loop on all the DetIds and extract the corresponding histogram
  for (const auto& detInfo : detInfo.getAllData()) {
    SiStripPedestals::InputVector theSiStripVector;

    // Take the path for each DetId and build the complete path + histogram name

    // MonitorElement * mE = getModuleHistogram(detInfo.first, "PedsPerStrip");
    const MonitorElement* mE{nullptr};
    std::string MEname("PedsPerStrip__det__" + std::to_string(detInfo.first));
    for (const MonitorElement* ime : MEs) {
      if (ime->getName() == MEname) {
        mE = ime;
        break;
      }
    }

    // find( MEs.begin(), MEs.end(), "PedsPerStrip__det__"+std::to_string(detInfo.first), findMEbyName() );
    // MonitorElement * mE = *(find( MEs.begin(), MEs.end(), findMEbyName("PedsPerStrip__det__"+std::to_string(detInfo.first)) ));

    if (mE) {
      TH1F* histo = mE->getTH1F();
      if (histo) {
        // Read the pedestals from the histograms
        uint32_t nBinsX = histo->GetXaxis()->GetNbins();

        if (nBinsX != stripsPerApv * (detInfo.second.nApvs)) {
          std::cout << "ERROR: number of bin = " << nBinsX
                    << " != number of strips = " << stripsPerApv * (detInfo.second.nApvs) << std::endl;
        }

        // std::cout << "Bin 0 = " << histo->GetBinContent(0) << std::endl;
        // TH1 bins start from 1, 0 is the underflow, nBinsX+1 the overflow.
        for (uint32_t iBin = 1; iBin <= nBinsX; ++iBin) {
          // encode the pedestal value and put it in the vector (push_back)
          m_obj.setData(histo->GetBinContent(iBin), theSiStripVector);
        }
      } else {
        std::cout << "ERROR: histo = " << histo << std::endl;
      }
    } else {
      std::cout << "ERROR: ME = " << mE << std::endl;
    }
    // If the ME was absent fill the vector with 0
    if (theSiStripVector.empty()) {
      for (unsigned short j = 0; j < 128 * detInfo.second.nApvs; ++j) {
        m_obj.setData(0, theSiStripVector);
      }
    }

    if (!m_obj.put(detInfo.first, theSiStripVector))
      edm::LogError("SiStripPedestalsFakeESSource::produce ") << " detid already exists" << std::endl;
  }
  getter.cd();
}

SiStripPedestals* SiStripPopConPedestalsHandlerFromDQM::getObj() const { return new SiStripPedestals(m_obj); }

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