SiStripApvGainFakeESSource.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    CalibTracker/SiStripESProducers
// Class:      SiStripApvGainFakeESSource
//
// system include files
#include <memory>
 
// user include files
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
 
#include "CondFormats/SiStripObjects/interface/SiStripApvGain.h"
#include "CondFormats/DataRecord/interface/SiStripCondDataRecords.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
 
class SiStripApvGainFakeESSource : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
public:
  SiStripApvGainFakeESSource(const edm::ParameterSet&);
  ~SiStripApvGainFakeESSource() override;
 
  void setIntervalFor(const edm::eventsetup::EventSetupRecordKey&,
                      const edm::IOVSyncValue& iov,
                      edm::ValidityInterval& iValidity) override;
 
  typedef std::unique_ptr<SiStripApvGain> ReturnType;
  ReturnType produce(const SiStripApvGainRcd&);
 
private:
  std::string m_genMode;
  double m_meanGain;
  double m_sigmaGain;
  double m_minimumPosValue;
  uint32_t m_printDebug;
  SiStripDetInfo m_detInfo;
};
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CLHEP/Random/RandGauss.h"
 
SiStripApvGainFakeESSource::SiStripApvGainFakeESSource(const edm::ParameterSet& iConfig) {
  setWhatProduced(this);
  findingRecord<SiStripApvGainRcd>();
 
  m_genMode = iConfig.getParameter<std::string>("genMode");
  m_meanGain = iConfig.getParameter<double>("MeanGain");
  m_sigmaGain = iConfig.getParameter<double>("SigmaGain");
  m_minimumPosValue = iConfig.getParameter<double>("MinPositiveGain");
  m_printDebug = iConfig.getUntrackedParameter<uint32_t>("printDebug", 5);
  m_detInfo = SiStripDetInfoFileReader::read(iConfig.getParameter<edm::FileInPath>("SiStripDetInfoFile").fullPath());
}
 
SiStripApvGainFakeESSource::~SiStripApvGainFakeESSource() {}
 
void SiStripApvGainFakeESSource::setIntervalFor(const edm::eventsetup::EventSetupRecordKey&,
                                                const edm::IOVSyncValue& iov,
                                                edm::ValidityInterval& iValidity) {
  iValidity = edm::ValidityInterval{iov.beginOfTime(), iov.endOfTime()};
}
 
// ------------ method called to produce the data  ------------
SiStripApvGainFakeESSource::ReturnType SiStripApvGainFakeESSource::produce(const SiStripApvGainRcd& iRecord) {
  using namespace edm::es;
 
  auto apvGain = std::make_unique<SiStripApvGain>();
 
  uint32_t count{0};
  for (const auto& elm : m_detInfo.getAllData()) {
    std::vector<float> theSiStripVector;
    for (unsigned short j = 0; j < elm.second.nApvs; ++j) {
      float gainValue;
      if (m_genMode == "default") {
        gainValue = m_meanGain;
      } else if (m_genMode == "gaussian") {
        gainValue = CLHEP::RandGauss::shoot(m_meanGain, m_sigmaGain);
        if (gainValue <= m_minimumPosValue) {
          gainValue = m_minimumPosValue;
        }
      } else {
        LogDebug("SiStripApvGain") << "ERROR: wrong genMode specifier : " << m_genMode
                                   << ", please select one of \"default\" or \"gaussian\"";
        exit(1);
      }
 
      if (count < m_printDebug) {
        edm::LogInfo("SiStripApvGainGenerator") << "detid: " << elm.first << " Apv: " << j << " gain: " << gainValue;
      }
      theSiStripVector.push_back(gainValue);
    }
    ++count;
 
    if (!apvGain->put(elm.first, SiStripApvGain::Range{theSiStripVector.begin(), theSiStripVector.end()})) {
      edm::LogError("SiStripApvGainGenerator") << " detid already exists";
    }
  }
 
  return apvGain;
}
 
//define this as a plug-in
#include "FWCore/Framework/interface/SourceFactory.h"
DEFINE_FWK_EVENTSETUP_SOURCE(SiStripApvGainFakeESSource);