#include <TrajSeedMatcher.h>

Classes
struct	Configuration

struct	MatchInfo

class	MatchingCuts

class	MatchingCutsV1

class	MatchingCutsV2

struct	SCHitMatch

struct	SeedWithInfo

Public Member Functions
std::vector< TrajSeedMatcher::SeedWithInfo >	operator() (const GlobalPoint &candPos, const float energy)

	TrajSeedMatcher (TrajectorySeedCollection const &seeds, math::XYZPoint const &vprim, Configuration const &cfg, edm::EventSetup const &iSetup, MeasurementTrackerEvent const &measTkEvt)

	~TrajSeedMatcher ()=default

Static Public Member Functions
static edm::ParameterSetDescription	makePSetDescription ()

Private Member Functions
void	clearCache ()

auto	ftsFromVertexToPoint (GlobalPoint const &point, GlobalPoint const &vertex, float energy, int charge) const

size_t	getNrHitsRequired (const int nrValidLayers) const

int	getNrValidLayersAlongTraj (const SCHitMatch &hit1, const SCHitMatch &hit2, const GlobalPoint &candPos, const float energy, const int charge)

int	getNrValidLayersAlongTraj (const DetId &hitId, const TrajectoryStateOnSurface &hitTrajState) const

const TrajectoryStateOnSurface &	getTrajStateFromPoint (const TrackingRecHit &hit, const FreeTrajectoryState &initialState, const GlobalPoint &point, const PropagatorWithMaterial &propagator)

const TrajectoryStateOnSurface &	getTrajStateFromVtx (const TrackingRecHit &hit, const TrajectoryStateOnSurface &initialState, const PropagatorWithMaterial &propagator)

bool	layerHasValidHits (const DetLayer &layer, const TrajectoryStateOnSurface &hitSurState, const Propagator &propToLayerFromState) const

TrajectoryStateOnSurface	makeTrajStateOnSurface (const GlobalPoint &pos, const float energy, const int charge) const

std::vector< SCHitMatch >	processSeed (const TrajectorySeed &seed, const GlobalPoint &candPos, const float energy, const TrajectoryStateOnSurface &initialTrajState)

Static Private Member Functions
static float	getZVtxFromExtrapolation (const GlobalPoint &primeVtxPos, const GlobalPoint &hitPos, const GlobalPoint &candPos)

Private Attributes
PropagatorWithMaterial	backwardPropagator_

Configuration const &	cfg_

DetLayerGeometry const &	detLayerGeom_

PropagatorWithMaterial	forwardPropagator_

MagneticField const &	magField_

MagneticField const &	magFieldParam_

MeasurementTrackerEvent const &	measTkEvt_

NavigationSchool const &	navSchool_

TrajectorySeedCollection const &	seeds_

IntGlobalPointPairUnorderedMap< TrajectoryStateOnSurface >	trajStateFromPointNegChargeCache_

IntGlobalPointPairUnorderedMap< TrajectoryStateOnSurface >	trajStateFromPointPosChargeCache_

std::unordered_map< int, TrajectoryStateOnSurface >	trajStateFromVtxNegChargeCache_

std::unordered_map< int, TrajectoryStateOnSurface >	trajStateFromVtxPosChargeCache_

const GlobalPoint	vprim_

Static Private Attributes
static constexpr float	kElectronMass_ = 0.000511

Detailed Description

Definition at line 54 of file TrajSeedMatcher.h.

Constructor & Destructor Documentation

◆ TrajSeedMatcher()

TrajSeedMatcher::TrajSeedMatcher	(	TrajectorySeedCollection const &	seeds,
		math::XYZPoint const &	vprim,
		TrajSeedMatcher::Configuration const &	cfg,
		edm::EventSetup const &	iSetup,
		MeasurementTrackerEvent const &	measTkEvt
	)

explicit

Definition at line 102 of file TrajSeedMatcher.cc.

References HLT_FULL_cff::seeds.

     : seeds_{seeds},
       vprim_(vprim.x(), vprim.y(), vprim.z()),
       cfg_{cfg},
       magField_{iSetup.getData(cfg_.magFieldToken)},
       magFieldParam_{iSetup.getData(cfg_.paramMagFieldToken)},
       measTkEvt_{measTkEvt},
       navSchool_{iSetup.getData(cfg_.navSchoolToken)},
       detLayerGeom_{iSetup.getData(cfg_.detLayerGeomToken)},
       forwardPropagator_(alongMomentum, kElectronMass_, &magField_),
       backwardPropagator_(oppositeToMomentum, kElectronMass_, &magField_) {}

