CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
CalibratedPhotonProducerT< T > Class Template Reference
Inheritance diagram for CalibratedPhotonProducerT< T >:
edm::stream::EDProducer<>

Public Member Functions

 CalibratedPhotonProducerT (const edm::ParameterSet &)
 
void produce (edm::Event &, const edm::EventSetup &) override
 
 ~CalibratedPhotonProducerT () 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

PhotonEnergyCalibrator energyCorrector_
 
edm::EDGetTokenT< edm::View< T > > photonToken_
 
bool produceCalibratedObjs_
 
edm::EDGetTokenT< EcalRecHitCollectionrecHitCollectionEBToken_
 
edm::EDGetTokenT< EcalRecHitCollectionrecHitCollectionEEToken_
 
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, GlobalCacheLuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCacheRunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 

Detailed Description

template<typename T>
class CalibratedPhotonProducerT< T >

Definition at line 32 of file CalibratedPhotonProducers.cc.

Constructor & Destructor Documentation

template<typename T >
CalibratedPhotonProducerT< T >::CalibratedPhotonProducerT ( const edm::ParameterSet conf)
explicit

Definition at line 92 of file CalibratedPhotonProducers.cc.

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

92  :
93  photonToken_(consumes<edm::View<T> >(conf.getParameter<edm::InputTag>("src"))),
94  energyCorrector_(conf.getParameter<std::string >("correctionFile")),
95  recHitCollectionEBToken_(consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitCollectionEB"))),
96  recHitCollectionEEToken_(consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitCollectionEE"))),
97  produceCalibratedObjs_(conf.getParameter<bool>("produceCalibratedObjs"))
98 {
99 
100  energyCorrector_.setMinEt(conf.getParameter<double>("minEtToCalibrate"));
101 
102  if (conf.getParameter<bool>("semiDeterministic")) {
103  semiDeterministicRng_.reset(new TRandom2());
105  }
106 
107  if(produceCalibratedObjs_) produces<std::vector<T>>();
108 
109  for(const auto& toStore : valMapsToStore_){
110  produces<edm::ValueMap<float>>(EGEnergySysIndex::name(toStore));
111  }
112 }
T getParameter(std::string const &) const
std::unique_ptr< TRandom > semiDeterministicRng_
void initPrivateRng(TRandom *rnd)
static const std::string & name(size_t index)
PhotonEnergyCalibrator energyCorrector_
edm::EDGetTokenT< EcalRecHitCollection > recHitCollectionEEToken_
edm::EDGetTokenT< EcalRecHitCollection > recHitCollectionEBToken_
edm::EDGetTokenT< edm::View< T > > photonToken_
static const std::vector< int > valMapsToStore_
template<typename T >
CalibratedPhotonProducerT< T >::~CalibratedPhotonProducerT ( )
inlineoverride

Member Function Documentation

template<typename T >
void CalibratedPhotonProducerT< T >::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 115 of file CalibratedPhotonProducers.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), EGEnergySysIndex::name(), AlCaHLTBitMon_QueryRunRegistry::string, and CalibratedPhotonProducerT< T >::valMapsToStore_.

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

116 {
118  desc.add<edm::InputTag>("src",edm::InputTag("gedGsfElectrons"));
119  desc.add<edm::InputTag>("recHitCollectionEB",edm::InputTag("reducedEcalRecHitsEB"));
120  desc.add<edm::InputTag>("recHitCollectionEE",edm::InputTag("reducedEcalRecHitsEE"));
121  desc.add<std::string>("correctionFile",std::string());
122  desc.add<double>("minEtToCalibrate",5.0);
123  desc.add<bool>("produceCalibratedObjs",true);
124  desc.add<bool>("semiDeterministic",true);
125  std::vector<std::string> valMapsProduced;
126  for(auto varToStore : valMapsToStore_) valMapsProduced.push_back(EGEnergySysIndex::name(varToStore));
127  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");
128  descriptions.addWithDefaultLabel(desc);
129 }
void addWithDefaultLabel(ParameterSetDescription const &psetDescription)
static const std::string & name(size_t index)
ParameterDescriptionBase * add(U const &iLabel, T const &value)
static const std::vector< int > valMapsToStore_
template<typename T >
void CalibratedPhotonProducerT< T >::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
override

Definition at line 133 of file CalibratedPhotonProducers.cc.

References PhotonEnergyCalibrator::calibrate(), PhotonEnergyCalibrator::DATA, CalibratedPhotonProducerT< T >::energyCorrector_, edm::Event::getByToken(), cmsBatch::handle, edm::EventBase::id(), edm::EventBase::isRealData(), EGEnergySysIndex::kNrSysErrs, PhotonEnergyCalibrator::MC, eostools::move(), EGEnergySysIndex::name(), MillePedeFileConverter_cfg::out, CalibratedPhotonProducerT< T >::photonToken_, CalibratedPhotonProducerT< T >::produceCalibratedObjs_, edm::Handle< T >::product(), edm::Event::put(), CalibratedPhotonProducerT< T >::recHitCollectionEBToken_, CalibratedPhotonProducerT< T >::recHitCollectionEEToken_, edm::EventID::run(), CalibratedPhotonProducerT< T >::semiDeterministicRng_, CalibratedPhotonProducerT< T >::setSemiDetRandomSeed(), edm::Event::size(), edm::Event::streamID(), and CalibratedPhotonProducerT< T >::valMapsToStore_.

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

133  {
134 
135  edm::Handle<edm::View<T>> inHandle;
136  iEvent.getByToken(photonToken_, inHandle);
137 
138  edm::Handle<EcalRecHitCollection> recHitCollectionEBHandle;
139  edm::Handle<EcalRecHitCollection> recHitCollectionEEHandle;
140 
141  iEvent.getByToken(recHitCollectionEBToken_, recHitCollectionEBHandle);
142  iEvent.getByToken(recHitCollectionEEToken_, recHitCollectionEEHandle);
143 
144  std::unique_ptr<std::vector<T>> out = std::make_unique<std::vector<T>>();
145 
146  size_t nrObj = inHandle->size();
147  std::array<std::vector<float>,EGEnergySysIndex::kNrSysErrs> results;
148  for(auto& res : results) res.reserve(nrObj);
149 
150  const PhotonEnergyCalibrator::EventType evtType = iEvent.isRealData() ?
152 
153  for (const auto& pho : *inHandle) {
154  out->emplace_back(pho);
155 
156  if(semiDeterministicRng_) setSemiDetRandomSeed(iEvent,pho,nrObj,out->size());
157 
158  const EcalRecHitCollection* recHits = (pho.isEB()) ? recHitCollectionEBHandle.product() : recHitCollectionEEHandle.product();
159  std::array<float,EGEnergySysIndex::kNrSysErrs> uncertainties = energyCorrector_.calibrate(out->back(), iEvent.id().run(), recHits, iEvent.streamID(), evtType);
160 
162  results[index].push_back(uncertainties[index]);
163  }
164  }
165 
166  auto fillAndStore = [&](auto handle){
167  for(const auto& mapToStore : valMapsToStore_){
168  fillAndStoreValueMap(iEvent,handle,results[mapToStore],EGEnergySysIndex::name(mapToStore));
169  }
170  };
171 
173  fillAndStore(iEvent.put(std::move(out)));
174  }else{
175  fillAndStore(inHandle);
176  }
177 
178 }
RunNumber_t run() const
Definition: EventID.h:39
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:125
size_t size() const
Definition: Event.cc:233
std::unique_ptr< TRandom > semiDeterministicRng_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
static const std::string & name(size_t index)
bool isRealData() const
Definition: EventBase.h:62
Definition: Electron.h:6
PhotonEnergyCalibrator energyCorrector_
void setSemiDetRandomSeed(const edm::Event &iEvent, const T &obj, size_t nrObjs, size_t objNr)
edm::EDGetTokenT< EcalRecHitCollection > recHitCollectionEEToken_
static constexpr size_t kNrSysErrs
edm::EDGetTokenT< EcalRecHitCollection > recHitCollectionEBToken_
T const * product() const
Definition: Handle.h:74
edm::EDGetTokenT< edm::View< T > > photonToken_
edm::EventID id() const
Definition: EventBase.h:59
StreamID streamID() const
Definition: Event.h:95
static const std::vector< int > valMapsToStore_
def move(src, dest)
Definition: eostools.py:511
std::array< float, EGEnergySysIndex::kNrSysErrs > calibrate(reco::Photon &photon, const unsigned int runNumber, const EcalRecHitCollection *recHits, edm::StreamID const &id, const EventType eventType) const
template<typename T >
void CalibratedPhotonProducerT< T >::setSemiDetRandomSeed ( const edm::Event iEvent,
const T obj,
size_t  nrObjs,
size_t  objNr 
)
private

