CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
EgammaHLTPixelMatchVarProducer Class Reference
Inheritance diagram for EgammaHLTPixelMatchVarProducer:
edm::stream::EDProducer<>

Public Member Functions

std::array< float, 4 > calS2 (const reco::ElectronSeed &seed, int charge) const
 
 EgammaHLTPixelMatchVarProducer (const edm::ParameterSet &)
 
void produce (edm::Event &, const edm::EventSetup &) override
 
 ~EgammaHLTPixelMatchVarProducer () 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 Attributes

egPM::Param< reco::ElectronSeeddPhi1Para_
 
egPM::Param< reco::ElectronSeeddPhi2Para_
 
egPM::Param< reco::ElectronSeeddRZ2Para_
 
size_t nrHits_
 
const edm::EDGetTokenT< reco::ElectronSeedCollectionpixelSeedsToken_
 
int productsToWrite_
 
const edm::EDGetTokenT< reco::RecoEcalCandidateCollectionrecoEcalCandidateToken_
 

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

Definition at line 102 of file EgammaHLTPixelMatchVarProducer.cc.

Constructor & Destructor Documentation

EgammaHLTPixelMatchVarProducer::EgammaHLTPixelMatchVarProducer ( const edm::ParameterSet config)
explicit

Definition at line 125 of file EgammaHLTPixelMatchVarProducer.cc.

References nrHits_, and productsToWrite_.

127  consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"))),
129  consumes<reco::ElectronSeedCollection>(config.getParameter<edm::InputTag>("pixelSeedsProducer"))),
130  dPhi1Para_(config.getParameter<edm::ParameterSet>("dPhi1SParams")),
131  dPhi2Para_(config.getParameter<edm::ParameterSet>("dPhi2SParams")),
132  dRZ2Para_(config.getParameter<edm::ParameterSet>("dRZ2SParams")),
133  productsToWrite_(config.getParameter<int>("productsToWrite")),
134  nrHits_(4)
135 
136 {
137  //register your products
138  produces<reco::RecoEcalCandidateIsolationMap>("s2");
139  if (productsToWrite_ >= 1) {
140  produces<reco::RecoEcalCandidateIsolationMap>("dPhi1BestS2");
141  produces<reco::RecoEcalCandidateIsolationMap>("dPhi2BestS2");
142  produces<reco::RecoEcalCandidateIsolationMap>("dzBestS2");
143  }
144  if (productsToWrite_ >= 2) {
145  //note for product names we start from index 1
146  for (size_t hitNr = 1; hitNr <= nrHits_; hitNr++) {
147  produces<reco::RecoEcalCandidateIsolationMap>("dPhi" + std::to_string(hitNr));
148  produces<reco::RecoEcalCandidateIsolationMap>("dPhi" + std::to_string(hitNr) + "Info");
149  produces<reco::RecoEcalCandidateIsolationMap>("dRZ" + std::to_string(hitNr));
150  produces<reco::RecoEcalCandidateIsolationMap>("dRZ" + std::to_string(hitNr) + "Info");
151  }
152  produces<reco::RecoEcalCandidateIsolationMap>("nrClus");
153  produces<reco::RecoEcalCandidateIsolationMap>("seedClusEFrac");
154  produces<reco::RecoEcalCandidateIsolationMap>("phiWidth");
155  produces<reco::RecoEcalCandidateIsolationMap>("etaWidth");
156  }
157 }
T getParameter(std::string const &) const
egPM::Param< reco::ElectronSeed > dPhi2Para_
egPM::Param< reco::ElectronSeed > dPhi1Para_
egPM::Param< reco::ElectronSeed > dRZ2Para_
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandidateToken_
const edm::EDGetTokenT< reco::ElectronSeedCollection > pixelSeedsToken_
EgammaHLTPixelMatchVarProducer::~EgammaHLTPixelMatchVarProducer ( )
override

Definition at line 159 of file EgammaHLTPixelMatchVarProducer.cc.

159 {}

Member Function Documentation

std::array< float, 4 > EgammaHLTPixelMatchVarProducer::calS2 ( const reco::ElectronSeed seed,
int  charge 
) const