◆ ~TrajSeedMatcher()

TrajSeedMatcher::~TrajSeedMatcher ( )

default

Member Function Documentation

◆ clearCache()

void TrajSeedMatcher::clearCache ( )

private

Definition at line 299 of file TrajSeedMatcher.cc.

References trajStateFromPointNegChargeCache_, trajStateFromPointPosChargeCache_, trajStateFromVtxNegChargeCache_, and trajStateFromVtxPosChargeCache_.

Referenced by operator()().

                                  {
   trajStateFromVtxPosChargeCache_.clear();
   trajStateFromVtxNegChargeCache_.clear();
   trajStateFromPointPosChargeCache_.clear();
   trajStateFromPointNegChargeCache_.clear();
 }

◆ ftsFromVertexToPoint()

auto TrajSeedMatcher::ftsFromVertexToPoint	(	GlobalPoint const &	point,
		GlobalPoint const &	vertex,
		float	energy,
		int	charge
	)		const

inlineprivate

Definition at line 192 of file TrajSeedMatcher.h.

References cfg_, ALCARECOTkAlJpsiMuMu_cff::charge, hcalRecHitTable_cff::energy, trackingTools::ftsFromVertexToPoint(), MagneticField::isDefined(), magField_, magFieldParam_, HLT_2023v12_cff::magneticField, point, TrajSeedMatcher::Configuration::useParamMagFieldIfDefined, and bphysicsOniaDQM_cfi::vertex.

Referenced by getNrValidLayersAlongTraj(), makeTrajStateOnSurface(), and processSeed().

                                                                                                                         {
     //parameterised b-fields may not be valid for entire detector, just tracker volume
     //however need we ecal so we auto select based on the position
     bool useMagFieldParam = cfg_.useParamMagFieldIfDefined && magFieldParam_.isDefined(point);
     auto const& magneticField = useMagFieldParam ? magFieldParam_ : magField_;
     return trackingTools::ftsFromVertexToPoint(magneticField, point, vertex, energy, charge);
   }

◆ getNrHitsRequired()

size_t TrajSeedMatcher::getNrHitsRequired ( const int nrValidLayers ) const

private

Definition at line 364 of file TrajSeedMatcher.cc.

References cfg_, TrajSeedMatcher::Configuration::minNrHits, and TrajSeedMatcher::Configuration::minNrHitsValidLayerBins.

Referenced by operator()().

                                                                        {
   for (size_t binNr = 0; binNr < cfg_.minNrHitsValidLayerBins.size(); binNr++) {
     if (nrValidLayers < cfg_.minNrHitsValidLayerBins[binNr])
       return cfg_.minNrHits[binNr];
   }
   return cfg_.minNrHits.back();
 }

◆ getNrValidLayersAlongTraj() [1/2]

int TrajSeedMatcher::getNrValidLayersAlongTraj	(	const SCHitMatch &	hit1,
		const SCHitMatch &	hit2,
		const GlobalPoint &	candPos,
		const float	energy,
		const int	charge
	)

private

Definition at line 306 of file TrajSeedMatcher.cc.

References cfg_, ALCARECOTkAlJpsiMuMu_cff::charge, hcalRecHitTable_cff::energy, forwardPropagator_, ftsFromVertexToPoint(), TrackingRecHit::geographicalId(), getTrajStateFromPoint(), getZVtxFromExtrapolation(), TrajSeedMatcher::SCHitMatch::hit, TrajSeedMatcher::SCHitMatch::hitPos, TrajSeedMatcher::Configuration::useRecoVertex, bphysicsOniaDQM_cfi::vertex, vprim_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and HLT_FULL_cff::zVertex.

Referenced by operator()().

                                                                                                                       {
   double zVertex = cfg_.useRecoVertex ? vprim_.z() : getZVtxFromExtrapolation(vprim_, hit1.hitPos, candPos);
   GlobalPoint vertex(vprim_.x(), vprim_.y(), zVertex);
 
   auto firstMatchFreeTraj = ftsFromVertexToPoint(hit1.hitPos, vertex, energy, charge);
   auto const& secondHitTraj = getTrajStateFromPoint(hit2.hit, firstMatchFreeTraj, hit1.hitPos, forwardPropagator_);
   return getNrValidLayersAlongTraj(hit2.hit.geographicalId(), secondHitTraj);
 }

