CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
LowPtGsfElectronSeedProducer Class Referencefinal

#include <RecoEgamma/EgammaElectronProducers/plugins/LowPtGsfElectronSeedProducer.cc>

Inheritance diagram for LowPtGsfElectronSeedProducer:
edm::stream::EDProducer< edm::GlobalCache< lowptgsfeleseed::HeavyObjectCache > >

Public Types

using TrackIndxMap = std::unordered_map< reco::TrackRef::key_type, size_t >
 
- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< lowptgsfeleseed::HeavyObjectCache > >
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
 

Public Member Functions

void beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override
 
 LowPtGsfElectronSeedProducer (const edm::ParameterSet &, const lowptgsfeleseed::HeavyObjectCache *)
 
void produce (edm::Event &, const edm::EventSetup &) override
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< lowptgsfeleseed::HeavyObjectCache > >
 EDProducer ()=default
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &)
 
static void globalEndJob (lowptgsfeleseed::HeavyObjectCache const *)
 
static std::unique_ptr< lowptgsfeleseed::HeavyObjectCache > initializeGlobalCache (const edm::ParameterSet &conf)
 

Private Member Functions

bool decision (const reco::PFRecTrackRef &pfTrackRef, reco::PreId &ecal, reco::PreId &hcal, double rho, const reco::BeamSpot &spot, noZS::EcalClusterLazyTools &ecalTools)
 
bool decision (const reco::TrackRef &kfTrackRef, reco::PreId &ecal, reco::PreId &hcal, double rho, const reco::BeamSpot &spot, noZS::EcalClusterLazyTools &ecalTools)
 
template<typename CollType >
void fillPreIdRefValueMap (edm::Handle< CollType > tracksHandle, const TrackIndxMap &trksToPreIdIndx, const edm::OrphanHandle< reco::PreIdCollection > &preIdHandle, edm::ValueMap< reco::PreIdRef >::Filler &filler)
 
reco::TrackRef getBaseRef (edm::Handle< std::vector< reco::PFRecTrack > > handle, int idx) const
 
reco::TrackRef getBaseRef (edm::Handle< std::vector< reco::Track > > handle, int idx) const
 
bool lightGsfTracking (reco::PreId &, const reco::TrackRef &, const reco::ElectronSeed &)
 
template<typename T >
void loop (const edm::Handle< std::vector< T > > &handle, edm::Handle< reco::PFClusterCollection > &hcalClusters, reco::ElectronSeedCollection &seeds, reco::PreIdCollection &ecalPreIds, reco::PreIdCollection &hcalPreIds, TrackIndxMap &trksToPreIdIndx, edm::Event &, const edm::EventSetup &)
 
void propagateTrackToCalo (const reco::PFRecTrackRef &pfTrackRef, const edm::Handle< reco::PFClusterCollection > &clusters, std::vector< int > &matchedClusters, reco::PreId &preId, bool ecal)
 
void propagateTrackToCalo (const reco::PFRecTrackRef &pfTrackRef, const edm::Handle< reco::PFClusterCollection > &ecalClusters, const edm::Handle< reco::PFClusterCollection > &hcalClusters, std::vector< int > &matchedEcalClusters, std::vector< int > &matchedHcalClusters, reco::PreId &ecalPreId, reco::PreId &hcalPreId)
 
void propagateTrackToCalo (const reco::TrackRef &pfTrack, const edm::Handle< reco::PFClusterCollection > &ecalClusters, const edm::Handle< reco::PFClusterCollection > &hcalClusters, std::vector< int > &matchedEcalClusters, std::vector< int > &matchedHcalClusters, reco::PreId &ecalPreId, reco::PreId &hcalPreId)
 

Private Attributes

const edm::EDGetTokenT< reco::BeamSpot > beamSpot_
 
const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord > builderToken_
 
const edm::EDGetTokenT< EcalRecHitCollection > ebRecHits_
 
const edm::EDGetTokenT< reco::PFClusterCollection > ecalClusters_
 
const noZS::EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
 
const edm::EDGetTokenT< EcalRecHitCollection > eeRecHits_
 
edm::ESHandle< MagneticField > field_
 
std::unique_ptr< TrajectoryFitter > fitterPtr_
 
edm::EDGetTokenT< reco::PFClusterCollection > hcalClusters_
 
edm::EDGetTokenT< reco::TrackCollection > kfTracks_
 
const double maxPtThreshold_
 
const double minPtThreshold_
 
const bool passThrough_
 
edm::EDGetTokenT< reco::PFRecTrackCollection > pfTracks_
 
const edm::EDGetTokenT< double > rho_
 
std::unique_ptr< TrajectorySmoother > smootherPtr_
 
const edm::ESGetToken< TrajectoryFitter, TrajectoryFitter::Record > trajectoryFitterToken_
 
const edm::ESGetToken< TrajectorySmoother, TrajectoryFitter::Record > trajectorySmootherToken_
 
const bool usePfTracks_
 

Static Private Attributes

static constexpr double boundary_ = 2.50746495928 * 2.50746495928
 
static constexpr double mass_ = 0.000511 * 0.000511
 

Detailed Description

Description: EDProducer of ElectronSeed objects Implementation: <Notes on="" implementation>="">

Definition at line 57 of file LowPtGsfElectronSeedProducer.cc.

Member Typedef Documentation

◆ TrackIndxMap

using LowPtGsfElectronSeedProducer::TrackIndxMap = std::unordered_map<reco::TrackRef::key_type, size_t>

Definition at line 60 of file LowPtGsfElectronSeedProducer.cc.

Constructor & Destructor Documentation

◆ LowPtGsfElectronSeedProducer()

LowPtGsfElectronSeedProducer::LowPtGsfElectronSeedProducer ( const edm::ParameterSet &  conf,
const lowptgsfeleseed::HeavyObjectCache *   
)
explicit

Definition at line 170 of file LowPtGsfElectronSeedProducer.cc.

