SiStripGainESProducer.cc

Go to the documentation of this file.

// system include files
#include <memory>
#include <utility>
 
// user include files
#include "FWCore/Framework/interface/ModuleFactory.h"
#include "FWCore/Framework/interface/ESProducer.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "CalibFormats/SiStripObjects/interface/SiStripGain.h"
#include "CondFormats/SiStripObjects/interface/SiStripApvGain.h"
#include "CondFormats/DataRecord/interface/SiStripCondDataRecords.h"
#include "CalibTracker/Records/interface/SiStripDependentRecords.h"
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
 
#include "SiStripGainFactor.h"
 
//
// class declaration
//
class SiStripGainESProducer : public edm::ESProducer {
public:
  SiStripGainESProducer(const edm::ParameterSet&);
  ~SiStripGainESProducer() override{};
 
  std::unique_ptr<SiStripGain> produce(const SiStripGainRcd&);
 
private:
  class GainGetter {
  public:
    GainGetter(std::string record, std::string label) : recordLabel_{std::move(record), std::move(label)} {}
    virtual ~GainGetter() = default;
    virtual const SiStripApvGain& gain(const SiStripGainRcd& rcd) const = 0;
 
    const auto& recordLabel() const { return recordLabel_; }
 
  private:
    std::pair<std::string, std::string> recordLabel_;
  };
 
  template <typename Record>
  class GainGetterT : public GainGetter {
  public:
    GainGetterT(edm::ESConsumesCollector& cc, std::string record, std::string label)
        : GainGetter(std::move(record), std::move(label)),
          token_{cc.consumesFrom<SiStripApvGain, Record>(edm::ESInputTag{"", recordLabel().second})} {}
 
    const SiStripApvGain& gain(const SiStripGainRcd& rcd) const override { return rcd.get(token_); }
 
  private:
    edm::ESGetToken<SiStripApvGain, Record> token_;
  };
 
  template <typename Record>
  auto make_GainGetter(edm::ESConsumesCollector& cc, std::string record, std::string label) {
    return std::make_unique<GainGetterT<Record>>(cc, std::move(record), std::move(label));
  }
 
  std::vector<std::unique_ptr<GainGetter>> gainGetters_;
 
  SiStripGainFactor factor_;
};
 
SiStripGainESProducer::SiStripGainESProducer(const edm::ParameterSet& iConfig) : factor_{iConfig} {
  auto cc = setWhatProduced(this);
 
  auto apvGainLabels = iConfig.getParameter<std::vector<edm::ParameterSet>>("APVGain");
  if (apvGainLabels.empty()) {
    throw cms::Exception("Configuration") << "Got empty APVGain vector, but need at least one entry";
  }
 
  // Fill the vector of apv labels
  for (const auto& gainPSet : apvGainLabels) {
    // Shouldn't all these parameters be tracked?
    auto record = gainPSet.getParameter<std::string>("Record");
    auto label = gainPSet.getUntrackedParameter<std::string>("Label", "");
    if (record == "SiStripApvGainRcd")
      gainGetters_.emplace_back(make_GainGetter<SiStripApvGainRcd>(cc, record, label));
    else if (record == "SiStripApvGain2Rcd")
      gainGetters_.emplace_back(make_GainGetter<SiStripApvGain2Rcd>(cc, record, label));
    else if (record == "SiStripApvGain3Rcd")
      gainGetters_.emplace_back(make_GainGetter<SiStripApvGain3Rcd>(cc, record, label));
    else
      throw cms::Exception("Configuration")
          << "SiStripGainESProducer::ctor ERROR: unrecognized record name " << record << std::endl
          << "please specify one of: SiStripApvGainRcd, SiStripApvGain2Rcd, SiStripApvGain3Rcd";
    factor_.push_back_norm(gainPSet.getUntrackedParameter<double>("NormalizationFactor", 1.));
  }
  factor_.resetIfBadNorm();
}
 
std::unique_ptr<SiStripGain> SiStripGainESProducer::produce(const SiStripGainRcd& iRecord) {
  const auto detInfo =
      SiStripDetInfoFileReader::read(edm::FileInPath{SiStripDetInfoFileReader::kDefaultFile}.fullPath());
 
  const auto& apvGain = gainGetters_[0]->gain(iRecord);
  // Create a new gain object and insert the ApvGain
  auto gain = std::make_unique<SiStripGain>(apvGain, factor_.get(apvGain, 0), gainGetters_[0]->recordLabel(), detInfo);
 
  for (unsigned int i = 1; i < gainGetters_.size(); ++i) {
    const auto& apvGain = gainGetters_[i]->gain(iRecord);
    // Add the new ApvGain to the gain object
    gain->multiply(apvGain, factor_.get(apvGain, i), gainGetters_[i]->recordLabel(), detInfo);
  }
 
  return gain;
}
 
#include "FWCore/Framework/interface/ModuleFactory.h"
 
DEFINE_FWK_EVENTSETUP_MODULE(SiStripGainESProducer);