Definition at line 304 of file EgammaHLTPixelMatchVarProducer.cc.

References DEFINE_FWK_MODULE, dPhi1Para_, dPhi2Para_, reco::ElectronSeed::dPhiNeg(), reco::ElectronSeed::dPhiPos(), dRZ2Para_, reco::ElectronSeed::dRZNeg(), reco::ElectronSeed::dRZPos(), and indexGen::s2.

Referenced by produce().

304  {
305  const float dPhi1Const = dPhi1Para_(seed);
306  const float dPhi2Const = dPhi2Para_(seed);
307  const float dRZ2Const = dRZ2Para_(seed);
308 
309  float dPhi1 = (charge < 0 ? seed.dPhiNeg(0) : seed.dPhiPos(0)) / dPhi1Const;
310  float dPhi2 = (charge < 0 ? seed.dPhiNeg(1) : seed.dPhiPos(1)) / dPhi2Const;
311  float dRz2 = (charge < 0 ? seed.dRZNeg(1) : seed.dRZPos(1)) / dRZ2Const;
312 
313  float s2 = dPhi1 * dPhi1 + dPhi2 * dPhi2 + dRz2 * dRz2;
314  return std::array<float, 4>{{s2, dPhi1, dPhi2, dRz2}};
315 }
egPM::Param< reco::ElectronSeed > dPhi2Para_
egPM::Param< reco::ElectronSeed > dPhi1Para_
float dRZPos(size_t hitNr) const
Definition: ElectronSeed.h:106
egPM::Param< reco::ElectronSeed > dRZ2Para_
float dRZNeg(size_t hitNr) const
Definition: ElectronSeed.h:107
float dPhiNeg(size_t hitNr) const
Definition: ElectronSeed.h:103
float dPhiPos(size_t hitNr) const
Definition: ElectronSeed.h:104
void EgammaHLTPixelMatchVarProducer::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 161 of file EgammaHLTPixelMatchVarProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::ifValue(), eostools::move(), and or.

161  {
163  desc.add<edm::InputTag>(("recoEcalCandidateProducer"), edm::InputTag("hltL1SeededRecoEcalCandidate"));
164  desc.add<edm::InputTag>(("pixelSeedsProducer"), edm::InputTag("electronPixelSeeds"));
165 
166  edm::ParameterSetDescription varParamDesc;
167  edm::ParameterSetDescription binParamDesc;
168 
169  auto binDescCases = "AbsEtaClus" >> (edm::ParameterDescription<double>("xMin", 0.0, true) and
170  edm::ParameterDescription<double>("xMax", 3.0, true) and
171  edm::ParameterDescription<int>("yMin", 0, true) and
172  edm::ParameterDescription<int>("yMax", 99999, true) and
173  edm::ParameterDescription<std::string>("funcType", "pol0", true) and
174  edm::ParameterDescription<std::vector<double>>("funcParams", {0.}, true)) or
175  "AbsEtaClusPhi" >> (edm::ParameterDescription<double>("xMin", 0.0, true) and
176  edm::ParameterDescription<double>("xMax", 3.0, true) and
177  edm::ParameterDescription<int>("yMin", 0, true) and
178  edm::ParameterDescription<int>("yMax", 99999, true) and
179  edm::ParameterDescription<std::string>("funcType", "pol0", true) and
180  edm::ParameterDescription<std::vector<double>>("funcParams", {0.}, true)) or
181  "AbsEtaClusEt" >> (edm::ParameterDescription<double>("xMin", 0.0, true) and
182  edm::ParameterDescription<double>("xMax", 3.0, true) and
183  edm::ParameterDescription<int>("yMin", 0, true) and
184  edm::ParameterDescription<int>("yMax", 99999, true) and
185  edm::ParameterDescription<std::string>("funcType", "pol0", true) and
186  edm::ParameterDescription<std::vector<double>>("funcParams", {0.}, true));
187 
188  binParamDesc.ifValue(edm::ParameterDescription<std::string>("binType", "AbsEtaClus", true), std::move(binDescCases));
189 
190  varParamDesc.addVPSet("bins", binParamDesc);
191  desc.add("dPhi1SParams", varParamDesc);
192  desc.add("dPhi2SParams", varParamDesc);
193  desc.add("dRZ2SParams", varParamDesc);
194  desc.add<int>("productsToWrite", 0);
195  descriptions.add(("hltEgammaHLTPixelMatchVarProducer"), desc);
196 }
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
ParameterDescriptionBase * add(U const &iLabel, T const &value)
void add(std::string const &label, ParameterSetDescription const &psetDescription)
HLT enums.
def move(src, dest)
Definition: eostools.py:511
void EgammaHLTPixelMatchVarProducer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
override