◆ getNrValidLayersAlongTraj() [2/2]

int TrajSeedMatcher::getNrValidLayersAlongTraj	(	const DetId &	hitId,
		const TrajectoryStateOnSurface &	hitTrajState
	)		const

private

Definition at line 316 of file TrajSeedMatcher.cc.

References alongMomentum, backwardPropagator_, NavigationSchool::compatibleLayers(), detLayerGeom_, forwardPropagator_, TrajectoryStateOnSurface::freeState(), DetLayerGeometry::idToLayer(), GeomDetEnumerators::isTrackerPixel(), nano_mu_digi_cff::layer, layerHasValidHits(), navSchool_, and oppositeToMomentum.

                                                                                                                      {
   const DetLayer* detLayer = detLayerGeom_.idToLayer(hitId);
   if (detLayer == nullptr)
     return 0;
 
   const FreeTrajectoryState& hitFreeState = *hitTrajState.freeState();
   auto const inLayers = navSchool_.compatibleLayers(*detLayer, hitFreeState, oppositeToMomentum);
   const auto outLayers = navSchool_.compatibleLayers(*detLayer, hitFreeState, alongMomentum);
 
   int nrValidLayers = 1;  //because our current hit is also valid and wont be included in the count otherwise
   for (auto layer : inLayers) {
     if (GeomDetEnumerators::isTrackerPixel(layer->subDetector())) {
       if (layerHasValidHits(*layer, hitTrajState, backwardPropagator_))
         nrValidLayers++;
     }
   }
   for (auto layer : outLayers) {
     if (GeomDetEnumerators::isTrackerPixel(layer->subDetector())) {
       if (layerHasValidHits(*layer, hitTrajState, forwardPropagator_))
         nrValidLayers++;
     }
   }
   return nrValidLayers;
 }

◆ getTrajStateFromPoint()

const TrajectoryStateOnSurface & TrajSeedMatcher::getTrajStateFromPoint	(	const TrackingRecHit &	hit,
		const FreeTrajectoryState &	initialState,
		const GlobalPoint &	point,
		const PropagatorWithMaterial &	propagator
	)

private

Definition at line 274 of file TrajSeedMatcher.cc.

References FreeTrajectoryState::charge(), submitPVResolutionJobs::key, point, TrackCandidateProducer_cfi::propagator, trajStateFromPointNegChargeCache_, trajStateFromPointPosChargeCache_, and heppy_batch::val.

Referenced by getNrValidLayersAlongTraj(), and processSeed().

                                                                                                                  {
   auto& trajStateFromPointCache =
       initialState.charge() == 1 ? trajStateFromPointPosChargeCache_ : trajStateFromPointNegChargeCache_;
 
   auto key = std::make_pair(hit.det()->gdetIndex(), point);
   auto res = trajStateFromPointCache.find(key);
   if (res != trajStateFromPointCache.end())
     return res->second;
   else {  //doesnt exist, need to make it
     //FIXME: check for efficiency
     auto val = trajStateFromPointCache.emplace(key, propagator.propagate(initialState, hit.det()->surface()));
     return val.first->second;
   }
 }

◆ getTrajStateFromVtx()

const TrajectoryStateOnSurface & TrajSeedMatcher::getTrajStateFromVtx	(	const TrackingRecHit &	hit,
		const TrajectoryStateOnSurface &	initialState,
		const PropagatorWithMaterial &	propagator
	)

private

Definition at line 257 of file TrajSeedMatcher.cc.

References TrajectoryStateOnSurface::charge(), submitPVResolutionJobs::key, TrackCandidateProducer_cfi::propagator, trajStateFromVtxNegChargeCache_, trajStateFromVtxPosChargeCache_, and heppy_batch::val.

Referenced by processSeed().

                                                                                                                {
   auto& trajStateFromVtxCache =
       initialState.charge() == 1 ? trajStateFromVtxPosChargeCache_ : trajStateFromVtxNegChargeCache_;
 
   auto key = hit.det()->gdetIndex();
   auto res = trajStateFromVtxCache.find(key);
   if (res != trajStateFromVtxCache.end())
     return res->second;
   else {  //doesnt exist, need to make it
     //FIXME: check for efficiency
     auto val = trajStateFromVtxCache.emplace(key, propagator.propagate(initialState, hit.det()->surface()));
     return val.first->second;
   }
 }

