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 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

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 77 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_.

78  : photonToken_(consumes<edm::View<T>>(conf.getParameter<edm::InputTag>("src"))),
79  energyCorrector_(conf.getParameter<std::string>("correctionFile")),
80  recHitCollectionEBToken_(consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitCollectionEB"))),
81  recHitCollectionEEToken_(consumes<EcalRecHitCollection>(conf.getParameter<edm::InputTag>("recHitCollectionEE"))),
82  produceCalibratedObjs_(conf.getParameter<bool>("produceCalibratedObjs")) {
83  energyCorrector_.setMinEt(conf.getParameter<double>("minEtToCalibrate"));
84 
85  if (conf.getParameter<bool>("semiDeterministic")) {
86  semiDeterministicRng_.reset(new TRandom2());
88  }
89 
91  produces<std::vector<T>>();
92 
93  for (const auto& toStore : valMapsToStore_) {
94  produces<edm::ValueMap<float>>(EGEnergySysIndex::name(toStore));
95  }
96 }
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 99 of file CalibratedPhotonProducers.cc.

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

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

99  {
101  desc.add<edm::InputTag>("src", edm::InputTag("gedGsfElectrons"));
102  desc.add<edm::InputTag>("recHitCollectionEB", edm::InputTag("reducedEcalRecHitsEB"));
103  desc.add<edm::InputTag>("recHitCollectionEE", edm::InputTag("reducedEcalRecHitsEE"));
104  desc.add<std::string>("correctionFile", std::string());
105  desc.add<double>("minEtToCalibrate", 5.0);
106  desc.add<bool>("produceCalibratedObjs", true);
107  desc.add<bool>("semiDeterministic", true);
108  std::vector<std::string> valMapsProduced;
109  for (auto varToStore : valMapsToStore_)
110  valMapsProduced.push_back(EGEnergySysIndex::name(varToStore));
111  desc.add<std::vector<std::string>>("valueMapsStored", valMapsProduced)
112  ->setComment(
113  "provides to python configs the list of valuemaps stored, can not be overriden in the python config");
114  descriptions.addWithDefaultLabel(desc);
115 }
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 118 of file CalibratedPhotonProducers.cc.

References PhotonEnergyCalibrator::calibrate(), PhotonEnergyCalibrator::DATA, CalibratedPhotonProducerT< T >::energyCorrector_, edm::Event::getByToken(), patZpeak::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_, FastTrackerRecHitMaskProducer_cfi::recHits, bookConverter::results, edm::EventID::run(), CalibratedPhotonProducerT< T >::semiDeterministicRng_, CalibratedPhotonProducerT< T >::setSemiDetRandomSeed(), edm::Event::size(), edm::Event::streamID(), and CalibratedPhotonProducerT< T >::valMapsToStore_.

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

118  {
119  edm::Handle<edm::View<T>> inHandle;
120  iEvent.getByToken(photonToken_, inHandle);
121 
122  edm::Handle<EcalRecHitCollection> recHitCollectionEBHandle;
123  edm::Handle<EcalRecHitCollection> recHitCollectionEEHandle;
124 
125  iEvent.getByToken(recHitCollectionEBToken_, recHitCollectionEBHandle);
126  iEvent.getByToken(recHitCollectionEEToken_, recHitCollectionEEHandle);
127 
128  std::unique_ptr<std::vector<T>> out = std::make_unique<std::vector<T>>();
129 
130  size_t nrObj = inHandle->size();
131  std::array<std::vector<float>, EGEnergySysIndex::kNrSysErrs> results;
132  for (auto& res : results)
133  res.reserve(nrObj);
134 
135  const PhotonEnergyCalibrator::EventType evtType =
137 
138  for (const auto& pho : *inHandle) {
139  out->emplace_back(pho);
140 
142  setSemiDetRandomSeed(iEvent, pho, nrObj, out->size());
143 
145  (pho.isEB()) ? recHitCollectionEBHandle.product() : recHitCollectionEEHandle.product();
146  std::array<float, EGEnergySysIndex::kNrSysErrs> uncertainties =
147  energyCorrector_.calibrate(out->back(), iEvent.id().run(), recHits, iEvent.streamID(), evtType);
148 
149  for (size_t index = 0; index < EGEnergySysIndex::kNrSysErrs; index++) {
150  results[index].push_back(uncertainties[index]);
151  }
152  }
153 
154  auto fillAndStore = [&](auto handle) {
155  for (const auto& mapToStore : valMapsToStore_) {
156  fillAndStoreValueMap(iEvent, handle, results[mapToStore], EGEnergySysIndex::name(mapToStore));
157  }
158  };
159 
161  fillAndStore(iEvent.put(std::move(out)));
162  } else {
163  fillAndStore(inHandle);
164  }
165 }
RunNumber_t run() const
Definition: EventID.h:38
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:131
size_t size() const
Definition: Event.cc:242
std::unique_ptr< TRandom > semiDeterministicRng_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
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:69
edm::EDGetTokenT< edm::View< T > > photonToken_
edm::EventID id() const
Definition: EventBase.h:59
StreamID streamID() const
Definition: Event.h:96
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 169 of file CalibratedPhotonProducers.cc.

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

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

172  {
173  if (obj.superCluster().isNonnull()) {
174  semiDeterministicRng_->SetSeed(egamma::getRandomSeedFromSC(iEvent, obj.superCluster()));
175  } else {
176  semiDeterministicRng_->SetSeed(egamma::getRandomSeedFromObj(iEvent, obj, nrObjs, objNr));
177  }
178 }
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