Definition at line 198 of file EgammaHLTPixelMatchVarProducer.cc.

References reco::ElectronSeed::caloCluster(), calS2(), edm::RefToBase< T >::castTo(), reco::ElectronSeed::dPhiBest(), reco::ElectronSeed::dRZBest(), edm::EventSetup::get(), edm::Event::getByToken(), edm::HandleBase::isValid(), SiStripPI::max, eostools::move(), nrHits_, pixelSeedsToken_, productsToWrite_, edm::Event::put(), and recoEcalCandidateToken_.

198  {
199  // Get the HLT filtered objects
201  iEvent.getByToken(recoEcalCandidateToken_, recoEcalCandHandle);
202 
204  iEvent.getByToken(pixelSeedsToken_, pixelSeedsHandle);
205 
206  if (!recoEcalCandHandle.isValid())
207  return;
208  else if (!pixelSeedsHandle.isValid()) {
209  auto s2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
210  for (unsigned int candNr = 0; candNr < recoEcalCandHandle->size(); candNr++) {
211  reco::RecoEcalCandidateRef candRef(recoEcalCandHandle, candNr);
212  s2Map->insert(candRef, 0);
213  }
214  iEvent.put(std::move(s2Map), "s2");
215  return;
216  }
217 
218  edm::ESHandle<TrackerTopology> trackerTopoHandle;
219  iSetup.get<TrackerTopologyRcd>().get(trackerTopoHandle);
220 
221  auto dPhi1BestS2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
222  auto dPhi2BestS2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
223  auto dzBestS2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
224  auto s2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
225 
226  auto nrClusMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
227  auto seedClusEFracMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
228  auto phiWidthMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
229  auto etaWidthMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
230 
231  std::vector<PixelData> pixelData;
232  for (size_t hitNr = 0; hitNr < nrHits_; hitNr++) {
233  pixelData.emplace_back(PixelData("dPhi", hitNr, &reco::ElectronSeed::dPhiBest, recoEcalCandHandle));
234  pixelData.emplace_back(PixelData("dRZ", hitNr, &reco::ElectronSeed::dRZBest, recoEcalCandHandle));
235  }
236 
237  for (unsigned int candNr = 0; candNr < recoEcalCandHandle->size(); candNr++) {
238  reco::RecoEcalCandidateRef candRef(recoEcalCandHandle, candNr);
239  reco::SuperClusterRef candSCRef = candRef->superCluster();
240 
241  std::array<float, 4> bestS2{{std::numeric_limits<float>::max(),
245  for (auto& seed : *pixelSeedsHandle) {
246  const edm::RefToBase<reco::CaloCluster>& pixelClusterRef = seed.caloCluster();
247  reco::SuperClusterRef pixelSCRef = pixelClusterRef.castTo<reco::SuperClusterRef>();
248  if (&(*candSCRef) == &(*pixelSCRef)) {
249  std::array<float, 4> s2Data = calS2(seed, -1);
250  std::array<float, 4> s2DataPos = calS2(seed, +1);
251  if (s2Data[0] < bestS2[0])
252  bestS2 = s2Data;
253  if (s2DataPos[0] < bestS2[0])
254  bestS2 = s2DataPos;
255 
256  if (productsToWrite_ >= 2) {
257  for (auto& pixelDatum : pixelData) {
258  pixelDatum.fill(seed);
259  }
260  }
261  }
262  }
263 
264  s2Map->insert(candRef, bestS2[0]);
265  if (productsToWrite_ >= 1) {
266  dPhi1BestS2Map->insert(candRef, bestS2[1]);
267  dPhi2BestS2Map->insert(candRef, bestS2[2]);
268  dzBestS2Map->insert(candRef, bestS2[3]);
269  }
270  if (productsToWrite_ >= 2) {
271  nrClusMap->insert(candRef, candSCRef->clustersSize());
272  float seedClusEFrac = candSCRef->rawEnergy() > 0 ? candSCRef->seed()->energy() / candSCRef->rawEnergy() : 0.;
273  // std::cout <<"cand "<<candSCRef->energy()<<" E Corr "<<candSCRef->correctedEnergyUncertainty()<<" "<<candSCRef->correctedEnergy()<<" width "<<candSCRef->phiWidth()<<std::endl;
274  // float seedClusEFrac = candSCRef->phiWidth();
275  seedClusEFracMap->insert(candRef, seedClusEFrac);
276  float phiWidth = candSCRef->phiWidth();
277  float etaWidth = candSCRef->etaWidth();
278  phiWidthMap->insert(candRef, phiWidth);
279  etaWidthMap->insert(candRef, etaWidth);
280 
281  for (auto& pixelDatum : pixelData) {
282  pixelDatum.fill(candRef);
283  }
284  }
285  }
286 
287  iEvent.put(std::move(s2Map), "s2");
288  if (productsToWrite_ >= 1) {
289  iEvent.put(std::move(dPhi1BestS2Map), "dPhi1BestS2");
290  iEvent.put(std::move(dPhi2BestS2Map), "dPhi2BestS2");
291  iEvent.put(std::move(dzBestS2Map), "dzBestS2");
292  }
293  if (productsToWrite_ >= 2) {
294  for (auto& pixelDatum : pixelData) {
295  pixelDatum.putInto(iEvent);
296  }
297  iEvent.put(std::move(nrClusMap), "nrClus");
298  iEvent.put(std::move(seedClusEFracMap), "seedClusEFrac");
299  iEvent.put(std::move(phiWidthMap), "phiWidth");
300  iEvent.put(std::move(etaWidthMap), "etaWidth");
301  }
302 }
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:131
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandidateToken_
bool isValid() const
Definition: HandleBase.h:70
float dPhiBest(size_t hitNr) const
Definition: ElectronSeed.h:105
std::array< float, 4 > calS2(const reco::ElectronSeed &seed, int charge) const
REF castTo() const
Definition: RefToBase.h:257
float dRZBest(size_t hitNr) const
Definition: ElectronSeed.h:108
T get() const
Definition: EventSetup.h:73
const edm::EDGetTokenT< reco::ElectronSeedCollection > pixelSeedsToken_
def move(src, dest)
Definition: eostools.py:511

Member Data Documentation

egPM::Param<reco::ElectronSeed> EgammaHLTPixelMatchVarProducer::dPhi1Para_
private

Definition at line 117 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

egPM::Param<reco::ElectronSeed> EgammaHLTPixelMatchVarProducer::dPhi2Para_
private

Definition at line 118 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

egPM::Param<reco::ElectronSeed> EgammaHLTPixelMatchVarProducer::dRZ2Para_
private

Definition at line 119 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

size_t EgammaHLTPixelMatchVarProducer::nrHits_
private

Definition at line 122 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by EgammaHLTPixelMatchVarProducer(), and produce().

const edm::EDGetTokenT<reco::ElectronSeedCollection> EgammaHLTPixelMatchVarProducer::pixelSeedsToken_
private

Definition at line 115 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by produce().

int EgammaHLTPixelMatchVarProducer::productsToWrite_
private

Definition at line 121 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by EgammaHLTPixelMatchVarProducer(), and produce().

const edm::EDGetTokenT<reco::RecoEcalCandidateCollection> EgammaHLTPixelMatchVarProducer::recoEcalCandidateToken_
private

Definition at line 114 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by produce().