CMS 3D CMS Logo
PATElectronProducer.cc
Go to the documentation of this file.
1 //
2 #include "PhysicsTools/PatAlgos/plugins/PATElectronProducer.h"
3 
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 #include "FWCore/ParameterSet/interface/FileInPath.h"
6 #include "FWCore/Utilities/interface/isFinite.h"
7 
8 #include "DataFormats/Common/interface/Association.h"
9 #include "DataFormats/Common/interface/ValueMap.h"
10 #include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
11 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
12 #include "DataFormats/PatCandidates/interface/PFIsolation.h"
13 
14 #include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"
15 #include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
16 
17 #include "PhysicsTools/PatUtils/interface/TrackerIsolationPt.h"
18 #include "PhysicsTools/PatUtils/interface/CaloIsolationEnergy.h"
19 
20 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
21 #include "DataFormats/VertexReco/interface/Vertex.h"
22 
23 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
24 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
25 
26 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
27 #include "TrackingTools/Records/interface/TransientTrackRecord.h"
28 #include "TrackingTools/TransientTrack/interface/TransientTrack.h"
29 #include "TrackingTools/IPTools/interface/IPTools.h"
30 
31 #include "DataFormats/GsfTrackReco/interface/GsfTrackFwd.h"
32 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
33 #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
34 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
35 #include "CommonTools/Egamma/interface/ConversionTools.h"
36 #include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
37 #include "Geometry/CaloTopology/interface/CaloTopology.h"
38 #include "Geometry/Records/interface/CaloGeometryRecord.h"
39 
40 #include "FWCore/ParameterSet/interface/EmptyGroupDescription.h"
41 
42 #include "FWCore/Utilities/interface/transform.h"
43 
44 #include "PhysicsTools/PatUtils/interface/MiniIsolation.h"
45 
46 #include <vector>
47 #include <memory>
48 
49 using namespace pat;
50 using namespace std;
51 
52 PATElectronProducer::PATElectronProducer(const edm::ParameterSet& iConfig)
53  : // general configurables
54  electronToken_(consumes<edm::View<reco::GsfElectron>>(iConfig.getParameter<edm::InputTag>("electronSource"))),
55  hConversionsToken_(consumes<reco::ConversionCollection>(edm::InputTag("allConversions"))),
56  embedGsfElectronCore_(iConfig.getParameter<bool>("embedGsfElectronCore")),
57  embedGsfTrack_(iConfig.getParameter<bool>("embedGsfTrack")),
58  embedSuperCluster_(iConfig.getParameter<bool>("embedSuperCluster")),
59  embedPflowSuperCluster_(iConfig.getParameter<bool>("embedPflowSuperCluster")),
60  embedSeedCluster_(iConfig.getParameter<bool>("embedSeedCluster")),
61  embedBasicClusters_(iConfig.getParameter<bool>("embedBasicClusters")),
62  embedPreshowerClusters_(iConfig.getParameter<bool>("embedPreshowerClusters")),
63  embedPflowBasicClusters_(iConfig.getParameter<bool>("embedPflowBasicClusters")),
64  embedPflowPreshowerClusters_(iConfig.getParameter<bool>("embedPflowPreshowerClusters")),
65  embedTrack_(iConfig.getParameter<bool>("embedTrack")),
66  addGenMatch_(iConfig.getParameter<bool>("addGenMatch")),
67  embedGenMatch_(addGenMatch_ ? iConfig.getParameter<bool>("embedGenMatch") : false),
68  embedRecHits_(iConfig.getParameter<bool>("embedRecHits")),
69  // pflow configurables
70  useParticleFlow_(iConfig.getParameter<bool>("useParticleFlow")),
71  usePfCandidateMultiMap_(iConfig.getParameter<bool>("usePfCandidateMultiMap")),
72  pfElecToken_(!usePfCandidateMultiMap_
73  ? consumes<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfElectronSource"))
74  : edm::EDGetTokenT<reco::PFCandidateCollection>()),
75  pfCandidateMapToken_(!usePfCandidateMultiMap_ ? mayConsume<edm::ValueMap<reco::PFCandidatePtr>>(
76  iConfig.getParameter<edm::InputTag>("pfCandidateMap"))
77  : edm::EDGetTokenT<edm::ValueMap<reco::PFCandidatePtr>>()),
78  pfCandidateMultiMapToken_(usePfCandidateMultiMap_
79  ? consumes<edm::ValueMap<std::vector<reco::PFCandidateRef>>>(
80  iConfig.getParameter<edm::InputTag>("pfCandidateMultiMap"))
81  : edm::EDGetTokenT<edm::ValueMap<std::vector<reco::PFCandidateRef>>>()),
82  embedPFCandidate_(iConfig.getParameter<bool>("embedPFCandidate")),
83  // mva input variables
84  addMVAVariables_(iConfig.getParameter<bool>("addMVAVariables")),
85  reducedBarrelRecHitCollection_(iConfig.getParameter<edm::InputTag>("reducedBarrelRecHitCollection")),
86  reducedBarrelRecHitCollectionToken_(mayConsume<EcalRecHitCollection>(reducedBarrelRecHitCollection_)),
87  reducedEndcapRecHitCollection_(iConfig.getParameter<edm::InputTag>("reducedEndcapRecHitCollection")),
88  reducedEndcapRecHitCollectionToken_(mayConsume<EcalRecHitCollection>(reducedEndcapRecHitCollection_)),
89  ecalClusterToolsESGetTokens_{consumesCollector()},
90  // PFCluster Isolation maps
91  addPFClusterIso_(iConfig.getParameter<bool>("addPFClusterIso")),
92  addPuppiIsolation_(iConfig.getParameter<bool>("addPuppiIsolation")),
93  ecalPFClusterIsoT_(consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("ecalPFClusterIsoMap"))),
94  hcalPFClusterIsoT_(consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("hcalPFClusterIsoMap"))),
95  // embed high level selection variables?
96  embedHighLevelSelection_(iConfig.getParameter<bool>("embedHighLevelSelection")),
97  beamLineToken_(consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamLineSrc"))),
98  pvToken_(mayConsume<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("pvSrc"))),
99  addElecID_(iConfig.getParameter<bool>("addElectronID")),
100  pTComparator_(),
101  isolator_(iConfig.exists("userIsolation") ? iConfig.getParameter<edm::ParameterSet>("userIsolation")
102  : edm::ParameterSet(),
103  consumesCollector(),
104  false),
105  addEfficiencies_(iConfig.getParameter<bool>("addEfficiencies")),
106  addResolutions_(iConfig.getParameter<bool>("addResolutions")),
107  useUserData_(iConfig.exists("userData"))
108 
109 {
110  // MC matching configurables (scheduled mode)
111 
112  if (addGenMatch_) {
113  if (iConfig.existsAs<edm::InputTag>("genParticleMatch")) {
114  genMatchTokens_.push_back(consumes<edm::Association<reco::GenParticleCollection>>(
115  iConfig.getParameter<edm::InputTag>("genParticleMatch")));
116  } else {
117  genMatchTokens_ = edm::vector_transform(
118  iConfig.getParameter<std::vector<edm::InputTag>>("genParticleMatch"),
119  [this](edm::InputTag const& tag) { return consumes<edm::Association<reco::GenParticleCollection>>(tag); });
120  }
121  }
122  // resolution configurables
123  if (addResolutions_) {
124  resolutionLoader_ = pat::helper::KinResolutionsLoader(iConfig.getParameter<edm::ParameterSet>("resolutions"));
125  }
126  if (addPuppiIsolation_) {
127  //puppi
128  PUPPIIsolation_charged_hadrons_ =
129  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationChargedHadrons"));
130  PUPPIIsolation_neutral_hadrons_ =
131  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationNeutralHadrons"));
132  PUPPIIsolation_photons_ =
133  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationPhotons"));
134  //puppiNoLeptons
135  PUPPINoLeptonsIsolation_charged_hadrons_ =
136  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiNoLeptonsIsolationChargedHadrons"));
137  PUPPINoLeptonsIsolation_neutral_hadrons_ =
138  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiNoLeptonsIsolationNeutralHadrons"));
139  PUPPINoLeptonsIsolation_photons_ =
140  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiNoLeptonsIsolationPhotons"));
141  }
142  // electron ID configurables
143  if (addElecID_) {
144  // it might be a single electron ID
145  if (iConfig.existsAs<edm::InputTag>("electronIDSource")) {
146  elecIDSrcs_.push_back(NameTag("", iConfig.getParameter<edm::InputTag>("electronIDSource")));
147  }
148  // or there might be many of them
149  if (iConfig.existsAs<edm::ParameterSet>("electronIDSources")) {
150  // please don't configure me twice
151  if (!elecIDSrcs_.empty()) {
152  throw cms::Exception("Configuration")
153  << "PATElectronProducer: you can't specify both 'electronIDSource' and 'electronIDSources'\n";
154  }
155  // read the different electron ID names
156  edm::ParameterSet idps = iConfig.getParameter<edm::ParameterSet>("electronIDSources");
157  std::vector<std::string> names = idps.getParameterNamesForType<edm::InputTag>();
158  for (std::vector<std::string>::const_iterator it = names.begin(), ed = names.end(); it != ed; ++it) {
159  elecIDSrcs_.push_back(NameTag(*it, idps.getParameter<edm::InputTag>(*it)));
160  }
161  }
162  // but in any case at least once
163  if (elecIDSrcs_.empty()) {
164  throw cms::Exception("Configuration")
165  << "PATElectronProducer: id addElectronID is true, you must specify either:\n"
166  << "\tInputTag electronIDSource = <someTag>\n"
167  << "or\n"
168  << "\tPSet electronIDSources = { \n"
169  << "\t\tInputTag <someName> = <someTag> // as many as you want \n "
170  << "\t}\n";
171  }
172  }
173  elecIDTokens_ = edm::vector_transform(
174  elecIDSrcs_, [this](NameTag const& tag) { return mayConsume<edm::ValueMap<float>>(tag.second); });
175  // construct resolution calculator
176 
177  // // IsoDeposit configurables
178  // if (iConfig.exists("isoDeposits")) {
179  // edm::ParameterSet depconf = iConfig.getParameter<edm::ParameterSet>("isoDeposits");
180  // if (depconf.exists("tracker")) isoDepositLabels_.push_back(std::make_pair(TrackerIso, depconf.getParameter<edm::InputTag>("tracker")));
181  // if (depconf.exists("ecal")) isoDepositLabels_.push_back(std::make_pair(ECalIso, depconf.getParameter<edm::InputTag>("ecal")));
182  // if (depconf.exists("hcal")) isoDepositLabels_.push_back(std::make_pair(HCalIso, depconf.getParameter<edm::InputTag>("hcal")));
183 
184  // if (depconf.exists("user")) {
185  // std::vector<edm::InputTag> userdeps = depconf.getParameter<std::vector<edm::InputTag> >("user");
186  // std::vector<edm::InputTag>::const_iterator it = userdeps.begin(), ed = userdeps.end();
187  // int key = UserBaseIso;
188  // for ( ; it != ed; ++it, ++key) {
189  // isoDepositLabels_.push_back(std::make_pair(IsolationKeys(key), *it));
190  // }
191  // }
192  // }
193  // isoDepositTokens_ = edm::vector_transform(isoDepositLabels_, [this](std::pair<IsolationKeys,edm::InputTag> const & label){return consumes<edm::ValueMap<IsoDeposit> >(label.second);});
194 
195  // for mini-iso
196  computeMiniIso_ = iConfig.getParameter<bool>("computeMiniIso");
197  miniIsoParamsE_ = iConfig.getParameter<std::vector<double>>("miniIsoParamsE");
198  miniIsoParamsB_ = iConfig.getParameter<std::vector<double>>("miniIsoParamsB");
199  if (computeMiniIso_ && (miniIsoParamsE_.size() != 9 || miniIsoParamsB_.size() != 9)) {
200  throw cms::Exception("ParameterError") << "miniIsoParams must have exactly 9 elements.\n";
201  }
202  if (computeMiniIso_)
203  pcToken_ = consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("pfCandsForMiniIso"));
204 
205  // read isoDeposit labels, for direct embedding
206  readIsolationLabels(iConfig, "isoDeposits", isoDepositLabels_, isoDepositTokens_);
207  // read isolation value labels, for direct embedding
208  readIsolationLabels(iConfig, "isolationValues", isolationValueLabels_, isolationValueTokens_);
209  // read isolation value labels for non PF identified electron, for direct embedding
210  readIsolationLabels(iConfig, "isolationValuesNoPFId", isolationValueLabelsNoPFId_, isolationValueNoPFIdTokens_);
211  // Efficiency configurables
212  if (addEfficiencies_) {
213  efficiencyLoader_ =
214  pat::helper::EfficiencyLoader(iConfig.getParameter<edm::ParameterSet>("efficiencies"), consumesCollector());
215  }
216  // Check to see if the user wants to add user data
217  if (useUserData_) {
218  userDataHelper_ =
219  PATUserDataHelper<Electron>(iConfig.getParameter<edm::ParameterSet>("userData"), consumesCollector());
220  }
221 
222  // consistency check
223  if (useParticleFlow_ && usePfCandidateMultiMap_)
224  throw cms::Exception("Configuration", "usePfCandidateMultiMap not supported when useParticleFlow is set to true");
225 
226  // produces vector of muons
227  produces<std::vector<Electron>>();
228 }
229 
230 PATElectronProducer::~PATElectronProducer() {}
231 
232 void PATElectronProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
233  // switch off embedding (in unschedules mode)
234  if (iEvent.isRealData()) {
235  addGenMatch_ = false;
236  embedGenMatch_ = false;
237  }
238 
239  edm::ESHandle<CaloTopology> theCaloTopology;
240  iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
241  ecalTopology_ = &(*theCaloTopology);
242 
243  // Get the collection of electrons from the event
244  edm::Handle<edm::View<reco::GsfElectron>> electrons;
245  iEvent.getByToken(electronToken_, electrons);
246 
247  edm::Handle<PackedCandidateCollection> pc;
248  if (computeMiniIso_)
249  iEvent.getByToken(pcToken_, pc);
250 
251  // for additional mva variables
252  edm::InputTag reducedEBRecHitCollection(string("reducedEcalRecHitsEB"));
253  edm::InputTag reducedEERecHitCollection(string("reducedEcalRecHitsEE"));
254  //EcalClusterLazyTools lazyTools(iEvent, iSetup, reducedEBRecHitCollection, reducedEERecHitCollection);
255  EcalClusterLazyTools lazyTools(iEvent,
256  ecalClusterToolsESGetTokens_.get(iSetup),
257  reducedBarrelRecHitCollectionToken_,
258  reducedEndcapRecHitCollectionToken_);
259 
260  // for conversion veto selection
261  edm::Handle<reco::ConversionCollection> hConversions;
262  iEvent.getByToken(hConversionsToken_, hConversions);
263 
264  // Get the ESHandle for the transient track builder, if needed for
265  // high level selection embedding
266  edm::ESHandle<TransientTrackBuilder> trackBuilder;
267 
268  if (isolator_.enabled())
269  isolator_.beginEvent(iEvent, iSetup);
270 
271  if (efficiencyLoader_.enabled())
272  efficiencyLoader_.newEvent(iEvent);
273  if (resolutionLoader_.enabled())
274  resolutionLoader_.newEvent(iEvent, iSetup);
275 
276  IsoDepositMaps deposits(isoDepositTokens_.size());
277  for (size_t j = 0, nd = isoDepositTokens_.size(); j < nd; ++j) {
278  iEvent.getByToken(isoDepositTokens_[j], deposits[j]);
279  }
280 
281  IsolationValueMaps isolationValues(isolationValueTokens_.size());
282  for (size_t j = 0; j < isolationValueTokens_.size(); ++j) {
283  iEvent.getByToken(isolationValueTokens_[j], isolationValues[j]);
284  }
285 
286  IsolationValueMaps isolationValuesNoPFId(isolationValueNoPFIdTokens_.size());
287  for (size_t j = 0; j < isolationValueNoPFIdTokens_.size(); ++j) {
288  iEvent.getByToken(isolationValueNoPFIdTokens_[j], isolationValuesNoPFId[j]);
289  }
290 
291  // prepare the MC matching
292  GenAssociations genMatches(genMatchTokens_.size());
293  if (addGenMatch_) {
294  for (size_t j = 0, nd = genMatchTokens_.size(); j < nd; ++j) {
295  iEvent.getByToken(genMatchTokens_[j], genMatches[j]);
296  }
297  }
298 
299  // prepare ID extraction
300  std::vector<edm::Handle<edm::ValueMap<float>>> idhandles;
301  std::vector<pat::Electron::IdPair> ids;
302  if (addElecID_) {
303  idhandles.resize(elecIDSrcs_.size());
304  ids.resize(elecIDSrcs_.size());
305  for (size_t i = 0; i < elecIDSrcs_.size(); ++i) {
306  iEvent.getByToken(elecIDTokens_[i], idhandles[i]);
307  ids[i].first = elecIDSrcs_[i].first;
308  }
309  }
310 
311  // prepare the high level selection:
312  // needs beamline
313  reco::TrackBase::Point beamPoint(0, 0, 0);
314  reco::Vertex primaryVertex;
315  reco::BeamSpot beamSpot;
316  bool beamSpotIsValid = false;
317  bool primaryVertexIsValid = false;
318 
319  // Get the beamspot
320  edm::Handle<reco::BeamSpot> beamSpotHandle;
321  iEvent.getByToken(beamLineToken_, beamSpotHandle);
322 
323  if (embedHighLevelSelection_) {
324  // Get the primary vertex
325  edm::Handle<std::vector<reco::Vertex>> pvHandle;
326  iEvent.getByToken(pvToken_, pvHandle);
327 
328  // This is needed by the IPTools methods from the tracking group
329  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", trackBuilder);
330 
331  if (pvHandle.isValid() && !pvHandle->empty()) {
332  primaryVertex = pvHandle->at(0);
333  primaryVertexIsValid = true;
334  } else {
335  edm::LogError("DataNotAvailable")
336  << "No primary vertex available from EventSetup, not adding high level selection \n";
337  }
338  }
339  //value maps for puppi isolation
340  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_charged_hadrons;
341  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_neutral_hadrons;
342  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_photons;
343  //value maps for puppiNoLeptons isolation
344  edm::Handle<edm::ValueMap<float>> PUPPINoLeptonsIsolation_charged_hadrons;
345  edm::Handle<edm::ValueMap<float>> PUPPINoLeptonsIsolation_neutral_hadrons;
346  edm::Handle<edm::ValueMap<float>> PUPPINoLeptonsIsolation_photons;
347  if (addPuppiIsolation_) {
348  //puppi
349  iEvent.getByToken(PUPPIIsolation_charged_hadrons_, PUPPIIsolation_charged_hadrons);
350  iEvent.getByToken(PUPPIIsolation_neutral_hadrons_, PUPPIIsolation_neutral_hadrons);
351  iEvent.getByToken(PUPPIIsolation_photons_, PUPPIIsolation_photons);
352  //puppiNoLeptons
353  iEvent.getByToken(PUPPINoLeptonsIsolation_charged_hadrons_, PUPPINoLeptonsIsolation_charged_hadrons);
354  iEvent.getByToken(PUPPINoLeptonsIsolation_neutral_hadrons_, PUPPINoLeptonsIsolation_neutral_hadrons);
355  iEvent.getByToken(PUPPINoLeptonsIsolation_photons_, PUPPINoLeptonsIsolation_photons);
356  }
357 
358  std::vector<Electron>* patElectrons = new std::vector<Electron>();
359 
360  if (useParticleFlow_) {
361  edm::Handle<reco::PFCandidateCollection> pfElectrons;
362  iEvent.getByToken(pfElecToken_, pfElectrons);
363  unsigned index = 0;
364 
365  for (reco::PFCandidateConstIterator i = pfElectrons->begin(); i != pfElectrons->end(); ++i, ++index) {
366  reco::PFCandidateRef pfRef(pfElectrons, index);
367  reco::PFCandidatePtr ptrToPFElectron(pfElectrons, index);
368  // reco::CandidateBaseRef pfBaseRef( pfRef );
369 
370  reco::GsfTrackRef PfTk = i->gsfTrackRef();
371 
372  bool Matched = false;
373  bool MatchedToAmbiguousGsfTrack = false;
374  for (edm::View<reco::GsfElectron>::const_iterator itElectron = electrons->begin(); itElectron != electrons->end();
375  ++itElectron) {
376  unsigned int idx = itElectron - electrons->begin();
377  auto elePtr = electrons->ptrAt(idx);
378  if (Matched || MatchedToAmbiguousGsfTrack)
379  continue;
380 
381  reco::GsfTrackRef EgTk = itElectron->gsfTrack();
382 
383  if (itElectron->gsfTrack() == i->gsfTrackRef()) {
384  Matched = true;
385  } else {
386  for (reco::GsfTrackRefVector::const_iterator it = itElectron->ambiguousGsfTracksBegin();
387  it != itElectron->ambiguousGsfTracksEnd();
388  it++) {
389  MatchedToAmbiguousGsfTrack |= (bool)(i->gsfTrackRef() == (*it));
390  }
391  }
392 
393  if (Matched || MatchedToAmbiguousGsfTrack) {
394  // ptr needed for finding the matched gen particle
395  reco::CandidatePtr ptrToGsfElectron(electrons, idx);
396 
397  // ref to base needed for the construction of the pat object
398  const edm::RefToBase<reco::GsfElectron>& elecsRef = electrons->refAt(idx);
399  Electron anElectron(elecsRef);
400  anElectron.setPFCandidateRef(pfRef);
401  if (addPuppiIsolation_) {
402  anElectron.setIsolationPUPPI((*PUPPIIsolation_charged_hadrons)[elePtr],
403  (*PUPPIIsolation_neutral_hadrons)[elePtr],
404  (*PUPPIIsolation_photons)[elePtr]);
405  anElectron.setIsolationPUPPINoLeptons((*PUPPINoLeptonsIsolation_charged_hadrons)[elePtr],
406  (*PUPPINoLeptonsIsolation_neutral_hadrons)[elePtr],
407  (*PUPPINoLeptonsIsolation_photons)[elePtr]);
408  } else {
409  anElectron.setIsolationPUPPI(-999., -999., -999.);
410  anElectron.setIsolationPUPPINoLeptons(-999., -999., -999.);
411  }
412 
413  //it should be always true when particleFlow electrons are used.
414  anElectron.setIsPF(true);
415 
416  if (embedPFCandidate_)
417  anElectron.embedPFCandidate();
418 
419  if (useUserData_) {
420  userDataHelper_.add(anElectron, iEvent, iSetup);
421  }
422 
423  double ip3d = -999; // for mva variable
424 
425  // embed high level selection
426  if (embedHighLevelSelection_) {
427  // get the global track
428  const reco::GsfTrackRef& track = PfTk;
429 
430  // Make sure the collection it points to is there
431  if (track.isNonnull() && track.isAvailable()) {
432  reco::TransientTrack tt = trackBuilder->build(track);
433  embedHighLevel(anElectron, track, tt, primaryVertex, primaryVertexIsValid, beamSpot, beamSpotIsValid);
434 
435  std::pair<bool, Measurement1D> ip3dpv = IPTools::absoluteImpactParameter3D(tt, primaryVertex);
436  ip3d = ip3dpv.second.value(); // for mva variable
437  }
438  }
439 
440  //Electron Id
441 
442  if (addElecID_) {
443  //STANDARD EL ID
444  for (size_t i = 0; i < elecIDSrcs_.size(); ++i) {
445  ids[i].second = (*idhandles[i])[elecsRef];
446  }
447  //SPECIFIC PF ID
448  ids.push_back(std::make_pair("pf_evspi", pfRef->mva_e_pi()));
449  ids.push_back(std::make_pair("pf_evsmu", pfRef->mva_e_mu()));
450  anElectron.setElectronIDs(ids);
451  }
452 
453  if (addMVAVariables_) {
454  // add missing mva variables
455  std::vector<float> vCov = lazyTools.localCovariances(*(itElectron->superCluster()->seed()));
456  anElectron.setMvaVariables(vCov[1], ip3d);
457  }
458  // PFClusterIso
459  if (addPFClusterIso_) {
460  // Get PFCluster Isolation
461  edm::Handle<edm::ValueMap<float>> ecalPFClusterIsoMapH;
462  iEvent.getByToken(ecalPFClusterIsoT_, ecalPFClusterIsoMapH);
463  edm::Handle<edm::ValueMap<float>> hcalPFClusterIsoMapH;
464  iEvent.getByToken(hcalPFClusterIsoT_, hcalPFClusterIsoMapH);
465  reco::GsfElectron::PflowIsolationVariables newPFIsol = anElectron.pfIsolationVariables();
466  newPFIsol.sumEcalClusterEt = (*ecalPFClusterIsoMapH)[elecsRef];
467  newPFIsol.sumHcalClusterEt = (*hcalPFClusterIsoMapH)[elecsRef];
468  anElectron.setPfIsolationVariables(newPFIsol);
469  }
470 
471  std::vector<DetId> selectedCells;
472  bool barrel = itElectron->isEB();
473  //loop over sub clusters
474  if (embedBasicClusters_) {
475  for (reco::CaloCluster_iterator clusIt = itElectron->superCluster()->clustersBegin();
476  clusIt != itElectron->superCluster()->clustersEnd();
477  ++clusIt) {
478  //get seed (max energy xtal)
479  DetId seed = lazyTools.getMaximum(**clusIt).first;
480  //get all xtals in 5x5 window around the seed
481  std::vector<DetId> dets5x5 =
482  (barrel) ? ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalBarrel)->getWindow(seed, 5, 5)
483  : ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalEndcap)->getWindow(seed, 5, 5);
484  selectedCells.insert(selectedCells.end(), dets5x5.begin(), dets5x5.end());
485 
486  //get all xtals belonging to cluster
487  for (const std::pair<DetId, float>& hit : (*clusIt)->hitsAndFractions()) {
488  selectedCells.push_back(hit.first);
489  }
490  }
491  }
492 
493  if (embedPflowBasicClusters_ && itElectron->parentSuperCluster().isNonnull()) {
494  for (reco::CaloCluster_iterator clusIt = itElectron->parentSuperCluster()->clustersBegin();
495  clusIt != itElectron->parentSuperCluster()->clustersEnd();
496  ++clusIt) {
497  //get seed (max energy xtal)
498  DetId seed = lazyTools.getMaximum(**clusIt).first;
499  //get all xtals in 5x5 window around the seed
500  std::vector<DetId> dets5x5 =
501  (barrel) ? ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalBarrel)->getWindow(seed, 5, 5)
502  : ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalEndcap)->getWindow(seed, 5, 5);
503  selectedCells.insert(selectedCells.end(), dets5x5.begin(), dets5x5.end());
504 
505  //get all xtals belonging to cluster
506  for (const std::pair<DetId, float>& hit : (*clusIt)->hitsAndFractions()) {
507  selectedCells.push_back(hit.first);
508  }
509  }
510  }
511 
512  //remove duplicates
513  std::sort(selectedCells.begin(), selectedCells.end());
514  std::unique(selectedCells.begin(), selectedCells.end());
515 
516  // Retrieve the corresponding RecHits
517 
518  edm::Handle<EcalRecHitCollection> rechitsH;
519  if (barrel)
520  iEvent.getByToken(reducedBarrelRecHitCollectionToken_, rechitsH);
521  else
522  iEvent.getByToken(reducedEndcapRecHitCollectionToken_, rechitsH);
523 
524  EcalRecHitCollection selectedRecHits;
525  const EcalRecHitCollection* recHits = rechitsH.product();
526 
527  unsigned nSelectedCells = selectedCells.size();
528  for (unsigned icell = 0; icell < nSelectedCells; ++icell) {
529  EcalRecHitCollection::const_iterator it = recHits->find(selectedCells[icell]);
530  if (it != recHits->end()) {
531  selectedRecHits.push_back(*it);
532  }
533  }
534  selectedRecHits.sort();
535  if (embedRecHits_)
536  anElectron.embedRecHits(&selectedRecHits);
537 
538  // set conversion veto selection
539  bool passconversionveto = false;
540  if (hConversions.isValid()) {
541  // this is recommended method
542  passconversionveto =
543  !ConversionTools::hasMatchedConversion(*itElectron, *hConversions, beamSpotHandle->position());
544  } else {
545  // use missing hits without vertex fit method
546  passconversionveto =
547  itElectron->gsfTrack()->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS) < 1;
548  }
549 
550  anElectron.setPassConversionVeto(passconversionveto);
551 
552  // fillElectron(anElectron,elecsRef,pfBaseRef,
553  // genMatches, deposits, isolationValues);
554 
555  //COLIN small warning !
556  // we are currently choosing to take the 4-momentum of the PFCandidate;
557  // the momentum of the GsfElectron is saved though
558  // we must therefore match the GsfElectron.
559  // because of this, we should not change the source of the electron matcher
560  // to the collection of PFElectrons in the python configuration
561  // I don't know what to do with the efficiencyLoader, since I don't know
562  // what this class is for.
563  fillElectron2(
564  anElectron, ptrToPFElectron, ptrToGsfElectron, ptrToGsfElectron, genMatches, deposits, isolationValues);
565 
566  //COLIN need to use fillElectron2 in the non-pflow case as well, and to test it.
567 
568  if (computeMiniIso_)
569  setElectronMiniIso(anElectron, pc.product());
570 
571  patElectrons->push_back(anElectron);
572  }
573  }
574  //if( !Matched && !MatchedToAmbiguousGsfTrack) std::cout << "!!!!A pf electron could not be matched to a gsf!!!!" << std::endl;
575  }
576  }
577 
578  else {
579  edm::Handle<reco::PFCandidateCollection> pfElectrons;
580  edm::Handle<edm::ValueMap<reco::PFCandidatePtr>> ValMapH;
581  edm::Handle<edm::ValueMap<std::vector<reco::PFCandidateRef>>> ValMultiMapH;
582  bool pfCandsPresent = false, valMapPresent = false;
583  if (usePfCandidateMultiMap_) {
584  iEvent.getByToken(pfCandidateMultiMapToken_, ValMultiMapH);
585  } else {
586  pfCandsPresent = iEvent.getByToken(pfElecToken_, pfElectrons);
587  valMapPresent = iEvent.getByToken(pfCandidateMapToken_, ValMapH);
588  }
589 
590  for (edm::View<reco::GsfElectron>::const_iterator itElectron = electrons->begin(); itElectron != electrons->end();
591  ++itElectron) {
592  // construct the Electron from the ref -> save ref to original object
593  //FIXME: looks like a lot of instances could be turned into const refs
594  unsigned int idx = itElectron - electrons->begin();
595  edm::RefToBase<reco::GsfElectron> elecsRef = electrons->refAt(idx);
596  reco::CandidateBaseRef elecBaseRef(elecsRef);
597  Electron anElectron(elecsRef);
598  auto elePtr = electrons->ptrAt(idx);
599 
600  // Is this GsfElectron also identified as an e- in the particle flow?
601  bool pfId = false;
602 
603  if (usePfCandidateMultiMap_) {
604  for (const reco::PFCandidateRef& pf : (*ValMultiMapH)[elePtr]) {
605  if (pf->particleId() == reco::PFCandidate::e) {
606  pfId = true;
607  anElectron.setPFCandidateRef(pf);
608  break;
609  }
610  }
611  } else if (pfCandsPresent) {
612  // PF electron collection not available.
613  const reco::GsfTrackRef& trkRef = itElectron->gsfTrack();
614  int index = 0;
615  for (reco::PFCandidateConstIterator ie = pfElectrons->begin(); ie != pfElectrons->end(); ++ie, ++index) {
616  if (ie->particleId() != reco::PFCandidate::e)
617  continue;
618  const reco::GsfTrackRef& pfTrkRef = ie->gsfTrackRef();
619  if (trkRef == pfTrkRef) {
620  pfId = true;
621  reco::PFCandidateRef pfRef(pfElectrons, index);
622  anElectron.setPFCandidateRef(pfRef);
623  break;
624  }
625  }
626  } else if (valMapPresent) {
627  // use value map if PF collection not available
628  const edm::ValueMap<reco::PFCandidatePtr>& myValMap(*ValMapH);
629  // Get the PFCandidate
630  const reco::PFCandidatePtr& pfElePtr(myValMap[elecsRef]);
631  pfId = pfElePtr.isNonnull();
632  }
633  // set PFId function
634  anElectron.setIsPF(pfId);
635 
636  // add resolution info
637 
638  // Isolation
639  if (isolator_.enabled()) {
640  isolator_.fill(*electrons, idx, isolatorTmpStorage_);
641  typedef pat::helper::MultiIsolator::IsolationValuePairs IsolationValuePairs;
642  // better to loop backwards, so the vector is resized less times
643  for (IsolationValuePairs::const_reverse_iterator it = isolatorTmpStorage_.rbegin(),
644  ed = isolatorTmpStorage_.rend();
645  it != ed;
646  ++it) {
647  anElectron.setIsolation(it->first, it->second);
648  }
649  }
650 
651  for (size_t j = 0, nd = deposits.size(); j < nd; ++j) {
652  anElectron.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[elecsRef]);
653  }
654 
655  // add electron ID info
656  if (addElecID_) {
657  for (size_t i = 0; i < elecIDSrcs_.size(); ++i) {
658  ids[i].second = (*idhandles[i])[elecsRef];
659  }
660  anElectron.setElectronIDs(ids);
661  }
662 
663  if (useUserData_) {
664  userDataHelper_.add(anElectron, iEvent, iSetup);
665  }
666 
667  double ip3d = -999; //for mva variable
668 
669  // embed high level selection
670  if (embedHighLevelSelection_) {
671  // get the global track
672  reco::GsfTrackRef track = itElectron->gsfTrack();
673 
674  // Make sure the collection it points to is there
675  if (track.isNonnull() && track.isAvailable()) {
676  reco::TransientTrack tt = trackBuilder->build(track);
677  embedHighLevel(anElectron, track, tt, primaryVertex, primaryVertexIsValid, beamSpot, beamSpotIsValid);
678 
679  std::pair<bool, Measurement1D> ip3dpv = IPTools::absoluteImpactParameter3D(tt, primaryVertex);
680  ip3d = ip3dpv.second.value(); // for mva variable
681  }
682  }
683 
684  if (addMVAVariables_) {
685  // add mva variables
686  std::vector<float> vCov = lazyTools.localCovariances(*(itElectron->superCluster()->seed()));
687  anElectron.setMvaVariables(vCov[1], ip3d);
688  }
689 
690  // PFCluster Isolation
691  if (addPFClusterIso_) {
692  // Get PFCluster Isolation
693  edm::Handle<edm::ValueMap<float>> ecalPFClusterIsoMapH;
694  iEvent.getByToken(ecalPFClusterIsoT_, ecalPFClusterIsoMapH);
695  edm::Handle<edm::ValueMap<float>> hcalPFClusterIsoMapH;
696  iEvent.getByToken(hcalPFClusterIsoT_, hcalPFClusterIsoMapH);
697  reco::GsfElectron::PflowIsolationVariables newPFIsol = anElectron.pfIsolationVariables();
698  newPFIsol.sumEcalClusterEt = (*ecalPFClusterIsoMapH)[elecsRef];
699  newPFIsol.sumHcalClusterEt = (*hcalPFClusterIsoMapH)[elecsRef];
700  anElectron.setPfIsolationVariables(newPFIsol);
701  }
702 
703  if (addPuppiIsolation_) {
704  anElectron.setIsolationPUPPI((*PUPPIIsolation_charged_hadrons)[elePtr],
705  (*PUPPIIsolation_neutral_hadrons)[elePtr],
706  (*PUPPIIsolation_photons)[elePtr]);
707  anElectron.setIsolationPUPPINoLeptons((*PUPPINoLeptonsIsolation_charged_hadrons)[elePtr],
708  (*PUPPINoLeptonsIsolation_neutral_hadrons)[elePtr],
709  (*PUPPINoLeptonsIsolation_photons)[elePtr]);
710  } else {
711  anElectron.setIsolationPUPPI(-999., -999., -999.);
712  anElectron.setIsolationPUPPINoLeptons(-999., -999., -999.);
713  }
714 
715  std::vector<DetId> selectedCells;
716  bool barrel = itElectron->isEB();
717  //loop over sub clusters
718  if (embedBasicClusters_) {
719  for (reco::CaloCluster_iterator clusIt = itElectron->superCluster()->clustersBegin();
720  clusIt != itElectron->superCluster()->clustersEnd();
721  ++clusIt) {
722  //get seed (max energy xtal)
723  DetId seed = lazyTools.getMaximum(**clusIt).first;
724  //get all xtals in 5x5 window around the seed
725  std::vector<DetId> dets5x5 =
726  (barrel) ? ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalBarrel)->getWindow(seed, 5, 5)
727  : ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalEndcap)->getWindow(seed, 5, 5);
728  selectedCells.insert(selectedCells.end(), dets5x5.begin(), dets5x5.end());
729 
730  //get all xtals belonging to cluster
731  for (const std::pair<DetId, float>& hit : (*clusIt)->hitsAndFractions()) {
732  selectedCells.push_back(hit.first);
733  }
734  }
735  }
736 
737  if (embedPflowBasicClusters_ && itElectron->parentSuperCluster().isNonnull()) {
738  for (reco::CaloCluster_iterator clusIt = itElectron->parentSuperCluster()->clustersBegin();
739  clusIt != itElectron->parentSuperCluster()->clustersEnd();
740  ++clusIt) {
741  //get seed (max energy xtal)
742  DetId seed = lazyTools.getMaximum(**clusIt).first;
743  //get all xtals in 5x5 window around the seed
744  std::vector<DetId> dets5x5 =
745  (barrel) ? ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalBarrel)->getWindow(seed, 5, 5)
746  : ecalTopology_->getSubdetectorTopology(DetId::Ecal, EcalEndcap)->getWindow(seed, 5, 5);
747  selectedCells.insert(selectedCells.end(), dets5x5.begin(), dets5x5.end());
748 
749  //get all xtals belonging to cluster
750  for (const std::pair<DetId, float>& hit : (*clusIt)->hitsAndFractions()) {
751  selectedCells.push_back(hit.first);
752  }
753  }
754  }
755 
756  //remove duplicates
757  std::sort(selectedCells.begin(), selectedCells.end());
758  std::unique(selectedCells.begin(), selectedCells.end());
759 
760  // Retrieve the corresponding RecHits
761 
762  edm::Handle<EcalRecHitCollection> rechitsH;
763  if (barrel)
764  iEvent.getByToken(reducedBarrelRecHitCollectionToken_, rechitsH);
765  else
766  iEvent.getByToken(reducedEndcapRecHitCollectionToken_, rechitsH);
767 
768  EcalRecHitCollection selectedRecHits;
769  const EcalRecHitCollection* recHits = rechitsH.product();
770 
771  unsigned nSelectedCells = selectedCells.size();
772  for (unsigned icell = 0; icell < nSelectedCells; ++icell) {
773  EcalRecHitCollection::const_iterator it = recHits->find(selectedCells[icell]);
774  if (it != recHits->end()) {
775  selectedRecHits.push_back(*it);
776  }
777  }
778  selectedRecHits.sort();
779  if (embedRecHits_)
780  anElectron.embedRecHits(&selectedRecHits);
781 
782  // set conversion veto selection
783  bool passconversionveto = false;
784  if (hConversions.isValid()) {
785  // this is recommended method
786  passconversionveto =
787  !ConversionTools::hasMatchedConversion(*itElectron, *hConversions, beamSpotHandle->position());
788  } else {
789  // use missing hits without vertex fit method
790  passconversionveto =
791  itElectron->gsfTrack()->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS) < 1;
792  }
793  anElectron.setPassConversionVeto(passconversionveto);
794 
795  // add sel to selected
796  fillElectron(
797  anElectron, elecsRef, elecBaseRef, genMatches, deposits, pfId, isolationValues, isolationValuesNoPFId);
798 
799  if (computeMiniIso_)
800  setElectronMiniIso(anElectron, pc.product());
801 
802  patElectrons->push_back(anElectron);
803  }
804  }
805 
806  // sort electrons in pt
807  std::sort(patElectrons->begin(), patElectrons->end(), pTComparator_);
808 
809  // add the electrons to the event output
810  std::unique_ptr<std::vector<Electron>> ptr(patElectrons);
811  iEvent.put(std::move(ptr));
812 
813  // clean up
814  if (isolator_.enabled())
815  isolator_.endEvent();
816 }
817 
818 void PATElectronProducer::fillElectron(Electron& anElectron,
819  const edm::RefToBase<reco::GsfElectron>& elecRef,
820  const reco::CandidateBaseRef& baseRef,
821  const GenAssociations& genMatches,
822  const IsoDepositMaps& deposits,
823  const bool pfId,
824  const IsolationValueMaps& isolationValues,
825  const IsolationValueMaps& isolationValuesNoPFId) const {
826  //COLIN: might want to use the PFCandidate 4-mom. Which one is in use now?
827  // if (useParticleFlow_)
828  // aMuon.setP4( aMuon.pfCandidateRef()->p4() );
829 
830  //COLIN:
831  //In the embedding case, the reference cannot be used to look into a value map.
832  //therefore, one has to had the PFCandidateRef to this function, which becomes a bit
833  //too much specific.
834 
835  // in fact, this function needs a baseref or ptr for genmatch
836  // and a baseref or ptr for isodeposits and isolationvalues.
837  // baseref is not needed
838  // the ptrForIsolation and ptrForMatching should be defined upstream.
839 
840  // is the concrete elecRef needed for the efficiency loader? what is this loader?
841  // how can we make it compatible with the particle flow electrons?
842 
843  if (embedGsfElectronCore_)
844  anElectron.embedGsfElectronCore();
845  if (embedGsfTrack_)
846  anElectron.embedGsfTrack();
847  if (embedSuperCluster_)
848  anElectron.embedSuperCluster();
849  if (embedPflowSuperCluster_)
850  anElectron.embedPflowSuperCluster();
851  if (embedSeedCluster_)
852  anElectron.embedSeedCluster();
853  if (embedBasicClusters_)
854  anElectron.embedBasicClusters();
855  if (embedPreshowerClusters_)
856  anElectron.embedPreshowerClusters();
857  if (embedPflowBasicClusters_)
858  anElectron.embedPflowBasicClusters();
859  if (embedPflowPreshowerClusters_)
860  anElectron.embedPflowPreshowerClusters();
861  if (embedTrack_)
862  anElectron.embedTrack();
863 
864  // store the match to the generated final state muons
865  if (addGenMatch_) {
866  for (size_t i = 0, n = genMatches.size(); i < n; ++i) {
867  if (useParticleFlow_) {
868  reco::GenParticleRef genElectron = (*genMatches[i])[anElectron.pfCandidateRef()];
869  anElectron.addGenParticleRef(genElectron);
870  } else {
871  reco::GenParticleRef genElectron = (*genMatches[i])[elecRef];
872  anElectron.addGenParticleRef(genElectron);
873  }
874  }
875  if (embedGenMatch_)
876  anElectron.embedGenParticle();
877  }
878 
879  if (efficiencyLoader_.enabled()) {
880  efficiencyLoader_.setEfficiencies(anElectron, elecRef);
881  }
882 
883  if (resolutionLoader_.enabled()) {
884  resolutionLoader_.setResolutions(anElectron);
885  }
886 
887  for (size_t j = 0, nd = deposits.size(); j < nd; ++j) {
888  if (useParticleFlow_) {
889  reco::PFCandidateRef pfcandref = anElectron.pfCandidateRef();
890  assert(!pfcandref.isNull());
891  reco::CandidatePtr source = pfcandref->sourceCandidatePtr(0);
892  anElectron.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[source]);
893  } else
894  anElectron.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[elecRef]);
895  }
896 
897  for (size_t j = 0; j < isolationValues.size(); ++j) {
898  if (useParticleFlow_) {
899  reco::CandidatePtr source = anElectron.pfCandidateRef()->sourceCandidatePtr(0);
900  anElectron.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[source]);
901  } else if (pfId) {
902  anElectron.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[elecRef]);
903  }
904  }
905 
906  //for electrons not identified as PF electrons
907  for (size_t j = 0; j < isolationValuesNoPFId.size(); ++j) {
908  if (!pfId) {
909  anElectron.setIsolation(isolationValueLabelsNoPFId_[j].first, (*isolationValuesNoPFId[j])[elecRef]);
910  }
911  }
912 }
913 
914 void PATElectronProducer::fillElectron2(Electron& anElectron,
915  const reco::CandidatePtr& candPtrForIsolation,
916  const reco::CandidatePtr& candPtrForGenMatch,
917  const reco::CandidatePtr& candPtrForLoader,
918  const GenAssociations& genMatches,
919  const IsoDepositMaps& deposits,
920  const IsolationValueMaps& isolationValues) const {
921  //COLIN/Florian: use the PFCandidate 4-mom.
922  anElectron.setEcalDrivenMomentum(anElectron.p4());
923  anElectron.setP4(anElectron.pfCandidateRef()->p4());
924 
925  // is the concrete elecRef needed for the efficiency loader? what is this loader?
926  // how can we make it compatible with the particle flow electrons?
927 
928  if (embedGsfElectronCore_)
929  anElectron.embedGsfElectronCore();
930  if (embedGsfTrack_)
931  anElectron.embedGsfTrack();
932  if (embedSuperCluster_)
933  anElectron.embedSuperCluster();
934  if (embedPflowSuperCluster_)
935  anElectron.embedPflowSuperCluster();
936  if (embedSeedCluster_)
937  anElectron.embedSeedCluster();
938  if (embedBasicClusters_)
939  anElectron.embedBasicClusters();
940  if (embedPreshowerClusters_)
941  anElectron.embedPreshowerClusters();
942  if (embedPflowBasicClusters_)
943  anElectron.embedPflowBasicClusters();
944  if (embedPflowPreshowerClusters_)
945  anElectron.embedPflowPreshowerClusters();
946  if (embedTrack_)
947  anElectron.embedTrack();
948 
949  // store the match to the generated final state muons
950 
951  if (addGenMatch_) {
952  for (size_t i = 0, n = genMatches.size(); i < n; ++i) {
953  reco::GenParticleRef genElectron = (*genMatches[i])[candPtrForGenMatch];
954  anElectron.addGenParticleRef(genElectron);
955  }
956  if (embedGenMatch_)
957  anElectron.embedGenParticle();
958  }
959 
960  //COLIN what's this? does it have to be GsfElectron specific?
961  if (efficiencyLoader_.enabled()) {
962  efficiencyLoader_.setEfficiencies(anElectron, candPtrForLoader);
963  }
964 
965  if (resolutionLoader_.enabled()) {
966  resolutionLoader_.setResolutions(anElectron);
967  }
968 
969  for (size_t j = 0, nd = deposits.size(); j < nd; ++j) {
970  if (isoDepositLabels_[j].first == pat::TrackIso || isoDepositLabels_[j].first == pat::EcalIso ||
971  isoDepositLabels_[j].first == pat::HcalIso || deposits[j]->contains(candPtrForGenMatch.id())) {
972  anElectron.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[candPtrForGenMatch]);
973  } else if (deposits[j]->contains(candPtrForIsolation.id())) {
974  anElectron.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[candPtrForIsolation]);
975  } else {
976  anElectron.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[candPtrForIsolation->sourceCandidatePtr(0)]);
977  }
978  }
979 
980  for (size_t j = 0; j < isolationValues.size(); ++j) {
981  if (isolationValueLabels_[j].first == pat::TrackIso || isolationValueLabels_[j].first == pat::EcalIso ||
982  isolationValueLabels_[j].first == pat::HcalIso || isolationValues[j]->contains(candPtrForGenMatch.id())) {
983  anElectron.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[candPtrForGenMatch]);
984  } else if (isolationValues[j]->contains(candPtrForIsolation.id())) {
985  anElectron.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[candPtrForIsolation]);
986  } else {
987  anElectron.setIsolation(isolationValueLabels_[j].first,
988  (*isolationValues[j])[candPtrForIsolation->sourceCandidatePtr(0)]);
989  }
990  }
991 }
992 
993 void PATElectronProducer::setElectronMiniIso(Electron& anElectron, const PackedCandidateCollection* pc) {
994  pat::PFIsolation miniiso;
995  if (anElectron.isEE())
996  miniiso = pat::getMiniPFIsolation(pc,
997  anElectron.polarP4(),
998  miniIsoParamsE_[0],
999  miniIsoParamsE_[1],
1000  miniIsoParamsE_[2],
1001  miniIsoParamsE_[3],
1002  miniIsoParamsE_[4],
1003  miniIsoParamsE_[5],
1004  miniIsoParamsE_[6],
1005  miniIsoParamsE_[7],
1006  miniIsoParamsE_[8]);
1007  else
1008  miniiso = pat::getMiniPFIsolation(pc,
1009  anElectron.polarP4(),
1010  miniIsoParamsB_[0],
1011  miniIsoParamsB_[1],
1012  miniIsoParamsB_[2],
1013  miniIsoParamsB_[3],
1014  miniIsoParamsB_[4],
1015  miniIsoParamsB_[5],
1016  miniIsoParamsB_[6],
1017  miniIsoParamsB_[7],
1018  miniIsoParamsB_[8]);
1019  anElectron.setMiniPFIsolation(miniiso);
1020 }
1021 
1022 // ParameterSet description for module
1023 void PATElectronProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
1024  edm::ParameterSetDescription iDesc;
1025  iDesc.setComment("PAT electron producer module");
1026 
1027  // input source
1028  iDesc.add<edm::InputTag>("pfCandidateMap", edm::InputTag("no default"))->setComment("input collection");
1029  iDesc.add<edm::InputTag>("electronSource", edm::InputTag("no default"))->setComment("input collection");
1030 
1031  iDesc.ifValue(
1032  edm::ParameterDescription<bool>("addPFClusterIso", false, true),
1033  true >> (edm::ParameterDescription<edm::InputTag>(
1034  "ecalPFClusterIsoMap", edm::InputTag("electronEcalPFClusterIsolationProducer"), true) and
1035  edm::ParameterDescription<edm::InputTag>(
1036  "hcalPFClusterIsoMap", edm::InputTag("electronHcalPFClusterIsolationProducer"), true)) or
1037  false >> (edm::ParameterDescription<edm::InputTag>("ecalPFClusterIsoMap", edm::InputTag(""), true) and
1038  edm::ParameterDescription<edm::InputTag>("hcalPFClusterIsoMap", edm::InputTag(""), true)));
1039 
1040  iDesc.ifValue(edm::ParameterDescription<bool>("addPuppiIsolation", false, true),
1041  true >> (edm::ParameterDescription<edm::InputTag>(
1042  "puppiIsolationChargedHadrons",
1043  edm::InputTag("egmElectronPUPPIIsolation", "h+-DR030-BarVeto000-EndVeto001"),
1044  true) and
1045  edm::ParameterDescription<edm::InputTag>(
1046  "puppiIsolationNeutralHadrons",
1047  edm::InputTag("egmElectronPUPPIIsolation", "h0-DR030-BarVeto000-EndVeto000"),
1048  true) and
1049  edm::ParameterDescription<edm::InputTag>(
1050  "puppiIsolationPhotons",
1051  edm::InputTag("egmElectronPUPPIIsolation", "gamma-DR030-BarVeto000-EndVeto008"),
1052  true) and
1053  edm::ParameterDescription<edm::InputTag>(
1054  "puppiNoLeptonsIsolationChargedHadrons",
1055  edm::InputTag("egmElectronPUPPINoLeptonsIsolation", "gamma-DR030-BarVeto000-EndVeto008"),
1056  true) and
1057  edm::ParameterDescription<edm::InputTag>(
1058  "puppiNoLeptonsIsolationNeutralHadrons",
1059  edm::InputTag("egmElectronPUPPINoLeptonsIsolation", "gamma-DR030-BarVeto000-EndVeto008"),
1060  true) and
1061  edm::ParameterDescription<edm::InputTag>(
1062  "puppiNoLeptonsIsolationPhotons",
1063  edm::InputTag("egmElectronPUPPINoLeptonsIsolation", "gamma-DR030-BarVeto000-EndVeto008"),
1064  true)) or
1065  false >> edm::EmptyGroupDescription());
1066 
1067  // embedding
1068  iDesc.add<bool>("embedGsfElectronCore", true)->setComment("embed external gsf electron core");
1069  iDesc.add<bool>("embedGsfTrack", true)->setComment("embed external gsf track");
1070  iDesc.add<bool>("embedSuperCluster", true)->setComment("embed external super cluster");
1071  iDesc.add<bool>("embedPflowSuperCluster", true)->setComment("embed external super cluster");
1072  iDesc.add<bool>("embedSeedCluster", true)->setComment("embed external seed cluster");
1073  iDesc.add<bool>("embedBasicClusters", true)->setComment("embed external basic clusters");
1074  iDesc.add<bool>("embedPreshowerClusters", true)->setComment("embed external preshower clusters");
1075  iDesc.add<bool>("embedPflowBasicClusters", true)->setComment("embed external pflow basic clusters");
1076  iDesc.add<bool>("embedPflowPreshowerClusters", true)->setComment("embed external pflow preshower clusters");
1077  iDesc.add<bool>("embedTrack", false)->setComment("embed external track");
1078  iDesc.add<bool>("embedRecHits", true)->setComment("embed external RecHits");
1079 
1080  // pf specific parameters
1081  iDesc.add<edm::InputTag>("pfElectronSource", edm::InputTag("pfElectrons"))
1082  ->setComment("particle flow input collection");
1083  auto&& usePfCandidateMultiMap = edm::ParameterDescription<bool>("usePfCandidateMultiMap", false, true);
1084  usePfCandidateMultiMap.setComment(
1085  "take ParticleFlow candidates from pfCandidateMultiMap instead of matching to pfElectrons by Gsf track "
1086  "reference");
1087  iDesc.ifValue(usePfCandidateMultiMap,
1088  true >> edm::ParameterDescription<edm::InputTag>("pfCandidateMultiMap", true) or
1089  false >> edm::EmptyGroupDescription());
1090  iDesc.add<bool>("useParticleFlow", false)->setComment("whether to use particle flow or not");
1091  iDesc.add<bool>("embedPFCandidate", false)->setComment("embed external particle flow object");
1092 
1093  // MC matching configurables
1094  iDesc.add<bool>("addGenMatch", true)->setComment("add MC matching");
1095  iDesc.add<bool>("embedGenMatch", false)->setComment("embed MC matched MC information");
1096  std::vector<edm::InputTag> emptySourceVector;
1097  iDesc
1098  .addNode(edm::ParameterDescription<edm::InputTag>("genParticleMatch", edm::InputTag(), true) xor
1099  edm::ParameterDescription<std::vector<edm::InputTag>>("genParticleMatch", emptySourceVector, true))
1100  ->setComment("input with MC match information");
1101 
1102  // electron ID configurables
1103  iDesc.add<bool>("addElectronID", true)->setComment("add electron ID variables");
1104  edm::ParameterSetDescription electronIDSourcesPSet;
1105  electronIDSourcesPSet.setAllowAnything();
1106  iDesc
1107  .addNode(
1108  edm::ParameterDescription<edm::InputTag>("electronIDSource", edm::InputTag(), true) xor
1109  edm::ParameterDescription<edm::ParameterSetDescription>("electronIDSources", electronIDSourcesPSet, true))
1110  ->setComment("input with electron ID variables");
1111 
1112  // mini-iso
1113  iDesc.add<bool>("computeMiniIso", false)->setComment("whether or not to compute and store electron mini-isolation");
1114  iDesc.add<edm::InputTag>("pfCandsForMiniIso", edm::InputTag("packedPFCandidates"))
1115  ->setComment("collection to use to compute mini-iso");
1116  iDesc.add<std::vector<double>>("miniIsoParamsE", std::vector<double>())
1117  ->setComment("mini-iso parameters to use for endcap electrons");
1118  iDesc.add<std::vector<double>>("miniIsoParamsB", std::vector<double>())
1119  ->setComment("mini-iso parameters to use for barrel electrons");
1120 
1121  // IsoDeposit configurables
1122  edm::ParameterSetDescription isoDepositsPSet;
1123  isoDepositsPSet.addOptional<edm::InputTag>("tracker");
1124  isoDepositsPSet.addOptional<edm::InputTag>("ecal");
1125  isoDepositsPSet.addOptional<edm::InputTag>("hcal");
1126  isoDepositsPSet.addOptional<edm::InputTag>("pfAllParticles");
1127  isoDepositsPSet.addOptional<edm::InputTag>("pfChargedHadrons");
1128  isoDepositsPSet.addOptional<edm::InputTag>("pfChargedAll");
1129  isoDepositsPSet.addOptional<edm::InputTag>("pfPUChargedHadrons");
1130  isoDepositsPSet.addOptional<edm::InputTag>("pfNeutralHadrons");
1131  isoDepositsPSet.addOptional<edm::InputTag>("pfPhotons");
1132  isoDepositsPSet.addOptional<std::vector<edm::InputTag>>("user");
1133  iDesc.addOptional("isoDeposits", isoDepositsPSet);
1134 
1135  // isolation values configurables
1136  edm::ParameterSetDescription isolationValuesPSet;
1137  isolationValuesPSet.addOptional<edm::InputTag>("tracker");
1138  isolationValuesPSet.addOptional<edm::InputTag>("ecal");
1139  isolationValuesPSet.addOptional<edm::InputTag>("hcal");
1140  isolationValuesPSet.addOptional<edm::InputTag>("pfAllParticles");
1141  isolationValuesPSet.addOptional<edm::InputTag>("pfChargedHadrons");
1142  isolationValuesPSet.addOptional<edm::InputTag>("pfChargedAll");
1143  isolationValuesPSet.addOptional<edm::InputTag>("pfPUChargedHadrons");
1144  isolationValuesPSet.addOptional<edm::InputTag>("pfNeutralHadrons");
1145  isolationValuesPSet.addOptional<edm::InputTag>("pfPhotons");
1146  isolationValuesPSet.addOptional<std::vector<edm::InputTag>>("user");
1147  iDesc.addOptional("isolationValues", isolationValuesPSet);
1148 
1149  // isolation values configurables
1150  edm::ParameterSetDescription isolationValuesNoPFIdPSet;
1151  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("tracker");
1152  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("ecal");
1153  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("hcal");
1154  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("pfAllParticles");
1155  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("pfChargedHadrons");
1156  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("pfChargedAll");
1157  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("pfPUChargedHadrons");
1158  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("pfNeutralHadrons");
1159  isolationValuesNoPFIdPSet.addOptional<edm::InputTag>("pfPhotons");
1160  isolationValuesNoPFIdPSet.addOptional<std::vector<edm::InputTag>>("user");
1161  iDesc.addOptional("isolationValuesNoPFId", isolationValuesNoPFIdPSet);
1162 
1163  // Efficiency configurables
1164  edm::ParameterSetDescription efficienciesPSet;
1165  efficienciesPSet.setAllowAnything(); // TODO: the pat helper needs to implement a description.
1166  iDesc.add("efficiencies", efficienciesPSet);
1167  iDesc.add<bool>("addEfficiencies", false);
1168 
1169  // Check to see if the user wants to add user data
1170  edm::ParameterSetDescription userDataPSet;
1171  PATUserDataHelper<Electron>::fillDescription(userDataPSet);
1172  iDesc.addOptional("userData", userDataPSet);
1173 
1174  // electron shapes
1175  iDesc.add<bool>("addMVAVariables", true)->setComment("embed extra variables in pat::Electron : sip3d, sigmaIEtaIPhi");
1176  iDesc.add<edm::InputTag>("reducedBarrelRecHitCollection", edm::InputTag("reducedEcalRecHitsEB"));
1177  iDesc.add<edm::InputTag>("reducedEndcapRecHitCollection", edm::InputTag("reducedEcalRecHitsEE"));
1178 
1179  edm::ParameterSetDescription isolationPSet;
1180  isolationPSet.setAllowAnything(); // TODO: the pat helper needs to implement a description.
1181  iDesc.add("userIsolation", isolationPSet);
1182 
1183  // Resolution configurables
1184  pat::helper::KinResolutionsLoader::fillDescription(iDesc);
1185 
1186  iDesc.add<bool>("embedHighLevelSelection", true)->setComment("embed high level selection");
1187  edm::ParameterSetDescription highLevelPSet;
1188  highLevelPSet.setAllowAnything();
1189  iDesc.addNode(edm::ParameterDescription<edm::InputTag>("beamLineSrc", edm::InputTag(), true))
1190  ->setComment("input with high level selection");
1191  iDesc.addNode(edm::ParameterDescription<edm::InputTag>("pvSrc", edm::InputTag(), true))
1192  ->setComment("input with high level selection");
1193 
1194  descriptions.add("PATElectronProducer", iDesc);
1195 }
1196 
1197 // embed various impact parameters with errors
1198 // embed high level selection
1199 void PATElectronProducer::embedHighLevel(pat::Electron& anElectron,
1200  reco::GsfTrackRef track,
1201  reco::TransientTrack& tt,
1202  reco::Vertex& primaryVertex,
1203  bool primaryVertexIsValid,
1204  reco::BeamSpot& beamspot,
1205  bool beamspotIsValid) {
1206  // Correct to PV
1207  // PV2D
1208  anElectron.setDB(track->dxy(primaryVertex.position()),
1209  track->dxyError(primaryVertex.position(), primaryVertex.covariance()),
1210  pat::Electron::PV2D);
1211 
1212  // PV3D
1213  std::pair<bool, Measurement1D> result =
1214  IPTools::signedImpactParameter3D(tt, GlobalVector(track->px(), track->py(), track->pz()), primaryVertex);
1215  double d0_corr = result.second.value();
1216  double d0_err = primaryVertexIsValid ? result.second.error() : -1.0;
1217  anElectron.setDB(d0_corr, d0_err, pat::Electron::PV3D);
1218 
1219  // Correct to beam spot
1220  // BS2D
1221  anElectron.setDB(track->dxy(beamspot), track->dxyError(beamspot), pat::Electron::BS2D);
1222 
1223  // make a fake vertex out of beam spot
1224  reco::Vertex vBeamspot(beamspot.position(), beamspot.covariance3D());
1225 
1226  // BS3D
1227  result = IPTools::signedImpactParameter3D(tt, GlobalVector(track->px(), track->py(), track->pz()), vBeamspot);
1228  d0_corr = result.second.value();
1229  d0_err = beamspotIsValid ? result.second.error() : -1.0;
1230  anElectron.setDB(d0_corr, d0_err, pat::Electron::BS3D);
1231 
1232  // PVDZ
1233  anElectron.setDB(
1234  track->dz(primaryVertex.position()), std::hypot(track->dzError(), primaryVertex.zError()), pat::Electron::PVDZ);
1235 }
1236 
1237 #include "FWCore/Framework/interface/MakerMacros.h"
1238 
1239 DEFINE_FWK_MODULE(PATElectronProducer);
ConfigurationDescriptions.h
reco::GsfElectron::isEE
bool isEE() const
Definition: GsfElectron.h:337
pat::PATElectronProducer::isolatorTmpStorage_
pat::helper::MultiIsolator::IsolationValuePairs isolatorTmpStorage_
Definition: PATElectronProducer.h:167
pat::Electron::embedTrack
void embedTrack()
method to store the electron's Track internally
pat::PATElectronProducer::isolationValueNoPFIdTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueNoPFIdTokens_
Definition: PATElectronProducer.h:173
electrons_cff.bool
bool
Definition: electrons_cff.py:393
pat::PATUserDataHelper
Assists in assimilating all pat::UserData into pat objects.
Definition: PATUserDataHelper.h:49
edm::SortedCollection::sort
void sort()
Definition: SortedCollection.h:302
mps_fire.i
i
Definition: mps_fire.py:428
pat::EcalIso
Definition: Isolation.h:11
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
pat::PATElectronProducer::embedGsfTrack_
const bool embedGsfTrack_
Definition: PATElectronProducer.h:64
edm::SortedCollection< EcalRecHit >::const_iterator
std::vector< EcalRecHit >::const_iterator const_iterator
Definition: SortedCollection.h:80
pwdgSkimBPark_cfi.beamSpot
beamSpot
Definition: pwdgSkimBPark_cfi.py:5
pat::PATElectronProducer::elecIDSrcs_
std::vector< NameTag > elecIDSrcs_
Definition: PATElectronProducer.h:160
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11713
pat::PATElectronProducer::genMatchTokens_
std::vector< edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > > > genMatchTokens_
Definition: PATElectronProducer.h:84
Reference_intrackfit_cff.barrel
list barrel
Definition: Reference_intrackfit_cff.py:37
PFIsolation.h
MessageLogger.h
pat::PATElectronProducer::addPFClusterIso_
const bool addPFClusterIso_
Definition: PATElectronProducer.h:102
funct::false
false
Definition: Factorize.h:29
pat::Electron::PV2D
Definition: Electron.h:190
dqmiodumpmetadata.n
n
Definition: dqmiodumpmetadata.py:28
reco::PFCandidatePtr
edm::Ptr< PFCandidate > PFCandidatePtr
persistent Ptr to a PFCandidate
Definition: PFCandidateFwd.h:27
edm::Handle::product
T const * product() const
Definition: Handle.h:70
electronProducer_cff.isolationValues
isolationValues
Definition: electronProducer_cff.py:26
pat::PATElectronProducer::embedHighLevel
void embedHighLevel(pat::Electron &anElectron, reco::GsfTrackRef track, reco::TransientTrack &tt, reco::Vertex &primaryVertex, bool primaryVertexIsValid, reco::BeamSpot &beamspot, bool beamspotIsValid)
Definition: PATElectronProducer.cc:1199
pat::PATElectronProducer::fillElectron2
void fillElectron2(Electron &anElectron, const reco::CandidatePtr &candPtrForIsolation, const reco::CandidatePtr &candPtrForGenMatch, const reco::CandidatePtr &candPtrForLoader, const GenAssociations &genMatches, const IsoDepositMaps &deposits, const IsolationValueMaps &isolationValues) const
Definition: PATElectronProducer.cc:914
pat::PATElectronProducer::PUPPINoLeptonsIsolation_photons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_photons_
Definition: PATElectronProducer.h:189
reco::PFCandidate::e
Definition: PFCandidate.h:47
groupFilesInBlocks.tt
int tt
Definition: groupFilesInBlocks.py:144
sistrip::View
View
Definition: ConstantsForView.h:26
pat::Electron::embedSeedCluster
void embedSeedCluster()
method to store the electron's seedcluster internally
PFCandidate.h
pat::Electron::embedPreshowerClusters
void embedPreshowerClusters()
method to store the electron's preshower clusters
pat::helper::EfficiencyLoader
Definition: EfficiencyLoader.h:16
pat::PATElectronProducer::setElectronMiniIso
void setElectronMiniIso(pat::Electron &anElectron, const pat::PackedCandidateCollection *pc)
Definition: PATElectronProducer.cc:993
pat::ConversionCollection
std::vector< Conversion > ConversionCollection
Definition: Conversion.h:13
pat::PATElectronProducer::resolutionLoader_
pat::helper::KinResolutionsLoader resolutionLoader_
Definition: PATElectronProducer.h:179
pat::PATElectronProducer::isolationValueLabels_
IsolationLabels isolationValueLabels_
Definition: PATElectronProducer.h:170
Electron
Definition: Electron.py:1
edm::Ref::isNull
bool isNull() const
Checks for null.
Definition: Ref.h:235
pat::helper::MultiIsolator::fill
void fill(const edm::View< T > &coll, int idx, IsolationValuePairs &isolations) const
Definition: MultiIsolator.h:84
edm
HLT enums.
Definition: AlignableModifier.h:19
pat::PATElectronProducer::embedRecHits_
const bool embedRecHits_
Definition: PATElectronProducer.h:75
IPTools::absoluteImpactParameter3D
std::pair< bool, Measurement1D > absoluteImpactParameter3D(const reco::TransientTrack &transientTrack, const reco::Vertex &vertex)
Definition: IPTools.cc:38
pat::PATUserDataHelper::fillDescription
static void fillDescription(edm::ParameterSetDescription &iDesc)
Definition: PATUserDataHelper.h:135
TrackerIsolationPt.h
pat::PATElectronProducer::pvToken_
const edm::EDGetTokenT< std::vector< reco::Vertex > > pvToken_
Definition: PATElectronProducer.h:110
reco::GsfElectron::PflowIsolationVariables
Definition: GsfElectron.h:606
pat::Electron::setEcalDrivenMomentum
void setEcalDrivenMomentum(const Candidate::LorentzVector &mom)
Definition: Electron.h:204
pat::Electron::embedGsfElectronCore
void embedGsfElectronCore()
method to store the electron's core internally
reco::PFCandidateConstIterator
PFCandidateCollection::const_iterator PFCandidateConstIterator
iterator
Definition: PFCandidateFwd.h:18
pat::Electron::setPassConversionVeto
void setPassConversionVeto(bool flag)
Definition: Electron.h:257
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89287
pat::PATElectronProducer::elecIDTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< float > > > elecIDTokens_
Definition: PATElectronProducer.h:161
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
pat::Electron::setIsolationPUPPINoLeptons
void setIsolationPUPPINoLeptons(float chargedhadrons_, float neutralhadrons_, float photons_)
sets PUPPINoLeptons isolations
Definition: Electron.h:167
pat::Lepton::setIsolation
void setIsolation(IsolationKeys key, float value)
Definition: Lepton.h:115
pat::helper::EfficiencyLoader::enabled
bool enabled() const
'true' if this there is at least one efficiency configured
Definition: EfficiencyLoader.h:25
pat::Electron::BS2D
Definition: Electron.h:190
cms::cuda::assert
assert(be >=bs)
pat::Electron::embedRecHits
void embedRecHits(const EcalRecHitCollection *rechits)
method to store the RecHits internally - can be called from the PATElectronProducer
edm::PtrVectorItr
Definition: PtrVector.h:51
pat::Electron::setIsolationPUPPI
void setIsolationPUPPI(float chargedhadrons_, float neutralhadrons_, float photons_)
sets PUPPI isolations
Definition: Electron.h:161
edm::SortedCollection< EcalRecHit >
pat::PATElectronProducer::embedPflowBasicClusters_
const bool embedPflowBasicClusters_
Definition: PATElectronProducer.h:70
electronProducer_cff.isolationValuesNoPFId
isolationValuesNoPFId
Definition: electronProducer_cff.py:34
pat::PATObject::addGenParticleRef
void addGenParticleRef(const reco::GenParticleRef &ref)
Definition: PATObject.h:753
pat::PATElectronProducer::embedHighLevelSelection_
const bool embedHighLevelSelection_
embed high level selection variables?
Definition: PATElectronProducer.h:108
pat::PATElectronProducer::IsoDepositMaps
std::vector< edm::Handle< edm::ValueMap< IsoDeposit > > > IsoDepositMaps
Definition: PATElectronProducer.h:113
pat::PATElectronProducer::produce
void produce(edm::Event &iEvent, const edm::EventSetup &iSetup) override
Definition: PATElectronProducer.cc:232
EmptyGroupDescription.h
EcalClusterLazyTools.h
CaloTopologyRecord
Definition: CaloTopologyRecord.h:10
pat::PATElectronProducer::userDataHelper_
pat::PATUserDataHelper< pat::Electron > userDataHelper_
Definition: PATElectronProducer.h:190
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
pat::PATElectronProducer::PUPPINoLeptonsIsolation_charged_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_charged_hadrons_
Definition: PATElectronProducer.h:187
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10
TransientTrack.h
pat::PATElectronProducer::isolationValueLabelsNoPFId_
IsolationLabels isolationValueLabelsNoPFId_
Definition: PATElectronProducer.h:172
Association.h
pat::Electron::embedPflowBasicClusters
void embedPflowBasicClusters()
method to store the electron's pflow basic clusters
pat::PATElectronProducer::embedGenMatch_
bool embedGenMatch_
Definition: PATElectronProducer.h:74
edm::Handle
Definition: AssociativeIterator.h:50
edm::ParameterSetDescription::addOptional
ParameterDescriptionBase * addOptional(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:105
pat::Electron::pfCandidateRef
reco::PFCandidateRef pfCandidateRef() const
reference to the source PFCandidates; null if this has been built from a standard electron
edm::SortedCollection::push_back
void push_back(T const &t)
Definition: SortedCollection.h:188
dqmdumpme.first
first
Definition: dqmdumpme.py:55
pat::Electron::embedPflowSuperCluster
void embedPflowSuperCluster()
method to store the electron's PflowSuperCluster internally
pat::helper::MultiIsolator::IsolationValuePairs
std::vector< std::pair< pat::IsolationKeys, float > > IsolationValuePairs
Definition: MultiIsolator.h:17
EcalBarrel
Definition: EcalSubdetector.h:10
pat::helper::KinResolutionsLoader::enabled
bool enabled() const
'true' if this there is at least one efficiency configured
Definition: KinResolutionsLoader.h:27
pat::Electron::embedGsfTrack
void embedGsfTrack()
method to store the electron's GsfTrack internally
pat::PATElectronProducer::useUserData_
const bool useUserData_
Definition: PATElectronProducer.h:181
pat::PATElectronProducer::embedPreshowerClusters_
const bool embedPreshowerClusters_
Definition: PATElectronProducer.h:69
pat::PATElectronProducer::pTComparator_
const GreaterByPt< Electron > pTComparator_
Definition: PATElectronProducer.h:164
pat::Electron::BS3D
Definition: Electron.h:190
pat::PATElectronProducer::pfCandidateMapToken_
const edm::EDGetTokenT< edm::ValueMap< reco::PFCandidatePtr > > pfCandidateMapToken_
Definition: PATElectronProducer.h:90
pat::helper::MultiIsolator::enabled
bool enabled() const
True if it has a non null configuration.
Definition: MultiIsolator.h:55
edm::Ref< PFCandidateCollection >
FileInPath.h
pat::PATElectronProducer::pfElecToken_
const edm::EDGetTokenT< reco::PFCandidateCollection > pfElecToken_
Definition: PATElectronProducer.h:89
pat::PATElectronProducer::embedPflowSuperCluster_
const bool embedPflowSuperCluster_
Definition: PATElectronProducer.h:66
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
pat::helper::MultiIsolator::endEvent
void endEvent()
Definition: MultiIsolator.cc:85
GenParticle.h
pat::PATElectronProducer::embedSeedCluster_
const bool embedSeedCluster_
Definition: PATElectronProducer.h:67
pat::Electron::setDB
void setDB(double dB, double edB, IPTYPE type)
Set impact parameter of a certain type and its uncertainty.
fileCollector.seed
seed
Definition: fileCollector.py:127
BeamMonitor_cff.primaryVertex
primaryVertex
hltOfflineBeamSpot for HLTMON
Definition: BeamMonitor_cff.py:7
DetId
Definition: DetId.h:17
electronIdMVAProducer_cfi.reducedEBRecHitCollection
reducedEBRecHitCollection
Definition: electronIdMVAProducer_cfi.py:7
MakerMacros.h
reco::GsfElectron::setP4
void setP4(P4Kind kind, const LorentzVector &p4, float p4Error, bool setCandidate)
Definition: GsfElectron.cc:188
pat::PATElectronProducer::usePfCandidateMultiMap_
const bool usePfCandidateMultiMap_
Definition: PATElectronProducer.h:88
edm::EventSetup::get
T get() const
Definition: EventSetup.h:80
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
BeamSpot.h
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
pat::PATElectronProducer::pfCandidateMultiMapToken_
const edm::EDGetTokenT< edm::ValueMap< std::vector< reco::PFCandidateRef > > > pfCandidateMultiMapToken_
Definition: PATElectronProducer.h:91
pat::getMiniPFIsolation
PFIsolation getMiniPFIsolation(const pat::PackedCandidateCollection *pfcands, const reco::Candidate::PolarLorentzVector &p4, float mindr=0.05, float maxdr=0.2, float kt_scale=10.0, float ptthresh=0.5, float deadcone_ch=0.0001, float deadcone_pu=0.01, float deadcone_ph=0.01, float deadcone_nh=0.01, float dZ_cut=0.0)
Definition: MiniIsolation.cc:19
l1t::PFCandidateRef
edm::Ref< l1t::PFCandidateCollection > PFCandidateRef
Definition: PFCandidate.h:58
pat::PATElectronProducer::isoDepositTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< IsoDeposit > > > isoDepositTokens_
Definition: PATElectronProducer.h:169
pat::PATElectronProducer::isolationValueTokens_
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueTokens_
Definition: PATElectronProducer.h:171
pat::Electron::embedPflowPreshowerClusters
void embedPflowPreshowerClusters()
method to store the electron's pflow preshower clusters
names
const std::string names[nVars_]
Definition: PhotonIDValueMapProducer.cc:124
GlobalPosition_Frontier_DevDB_cff.tag
tag
Definition: GlobalPosition_Frontier_DevDB_cff.py:11
IPTools::signedImpactParameter3D
std::pair< bool, Measurement1D > signedImpactParameter3D(const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)
Definition: IPTools.cc:81
pat::PATUserDataHelper::add
void add(ObjectType &patObject, edm::Event const &iEvent, edm::EventSetup const &iSetup)
Definition: PATUserDataHelper.h:114
electronIdMVAProducer_cfi.reducedEERecHitCollection
reducedEERecHitCollection
Definition: electronIdMVAProducer_cfi.py:8
pat::PATElectronProducer::embedBasicClusters_
const bool embedBasicClusters_
Definition: PATElectronProducer.h:68
pat::helper::EfficiencyLoader::newEvent
void newEvent(const edm::Event &event)
To be called for each new event, reads in the ValueMaps for efficiencies.
Definition: EfficiencyLoader.cc:21
pat::PATElectronProducer::embedGsfElectronCore_
const bool embedGsfElectronCore_
Definition: PATElectronProducer.h:63
pat::Electron::PV3D
Definition: Electron.h:190
GenParticleFwd.h
ConversionTools::hasMatchedConversion
static bool hasMatchedConversion(const reco::GsfElectron &ele, const reco::ConversionCollection &convCol, const math::XYZPoint &beamspot, bool allowCkfMatch=true, float lxyMin=2.0, float probMin=1e-6, unsigned int nHitsBeforeVtxMax=0)
Definition: ConversionTools.cc:183
reco::BeamSpot
Definition: BeamSpot.h:21
TransientTrackRecord
Definition: TransientTrackRecord.h:11
source
static const std::string source
Definition: EdmProvDump.cc:47
pat::helper::KinResolutionsLoader
Definition: KinResolutionsLoader.h:18
pat::Lepton::setIsoDeposit
void setIsoDeposit(IsolationKeys key, const IsoDeposit &dep)
Sets the IsoDeposit associated with some key; if it is already existent, it is overwritten.
Definition: Lepton.h:191
edm::ESHandle< CaloTopology >
PATElectronProducer.h
beamspot
Definition: BeamSpotWrite2Txt.h:8
pat::Electron::embedBasicClusters
void embedBasicClusters()
method to store the electron's basic clusters
MiniIsolation.h
CaloIsolationEnergy.h
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:59
ParameterSetDescription.h
EcalClusterLazyTools
pat::Electron::embedSuperCluster
void embedSuperCluster()
method to store the electron's SuperCluster internally
CaloGeometryRecord.h
edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:28
edm::vector_transform
auto vector_transform(std::vector< InputType > const &input, Function predicate) -> std::vector< typename std::remove_cv< typename std::remove_reference< decltype(predicate(input.front()))>::type >::type >
Definition: transform.h:11
Vertex.h
EcalSubdetector.h
EcalEndcap
Definition: EcalSubdetector.h:10
FastTrackerRecHitMaskProducer_cfi.recHits
recHits
Definition: FastTrackerRecHitMaskProducer_cfi.py:8
pat::Electron::setMvaVariables
void setMvaVariables(double sigmaIetaIphi, double ip3d)
set missing mva input variables
pat::helper::KinResolutionsLoader::setResolutions
void setResolutions(pat::PATObject< T > &obj) const
Sets the efficiencies for this object, using the reference to the original objects.
Definition: KinResolutionsLoader.h:49
pat::PATElectronProducer::PUPPIIsolation_charged_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_charged_hadrons_
Definition: PATElectronProducer.h:183
pat::Electron::setIsPF
void setIsPF(bool hasPFCandidate)
Definition: Electron.h:174
TransientTrackBuilder.h
pat::PATElectronProducer::hcalPFClusterIsoT_
const edm::EDGetTokenT< edm::ValueMap< float > > hcalPFClusterIsoT_
Definition: PATElectronProducer.h:105
edm::ParameterSet
Definition: ParameterSet.h:47
edm::ParameterSetDescription::setComment
void setComment(std::string const &value)
Definition: ParameterSetDescription.cc:33
edm::Ptr::id
ProductID id() const
Accessor for product ID.
Definition: Ptr.h:158
electronIsolatorFromEffectiveArea_cfi.pfElectrons
pfElectrons
Definition: electronIsolatorFromEffectiveArea_cfi.py:5
pat::helper::KinResolutionsLoader::fillDescription
static void fillDescription(edm::ParameterSetDescription &iDesc)
Method for documentation and validation of PSet.
Definition: KinResolutionsLoader.cc:34
CaloTopology::getSubdetectorTopology
const CaloSubdetectorTopology * getSubdetectorTopology(const DetId &id) const
access the subdetector Topology for the given subdetector directly
Definition: CaloTopology.cc:17
pat::PATElectronProducer::embedSuperCluster_
const bool embedSuperCluster_
Definition: PATElectronProducer.h:65
reco::GsfElectron::setPfIsolationVariables
void setPfIsolationVariables(const PflowIsolationVariables &iso)
Definition: GsfElectron.h:667
CaloTopologyRecord.h
edm::ParameterSetDescription::ifValue
ParameterDescriptionNode * ifValue(ParameterDescription< T > const &switchParameter, std::unique_ptr< ParameterDescriptionCases< T >> cases)
Definition: ParameterSetDescription.h:220
pat::TrackIso
Definition: Isolation.h:10
electrons_cff.ip3d
ip3d
Definition: electrons_cff.py:384
pat::PATElectronProducer::reducedEndcapRecHitCollectionToken_
const edm::EDGetTokenT< EcalRecHitCollection > reducedEndcapRecHitCollectionToken_
Definition: PATElectronProducer.h:99
pat::PATElectronProducer::embedPflowPreshowerClusters_
const bool embedPflowPreshowerClusters_
Definition: PATElectronProducer.h:71
edm::ParameterSet::getParameterNamesForType
std::vector< std::string > getParameterNamesForType(bool trackiness=true) const
Definition: ParameterSet.h:179
iEvent
int iEvent
Definition: GenABIO.cc:224
pat::HcalIso
Definition: Isolation.h:12
pat::PATElectronProducer::addElecID_
const bool addElecID_
Definition: PATElectronProducer.h:158
GsfTrack.h
trackerHitRTTI::vector
Definition: trackerHitRTTI.h:21
pat::PATElectronProducer::electronToken_
const edm::EDGetTokenT< edm::View< reco::GsfElectron > > electronToken_
Definition: PATElectronProducer.h:61
pat::Lepton::setMiniPFIsolation
void setMiniPFIsolation(PFIsolation const &iso)
Definition: Lepton.h:217
pat::PATElectronProducer::ecalTopology_
const CaloTopology * ecalTopology_
Definition: PATElectronProducer.h:192
edm::Association
Definition: Association.h:18
edm::EventSetup
Definition: EventSetup.h:57
edm::ParameterSetDescription::addNode
ParameterDescriptionNode * addNode(ParameterDescriptionNode const &node)
Definition: ParameterSetDescription.cc:41
pat::PATElectronProducer::embedPFCandidate_
const bool embedPFCandidate_
Definition: PATElectronProducer.h:92
pat
Definition: HeavyIon.h:7
edm::contains
bool contains(EventRange const &lh, EventID const &rh)
Definition: EventRange.cc:37
edm::ParameterSetDescription::setAllowAnything
void setAllowAnything()
allow any parameter label/value pairs
Definition: ParameterSetDescription.cc:37
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
TransientTrackRecord.h
DetId::Ecal
Definition: DetId.h:27
pat::PATElectronProducer::hConversionsToken_
const edm::EDGetTokenT< reco::ConversionCollection > hConversionsToken_
Definition: PATElectronProducer.h:62
get
#define get
pat::helper::KinResolutionsLoader::newEvent
void newEvent(const edm::Event &event, const edm::EventSetup &setup)
To be called for each new event, reads in the EventSetup object.
Definition: KinResolutionsLoader.cc:27
reco::GsfElectron::p4
const LorentzVector & p4(P4Kind kind) const
Definition: GsfElectron.cc:211
pat::PATElectronProducer::~PATElectronProducer
~PATElectronProducer() override
Definition: PATElectronProducer.cc:230
l1t::PFCandidateCollection
std::vector< l1t::PFCandidate > PFCandidateCollection
Definition: PFCandidate.h:57
pat::Electron::embedPFCandidate
void embedPFCandidate()
embed the PFCandidate pointed to by pfCandidateRef_
EcalClusterLazyToolsBase::ESGetTokens::get
ESData get(edm::EventSetup const &eventSetup) const
Definition: EcalClusterLazyTools.h:62
pat::PATElectronProducer::PUPPIIsolation_photons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_photons_
Definition: PATElectronProducer.h:185
pat::PATElectronProducer::addMVAVariables_
const bool addMVAVariables_
mva input variables
Definition: PATElectronProducer.h:95
pat::PATElectronProducer::useParticleFlow_
const bool useParticleFlow_
pflow specific
Definition: PATElectronProducer.h:87
edm::Ptr< PFCandidate >
pat::PATElectronProducer::addPuppiIsolation_
const bool addPuppiIsolation_
Definition: PATElectronProducer.h:103
pat::PackedCandidateCollection
std::vector< pat::PackedCandidate > PackedCandidateCollection
Definition: PackedCandidate.h:1130
pat::PATElectronProducer::PATElectronProducer
PATElectronProducer(const edm::ParameterSet &iConfig)
Definition: PATElectronProducer.cc:52
pat::PATElectronProducer::isolator_
pat::helper::MultiIsolator isolator_
Definition: PATElectronProducer.h:166
CaloTopology.h
ValueMap.h
pat::PATElectronProducer::IsolationValueMaps
std::vector< edm::Handle< edm::ValueMap< double > > > IsolationValueMaps
Definition: PATElectronProducer.h:114
reco::TrackBase::Point
math::XYZPoint Point
point in the space
Definition: TrackBase.h:80
CaloSubdetectorTopology::getWindow
virtual std::vector< DetId > getWindow(const DetId &id, const int &northSouthSize, const int &eastWestSize) const
Definition: CaloSubdetectorTopology.cc:4
pat::PATElectronProducer::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: PATElectronProducer.cc:1023
pat::PFIsolation
Definition: PFIsolation.h:12
pat::PATElectronProducer::embedTrack_
const bool embedTrack_
Definition: PATElectronProducer.h:72
reco::Matched
Definition: TrackInfoEnum.h:16
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
IPTools.h
pat::PATElectronProducer::pcToken_
edm::EDGetTokenT< pat::PackedCandidateCollection > pcToken_
Definition: PATElectronProducer.h:77
GsfTrackFwd.h
reco::TransientTrack
Definition: TransientTrack.h:19
pat::helper::MultiIsolator::beginEvent
void beginEvent(const edm::Event &event, const edm::EventSetup &eventSetup)
Definition: MultiIsolator.cc:79
pat::PATElectronProducer::miniIsoParamsB_
std::vector< double > miniIsoParamsB_
Definition: PATElectronProducer.h:80
tier0.unique
def unique(seq, keepstr=True)
Definition: tier0.py:24
isFinite.h
pwdgSkimBPark_cfi.electrons
electrons
Definition: pwdgSkimBPark_cfi.py:6
pat::helper::EfficiencyLoader::setEfficiencies
void setEfficiencies(pat::PATObject< T > &obj, const R &originalRef) const
Sets the efficiencies for this object, using the reference to the original objects.
Definition: EfficiencyLoader.h:41
pat::PATElectronProducer::fillElectron
void fillElectron(Electron &aElectron, const ElectronBaseRef &electronRef, const reco::CandidateBaseRef &baseRef, const GenAssociations &genMatches, const IsoDepositMaps &deposits, const bool pfId, const IsolationValueMaps &isolationValues, const IsolationValueMaps &isolationValuesNoPFId) const
common electron filling, for both the standard and PF2PAT case
Definition: PATElectronProducer.cc:818
packedPFCandidateRefMixer_cfi.pf
pf
Definition: packedPFCandidateRefMixer_cfi.py:4
pat::PATElectronProducer::addGenMatch_
bool addGenMatch_
Definition: PATElectronProducer.h:73
transform.h
edm::ValueMap< float >
reco::GsfElectron::PflowIsolationVariables::sumHcalClusterEt
float sumHcalClusterEt
Definition: GsfElectron.h:618
Exception
Definition: hltDiff.cc:246
edm::RefVectorIterator
Definition: EDProductfwd.h:33
pat::PATElectronProducer::ecalPFClusterIsoT_
const edm::EDGetTokenT< edm::ValueMap< float > > ecalPFClusterIsoT_
Definition: PATElectronProducer.h:104
reco::HitPattern::MISSING_INNER_HITS
Definition: HitPattern.h:155
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
pat::PATElectronProducer::beamLineToken_
const edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
Definition: PATElectronProducer.h:109
pat::Electron::PVDZ
Definition: Electron.h:190
edm::Ptr::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ptr.h:146
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::RefToBase
Definition: AssociativeIterator.h:54
pat::PATElectronProducer::PUPPINoLeptonsIsolation_neutral_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPINoLeptonsIsolation_neutral_hadrons_
Definition: PATElectronProducer.h:188
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
patElectronEAIsoCorrectionProducer_cfi.patElectrons
patElectrons
Definition: patElectronEAIsoCorrectionProducer_cfi.py:4
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:20
pat::PATElectronProducer::GenAssociations
std::vector< edm::Handle< edm::Association< reco::GenParticleCollection > > > GenAssociations
Definition: PATElectronProducer.h:82
electronProducer_cfi.usePfCandidateMultiMap
usePfCandidateMultiMap
Definition: electronProducer_cfi.py:11
mps_fire.result
result
Definition: mps_fire.py:311
pat::Electron
Analysis-level electron class.
Definition: Electron.h:51
pat::PATElectronProducer::ecalClusterToolsESGetTokens_
const EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
Definition: PATElectronProducer.h:100
reco::LeafCandidate::polarP4
const PolarLorentzVector & polarP4() const final
four-momentum Lorentz vector
Definition: LeafCandidate.h:116
reco::GsfElectron::pfIsolationVariables
const PflowIsolationVariables & pfIsolationVariables() const
Definition: GsfElectron.h:658
edm::ParameterDescriptionNode::setComment
void setComment(std::string const &value)
Definition: ParameterDescriptionNode.cc:106
pat::PATElectronProducer::isoDepositLabels_
IsolationLabels isoDepositLabels_
Definition: PATElectronProducer.h:168
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
pat::Electron::setPFCandidateRef
void setPFCandidateRef(const reco::PFCandidateRef &ref)
add a reference to the source IsolatedPFCandidate
Definition: Electron.h:179
pat::PATObject::embedGenParticle
void embedGenParticle()
Definition: PATObject.h:770
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
edm::Event
Definition: Event.h:73
pat::PATElectronProducer::miniIsoParamsE_
std::vector< double > miniIsoParamsE_
Definition: PATElectronProducer.h:79
pat::Electron::setElectronIDs
void setElectronIDs(const std::vector< IdPair > &ids)
Store multiple electron ID values, discarding existing ones. The first one in the list becomes the 'd...
Definition: Electron.h:141
edm::EmptyGroupDescription
Definition: EmptyGroupDescription.h:15
edm::ParameterDescription
Definition: ParameterDescription.h:110
CandIsolatorFromDeposits_cfi.deposits
deposits
Definition: CandIsolatorFromDeposits_cfi.py:4
reco::GsfElectron::PflowIsolationVariables::sumEcalClusterEt
float sumEcalClusterEt
Definition: GsfElectron.h:617
edm::InputTag
Definition: InputTag.h:15
pat::PATElectronProducer::efficiencyLoader_
pat::helper::EfficiencyLoader efficiencyLoader_
Definition: PATElectronProducer.h:176
reco::Vertex
Definition: Vertex.h:35
pat::PATElectronProducer
Produces pat::Electron's.
Definition: PATElectronProducer.h:50
hit
Definition: SiStripHitEffFromCalibTree.cc:88
pat::PATElectronProducer::computeMiniIso_
bool computeMiniIso_
Definition: PATElectronProducer.h:78
PFCandidateFwd.h
ConversionTools.h
pat::PATElectronProducer::reducedBarrelRecHitCollectionToken_
const edm::EDGetTokenT< EcalRecHitCollection > reducedBarrelRecHitCollectionToken_
Definition: PATElectronProducer.h:97
pat::PATElectronProducer::PUPPIIsolation_neutral_hadrons_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_neutral_hadrons_
Definition: PATElectronProducer.h:184
Generated for CMSSW Reference Manual by doxygen 1.8.16