CalibratedElectronProducerT< T > Class Template Reference

Inheritance diagram for CalibratedElectronProducerT< T >:

Public Member Functions
	CalibratedElectronProducerT (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~CalibratedElectronProducerT () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	setSemiDetRandomSeed (const edm::Event &iEvent, const T &obj, size_t nrObjs, size_t objNr)

Private Attributes
edm::EDGetTokenT< edm::View< T > >	electronToken_
ElectronEnergyCalibrator	energyCorrector_
EpCombinationTool	epCombinationTool_
bool	produceCalibratedObjs_
edm::EDGetTokenT< EcalRecHitCollection >	recHitCollectionEBToken_
edm::EDGetTokenT< EcalRecHitCollection >	recHitCollectionEEToken_
std::unique_ptr< TRandom >	semiDeterministicRng_

Static Private Attributes
static const std::vector< int >	valMapsToStore_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

template<typename T>
class CalibratedElectronProducerT< T >

Definition at line 37 of file CalibratedElectronProducers.cc.

Constructor & Destructor Documentation

template<typename T >

CalibratedElectronProducerT< T >::CalibratedElectronProducerT ( const edm::ParameterSet &  conf )

explicit

Definition at line 85 of file CalibratedElectronProducers.cc.

References CalibratedElectronProducerT< T >::energyCorrector_, edm::ParameterSet::getParameter(), ElectronEnergyCalibrator::initPrivateRng(), EGEnergySysIndex::name(), CalibratedElectronProducerT< T >::produceCalibratedObjs_, CalibratedElectronProducerT< T >::semiDeterministicRng_, ElectronEnergyCalibrator::setMinEt(), and CalibratedElectronProducerT< T >::valMapsToStore_.

    : electronToken_(consumes<edm::View<T>>(conf.getParameter<edm::InputTag>("src"))),
      epCombinationTool_(conf.getParameter<edm::ParameterSet>("epCombConfig")),
      energyCorrector_(epCombinationTool_, conf.getParameter<std::string>("correctionFile")),
      recHitCollectionEBToken_(consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitCollectionEB"))),
      recHitCollectionEEToken_(consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitCollectionEE"))),
      produceCalibratedObjs_(conf.getParameter<bool>("produceCalibratedObjs")) {
  energyCorrector_.setMinEt(conf.getParameter<double>("minEtToCalibrate"));

  if (conf.getParameter<bool>("semiDeterministic")) {
    semiDeterministicRng_.reset(new TRandom2());
    energyCorrector_.initPrivateRng(semiDeterministicRng_.get());
  }

  if (produceCalibratedObjs_)
    produces<std::vector<T>>();

  for (const auto& toStore : valMapsToStore_) {
    produces<edm::ValueMap<float>>(EGEnergySysIndex::name(toStore));
  }
}

template<typename T >

CalibratedElectronProducerT< T >::~CalibratedElectronProducerT ( )

inlineoverride

Definition at line 40 of file CalibratedElectronProducers.cc.

References CalibratedElectronProducerT< T >::fillDescriptions(), iEvent, getGTfromDQMFile::obj, CalibratedElectronProducerT< T >::produce(), and CalibratedElectronProducerT< T >::setSemiDetRandomSeed().

{}

Member Function Documentation

template<typename T >

void CalibratedElectronProducerT< T >::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 108 of file CalibratedElectronProducers.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), HLT_2018_cff::InputTag, EpCombinationTool::makePSetDescription(), EGEnergySysIndex::name(), AlCaHLTBitMon_QueryRunRegistry::string, and CalibratedElectronProducerT< T >::valMapsToStore_.

Referenced by CalibratedElectronProducerT< T >::~CalibratedElectronProducerT().

                                                                                                {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("src", edm::InputTag("gedPhotons"));
  desc.add<edm::ParameterSetDescription>("epCombConfig", EpCombinationTool::makePSetDescription());
  desc.add<edm::InputTag>("recHitCollectionEB", edm::InputTag("reducedEcalRecHitsEB"));
  desc.add<edm::InputTag>("recHitCollectionEE", edm::InputTag("reducedEcalRecHitsEE"));
  desc.add<std::string>("correctionFile", std::string());
  desc.add<double>("minEtToCalibrate", 5.0);
  desc.add<bool>("produceCalibratedObjs", true);
  desc.add<bool>("semiDeterministic", true);
  std::vector<std::string> valMapsProduced;
  for (auto varToStore : valMapsToStore_)
    valMapsProduced.push_back(EGEnergySysIndex::name(varToStore));
  desc.add<std::vector<std::string>>("valueMapsStored", valMapsProduced)
      ->setComment(
          "provides to python configs the list of valuemaps stored, can not be overriden in the python config");
  descriptions.addWithDefaultLabel(desc);
}

template<typename T >

void CalibratedElectronProducerT< T >::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 128 of file CalibratedElectronProducers.cc.

References ElectronEnergyCalibrator::calibrate(), ElectronEnergyCalibrator::DATA, CalibratedElectronProducerT< T >::electronToken_, CalibratedElectronProducerT< T >::energyCorrector_, CalibratedElectronProducerT< T >::epCombinationTool_, edm::Event::getByToken(), patZpeak::handle, edm::EventBase::id(), edm::EventBase::isRealData(), EGEnergySysIndex::kNrSysErrs, ElectronEnergyCalibrator::MC, eostools::move(), EGEnergySysIndex::name(), MillePedeFileConverter_cfg::out, CalibratedElectronProducerT< T >::produceCalibratedObjs_, edm::Handle< T >::product(), edm::Event::put(), CalibratedElectronProducerT< T >::recHitCollectionEBToken_, CalibratedElectronProducerT< T >::recHitCollectionEEToken_, FastTrackerRecHitMaskProducer_cfi::recHits, bookConverter::results, edm::EventID::run(), CalibratedElectronProducerT< T >::semiDeterministicRng_, EpCombinationTool::setEventContent(), CalibratedElectronProducerT< T >::setSemiDetRandomSeed(), edm::Event::size(), edm::Event::streamID(), and CalibratedElectronProducerT< T >::valMapsToStore_.

Referenced by CalibratedElectronProducerT< T >::~CalibratedElectronProducerT().

                                                                                          {
  epCombinationTool_.setEventContent(iSetup);

  edm::Handle<edm::View<T>> inHandle;
  iEvent.getByToken(electronToken_, inHandle);

  edm::Handle<EcalRecHitCollection> recHitCollectionEBHandle;
  edm::Handle<EcalRecHitCollection> recHitCollectionEEHandle;

  iEvent.getByToken(recHitCollectionEBToken_, recHitCollectionEBHandle);
  iEvent.getByToken(recHitCollectionEEToken_, recHitCollectionEEHandle);

  std::unique_ptr<std::vector<T>> out = std::make_unique<std::vector<T>>();

  size_t nrObj = inHandle->size();
  std::array<std::vector<float>, EGEnergySysIndex::kNrSysErrs> results;
  for (auto& res : results)
    res.reserve(nrObj);

  const ElectronEnergyCalibrator::EventType evtType =
      iEvent.isRealData() ? ElectronEnergyCalibrator::EventType::DATA : ElectronEnergyCalibrator::EventType::MC;

  for (const auto& ele : *inHandle) {
    out->push_back(ele);

    if (semiDeterministicRng_)
      setSemiDetRandomSeed(iEvent, ele, nrObj, out->size());

    const EcalRecHitCollection* recHits =
        (ele.isEB()) ? recHitCollectionEBHandle.product() : recHitCollectionEEHandle.product();
    std::array<float, EGEnergySysIndex::kNrSysErrs> uncertainties =
        energyCorrector_.calibrate(out->back(), iEvent.id().run(), recHits, iEvent.streamID(), evtType);

    for (size_t index = 0; index < EGEnergySysIndex::kNrSysErrs; index++) {
      results[index].push_back(uncertainties[index]);
    }
  }

  auto fillAndStore = [&](auto handle) {
    for (const auto& mapToStore : valMapsToStore_) {
      fillAndStoreValueMap(iEvent, handle, results[mapToStore], EGEnergySysIndex::name(mapToStore));
    }
  };

  if (produceCalibratedObjs_) {
    fillAndStore(iEvent.put(std::move(out)));
  } else {
    fillAndStore(inHandle);
  }
}

template<typename T >

void CalibratedElectronProducerT< T >::setSemiDetRandomSeed ( const edm::Event &  iEvent, const T &  obj, size_t  nrObjs, size_t  objNr  )

private

Definition at line 180 of file CalibratedElectronProducers.cc.

References egamma::getRandomSeedFromObj(), egamma::getRandomSeedFromSC(), and CalibratedElectronProducerT< T >::semiDeterministicRng_.

Referenced by CalibratedElectronProducerT< T >::produce(), and CalibratedElectronProducerT< T >::~CalibratedElectronProducerT().

                                                                        {
  if (obj.superCluster().isNonnull()) {
    semiDeterministicRng_->SetSeed(egamma::getRandomSeedFromSC(iEvent, obj.superCluster()));
  } else {
    semiDeterministicRng_->SetSeed(egamma::getRandomSeedFromObj(iEvent, obj, nrObjs, objNr));
  }
}

Member Data Documentation

template<typename T >

edm::EDGetTokenT<edm::View<T> > CalibratedElectronProducerT< T >::electronToken_

private

Definition at line 47 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::produce().

template<typename T >

ElectronEnergyCalibrator CalibratedElectronProducerT< T >::energyCorrector_

private

Definition at line 50 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::CalibratedElectronProducerT(), and CalibratedElectronProducerT< T >::produce().

template<typename T >

EpCombinationTool CalibratedElectronProducerT< T >::epCombinationTool_

private

Definition at line 49 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::produce().

template<typename T >

bool CalibratedElectronProducerT< T >::produceCalibratedObjs_

private

Definition at line 54 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::CalibratedElectronProducerT(), and CalibratedElectronProducerT< T >::produce().

template<typename T >

edm::EDGetTokenT<EcalRecHitCollection> CalibratedElectronProducerT< T >::recHitCollectionEBToken_

private

Definition at line 52 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::produce().

template<typename T >

edm::EDGetTokenT<EcalRecHitCollection> CalibratedElectronProducerT< T >::recHitCollectionEEToken_

private

Definition at line 53 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::produce().

template<typename T >

std::unique_ptr<TRandom> CalibratedElectronProducerT< T >::semiDeterministicRng_

private

Definition at line 51 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::CalibratedElectronProducerT(), CalibratedElectronProducerT< T >::produce(), and CalibratedElectronProducerT< T >::setSemiDetRandomSeed().

template<typename T >

const std::vector< int > CalibratedElectronProducerT< T >::valMapsToStore_

staticprivate

Initial value:

= {
    EGEnergySysIndex::kScaleStatUp,       EGEnergySysIndex::kScaleStatDown,     EGEnergySysIndex::kScaleSystUp,
    EGEnergySysIndex::kScaleSystDown,     EGEnergySysIndex::kScaleGainUp,       EGEnergySysIndex::kScaleGainDown,
    EGEnergySysIndex::kSmearRhoUp,        EGEnergySysIndex::kSmearRhoDown,      EGEnergySysIndex::kSmearPhiUp,
    EGEnergySysIndex::kSmearPhiDown,      EGEnergySysIndex::kScaleUp,           EGEnergySysIndex::kScaleDown,
    EGEnergySysIndex::kSmearUp,           EGEnergySysIndex::kSmearDown,         EGEnergySysIndex::kScaleValue,
    EGEnergySysIndex::kSmearValue,        EGEnergySysIndex::kSmearNrSigma,      EGEnergySysIndex::kEcalPreCorr,
    EGEnergySysIndex::kEcalErrPreCorr,    EGEnergySysIndex::kEcalPostCorr,      EGEnergySysIndex::kEcalErrPostCorr,
    EGEnergySysIndex::kEcalTrkPreCorr,    EGEnergySysIndex::kEcalTrkErrPreCorr, EGEnergySysIndex::kEcalTrkPostCorr,
    EGEnergySysIndex::kEcalTrkErrPostCorr}

Definition at line 55 of file CalibratedElectronProducers.cc.

Referenced by CalibratedElectronProducerT< T >::CalibratedElectronProducerT(), CalibratedElectronProducerT< T >::fillDescriptions(), and CalibratedElectronProducerT< T >::produce().