◆ getZVtxFromExtrapolation()

float TrajSeedMatcher::getZVtxFromExtrapolation	(	const GlobalPoint &	primeVtxPos,
		const GlobalPoint &	hitPos,
		const GlobalPoint &	candPos
	)

staticprivate

Definition at line 245 of file TrajSeedMatcher.cc.

References LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, mathSSE::sqrt(), x, and PV3DBase< T, PVType, FrameType >::z().

Referenced by getNrValidLayersAlongTraj(), and processSeed().

                                                                             {
   auto sq = [](float x) { return x * x; };
   auto calRDiff = [sq](const GlobalPoint& p1, const GlobalPoint& p2) {
     return std::sqrt(sq(p2.x() - p1.x()) + sq(p2.y() - p1.y()));
   };
   const double r1Diff = calRDiff(primeVtxPos, hitPos);
   const double r2Diff = calRDiff(hitPos, candPos);
   return hitPos.z() - r1Diff * (candPos.z() - hitPos.z()) / r2Diff;
 }

◆ layerHasValidHits()

bool TrajSeedMatcher::layerHasValidHits	(	const DetLayer &	layer,
		const TrajectoryStateOnSurface &	hitSurState,
		const Propagator &	propToLayerFromState
	)		const

private

Definition at line 341 of file TrajSeedMatcher.cc.

References MeasurementDetWithData::hasBadComponents(), MeasurementTrackerEvent::idToDet(), MeasurementDetWithData::isActive(), nano_mu_digi_cff::layer, and measTkEvt_.

Referenced by getNrValidLayersAlongTraj().

                                                                                       {
   //FIXME: do not hardcode with weird magic numbers stolen from ancient tracking code
   //its taken from https://cmssdt.cern.ch/dxr/CMSSW/source/RecoTracker/TrackProducer/interface/TrackProducerBase.icc#165
   //which inspires this code
   Chi2MeasurementEstimator estimator(30., -3.0, 0.5, 2.0, 0.5, 1.e12);  // same as defauts....
 
   const std::vector<GeometricSearchDet::DetWithState>& detWithState =
       layer.compatibleDets(hitSurState, propToLayerFromState, estimator);
   if (detWithState.empty())
     return false;
   else {
     DetId id = detWithState.front().first->geographicalId();
     MeasurementDetWithData measDet = measTkEvt_.idToDet(id);
     //Below, measDet.hasBadComponents handles the check that a Pixel module has or not errors like FED25
     if (measDet.isActive() && !measDet.hasBadComponents(detWithState.front().second))
       return true;
     else
       return false;
   }
 }

◆ makePSetDescription()

edm::ParameterSetDescription TrajSeedMatcher::makePSetDescription ( )

static

Definition at line 118 of file TrajSeedMatcher.cc.

References common_cff::defaults, submitPVResolutionJobs::desc, eostools::move(), and or.

