SiStripPopConNoisesDQM.cc

Go to the documentation of this file.

#include "DQMOffline/CalibTracker/plugins/SiStripDQMPopConSourceHandler.h"
#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
class SiStripPopConNoisesHandlerFromDQM : public SiStripDQMPopConSourceHandler<SiStripNoises> {
public:
  typedef dqm::legacy::MonitorElement MonitorElement;
  typedef dqm::legacy::DQMStore DQMStore;

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

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

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

SiStripPopConNoisesHandlerFromDQM::SiStripPopConNoisesHandlerFromDQM(const edm::ParameterSet& iConfig)
    : SiStripDQMPopConSourceHandler<SiStripNoises>(iConfig),
      fp_{iConfig.getUntrackedParameter<edm::FileInPath>(
          "file", edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"))},
      MEDir_{iConfig.getUntrackedParameter<std::string>("ME_DIR", "DQMData")} {
  edm::LogInfo("SiStripNoisesDQMService") << "[SiStripNoisesDQMService::SiStripNoisesDQMService]";
}

SiStripPopConNoisesHandlerFromDQM::~SiStripPopConNoisesHandlerFromDQM() {
  edm::LogInfo("SiStripNoisesDQMService") << "[SiStripNoisesDQMService::~SiStripNoisesDQMService]";
}

void SiStripPopConNoisesHandlerFromDQM::dqmEndJob(DQMStore::IBooker&, DQMStore::IGetter& getter) {
  std::cout << "SiStripNoisesDQMService::readNoises" << std::endl;

  m_obj = SiStripNoises();

  SiStripDetInfoFileReader reader(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("CMSubNoisePerStrip__det__");
                           }),
            MEs.end());

  // The histograms are one per DetId, loop on all the DetIds and extract the corresponding histogram
  for (const auto& detInfo : reader.getAllData()) {
    SiStripNoises::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("CMSubNoisePerStrip__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 != nullptr) {
        // Read the noise 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 50 (we want a high noise to avoid these modules being considered good by mistake)
    if (theSiStripVector.empty()) {
      for (unsigned short j = 0; j < 128 * detInfo.second.nApvs; ++j) {
        m_obj.setData(50, theSiStripVector);
      }
    }

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

SiStripNoises* SiStripPopConNoisesHandlerFromDQM::getObj() const { return new SiStripNoises(m_obj); }

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