Definition at line 182 of file CalibratedPhotonProducers.cc.

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

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

183 {
184  if(obj.superCluster().isNonnull()){
185  semiDeterministicRng_->SetSeed(egamma::getRandomSeedFromSC(iEvent,obj.superCluster()));
186  }else{
187  semiDeterministicRng_->SetSeed(egamma::getRandomSeedFromObj(iEvent,obj,nrObjs,objNr));
188  }
189 }
std::unique_ptr< TRandom > semiDeterministicRng_
uint32_t getRandomSeedFromObj(const edm::Event &iEvent, const T &obj, size_t nrObjs, size_t objNr)
uint32_t getRandomSeedFromSC(const edm::Event &iEvent, const reco::SuperClusterRef scRef)

Member Data Documentation

template<typename T >
PhotonEnergyCalibrator CalibratedPhotonProducerT< T >::energyCorrector_
private
template<typename T >
edm::EDGetTokenT<edm::View<T> > CalibratedPhotonProducerT< T >::photonToken_
private
template<typename T >
bool CalibratedPhotonProducerT< T >::produceCalibratedObjs_
private
template<typename T >
edm::EDGetTokenT<EcalRecHitCollection> CalibratedPhotonProducerT< T >::recHitCollectionEBToken_
private
template<typename T >
edm::EDGetTokenT<EcalRecHitCollection> CalibratedPhotonProducerT< T >::recHitCollectionEEToken_
private
template<typename T >
std::unique_ptr<TRandom> CalibratedPhotonProducerT< T >::semiDeterministicRng_
private
template<typename T >
const std::vector< int > CalibratedPhotonProducerT< T >::valMapsToStore_
staticprivate