Referenced by ElectronNHitSeedProducer::fillDescriptions().

                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<bool>("useRecoVertex", false);
   desc.add<bool>("enableHitSkipping", false);
   desc.add<bool>("requireExactMatchCount", true);
   desc.add<bool>("useParamMagFieldIfDefined", true);
   desc.add<edm::ESInputTag>("paramMagField", edm::ESInputTag{"", "ParabolicMf"});
   desc.add<edm::ESInputTag>("navSchool", edm::ESInputTag{"", "SimpleNavigationSchool"});
   desc.add<edm::ESInputTag>("detLayerGeom", edm::ESInputTag{"", "hltESPGlobalDetLayerGeometry"});
   desc.add<std::vector<int> >("minNrHitsValidLayerBins", {4});
   desc.add<std::vector<unsigned int> >("minNrHits", {2, 3});
 
   edm::ParameterSetDescription cutsDesc;
   auto cutDescCases = 1 >> (edm::ParameterDescription<double>("dPhiMax", 0.04, true) and
                             edm::ParameterDescription<double>("dRZMax", 0.09, true) and
                             edm::ParameterDescription<double>("dRZMaxLowEtThres", 20., true) and
                             edm::ParameterDescription<std::vector<double> >("dRZMaxLowEtEtaBins", {1., 1.5}, true) and
                             edm::ParameterDescription<std::vector<double> >("dRZMaxLowEt", {0.09, 0.15, 0.09}, true)) or
                       2 >> (edm::ParameterDescription<std::vector<double> >("dPhiMaxHighEt", {0.003}, true) and
                             edm::ParameterDescription<std::vector<double> >("dPhiMaxHighEtThres", {0.0}, true) and
                             edm::ParameterDescription<std::vector<double> >("dPhiMaxLowEtGrad", {0.0}, true) and
                             edm::ParameterDescription<std::vector<double> >("dRZMaxHighEt", {0.005}, true) and
                             edm::ParameterDescription<std::vector<double> >("dRZMaxHighEtThres", {30}, true) and
                             edm::ParameterDescription<std::vector<double> >("dRZMaxLowEtGrad", {-0.002}, true) and
                             edm::ParameterDescription<std::vector<double> >("etaBins", {}, true));
   cutsDesc.ifValue(edm::ParameterDescription<int>("version", 1, true), std::move(cutDescCases));
 
   edm::ParameterSet defaults;
   defaults.addParameter<double>("dPhiMax", 0.04);
   defaults.addParameter<double>("dRZMax", 0.09);
   defaults.addParameter<double>("dRZMaxLowEtThres", 0.09);
   defaults.addParameter<std::vector<double> >("dRZMaxLowEtEtaBins", std::vector<double>{1., 1.5});
   defaults.addParameter<std::vector<double> >("dRZMaxLowEt", std::vector<double>{0.09, 0.09, 0.09});
   defaults.addParameter<int>("version", 1);
   desc.addVPSet("matchingCuts", cutsDesc, std::vector<edm::ParameterSet>{defaults, defaults, defaults});
   return desc;
 }

◆ makeTrajStateOnSurface()

TrajectoryStateOnSurface TrajSeedMatcher::makeTrajStateOnSurface	(	const GlobalPoint &	pos,
		const float	energy,
		const int	charge
	)		const

private

Definition at line 292 of file TrajSeedMatcher.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, hcalRecHitTable_cff::energy, ftsFromVertexToPoint(), and vprim_.

Referenced by operator()().

                                                                                          {
   auto freeTS = ftsFromVertexToPoint(pos, vprim_, energy, charge);
   return TrajectoryStateOnSurface(freeTS, *PerpendicularBoundPlaneBuilder{}(freeTS.position(), freeTS.momentum()));
 }

◆ operator()()

std::vector< TrajSeedMatcher::SeedWithInfo > TrajSeedMatcher::operator()	(	const GlobalPoint &	candPos,
		const float	energy
	)

Definition at line 156 of file TrajSeedMatcher.cc.

References cfg_, clearCache(), hcalRecHitTable_cff::energy, getNrHitsRequired(), getNrValidLayersAlongTraj(), makeTrajStateOnSurface(), SiStripPI::max, processSeed(), TrajSeedMatcher::Configuration::requireExactMatchCount, fileCollector::seed, and seeds_.

                                                                                                                  {
   clearCache();
 
   std::vector<SeedWithInfo> matchedSeeds;
 
   //these are super expensive functions
   TrajectoryStateOnSurface scTrajStateOnSurfNeg = makeTrajStateOnSurface(candPos, energy, -1);
   TrajectoryStateOnSurface scTrajStateOnSurfPos = makeTrajStateOnSurface(candPos, energy, 1);
 
   for (const auto& seed : seeds_) {
     std::vector<SCHitMatch> matchedHitsNeg = processSeed(seed, candPos, energy, scTrajStateOnSurfNeg);
     std::vector<SCHitMatch> matchedHitsPos = processSeed(seed, candPos, energy, scTrajStateOnSurfPos);
 
     int nrValidLayersPos = 0;
     int nrValidLayersNeg = 0;
     if (matchedHitsNeg.size() >= 2) {
       nrValidLayersNeg = getNrValidLayersAlongTraj(matchedHitsNeg[0], matchedHitsNeg[1], candPos, energy, -1);
     }
     if (matchedHitsPos.size() >= 2) {
       nrValidLayersPos = getNrValidLayersAlongTraj(matchedHitsPos[0], matchedHitsPos[1], candPos, energy, +1);
     }
 
     int nrValidLayers = std::max(nrValidLayersNeg, nrValidLayersPos);
     size_t nrHitsRequired = getNrHitsRequired(nrValidLayers);
     bool matchCountPasses;
     if (cfg_.requireExactMatchCount) {
       // If the input seed collection is not cross-cleaned, an exact match is necessary to
       // prevent redundant seeds.
       matchCountPasses = matchedHitsNeg.size() == nrHitsRequired || matchedHitsPos.size() == nrHitsRequired;
     } else {
       matchCountPasses = matchedHitsNeg.size() >= nrHitsRequired || matchedHitsPos.size() >= nrHitsRequired;
     }
     if (matchCountPasses) {
       matchedSeeds.push_back({seed, makeMatchInfoVector(matchedHitsPos, matchedHitsNeg), nrValidLayers});
     }
   }
   return matchedSeeds;
 }