172  : field_(),
173  fitterPtr_(),
174  smootherPtr_(),
175  kfTracks_(),
176  pfTracks_(),
177  ecalClusters_{consumes(conf.getParameter<edm::InputTag>("ecalClusters"))},
178  hcalClusters_(),
179  ebRecHits_{consumes(conf.getParameter<edm::InputTag>("EBRecHits"))},
180  eeRecHits_{consumes(conf.getParameter<edm::InputTag>("EERecHits"))},
181  rho_(consumes(conf.getParameter<edm::InputTag>("rho"))),
182  beamSpot_(consumes(conf.getParameter<edm::InputTag>("BeamSpot"))),
183  trajectoryFitterToken_{esConsumes(conf.getParameter<edm::ESInputTag>("Fitter"))},
184  trajectorySmootherToken_{esConsumes(conf.getParameter<edm::ESInputTag>("Smoother"))},
185  builderToken_{esConsumes(conf.getParameter<edm::ESInputTag>("TTRHBuilder"))},
186  ecalClusterToolsESGetTokens_{consumesCollector()},
187  passThrough_(conf.getParameter<bool>("PassThrough")),
188  usePfTracks_(conf.getParameter<bool>("UsePfTracks")),
189  minPtThreshold_(conf.getParameter<double>("MinPtThreshold")),
190  maxPtThreshold_(conf.getParameter<double>("MaxPtThreshold")) {
191  if (usePfTracks_) {
192  pfTracks_ = consumes(conf.getParameter<edm::InputTag>("pfTracks"));
193  hcalClusters_ = consumes(conf.getParameter<edm::InputTag>("hcalClusters"));
194  }
195  kfTracks_ = consumes(conf.getParameter<edm::InputTag>("tracks"));
196 
197  produces<reco::ElectronSeedCollection>();
198  produces<reco::PreIdCollection>();
199  produces<reco::PreIdCollection>("HCAL");
200  produces<edm::ValueMap<reco::PreIdRef> >(); // indexed by edm::Ref<ElectronSeed>.index()
201 }

References edm::ParameterSet::getParameter().

Member Function Documentation

◆ beginLuminosityBlock()

void LowPtGsfElectronSeedProducer::beginLuminosityBlock ( edm::LuminosityBlock const &  ,
edm::EventSetup const &  setup 
)
override

Definition at line 205 of file LowPtGsfElectronSeedProducer.cc.

205  {
206  setup.get<IdealMagneticFieldRecord>().get(field_);
207 }

References field_, get, and singleTopDQM_cfi::setup.

◆ decision() [1/2]

bool LowPtGsfElectronSeedProducer::decision ( const reco::PFRecTrackRef &  pfTrackRef,
reco::PreId &  ecal,
reco::PreId &  hcal,
double  rho,
const reco::BeamSpot &  spot,
noZS::EcalClusterLazyTools &  ecalTools 
)
private

Definition at line 595 of file LowPtGsfElectronSeedProducer.cc.

600  {
601  bool result = false;
602  for (auto& name : globalCache()->modelNames()) {
603  result |= globalCache()->eval(name, ecalPreId, hcalPreId, rho, spot, ecalTools);
604  }
605  return passThrough_ || (pfTrackRef->trackRef()->pt() > maxPtThreshold_) || result;
606 }

References maxPtThreshold_, Skims_PA_cff::name, passThrough_, mps_fire::result, and rho.

Referenced by loop().

◆ decision() [2/2]

bool LowPtGsfElectronSeedProducer::decision ( const reco::TrackRef &  kfTrackRef,
reco::PreId &  ecal,
reco::PreId &  hcal,
double  rho,
const reco::BeamSpot &  spot,
noZS::EcalClusterLazyTools &  ecalTools 
)
private

Definition at line 610 of file LowPtGsfElectronSeedProducer.cc.

615  {
616  // No implementation currently
617  return passThrough_;
618 }

References passThrough_.

◆ fillDescriptions()

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

Definition at line 643 of file LowPtGsfElectronSeedProducer.cc.

