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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

template<typename T>
class CalibratedElectronProducerT< T >

Definition at line 40 of file CalibratedElectronProducers.cc.

Constructor & Destructor Documentation

template<typename T >

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

explicit

Definition at line 88 of file CalibratedElectronProducers.cc.

References edm::ParameterSet::getParameterSet().

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

  if (conf.getParameter<bool>("semiDeterministic")) {
    semiDeterministicRng_ = std::make_unique<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 43 of file CalibratedElectronProducers.cc.

{}

Member Function Documentation

template<typename T >

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

static

Definition at line 111 of file CalibratedElectronProducers.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, EpCombinationTool::makePSetDescription(), EGEnergySysIndex::name(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                {
  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;
  valMapsProduced.reserve(valMapsToStore_.size());
  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 132 of file CalibratedElectronProducers.cc.

References ElectronEnergyCalibrator::DATA, edm::Event::getHandle(), patZpeak::handle, edm::EventBase::id(), edm::EventBase::isRealData(), EGEnergySysIndex::kNrSysErrs, ElectronEnergyCalibrator::MC, eostools::move(), EGEnergySysIndex::name(), submitPVResolutionJobs::out, edm::Event::put(), HLT_FULL_cff::recHits, bookConverter::results, edm::EventID::run(), edm::Event::size(), and edm::Event::streamID().

                                                                                          {
  epCombinationTool_.setEventContent(iSetup);

  auto inHandle = iEvent.getHandle(electronToken_);

  auto recHitCollectionEBHandle = iEvent.getHandle(recHitCollectionEBToken_);
  auto recHitCollectionEEHandle = iEvent.getHandle(recHitCollectionEEToken_);

  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(), and egamma::getRandomSeedFromSC().

                                                                        {
  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 50 of file CalibratedElectronProducers.cc.

template<typename T >

ElectronEnergyCalibrator CalibratedElectronProducerT< T >::energyCorrector_

private

Definition at line 53 of file CalibratedElectronProducers.cc.

template<typename T >

EpCombinationTool CalibratedElectronProducerT< T >::epCombinationTool_

private

Definition at line 52 of file CalibratedElectronProducers.cc.

template<typename T >

bool CalibratedElectronProducerT< T >::produceCalibratedObjs_

private

Definition at line 57 of file CalibratedElectronProducers.cc.

template<typename T >

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

private

Definition at line 55 of file CalibratedElectronProducers.cc.

template<typename T >

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

private

Definition at line 56 of file CalibratedElectronProducers.cc.

template<typename T >

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

private

Definition at line 54 of file CalibratedElectronProducers.cc.

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 58 of file CalibratedElectronProducers.cc.