◆ processSeed()

std::vector< TrajSeedMatcher::SCHitMatch > TrajSeedMatcher::processSeed	(	const TrajectorySeed &	seed,
		const GlobalPoint &	candPos,
		const float	energy,
		const TrajectoryStateOnSurface &	initialTrajState
	)

private

Definition at line 195 of file TrajSeedMatcher.cc.

References backwardPropagator_, cfg_, ALCARECOTkAlJpsiMuMu_cff::charge, TrajectoryStateOnSurface::charge(), TrajSeedMatcher::Configuration::enableHitSkipping, hcalRecHitTable_cff::energy, PV3DBase< T, PVType, FrameType >::eta(), forwardPropagator_, ftsFromVertexToPoint(), getTrajStateFromPoint(), getTrajStateFromVtx(), getZVtxFromExtrapolation(), TrajectoryStateOnSurface::isValid(), match(), oniaPATMuonsWithTrigger_cff::matches, TrajSeedMatcher::Configuration::matchingCuts, PV3DBase< T, PVType, FrameType >::phi(), rpcPointValidation_cfi::recHit, fileCollector::seed, funct::sin(), PV3DBase< T, PVType, FrameType >::theta(), TrajSeedMatcher::Configuration::useRecoVertex, bphysicsOniaDQM_cfi::vertex, vprim_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and HLT_FULL_cff::zVertex.

Referenced by operator()().

                                                                                                                         {
   //next try passing these variables in once...
   const float candEta = candPos.eta();
   const float candEt = energy * std::sin(candPos.theta());
   const int charge = initialTrajState.charge();
 
   std::vector<SCHitMatch> matches;
   FreeTrajectoryState firstMatchFreeTraj;
   GlobalPoint prevHitPos;
   GlobalPoint vertex;
   const auto nCuts = cfg_.matchingCuts.size();
   for (size_t iHit = 0;
        matches.size() < nCuts && iHit < seed.nHits() && (cfg_.enableHitSkipping || iHit == matches.size());
        iHit++) {
     auto const& recHit = *(seed.recHits().begin() + iHit);
 
     if (!recHit.isValid()) {
       continue;
     }
 
     const bool doFirstMatch = matches.empty();
 
     auto const& trajState = doFirstMatch
                                 ? getTrajStateFromVtx(recHit, initialTrajState, backwardPropagator_)
                                 : getTrajStateFromPoint(recHit, firstMatchFreeTraj, prevHitPos, forwardPropagator_);
     if (!trajState.isValid()) {
       continue;
     }
 
     auto const& vtxForMatchObject = doFirstMatch ? vprim_ : vertex;
     auto match = makeSCHitMatch(vtxForMatchObject, trajState, recHit, candEt, candEta, candPos.phi(), charge, 1);
 
     if ((*cfg_.matchingCuts[matches.size()])(match)) {
       matches.push_back(match);
       if (doFirstMatch) {
         //now we can figure out the z vertex
         double zVertex = cfg_.useRecoVertex ? vprim_.z() : getZVtxFromExtrapolation(vprim_, match.hitPos, candPos);
         vertex = GlobalPoint(vprim_.x(), vprim_.y(), zVertex);
         firstMatchFreeTraj = ftsFromVertexToPoint(match.hitPos, vertex, energy, charge);
       }
       prevHitPos = match.hitPos;
     }
   }
   return matches;
 }

Member Data Documentation

◆ backwardPropagator_

PropagatorWithMaterial TrajSeedMatcher::backwardPropagator_

private

Definition at line 215 of file TrajSeedMatcher.h.

Referenced by getNrValidLayersAlongTraj(), and processSeed().

◆ cfg_

Configuration const& TrajSeedMatcher::cfg_

private

Definition at line 206 of file TrajSeedMatcher.h.

