CMS 3D CMS Logo
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::ElectronSeed > dPhi1Para_
 
egPM::Param< reco::ElectronSeed > dPhi2Para_
 
egPM::Param< reco::ElectronSeed > dRZ2Para_
 
size_t nrHits_
 
const edm::EDGetTokenT< reco::ElectronSeedCollection > pixelSeedsToken_
 
int productsToWrite_
 
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandidateToken_
 

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

Definition at line 102 of file EgammaHLTPixelMatchVarProducer.cc.

Constructor & Destructor Documentation

◆ EgammaHLTPixelMatchVarProducer()

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

Definition at line 125 of file EgammaHLTPixelMatchVarProducer.cc.

126  : recoEcalCandidateToken_(
127  consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"))),
128  pixelSeedsToken_(
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 }

References nrHits_, and productsToWrite_.

◆ ~EgammaHLTPixelMatchVarProducer()

EgammaHLTPixelMatchVarProducer::~EgammaHLTPixelMatchVarProducer ( )
override

Definition at line 159 of file EgammaHLTPixelMatchVarProducer.cc.

159 {}

Member Function Documentation

◆ calS2()

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

Definition at line 304 of file EgammaHLTPixelMatchVarProducer.cc.

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 }

References ALCARECOTkAlJpsiMuMu_cff::charge, dPhi1Para_, dPhi2Para_, dRZ2Para_, indexGen::s2, and SurveyInfoScenario_cff::seed.

Referenced by produce().

◆ fillDescriptions()

void EgammaHLTPixelMatchVarProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 161 of file EgammaHLTPixelMatchVarProducer.cc.

161  {
162  edm::ParameterSetDescription desc;
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 }

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

◆ produce()

void EgammaHLTPixelMatchVarProducer::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
override

Definition at line 198 of file EgammaHLTPixelMatchVarProducer.cc.

198  {
199  // Get the HLT filtered objects
200  edm::Handle<reco::RecoEcalCandidateCollection> recoEcalCandHandle;
201  iEvent.getByToken(recoEcalCandidateToken_, recoEcalCandHandle);
202 
203  edm::Handle<reco::ElectronSeedCollection> pixelSeedsHandle;
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(),
242  std::numeric_limits<float>::max(),
243  std::numeric_limits<float>::max(),
244  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 }

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

Member Data Documentation

◆ dPhi1Para_

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

Definition at line 117 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

◆ dPhi2Para_

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

Definition at line 118 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

◆ dRZ2Para_

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

Definition at line 119 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

◆ nrHits_

size_t EgammaHLTPixelMatchVarProducer::nrHits_
private

Definition at line 122 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by EgammaHLTPixelMatchVarProducer(), and produce().

◆ pixelSeedsToken_

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

Definition at line 115 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by produce().

◆ productsToWrite_

int EgammaHLTPixelMatchVarProducer::productsToWrite_
private

Definition at line 121 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by EgammaHLTPixelMatchVarProducer(), and produce().

◆ recoEcalCandidateToken_

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

Definition at line 114 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by produce().

ParameterDescription
EgammaHLTPixelMatchVarProducer::pixelSeedsToken_
const edm::EDGetTokenT< reco::ElectronSeedCollection > pixelSeedsToken_
Definition: EgammaHLTPixelMatchVarProducer.cc:115
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
edm
HLT enums.
Definition: AlignableModifier.h:19
funcType
funcType
Definition: SiPixelActionExecutor.h:21
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
photonAnalyzer_cfi.yMin
yMin
Definition: photonAnalyzer_cfi.py:86
indexGen.s2
s2
Definition: indexGen.py:107
EgammaHLTPixelMatchVarProducer::productsToWrite_
int productsToWrite_
Definition: EgammaHLTPixelMatchVarProducer.cc:121
EgammaHLTPixelMatchVarProducer::nrHits_
size_t nrHits_
Definition: EgammaHLTPixelMatchVarProducer.cc:122
edm::Handle
Definition: AssociativeIterator.h:50
edm::Ref
Definition: AssociativeIterator.h:58
photonAnalyzer_cfi.xMin
xMin
Definition: photonAnalyzer_cfi.py:82
config
Definition: config.py:1
EgammaHLTPixelMatchVarProducer::recoEcalCandidateToken_
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection > recoEcalCandidateToken_
Definition: EgammaHLTPixelMatchVarProducer.cc:114
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
edm::ESHandle< TrackerTopology >
multiplicitycorr_cfi.yMax
yMax
Definition: multiplicitycorr_cfi.py:6
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
funct::true
true
Definition: Factorize.h:173
edm::ParameterSet
Definition: ParameterSet.h:36
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
PixelData
Definition: EgammaHLTPixelMatchVarProducer.cc:55
iEvent
int iEvent
Definition: GenABIO.cc:224
EgammaHLTPixelMatchVarProducer::dPhi1Para_
egPM::Param< reco::ElectronSeed > dPhi1Para_
Definition: EgammaHLTPixelMatchVarProducer.cc:117
reco::ElectronSeed::dPhiBest
float dPhiBest(size_t hitNr) const
Definition: ElectronSeed.h:105
edm::RefToBase::castTo
REF castTo() const
Definition: RefToBase.h:257
get
#define get
EgammaHLTPixelMatchVarProducer::dPhi2Para_
egPM::Param< reco::ElectronSeed > dPhi2Para_
Definition: EgammaHLTPixelMatchVarProducer.cc:118
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
or
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
multiplicitycorr_cfi.xMax
xMax
Definition: multiplicitycorr_cfi.py:5
HLT_2018_cff.funcParams
funcParams
Definition: HLT_2018_cff.py:13908
edm::RefToBase
Definition: AssociativeIterator.h:54
EgammaHLTPixelMatchVarProducer::calS2
std::array< float, 4 > calS2(const reco::ElectronSeed &seed, int charge) const
Definition: EgammaHLTPixelMatchVarProducer.cc:304
reco::ElectronSeed::dRZBest
float dRZBest(size_t hitNr) const
Definition: ElectronSeed.h:108
TrackerTopologyRcd
Definition: TrackerTopologyRcd.h:10
EgammaHLTPixelMatchVarProducer::dRZ2Para_
egPM::Param< reco::ElectronSeed > dRZ2Para_
Definition: EgammaHLTPixelMatchVarProducer.cc:119
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
edm::ParameterDescription
Definition: ParameterDescription.h:110
edm::InputTag
Definition: InputTag.h:15
SurveyInfoScenario_cff.seed
seed
Definition: SurveyInfoScenario_cff.py:295
Generated for CMSSW Reference Manual by doxygen 1.8.16