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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 38 of file CalibratedElectronProducers.cc.

Constructor & Destructor Documentation

template<typename T >

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

explicit

Definition at line 104 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 42 of file CalibratedElectronProducers.cc.

References CalibratedElectronProducerT< T >::fillDescriptions(), iEvent, hgcalPlots::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 127 of file CalibratedElectronProducers.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), 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 146 of file CalibratedElectronProducers.cc.

References ElectronEnergyCalibrator::calibrate(), ElectronEnergyCalibrator::DATA, CalibratedElectronProducerT< T >::electronToken_, CalibratedElectronProducerT< T >::energyCorrector_, CalibratedElectronProducerT< T >::epCombinationTool_, edm::Event::getByToken(), cmsBatch::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_, mps_update::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 198 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 49 of file CalibratedElectronProducers.cc.

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

template<typename T >

ElectronEnergyCalibrator CalibratedElectronProducerT< T >::energyCorrector_

private

Definition at line 52 of file CalibratedElectronProducers.cc.

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

template<typename T >

EpCombinationTool CalibratedElectronProducerT< T >::epCombinationTool_

private

Definition at line 51 of file CalibratedElectronProducers.cc.

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

template<typename T >

bool CalibratedElectronProducerT< T >::produceCalibratedObjs_

private

Definition at line 56 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 54 of file CalibratedElectronProducers.cc.

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

template<typename T >

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

private

Definition at line 55 of file CalibratedElectronProducers.cc.

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

template<typename T >

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

private

Definition at line 53 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 57 of file CalibratedElectronProducers.cc.

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