Referenced by ftsFromVertexToPoint(), getNrHitsRequired(), getNrValidLayersAlongTraj(), operator()(), and processSeed().

◆ detLayerGeom_

DetLayerGeometry const& TrajSeedMatcher::detLayerGeom_

private

Definition at line 212 of file TrajSeedMatcher.h.

Referenced by getNrValidLayersAlongTraj().

◆ forwardPropagator_

PropagatorWithMaterial TrajSeedMatcher::forwardPropagator_

private

Definition at line 214 of file TrajSeedMatcher.h.

Referenced by getNrValidLayersAlongTraj(), and processSeed().

◆ kElectronMass_

constexpr float TrajSeedMatcher::kElectronMass_ = 0.000511

staticprivate

Definition at line 201 of file TrajSeedMatcher.h.

◆ magField_

MagneticField const& TrajSeedMatcher::magField_

private

Definition at line 208 of file TrajSeedMatcher.h.

Referenced by ftsFromVertexToPoint().

◆ magFieldParam_

MagneticField const& TrajSeedMatcher::magFieldParam_

private

Definition at line 209 of file TrajSeedMatcher.h.

Referenced by ftsFromVertexToPoint().

◆ measTkEvt_

MeasurementTrackerEvent const& TrajSeedMatcher::measTkEvt_

private

Definition at line 210 of file TrajSeedMatcher.h.

Referenced by layerHasValidHits().

◆ navSchool_

NavigationSchool const& TrajSeedMatcher::navSchool_

private

Definition at line 211 of file TrajSeedMatcher.h.

Referenced by getNrValidLayersAlongTraj().

◆ seeds_

TrajectorySeedCollection const& TrajSeedMatcher::seeds_

private

Definition at line 203 of file TrajSeedMatcher.h.

Referenced by operator()().

◆ trajStateFromPointNegChargeCache_

IntGlobalPointPairUnorderedMap<TrajectoryStateOnSurface> TrajSeedMatcher::trajStateFromPointNegChargeCache_

private

Definition at line 221 of file TrajSeedMatcher.h.

Referenced by clearCache(), and getTrajStateFromPoint().

◆ trajStateFromPointPosChargeCache_

IntGlobalPointPairUnorderedMap<TrajectoryStateOnSurface> TrajSeedMatcher::trajStateFromPointPosChargeCache_

private

Definition at line 220 of file TrajSeedMatcher.h.

Referenced by clearCache(), and getTrajStateFromPoint().

◆ trajStateFromVtxNegChargeCache_

std::unordered_map<int, TrajectoryStateOnSurface> TrajSeedMatcher::trajStateFromVtxNegChargeCache_

private

Definition at line 218 of file TrajSeedMatcher.h.

Referenced by clearCache(), and getTrajStateFromVtx().

◆ trajStateFromVtxPosChargeCache_

std::unordered_map<int, TrajectoryStateOnSurface> TrajSeedMatcher::trajStateFromVtxPosChargeCache_

private

Definition at line 217 of file TrajSeedMatcher.h.

Referenced by clearCache(), and getTrajStateFromVtx().

◆ vprim_

const GlobalPoint TrajSeedMatcher::vprim_

private

Definition at line 204 of file TrajSeedMatcher.h.

Referenced by getNrValidLayersAlongTraj(), makeTrajStateOnSurface(), and processSeed().

Classes

Public Member Functions

Static Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TrajSeedMatcher()

◆ ~TrajSeedMatcher()

Member Function Documentation

◆ clearCache()

◆ ftsFromVertexToPoint()

◆ getNrHitsRequired()

◆ getNrValidLayersAlongTraj() [1/2]

◆ getNrValidLayersAlongTraj() [2/2]

◆ getTrajStateFromPoint()

◆ getTrajStateFromVtx()

◆ getZVtxFromExtrapolation()

◆ layerHasValidHits()

◆ makePSetDescription()

◆ makeTrajStateOnSurface()

◆ operator()()

◆ processSeed()

Member Data Documentation

◆ backwardPropagator_

◆ cfg_

◆ detLayerGeom_

◆ forwardPropagator_

◆ kElectronMass_

◆ magField_

◆ magFieldParam_

◆ measTkEvt_

◆ navSchool_

◆ seeds_

◆ trajStateFromPointNegChargeCache_

◆ trajStateFromPointPosChargeCache_

◆ trajStateFromVtxNegChargeCache_

◆ trajStateFromVtxPosChargeCache_

◆ vprim_