643  {
644  edm::ParameterSetDescription desc;
645  desc.add<edm::InputTag>("tracks", edm::InputTag("generalTracks"));
646  desc.add<edm::InputTag>("pfTracks", edm::InputTag("lowPtGsfElePfTracks"));
647  desc.add<edm::InputTag>("ecalClusters", edm::InputTag("particleFlowClusterECAL"));
648  desc.add<edm::InputTag>("hcalClusters", edm::InputTag("particleFlowClusterHCAL"));
649  desc.add<edm::InputTag>("EBRecHits", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
650  desc.add<edm::InputTag>("EERecHits", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
651  desc.add<edm::InputTag>("rho", edm::InputTag("fixedGridRhoFastjetAllTmp"));
652  desc.add<edm::InputTag>("BeamSpot", edm::InputTag("offlineBeamSpot"));
653  desc.add<edm::ESInputTag>("Fitter", edm::ESInputTag{"", "GsfTrajectoryFitter_forPreId"});
654  desc.add<edm::ESInputTag>("Smoother", edm::ESInputTag{"", "GsfTrajectorySmoother_forPreId"});
655  desc.add<edm::ESInputTag>("TTRHBuilder", edm::ESInputTag{"", "WithAngleAndTemplate"});
656  desc.add<std::vector<std::string> >("ModelNames", {});
657  desc.add<std::vector<std::string> >("ModelWeights", {});
658  desc.add<std::vector<double> >("ModelThresholds", {});
659  desc.add<bool>("PassThrough", false);
660  desc.add<bool>("UsePfTracks", true);
661  desc.add<double>("MinPtThreshold", 1.0);
662  desc.add<double>("MaxPtThreshold", 15.);
663  descriptions.add("lowPtGsfElectronSeeds", desc);
664 }

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

◆ fillPreIdRefValueMap()

template<typename CollType >
void LowPtGsfElectronSeedProducer::fillPreIdRefValueMap ( edm::Handle< CollType >  tracksHandle,
const TrackIndxMap &  trksToPreIdIndx,
const edm::OrphanHandle< reco::PreIdCollection > &  preIdHandle,
edm::ValueMap< reco::PreIdRef >::Filler &  filler 
)
private

Definition at line 621 of file LowPtGsfElectronSeedProducer.cc.

624  {
625  std::vector<reco::PreIdRef> values;
626 
627  unsigned ntracks = tracksHandle->size();
628  for (unsigned itrack = 0; itrack < ntracks; ++itrack) {
629  edm::Ref<CollType> trackRef(tracksHandle, itrack);
630  auto preIdRefIt = trksToPreIdIndx.find(trackRef.index());
631  if (preIdRefIt == trksToPreIdIndx.end()) {
632  values.push_back(reco::PreIdRef());
633  } else {
634  edm::Ref<reco::PreIdCollection> preIdRef(preIdHandle, preIdRefIt->second);
635  values.push_back(preIdRef);
636  }
637  }
638  filler.insert(tracksHandle, values.begin(), values.end());
639 }

References trigObjTnPSource_cfi::filler, edm::Ref< C, T, F >::index(), vertices_cff::ntracks, and contentValuesCheck::values.

Referenced by produce().

◆ getBaseRef() [1/2]

reco::TrackRef LowPtGsfElectronSeedProducer::getBaseRef ( edm::Handle< std::vector< reco::PFRecTrack > >  handle,
int  idx 
) const
private

Definition at line 281 of file LowPtGsfElectronSeedProducer.cc.

282  {
283  return reco::PFRecTrackRef(handle, idx)->trackRef();
284 }

References patZpeak::handle, and heavyIonCSV_trainingSettings::idx.

◆ getBaseRef() [2/2]

reco::TrackRef LowPtGsfElectronSeedProducer::getBaseRef ( edm::Handle< std::vector< reco::Track > >  handle,
int  idx 
) const
private

Definition at line 277 of file LowPtGsfElectronSeedProducer.cc.

277  {
278  return reco::TrackRef(handle, idx);
279 }

References patZpeak::handle, and heavyIonCSV_trainingSettings::idx.

Referenced by loop().

◆ globalEndJob()

static void LowPtGsfElectronSeedProducer::globalEndJob ( lowptgsfeleseed::HeavyObjectCache const *  )
inlinestatic

Definition at line 67 of file LowPtGsfElectronSeedProducer.cc.

67 {}

◆ initializeGlobalCache()

static std::unique_ptr<lowptgsfeleseed::HeavyObjectCache> LowPtGsfElectronSeedProducer::initializeGlobalCache ( const edm::ParameterSet &  conf)
inlinestatic

Definition at line 63 of file LowPtGsfElectronSeedProducer.cc.

63  {
64  return std::make_unique<lowptgsfeleseed::HeavyObjectCache>(lowptgsfeleseed::HeavyObjectCache(conf));
65  }

◆ lightGsfTracking()

bool LowPtGsfElectronSeedProducer::lightGsfTracking ( reco::PreId &  preId,
const reco::TrackRef &  trackRef,
const reco::ElectronSeed &  seed 
)
private

Definition at line 557 of file LowPtGsfElectronSeedProducer.cc.

559  {
560  Trajectory::ConstRecHitContainer hits;
561  for (unsigned int ihit = 0; ihit < trackRef->recHitsSize(); ++ihit) {
562  hits.push_back(trackRef->recHit(ihit)->cloneSH());
563  }
564 
565  GlobalVector gv(trackRef->innerMomentum().x(), trackRef->innerMomentum().y(), trackRef->innerMomentum().z());
566  GlobalPoint gp(trackRef->innerPosition().x(), trackRef->innerPosition().y(), trackRef->innerPosition().z());
567 
568  GlobalTrajectoryParameters gtps(gp, gv, trackRef->charge(), &*field_);
569 
570  TrajectoryStateOnSurface tsos(gtps, trackRef->innerStateCovariance(), *hits[0]->surface());
571 
572  // Track fitted and smoothed under electron hypothesis
573  Trajectory traj1 = fitterPtr_->fitOne(seed, hits, tsos);
574  if (!traj1.isValid()) {
575  return false;
576  }
577  Trajectory traj2 = smootherPtr_->trajectory(traj1);
578  if (!traj2.isValid()) {
579  return false;
580  }
581 
582  // Set PreId content
583  float chi2Ratio = trackRef->chi2() > 0. ? traj2.chiSquared() / trackRef->chi2() : -1.;
584  float gsfReducedChi2 = chi2Ratio > -1. ? chi2Ratio * trackRef->normalizedChi2() : -1.;
585  float ptOut = traj2.firstMeasurement().updatedState().globalMomentum().perp();
586  float ptIn = traj2.lastMeasurement().updatedState().globalMomentum().perp();
587  float gsfDpt = (ptIn > 0) ? fabs(ptOut - ptIn) / ptIn : 0.;
588  preId.setTrackProperties(gsfReducedChi2, chi2Ratio, gsfDpt);
589 
590  return true;
591 }

References Trajectory::chiSquared(), field_, Trajectory::firstMeasurement(), fitterPtr_, TrajectoryStateOnSurface::globalMomentum(), runTauDisplay::gp, hfClusterShapes_cfi::hits, Trajectory::isValid(), Trajectory::lastMeasurement(), PV3DBase< T, PVType, FrameType >::perp(), fileCollector::seed, reco::PreId::setTrackProperties(), smootherPtr_, and TrajectoryMeasurement::updatedState().

Referenced by loop().

◆ loop()

template<typename T >
void LowPtGsfElectronSeedProducer::loop ( const edm::Handle< std::vector< T > > &  handle,
edm::Handle< reco::PFClusterCollection > &  hcalClusters,
reco::ElectronSeedCollection &  seeds,
reco::PreIdCollection &  ecalPreIds,
reco::PreIdCollection &  hcalPreIds,
TrackIndxMap &  trksToPreIdIndx,
edm::Event &  event,
const edm::EventSetup &  setup 
)
private

Definition at line 289 of file LowPtGsfElectronSeedProducer.cc.

296  {
297  // Pileup
298  auto const& rho = event.get(rho_);
299 
300  // Beam spot
301  auto const& spot = event.get(beamSpot_);
302 
303  // Track fitter
304  fitterPtr_ = setup.getData(trajectoryFitterToken_).clone();
305 
306  // Track smoother
307  smootherPtr_.reset(setup.getData(trajectorySmootherToken_).clone());
308 
309  // RecHit cloner
310  TkClonerImpl hitCloner = static_cast<TkTransientTrackingRecHitBuilder const&>(setup.getData(builderToken_)).cloner();
311  fitterPtr_->setHitCloner(&hitCloner);
312  smootherPtr_->setHitCloner(&hitCloner);
313 
314  // ECAL clusters
315  auto ecalClusters = event.getHandle(ecalClusters_);
316 
317  // Utility to access to shower shape vars
318  noZS::EcalClusterLazyTools ecalTools(event, ecalClusterToolsESGetTokens_.get(setup), ebRecHits_, eeRecHits_);
319 
320  // Ensure each cluster is only matched once to a track
321  std::vector<int> matchedEcalClusters;
322  std::vector<int> matchedHcalClusters;
323 
324  // Reserve
325  seeds.reserve(handle->size());
326  ecalPreIds.reserve(handle->size());
327  hcalPreIds.reserve(handle->size());
328 
329  // Iterate through (PF or KF) tracks
330  for (unsigned int itrk = 0; itrk < handle.product()->size(); itrk++) {
331  edm::Ref<std::vector<T> > templatedRef(handle, itrk); // TrackRef or PFRecTrackRef
332  reco::TrackRef trackRef = getBaseRef(handle, itrk);
333 
334  if (!(trackRef->quality(reco::TrackBase::qualityByName("highPurity")))) {
335  continue;
336  }
337  if (!passThrough_ && (trackRef->pt() < minPtThreshold_)) {
338  continue;
339  }
340 
341  // Create ElectronSeed
342  reco::ElectronSeed seed(*(trackRef->seedRef()));
343  seed.setCtfTrack(trackRef);
344 
345  // Create PreIds
346  unsigned int nModels = globalCache()->modelNames().size();
347  reco::PreId ecalPreId(nModels);
348  reco::PreId hcalPreId(nModels);
349 
350  // Add track ref to PreId
351  ecalPreId.setTrack(trackRef);
352  hcalPreId.setTrack(trackRef);
353 
354  // Add Track-Calo matching variables to PreIds
355  propagateTrackToCalo(
356  templatedRef, ecalClusters, hcalClusters, matchedEcalClusters, matchedHcalClusters, ecalPreId, hcalPreId);
357 
358  // Add variables related to GSF tracks to PreId
359  lightGsfTracking(ecalPreId, trackRef, seed);
360 
361  // Decision based on BDT
362  bool result = decision(templatedRef, ecalPreId, hcalPreId, rho, spot, ecalTools);
363 
364  // If fails BDT, do not store seed
365  if (!result) {
366  continue;
367  }
368 
369  // Store PreId
370  ecalPreIds.push_back(ecalPreId);
371  hcalPreIds.push_back(hcalPreId);
372  trksToPreIdIndx[trackRef.index()] = ecalPreIds.size() - 1;
373 
374  // Store ElectronSeed and corresponding edm::Ref<PreId>.index()
375  seeds.push_back(seed);
376  }
377 }

References beamSpot_, builderToken_, decision(), ebRecHits_, RecoEcal_cff::ecalClusters, ecalClusters_, ecalClusterToolsESGetTokens_, eeRecHits_, fitterPtr_, EcalClusterLazyToolsBase::ESGetTokens::get(), getBaseRef(), patZpeak::handle, lowPtGsfElectronSeeds_cfi::hcalClusters, edm::Ref< C, T, F >::index(), lightGsfTracking(), minPtThreshold_, passThrough_, propagateTrackToCalo(), reco::TrackBase::qualityByName(), mps_fire::result, rho, rho_, fileCollector::seed, InitialStep_cff::seeds, reco::PreId::setTrack(), singleTopDQM_cfi::setup, smootherPtr_, trajectoryFitterToken_, and trajectorySmootherToken_.

Referenced by produce().

◆ produce()

void LowPtGsfElectronSeedProducer::produce ( edm::Event &  event,
const edm::EventSetup &  setup 
)
override

Definition at line 211 of file LowPtGsfElectronSeedProducer.cc.

211  {
212  // Products
213  auto seeds = std::make_unique<reco::ElectronSeedCollection>();
214  auto ecalPreIds = std::make_unique<reco::PreIdCollection>();
215  auto hcalPreIds = std::make_unique<reco::PreIdCollection>();
216 
217  const edm::RefProd<reco::PreIdCollection> preIdsRefProd = event.getRefBeforePut<reco::PreIdCollection>();
218 
219  // HCAL clusters (only used with PF tracks)
220  edm::Handle<reco::PFClusterCollection> hcalClusters;
221 
222  //we always need kftracks as we link the preids to them
223  edm::Handle<reco::TrackCollection> kfTracks;
224  event.getByToken(kfTracks_, kfTracks);
225 
226  TrackIndxMap trksToPreIdIndx;
227  if (usePfTracks_) {
228  edm::Handle<reco::PFRecTrackCollection> pfTracks;
229  event.getByToken(pfTracks_, pfTracks);
230  event.getByToken(hcalClusters_, hcalClusters);
231 
232  //check consistency between kfTracks and pfTracks collection
233  for (auto& trk : *pfTracks) {
234  if (trk.trackRef().isNonnull()) {
235  if (trk.trackRef().id() != kfTracks.id()) {
236  throw cms::Exception("ConfigError")
237  << "kfTracks is not the collection that pfTracks was built from, please fix this";
238  }
239  break; //we only need one valid track ref to check this so end the loop
240  }
241  }
242 
243  loop(pfTracks, // PF tracks
244  hcalClusters,
245  *seeds,
246  *ecalPreIds,
247  *hcalPreIds,
248  trksToPreIdIndx,
249  event,
250  setup);
251 
252  } else {
253  loop(kfTracks, // KF tracks
254  hcalClusters,
255  *seeds,
256  *ecalPreIds,
257  *hcalPreIds,
258  trksToPreIdIndx,
259  event,
260  setup);
261  }
262 
263  auto ecalPreIdsHandle = event.put(std::move(ecalPreIds));
264  event.put(std::move(hcalPreIds), "HCAL");
265  event.put(std::move(seeds));
266 
267  auto preIdVMOut = std::make_unique<edm::ValueMap<reco::PreIdRef> >();
268  edm::ValueMap<reco::PreIdRef>::Filler mapFiller(*preIdVMOut);
269  fillPreIdRefValueMap(kfTracks, trksToPreIdIndx, ecalPreIdsHandle, mapFiller);
270  mapFiller.fill();
271  event.put(std::move(preIdVMOut));
272 }

References Exception, fillPreIdRefValueMap(), lowPtGsfElectronSeeds_cfi::hcalClusters, hcalClusters_, edm::HandleBase::id(), kfTracks_, loop(), eostools::move(), lowPtGsfElectronSeeds_cfi::pfTracks, pfTracks_, InitialStep_cff::seeds, singleTopDQM_cfi::setup, and usePfTracks_.

◆ propagateTrackToCalo() [1/3]

void LowPtGsfElectronSeedProducer::propagateTrackToCalo ( const reco::PFRecTrackRef &  pfTrackRef,
const edm::Handle< reco::PFClusterCollection > &  clusters,
std::vector< int > &  matchedClusters,
reco::PreId &  preId,
bool  ecal 
)
private

Definition at line 416 of file LowPtGsfElectronSeedProducer.cc.

420  {
421  // Store info for PreId
422  struct Info {
423  reco::PFClusterRef cluRef = reco::PFClusterRef();
424  float dr2min = 1.e6;
425  float deta = 1.e6;
426  float dphi = 1.e6;
427  math::XYZPoint showerPos = math::XYZPoint(0., 0., 0.);
428  } info;
429 
430  // Find closest "seed cluster" to KF track extrapolated to ECAL (or HCAL)
431  reco::PFTrajectoryPoint point;
432  if (ecal) {
433  point = pfTrackRef->extrapolatedPoint(reco::PFTrajectoryPoint::LayerType::ECALShowerMax);
434  } else {
435  point = pfTrackRef->extrapolatedPoint(reco::PFTrajectoryPoint::LayerType::HCALEntrance);
436  }
437 
438  if (point.isValid()) {
439  Info info;
440  for (unsigned int iclu = 0; iclu < clusters.product()->size(); iclu++) {
441  if (std::find(matched.begin(), matched.end(), iclu) == matched.end()) {
442  reco::PFClusterRef cluRef(clusters, iclu);
443 
444  // Determine dR squared
445  float dr2 = reco::deltaR2(cluRef->positionREP(), point.positionREP());
446 
447  if (dr2 < info.dr2min) {
448  info.dr2min = dr2;
449  info.cluRef = cluRef;
450  info.deta = cluRef->positionREP().eta() - point.positionREP().eta();
451  info.dphi =
452  reco::deltaPhi(cluRef->positionREP().phi(), point.positionREP().phi()) * pfTrackRef->trackRef()->charge();
453  info.showerPos = point.position();
454  }
455  }
456  }
457 
458  // Set PreId content if match found
459  if (info.dr2min < 1.e5) {
460  float ep = info.cluRef->correctedEnergy() / std::sqrt(pfTrackRef->trackRef()->innerMomentum().mag2());
461  preId.setECALMatchingProperties(info.cluRef,
462  point.position(), // ECAL or HCAL surface
463  info.showerPos, //
464  info.deta,
465  info.dphi,
466  0.f, // chieta
467  0.f, // chiphi
468  pfTrackRef->trackRef()->normalizedChi2(), // chi2
469  ep);
470  }
471 
472  } // clusters
473 }

References bsc_activity_cfg::clusters, reco::deltaPhi(), reco::deltaR2(), bsc_activity_cfg::ecal, SiStripBadComponentsDQMServiceTemplate_cfg::ep, spr::find(), info(), muonTagProbeFilters_cff::matched, point, reco::PreId::setECALMatchingProperties(), and mathSSE::sqrt().

◆ propagateTrackToCalo() [2/3]

void LowPtGsfElectronSeedProducer::propagateTrackToCalo ( const reco::PFRecTrackRef &  pfTrackRef,
const edm::Handle< reco::PFClusterCollection > &  ecalClusters,
const edm::Handle< reco::PFClusterCollection > &  hcalClusters,
std::vector< int > &  matchedEcalClusters,
std::vector< int > &  matchedHcalClusters,
reco::PreId &  ecalPreId,
reco::PreId &  hcalPreId 
)
private

Definition at line 403 of file LowPtGsfElectronSeedProducer.cc.

409  {
410  propagateTrackToCalo(pfTrackRef, ecalClusters, matchedEcalClusters, ecalPreId, true);
411  propagateTrackToCalo(pfTrackRef, hcalClusters, matchedHcalClusters, hcalPreId, false);
412 }

References RecoEcal_cff::ecalClusters, and lowPtGsfElectronSeeds_cfi::hcalClusters.

Referenced by loop().

◆ propagateTrackToCalo() [3/3]

void LowPtGsfElectronSeedProducer::propagateTrackToCalo ( const reco::TrackRef &  pfTrack,
const edm::Handle< reco::PFClusterCollection > &  ecalClusters,
const edm::Handle< reco::PFClusterCollection > &  hcalClusters,
std::vector< int > &  matchedEcalClusters,
std::vector< int > &  matchedHcalClusters,
reco::PreId &  ecalPreId,
reco::PreId &  hcalPreId 
)
private

Definition at line 477 of file LowPtGsfElectronSeedProducer.cc.

484  {
485  // Store info for PreId
486  struct Info {
487  reco::PFClusterRef cluRef = reco::PFClusterRef();
488  float dr2min = 1.e6;
489  float deta = 1.e6;
490  float dphi = 1.e6;
491  math::XYZPoint showerPos = math::XYZPoint(0., 0., 0.);
492  } info;
493 
494  // Propagate 'electron' to ECAL surface
495  float energy = sqrt(mass_ + kfTrackRef->outerMomentum().Mag2());
496  XYZTLorentzVector mom = XYZTLorentzVector(
497  kfTrackRef->outerMomentum().x(), kfTrackRef->outerMomentum().y(), kfTrackRef->outerMomentum().z(), energy);
498  XYZTLorentzVector pos = XYZTLorentzVector(
499  kfTrackRef->outerPosition().x(), kfTrackRef->outerPosition().y(), kfTrackRef->outerPosition().z(), 0.);
500  math::XYZVector field(field_->inTesla(GlobalPoint(0, 0, 0)));
501  BaseParticlePropagator particle(RawParticle(mom, pos, kfTrackRef->charge()), 0, 0, field.z());
502  particle.propagateToEcalEntrance(false);
503  if (particle.getSuccess() == 0) {
504  return;
505  }
506 
507  // ECAL entry point for track
508  GlobalPoint ecal_pos(
509  particle.particle().vertex().x(), particle.particle().vertex().y(), particle.particle().vertex().z());
510 
511  // Preshower limit
512  bool below_ps = pow(ecal_pos.z(), 2.) > boundary_ * ecal_pos.perp2();
513 
514  // Iterate through ECAL clusters
515  for (unsigned int iclu = 0; iclu < ecalClusters.product()->size(); iclu++) {
516  reco::PFClusterRef cluRef(ecalClusters, iclu);
517 
518  // Correct ecal_pos for shower depth
519  double shower_depth = reco::PFCluster::getDepthCorrection(cluRef->correctedEnergy(), below_ps, false);
520  GlobalPoint showerPos = ecal_pos + GlobalVector(particle.particle().momentum().x(),
521  particle.particle().momentum().y(),
522  particle.particle().momentum().z())
523  .unit() *
524  shower_depth;
525 
526  // Determine dR squared
527  float dr2 = reco::deltaR2(cluRef->positionREP(), showerPos);
528 
529  // Find nearest ECAL cluster
530  if (dr2 < info.dr2min) {
531  info.dr2min = dr2;
532  info.cluRef = cluRef;
533  info.deta = std::abs(cluRef->positionREP().eta() - showerPos.eta());
534  info.dphi = std::abs(reco::deltaPhi(cluRef->positionREP().phi(), showerPos.phi())) * kfTrackRef->charge();
535  info.showerPos = showerPos;
536  }
537  }
538 
539  // Populate PreId object
540  math::XYZPoint point(ecal_pos.x(), ecal_pos.y(), ecal_pos.z());
541 
542  // Set PreId content
543  ecalPreId.setECALMatchingProperties(
544  info.cluRef,
545  point,
546  info.showerPos,
547  info.deta,
548  info.dphi,
549  0.f, // chieta
550  0.f, // chiphi
551  kfTrackRef->normalizedChi2(), // chi2
552  info.cluRef->correctedEnergy() / std::sqrt(kfTrackRef->innerMomentum().mag2())); // E/p
553 }

References funct::abs(), boundary_, reco::deltaPhi(), reco::deltaR2(), RecoEcal_cff::ecalClusters, HCALHighEnergyHPDFilter_cfi::energy, field_, reco::PFCluster::getDepthCorrection(), info(), MagneticField::inTesla(), mass_, point, funct::pow(), BaseParticlePropagator::propagateToEcalEntrance(), reco::PreId::setECALMatchingProperties(), mathSSE::sqrt(), unit(), and RawParticle::z().

Member Data Documentation

◆ beamSpot_

const edm::EDGetTokenT<reco::BeamSpot> LowPtGsfElectronSeedProducer::beamSpot_
private

Definition at line 150 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ boundary_

constexpr double LowPtGsfElectronSeedProducer::boundary_ = 2.50746495928 * 2.50746495928
staticconstexprprivate

Definition at line 163 of file LowPtGsfElectronSeedProducer.cc.

Referenced by propagateTrackToCalo().

◆ builderToken_

const edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> LowPtGsfElectronSeedProducer::builderToken_
private

Definition at line 154 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ ebRecHits_

const edm::EDGetTokenT<EcalRecHitCollection> LowPtGsfElectronSeedProducer::ebRecHits_
private

Definition at line 147 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ ecalClusters_

const edm::EDGetTokenT<reco::PFClusterCollection> LowPtGsfElectronSeedProducer::ecalClusters_
private

Definition at line 145 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ ecalClusterToolsESGetTokens_

const noZS::EcalClusterLazyTools::ESGetTokens LowPtGsfElectronSeedProducer::ecalClusterToolsESGetTokens_
private

Definition at line 155 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ eeRecHits_

const edm::EDGetTokenT<EcalRecHitCollection> LowPtGsfElectronSeedProducer::eeRecHits_
private

Definition at line 148 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ field_

edm::ESHandle<MagneticField> LowPtGsfElectronSeedProducer::field_
private

Definition at line 140 of file LowPtGsfElectronSeedProducer.cc.

Referenced by beginLuminosityBlock(), lightGsfTracking(), and propagateTrackToCalo().

◆ fitterPtr_

std::unique_ptr<TrajectoryFitter> LowPtGsfElectronSeedProducer::fitterPtr_
private

Definition at line 141 of file LowPtGsfElectronSeedProducer.cc.

Referenced by lightGsfTracking(), and loop().

◆ hcalClusters_

edm::EDGetTokenT<reco::PFClusterCollection> LowPtGsfElectronSeedProducer::hcalClusters_
private

Definition at line 146 of file LowPtGsfElectronSeedProducer.cc.

Referenced by produce().

◆ kfTracks_

edm::EDGetTokenT<reco::TrackCollection> LowPtGsfElectronSeedProducer::kfTracks_
private

Definition at line 143 of file LowPtGsfElectronSeedProducer.cc.

Referenced by produce().

◆ mass_

constexpr double LowPtGsfElectronSeedProducer::mass_ = 0.000511 * 0.000511
staticconstexprprivate

Definition at line 165 of file LowPtGsfElectronSeedProducer.cc.

Referenced by propagateTrackToCalo().

◆ maxPtThreshold_

const double LowPtGsfElectronSeedProducer::maxPtThreshold_
private

Definition at line 160 of file LowPtGsfElectronSeedProducer.cc.

Referenced by decision().

◆ minPtThreshold_

const double LowPtGsfElectronSeedProducer::minPtThreshold_
private

Definition at line 159 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ passThrough_

const bool LowPtGsfElectronSeedProducer::passThrough_
private

Definition at line 157 of file LowPtGsfElectronSeedProducer.cc.

Referenced by decision(), and loop().

◆ pfTracks_

edm::EDGetTokenT<reco::PFRecTrackCollection> LowPtGsfElectronSeedProducer::pfTracks_
private

Definition at line 144 of file LowPtGsfElectronSeedProducer.cc.

Referenced by produce().

◆ rho_

const edm::EDGetTokenT<double> LowPtGsfElectronSeedProducer::rho_
private

Definition at line 149 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ smootherPtr_

std::unique_ptr<TrajectorySmoother> LowPtGsfElectronSeedProducer::smootherPtr_
private

Definition at line 142 of file LowPtGsfElectronSeedProducer.cc.

Referenced by lightGsfTracking(), and loop().

◆ trajectoryFitterToken_

const edm::ESGetToken<TrajectoryFitter, TrajectoryFitter::Record> LowPtGsfElectronSeedProducer::trajectoryFitterToken_
private

Definition at line 152 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ trajectorySmootherToken_

const edm::ESGetToken<TrajectorySmoother, TrajectoryFitter::Record> LowPtGsfElectronSeedProducer::trajectorySmootherToken_
private

Definition at line 153 of file LowPtGsfElectronSeedProducer.cc.

Referenced by loop().

◆ usePfTracks_

const bool LowPtGsfElectronSeedProducer::usePfTracks_
private

Definition at line 158 of file LowPtGsfElectronSeedProducer.cc.

Referenced by produce().

Vector3DBase
Definition: Vector3DBase.h:8
muonTagProbeFilters_cff.matched
matched
Definition: muonTagProbeFilters_cff.py:62
edm::RefProd
Definition: EDProductfwd.h:25
LowPtGsfElectronSeedProducer::pfTracks_
edm::EDGetTokenT< reco::PFRecTrackCollection > pfTracks_
Definition: LowPtGsfElectronSeedProducer.cc:144
lowPtGsfElectronSeeds_cfi.pfTracks
pfTracks
Definition: lowPtGsfElectronSeeds_cfi.py:10
LowPtGsfElectronSeedProducer::rho_
const edm::EDGetTokenT< double > rho_
Definition: LowPtGsfElectronSeedProducer.cc:149
LowPtGsfElectronSeedProducer::usePfTracks_
const bool usePfTracks_
Definition: LowPtGsfElectronSeedProducer.cc:158
MagneticField::inTesla
virtual GlobalVector inTesla(const GlobalPoint &gp) const =0
Field value ad specified global point, in Tesla.
RawParticle
Definition: RawParticle.h:37
edm::ESInputTag
Definition: ESInputTag.h:87
BaseParticlePropagator::propagateToEcalEntrance
bool propagateToEcalEntrance(bool first=true)
Definition: BaseParticlePropagator.cc:450
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
lowPtGsfElectronSeeds_cfi.hcalClusters
hcalClusters
Definition: lowPtGsfElectronSeeds_cfi.py:12
Trajectory::chiSquared
float chiSquared() const
Definition: Trajectory.h:241
patZpeak.handle
handle
Definition: patZpeak.py:23
LowPtGsfElectronSeedProducer::eeRecHits_
const edm::EDGetTokenT< EcalRecHitCollection > eeRecHits_
Definition: LowPtGsfElectronSeedProducer.cc:148
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
pos
Definition: PixelAliasList.h:18
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89287
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:153
TrajectoryMeasurement::updatedState
TrajectoryStateOnSurface const & updatedState() const
Definition: TrajectoryMeasurement.h:184
XYZTLorentzVector
math::XYZTLorentzVector XYZTLorentzVector
Definition: RawParticle.h:25
LowPtGsfElectronSeedProducer::boundary_
static constexpr double boundary_
Definition: LowPtGsfElectronSeedProducer.cc:163
vertices_cff.ntracks
ntracks
Definition: vertices_cff.py:34
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
reco::PreId
Definition: PreId.h:13
edm::Ref::index
key_type index() const
Definition: Ref.h:253
edm::Handle
Definition: AssociativeIterator.h:50
singleTopDQM_cfi.setup
setup
Definition: singleTopDQM_cfi.py:37
BaseParticlePropagator
Definition: BaseParticlePropagator.h:82
edm::Ref
Definition: AssociativeIterator.h:58
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
reco::PFClusterRef
edm::Ref< PFClusterCollection > PFClusterRef
persistent reference to PFCluster objects
Definition: PFClusterFwd.h:15
fileCollector.seed
seed
Definition: fileCollector.py:127
IdealMagneticFieldRecord
Definition: IdealMagneticFieldRecord.h:11
TrajectoryStateOnSurface
Definition: TrajectoryStateOnSurface.h:16
reco::ElectronSeed
Definition: ElectronSeed.h:51
LowPtGsfElectronSeedProducer::maxPtThreshold_
const double maxPtThreshold_
Definition: LowPtGsfElectronSeedProducer.cc:160
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
LowPtGsfElectronSeedProducer::builderToken_
const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord > builderToken_
Definition: LowPtGsfElectronSeedProducer.cc:154
contentValuesCheck.values
values
Definition: contentValuesCheck.py:38
LowPtGsfElectronSeedProducer::getBaseRef
reco::TrackRef getBaseRef(edm::Handle< std::vector< reco::Track > > handle, int idx) const
Definition: LowPtGsfElectronSeedProducer.cc:277
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
reco::PFRecTrackRef
edm::Ref< PFRecTrackCollection > PFRecTrackRef
persistent reference to PFRecTrack objects
Definition: PFRecTrackFwd.h:15
HCALHighEnergyHPDFilter_cfi.energy
energy
Definition: HCALHighEnergyHPDFilter_cfi.py:5
GlobalTrajectoryParameters
Definition: GlobalTrajectoryParameters.h:15
Point3DBase< float, GlobalTag >
RawParticle::z
double z() const
z of vertex
Definition: RawParticle.h:284
LowPtGsfElectronSeedProducer::fitterPtr_
std::unique_ptr< TrajectoryFitter > fitterPtr_
Definition: LowPtGsfElectronSeedProducer.cc:141
reco::TrackRef
edm::Ref< TrackCollection > TrackRef
persistent reference to a Track
Definition: TrackFwd.h:20
DDAxes::rho
Trajectory::ConstRecHitContainer
TrackingRecHit::ConstRecHitContainer ConstRecHitContainer
Definition: Trajectory.h:41
bsc_activity_cfg.clusters
clusters
Definition: bsc_activity_cfg.py:36
InitialStep_cff.seeds
seeds
Definition: InitialStep_cff.py:230
runTauDisplay.gp
gp
Definition: runTauDisplay.py:431
EcalClusterLazyToolsT
Definition: EcalClusterLazyTools.h:135
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
RecoEcal_cff.ecalClusters
ecalClusters
Definition: RecoEcal_cff.py:26
Trajectory::lastMeasurement
TrajectoryMeasurement const & lastMeasurement() const
Definition: Trajectory.h:150
trigObjTnPSource_cfi.filler
filler
Definition: trigObjTnPSource_cfi.py:21
reco::deltaR2
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16
LowPtGsfElectronSeedProducer::field_
edm::ESHandle< MagneticField > field_
Definition: LowPtGsfElectronSeedProducer.cc:140
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:31
reco::PFCluster::getDepthCorrection
static double getDepthCorrection(double energy, bool isBelowPS=false, bool isHadron=false)
Definition: PFCluster.cc:39
LowPtGsfElectronSeedProducer::ecalClusters_
const edm::EDGetTokenT< reco::PFClusterCollection > ecalClusters_
Definition: LowPtGsfElectronSeedProducer.cc:145
LowPtGsfElectronSeedProducer::smootherPtr_
std::unique_ptr< TrajectorySmoother > smootherPtr_
Definition: LowPtGsfElectronSeedProducer.cc:142
LowPtGsfElectronSeedProducer::hcalClusters_
edm::EDGetTokenT< reco::PFClusterCollection > hcalClusters_
Definition: LowPtGsfElectronSeedProducer.cc:146
LowPtGsfElectronSeedProducer::decision
bool decision(const reco::PFRecTrackRef &pfTrackRef, reco::PreId &ecal, reco::PreId &hcal, double rho, const reco::BeamSpot &spot, noZS::EcalClusterLazyTools &ecalTools)
Definition: LowPtGsfElectronSeedProducer.cc:595
LowPtGsfElectronSeedProducer::propagateTrackToCalo
void propagateTrackToCalo(const reco::PFRecTrackRef &pfTrackRef, const edm::Handle< reco::PFClusterCollection > &ecalClusters, const edm::Handle< reco::PFClusterCollection > &hcalClusters, std::vector< int > &matchedEcalClusters, std::vector< int > &matchedHcalClusters, reco::PreId &ecalPreId, reco::PreId &hcalPreId)
Definition: LowPtGsfElectronSeedProducer.cc:403
get
#define get
LowPtGsfElectronSeedProducer::trajectorySmootherToken_
const edm::ESGetToken< TrajectorySmoother, TrajectoryFitter::Record > trajectorySmootherToken_
Definition: LowPtGsfElectronSeedProducer.cc:153
EcalClusterLazyToolsBase::ESGetTokens::get
ESData get(edm::EventSetup const &eventSetup) const
Definition: EcalClusterLazyTools.h:62
unit
Basic3DVector unit() const
Definition: Basic3DVectorLD.h:162
LowPtGsfElectronSeedProducer::lightGsfTracking
bool lightGsfTracking(reco::PreId &, const reco::TrackRef &, const reco::ElectronSeed &)
Definition: LowPtGsfElectronSeedProducer.cc:557
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
LowPtGsfElectronSeedProducer::beamSpot_
const edm::EDGetTokenT< reco::BeamSpot > beamSpot_
Definition: LowPtGsfElectronSeedProducer.cc:150
LowPtGsfElectronSeedProducer::passThrough_
const bool passThrough_
Definition: LowPtGsfElectronSeedProducer.cc:157
TrajectoryStateOnSurface::globalMomentum
GlobalVector globalMomentum() const
Definition: TrajectoryStateOnSurface.h:66
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
eostools.move
def move(src, dest)
Definition: eostools.py:511
TkClonerImpl
Definition: TkClonerImpl.h:12
Trajectory::firstMeasurement
TrajectoryMeasurement const & firstMeasurement() const
Definition: Trajectory.h:166
LowPtGsfElectronSeedProducer::TrackIndxMap
std::unordered_map< reco::TrackRef::key_type, size_t > TrackIndxMap
Definition: LowPtGsfElectronSeedProducer.cc:60
Trajectory
Definition: Trajectory.h:38
LowPtGsfElectronSeedProducer::trajectoryFitterToken_
const edm::ESGetToken< TrajectoryFitter, TrajectoryFitter::Record > trajectoryFitterToken_
Definition: LowPtGsfElectronSeedProducer.cc:152
LowPtGsfElectronSeedProducer::ecalClusterToolsESGetTokens_
const noZS::EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
Definition: LowPtGsfElectronSeedProducer.cc:155
lowptgsfeleseed::HeavyObjectCache
Definition: LowPtGsfElectronSeedHeavyObjectCache.h:16
reco::PFTrajectoryPoint
A PFTrack holds several trajectory points, which basically contain the position and momentum of a tra...
Definition: PFTrajectoryPoint.h:26
Exception
Definition: hltDiff.cc:246
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
bsc_activity_cfg.ecal
ecal
Definition: bsc_activity_cfg.py:25
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
mps_fire.result
result
Definition: mps_fire.py:311
edm::helper::Filler
Definition: ValueMap.h:22
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
LowPtGsfElectronSeedProducer::fillPreIdRefValueMap
void fillPreIdRefValueMap(edm::Handle< CollType > tracksHandle, const TrackIndxMap &trksToPreIdIndx, const edm::OrphanHandle< reco::PreIdCollection > &preIdHandle, edm::ValueMap< reco::PreIdRef >::Filler &filler)
Definition: LowPtGsfElectronSeedProducer.cc:621
reco::PreId::setTrackProperties
void setTrackProperties(float newchi2, float chi2ratio, float dpt)
Definition: PreId.h:53
point
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
LowPtGsfElectronSeedProducer::loop
void loop(const edm::Handle< std::vector< T > > &handle, edm::Handle< reco::PFClusterCollection > &hcalClusters, reco::ElectronSeedCollection &seeds, reco::PreIdCollection &ecalPreIds, reco::PreIdCollection &hcalPreIds, TrackIndxMap &trksToPreIdIndx, edm::Event &, const edm::EventSetup &)
Definition: LowPtGsfElectronSeedProducer.cc:289
event
Definition: event.py:1
LowPtGsfElectronSeedProducer::mass_
static constexpr double mass_
Definition: LowPtGsfElectronSeedProducer.cc:165
edm::HandleBase::id
ProductID id() const
Definition: HandleBase.cc:29
reco::PreId::setECALMatchingProperties
void setECALMatchingProperties(PFClusterRef clusterRef, const math::XYZPoint &ecalpos, const math::XYZPoint &meanShower, float deta, float dphi, float chieta, float chiphi, float chi2, float eop)
Definition: PreId.h:33
LowPtGsfElectronSeedProducer::kfTracks_
edm::EDGetTokenT< reco::TrackCollection > kfTracks_
Definition: LowPtGsfElectronSeedProducer.cc:143
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
SiStripBadComponentsDQMServiceTemplate_cfg.ep
ep
Definition: SiStripBadComponentsDQMServiceTemplate_cfg.py:86
LowPtGsfElectronSeedProducer::minPtThreshold_
const double minPtThreshold_
Definition: LowPtGsfElectronSeedProducer.cc:159
edm::InputTag
Definition: InputTag.h:15
reco::PreIdCollection
std::vector< reco::PreId > PreIdCollection
Definition: PreIdFwd.h:6
Trajectory::isValid
bool isValid() const
Definition: Trajectory.h:257
LowPtGsfElectronSeedProducer::ebRecHits_
const edm::EDGetTokenT< EcalRecHitCollection > ebRecHits_
Definition: LowPtGsfElectronSeedProducer.cc:147
Generated for CMSSW Reference Manual by doxygen 1.8.16