TrackToTrackingParticleAssociator that associates by hits a bit quicker than the normal TrackAssociatorByHitsImpl class. More...

#include <QuickTrackAssociatorByHitsImpl.h>

Inheritance diagram for QuickTrackAssociatorByHitsImpl:

Public Types
enum	SimToRecoDenomType { denomnone, denomsim, denomreco }

Public Member Functions
reco::RecoToSimCollection	associateRecoToSim (const edm::Handle< edm::View< reco::Track > > &trackCollectionHandle, const edm::Handle< TrackingParticleCollection > &trackingParticleCollectionHandle) const override

reco::RecoToSimCollection	associateRecoToSim (const edm::RefToBaseVector< reco::Track > &trackCollection, const edm::RefVector< TrackingParticleCollection > &trackingParticleCollection) const override
	Association Reco To Sim with Collections. More...

reco::RecoToSimCollectionSeed	associateRecoToSim (const edm::Handle< edm::View< TrajectorySeed > > &, const edm::Handle< TrackingParticleCollection > &) const override

reco::RecoToSimCollectionTCandidate	associateRecoToSim (const edm::Handle< TrackCandidateCollection > &, const edm::RefVector< TrackingParticleCollection > &) const override

reco::SimToRecoCollection	associateSimToReco (const edm::Handle< edm::View< reco::Track > > &trackCollectionHandle, const edm::Handle< TrackingParticleCollection > &trackingParticleCollectionHandle) const override

reco::SimToRecoCollection	associateSimToReco (const edm::RefToBaseVector< reco::Track > &trackCollection, const edm::RefVector< TrackingParticleCollection > &trackingParticleCollection) const override
	Association Sim To Reco with Collections. More...

reco::SimToRecoCollectionSeed	associateSimToReco (const edm::Handle< edm::View< TrajectorySeed > > &, const edm::Handle< TrackingParticleCollection > &) const override

reco::SimToRecoCollectionTCandidate	associateSimToReco (const edm::Handle< TrackCandidateCollection > &, const edm::RefVector< TrackingParticleCollection > &) const override

	QuickTrackAssociatorByHitsImpl (edm::EDProductGetter const &productGetter, std::unique_ptr< const TrackerHitAssociator > hitAssoc, const ClusterTPAssociation *clusterToTPMap, bool absoluteNumberOfHits, double qualitySimToReco, double puritySimToReco, double cutRecoToSim, double pixelHitWeight, bool threeHitTracksAreSpecial, SimToRecoDenomType simToRecoDenominator)

template<>
double	weightedNumberOfTrackClusters (const reco::Track &track, const TrackerHitAssociator &) const

template<>
double	weightedNumberOfTrackClusters (const reco::Track &track, const ClusterTPAssociation &) const

Public Member Functions inherited from reco::TrackToTrackingParticleAssociatorBaseImpl
virtual reco::RecoToSimCollection	associateRecoToSim (const edm::Handle< edm::View< reco::Track >> &tCH, const edm::Handle< TrackingParticleCollection > &tPCH) const

virtual reco::RecoToSimCollectionSeed	associateRecoToSim (const edm::Handle< edm::View< TrajectorySeed >> &, const edm::Handle< TrackingParticleCollection > &) const

virtual reco::SimToRecoCollection	associateSimToReco (const edm::Handle< edm::View< reco::Track >> &tCH, const edm::Handle< TrackingParticleCollection > &tPCH) const

virtual reco::SimToRecoCollectionSeed	associateSimToReco (const edm::Handle< edm::View< TrajectorySeed >> &, const edm::Handle< TrackingParticleCollection > &) const

	TrackToTrackingParticleAssociatorBaseImpl ()
	Constructor. More...

virtual	~TrackToTrackingParticleAssociatorBaseImpl ()
	Destructor. More...

Private Types
typedef std::pair< uint32_t, EncodedEventId >	SimTrackIdentifiers
	This is enough information to uniquely identify a sim track. More...

typedef edm::IndexSet	TrackingParticleRefKeySet
	Set for TrackingParticleRef keys. More...

Private Member Functions
template<class T_RecoToSimCollection , class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >
T_RecoToSimCollection	associateRecoToSimImplementation (const T_TrackCollection &trackCollection, const T_TrackingParticleCollection &trackingParticleCollection, const TrackingParticleRefKeySet *trackingParticleKeys, T_hitOrClusterAssociator hitOrClusterAssociator) const
	The method that does the work for both overloads of associateRecoToSim. More...

template<class T_SimToRecoCollection , class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >
T_SimToRecoCollection	associateSimToRecoImplementation (const T_TrackCollection &trackCollection, const T_TrackingParticleCollection &trackingParticleCollection, const TrackingParticleRefKeySet *trackingParticleKeys, T_hitOrClusterAssociator hitOrClusterAssociator) const
	The method that does the work for both overloads of associateSimToReco. More...

template<typename T_TPCollection , typename iter >
std::vector< std::pair< edm::Ref< TrackingParticleCollection >, double > >	associateTrack (const TrackerHitAssociator &hitAssociator, const T_TPCollection &trackingParticles, const TrackingParticleRefKeySet *trackingParticleKeys, iter begin, iter end) const
	Returns the TrackingParticle that has the most associated hits to the given track. More...

template<typename T_TPCollection , typename iter >
std::vector< std::pair< edm::Ref< TrackingParticleCollection >, double > >	associateTrack (const ClusterTPAssociation &clusterToTPMap, const T_TPCollection &trackingParticles, const TrackingParticleRefKeySet *trackingParticleKeys, iter begin, iter end) const
	Returns the TrackingParticle that has the most associated hits to the given track. More...

template<typename iter >
std::vector< std::pair< SimTrackIdentifiers, double > >	getAllSimTrackIdentifiers (const TrackerHitAssociator &hitAssociator, iter begin, iter end) const
	Returns a vector of pairs where first is a SimTrackIdentifiers (see typedef above) and second is the number of hits that came from that sim track. More...

template<typename iter >
double	getDoubleCount (const TrackerHitAssociator &hitAssociator, iter begin, iter end, TrackingParticleRef associatedTrackingParticle) const
	This method was copied almost verbatim from the standard TrackAssociatorByHits. More...

template<typename iter >
double	getDoubleCount (const ClusterTPAssociation &clusterToTPList, iter begin, iter end, TrackingParticleRef associatedTrackingParticle) const
	Overload for when using cluster to TrackingParticle association list. More...

bool	trackingParticleContainsIdentifier (const TrackingParticle *pTrackingParticle, const SimTrackIdentifiers &identifier) const
	Returns true if the supplied TrackingParticle has the supplied g4 track identifiers. More...

template<typename T_Track >
double	weightedNumberOfTrackClusters (const T_Track &track, const TrackerHitAssociator &) const

template<typename T_Track >
double	weightedNumberOfTrackClusters (const T_Track &track, const ClusterTPAssociation &) const

template<typename iter >
double	weightedNumberOfTrackClusters (iter begin, iter end) const

Static Private Member Functions
static bool	tpIntPairGreater (std::pair< edm::Ref< TrackingParticleCollection >, size_t > i, std::pair< edm::Ref< TrackingParticleCollection >, size_t > j)

Private Attributes
bool	absoluteNumberOfHits_

const ClusterTPAssociation *	clusterToTPMap_

double	cutRecoToSim_

std::unique_ptr< const TrackerHitAssociator >	hitAssociator_

double	pixelHitWeight_

edm::EDProductGetter const *	productGetter_
	creates either a ClusterTPAssociation OR a TrackerHitAssociator and stores it in the provided unique_ptr. The other will be null. More...

double	puritySimToReco_

double	qualitySimToReco_

SimToRecoDenomType	simToRecoDenominator_

bool	threeHitTracksAreSpecial_

Detailed Description

TrackToTrackingParticleAssociator that associates by hits a bit quicker than the normal TrackAssociatorByHitsImpl class.

NOTE - Doesn't implement the TrackCandidate association methods (from TrackAssociatorBase) so will always return empty associations for those.

This track associator (mostly) does the same as TrackAssociatorByHitsImpl, but faster. I've tested it a fair bit and can't find any differences between the results of this and the standard TrackAssociatorByHitsImpl.

Configuration parameters:

AbsoluteNumberOfHits - bool - if true, Quality_SimToReco and Cut_RecoToSim are the absolute number of shared hits required for association, not the percentage.

Quality_SimToReco - double - The minimum amount of shared hits required, as a percentage of either the reconstructed hits or simulated hits (see SimToRecoDenominator), for the track to be considered associated during a call to associateSimToReco. See also AbsoluteNumberOfHits.

Purity_SimToReco - double - The minimum amount of shared hits required, as a percentage of the reconstructed hits, for the track to be considered associated during a call to associateSimToReco. Has no effect if AbsoluteNumberOfHits is true.

Cut_RecoToSim - double - The minimum amount of shared hits required, as a percentage of the reconstructed hits, for the track to be considered associated during a call to associateRecoToSim. See also AbsoluteNumberOfHits.

ThreeHitTracksAreSpecial - bool - If true, tracks with 3 hits must have all their hits associated.

SimToRecoDenominator - string - Must be either "sim" or "reco". If "sim" Quality_SimToReco is the percentage of simulated hits that need to be shared. If "reco" then it's the percentage of reconstructed hits (i.e. same as Purity_SimToReco).

associatePixel - bool - Passed on to the hit associator.

associateStrip - bool - Passed on to the hit associator.

requireStoredHits - bool - Whether or not to insist all TrackingParticles have at least one PSimHit. The PSimHits are not required for the association, but the old TrackAssociatorByHitsImpl still had this requirement. Storing PSimHits in the TrackingParticle is now optional (see TrackingTruthAccumulator which replaces TrackingTruthProducer). Having requireStoredHits set to true will mean no TrackingParticles will be associated if you have chosen not to store the hits. The flag is only kept in order to retain the old behaviour which can give very slightly different results.

Note that the TrackAssociatorByHitsImpl parameters UseGrouped and UseSplitting are not used.

Author: Mark Grimes (mark..nosp@m.grim.nosp@m.es@ce.nosp@m.rn.c.nosp@m.h)

Date: 09/Nov/2010 Significant changes to remove any differences to the standard TrackAssociatorByHitsImpl results 07/Jul/2011. Association for TrajectorySeeds added by Giuseppe Cerati sometime between 2011 and 2013. Functionality to associate using pre calculated cluster to TrackingParticle maps added by Subir Sarker sometime in 2013. Overhauled to remove mutables to make it thread safe by Mark Grimes 01/May/2014.

Definition at line 66 of file QuickTrackAssociatorByHitsImpl.h.

Member Typedef Documentation

◆ SimTrackIdentifiers

typedef std::pair<uint32_t, EncodedEventId> QuickTrackAssociatorByHitsImpl::SimTrackIdentifiers

private

This is enough information to uniquely identify a sim track.

Definition at line 111 of file QuickTrackAssociatorByHitsImpl.h.

◆ TrackingParticleRefKeySet

typedef edm::IndexSet QuickTrackAssociatorByHitsImpl::TrackingParticleRefKeySet

private

Set for TrackingParticleRef keys.

Definition at line 113 of file QuickTrackAssociatorByHitsImpl.h.

Member Enumeration Documentation

◆ SimToRecoDenomType

enum QuickTrackAssociatorByHitsImpl::SimToRecoDenomType

Enumerator
denomnone
denomsim
denomreco

Definition at line 68 of file QuickTrackAssociatorByHitsImpl.h.

68 { denomnone, denomsim, denomreco };

QuickTrackAssociatorByHitsImpl::denomsim

Definition: QuickTrackAssociatorByHitsImpl.h:68

QuickTrackAssociatorByHitsImpl::denomnone

Definition: QuickTrackAssociatorByHitsImpl.h:68

QuickTrackAssociatorByHitsImpl::denomreco

Definition: QuickTrackAssociatorByHitsImpl.h:68

Constructor & Destructor Documentation

◆ QuickTrackAssociatorByHitsImpl()

QuickTrackAssociatorByHitsImpl::QuickTrackAssociatorByHitsImpl	(	edm::EDProductGetter const &	productGetter,
		std::unique_ptr< const TrackerHitAssociator >	hitAssoc,
		const ClusterTPAssociation *	clusterToTPMap,
		bool	absoluteNumberOfHits,
		double	qualitySimToReco,
		double	puritySimToReco,
		double	cutRecoToSim,
		double	pixelHitWeight,
		bool	threeHitTracksAreSpecial,
		SimToRecoDenomType	simToRecoDenominator
	)

Definition at line 111 of file QuickTrackAssociatorByHitsImpl.cc.

     : productGetter_(&productGetter),
       hitAssociator_(std::move(hitAssoc)),
       clusterToTPMap_(clusterToTPMap),
       qualitySimToReco_(qualitySimToReco),
       puritySimToReco_(puritySimToReco),
       pixelHitWeight_(pixelHitWeight),
       cutRecoToSim_(cutRecoToSim),
       simToRecoDenominator_(simToRecoDenominator),
       threeHitTracksAreSpecial_(threeHitTracksAreSpecial),
       absoluteNumberOfHits_(absoluteNumberOfHits) {}

Member Function Documentation

◆ associateRecoToSim() [1/4]

reco::RecoToSimCollection QuickTrackAssociatorByHitsImpl::associateRecoToSim	(	const edm::Handle< edm::View< reco::Track > > &	trackCollectionHandle,
		const edm::Handle< TrackingParticleCollection > &	trackingParticleCollectionHandle
	)		const

override

◆ associateRecoToSim() [2/4]

reco::RecoToSimCollection QuickTrackAssociatorByHitsImpl::associateRecoToSim	(	const edm::RefToBaseVector< reco::Track > &	tc,
		const edm::RefVector< TrackingParticleCollection > &	tpc
	)		const

overridevirtual

Association Reco To Sim with Collections.

Implements reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 157 of file QuickTrackAssociatorByHitsImpl.cc.

References clusterToTPMap_, hitAssociator_, JetHT_cfg::trackCollection, and PrimaryVertexAnalyzer4PUSlimmed_cfi::trackingParticleCollection.

                                                                                       {
   // Only pass the one that was successfully created to the templated method.
   if (not clusterToTPMap_)
     return associateRecoToSimImplementation<reco::RecoToSimCollection>(
         trackCollection, trackingParticleCollection, nullptr, *hitAssociator_);
   else {
     TrackingParticleRefKeySet tpKeys;
     fillKeys(tpKeys, trackingParticleCollection);
     return associateRecoToSimImplementation<reco::RecoToSimCollection>(
         trackCollection, trackingParticleCollection, &tpKeys, *clusterToTPMap_);
   }
 }

◆ associateRecoToSim() [3/4]

reco::RecoToSimCollectionSeed QuickTrackAssociatorByHitsImpl::associateRecoToSim	(	const edm::Handle< edm::View< TrajectorySeed > > &	,
		const edm::Handle< TrackingParticleCollection > &
	)		const

override

◆ associateRecoToSim() [4/4]

reco::RecoToSimCollectionTCandidate QuickTrackAssociatorByHitsImpl::associateRecoToSim	(	const edm::Handle< TrackCandidateCollection > &	trackCollectionHandle,
		const edm::RefVector< TrackingParticleCollection > &	trackingParticleCollectionHandle
	)		const

overridevirtual

Reimplemented from reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 781 of file QuickTrackAssociatorByHitsImpl.cc.

References clusterToTPMap_, and hitAssociator_.

                                                                                             {
   // Only pass the one that was successfully created to the templated method.
   if (not clusterToTPMap_)
     return associateRecoToSimImplementation<reco::RecoToSimCollectionTCandidate>(
         trackCollectionHandle, trackingParticleCollectionHandle, nullptr, *hitAssociator_);
   else
     return associateRecoToSimImplementation<reco::RecoToSimCollectionTCandidate>(
         trackCollectionHandle, trackingParticleCollectionHandle, nullptr, *clusterToTPMap_);
 }

◆ associateRecoToSimImplementation()

template<class T_RecoToSimCollection , class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >

T_RecoToSimCollection QuickTrackAssociatorByHitsImpl::associateRecoToSimImplementation	(	const T_TrackCollection &	trackCollection,
		const T_TrackingParticleCollection &	trackingParticleCollection,
		const TrackingParticleRefKeySet *	trackingParticleKeys,
		T_hitOrClusterAssociator	hitOrClusterAssociator
	)		const

private

The method that does the work for both overloads of associateRecoToSim.

Parts that actually rely on the type of the collections are delegated out to overloaded functions in the unnamed namespace of the .cc file. Parts that rely on the type of T_hitOrClusterAssociator are delegated out to overloaded methods.

Definition at line 191 of file QuickTrackAssociatorByHitsImpl.cc.

References funct::abs(), absoluteNumberOfHits_, associateTrack(), cutRecoToSim_, getDoubleCount(), mps_fire::i, TrackCollections2monitor_cff::numberOfValidHits, productGetter_, quality, threeHitTracksAreSpecial_, JetHT_cfg::trackCollection, PrimaryVertexAnalyzer4PUSlimmed_cfi::trackingParticleCollection, and weightedNumberOfTrackClusters().

                                                            {
   T_RecoToSimCollection returnValue(productGetter_);
   if (::collectionSize(trackingParticleCollection) == 0)
     return returnValue;
 
   checkClusterMapProductID(hitOrClusterAssociator, trackingParticleCollection);
 
   size_t collectionSize = ::collectionSize(trackCollection);  // Delegate away type specific part
 
   for (size_t i = 0; i < collectionSize; ++i) {
     // Get a normal pointer for ease of use. This part is type specific so delegate.
     const auto* pTrack = ::getTrackAt(trackCollection, i);
 
     // The return of this function has first as the index and second as the number of associated hits
     std::vector<std::pair<edm::Ref<TrackingParticleCollection>, double>> trackingParticleQualityPairs =
         associateTrack(hitOrClusterAssociator,
                        trackingParticleCollection,
                        trackingParticleKeys,
                        pTrack->recHits().begin(),
                        pTrack->recHits().end());
 
     // int nt = 0;
     for (auto iTrackingParticleQualityPair = trackingParticleQualityPairs.begin();
          iTrackingParticleQualityPair != trackingParticleQualityPairs.end();
          ++iTrackingParticleQualityPair) {
       const edm::Ref<TrackingParticleCollection>& trackingParticleRef = iTrackingParticleQualityPair->first;
       double numberOfSharedClusters = iTrackingParticleQualityPair->second;
       double numberOfValidTrackClusters = weightedNumberOfTrackClusters(*pTrack, hitOrClusterAssociator);
 
       if (numberOfSharedClusters == 0.0)
         continue;  // No point in continuing if there was no association
 
       //if electron subtract double counting
       if (abs(trackingParticleRef->pdgId()) == 11 &&
           (trackingParticleRef->g4Track_end() - trackingParticleRef->g4Track_begin()) > 1) {
         numberOfSharedClusters -= getDoubleCount(
             hitOrClusterAssociator, pTrack->recHits().begin(), pTrack->recHits().end(), trackingParticleRef);
       }
 
       double quality;
       if (absoluteNumberOfHits_)
         quality = numberOfSharedClusters;
       else if (numberOfValidTrackClusters != 0.0)
         quality = numberOfSharedClusters / numberOfValidTrackClusters;
       else
         quality = 0;
       if (quality > cutRecoToSim_ &&
           !(threeHitTracksAreSpecial_ && ::numberOfValidHits(*pTrack) == 3 && numberOfSharedClusters < 3.0)) {
         // Getting the RefToBase is dependent on the type of trackCollection, so delegate that to an overload.
         returnValue.insert(::getRefToTrackAt(trackCollection, i), std::make_pair(trackingParticleRef, quality));
       }
     }
   }
   returnValue.post_insert();
   return returnValue;
 }

◆ associateSimToReco() [1/4]

reco::SimToRecoCollection QuickTrackAssociatorByHitsImpl::associateSimToReco	(	const edm::Handle< edm::View< reco::Track > > &	trackCollectionHandle,
		const edm::Handle< TrackingParticleCollection > &	trackingParticleCollectionHandle
	)		const

override

◆ associateSimToReco() [2/4]

reco::SimToRecoCollection QuickTrackAssociatorByHitsImpl::associateSimToReco	(	const edm::RefToBaseVector< reco::Track > &	tc,
		const edm::RefVector< TrackingParticleCollection > &	tpc
	)		const

overridevirtual

Association Sim To Reco with Collections.

Implements reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 172 of file QuickTrackAssociatorByHitsImpl.cc.

References clusterToTPMap_, hitAssociator_, JetHT_cfg::trackCollection, and PrimaryVertexAnalyzer4PUSlimmed_cfi::trackingParticleCollection.

                                                                                       {
   // Only pass the one that was successfully created to the templated method.
   if (not clusterToTPMap_)
     return associateSimToRecoImplementation<reco::SimToRecoCollection>(
         trackCollection, trackingParticleCollection, nullptr, *hitAssociator_);
   else {
     TrackingParticleRefKeySet tpKeys;
     fillKeys(tpKeys, trackingParticleCollection);
     return associateSimToRecoImplementation<reco::SimToRecoCollection>(
         trackCollection, trackingParticleCollection, &tpKeys, *clusterToTPMap_);
   }
 }

◆ associateSimToReco() [3/4]

reco::SimToRecoCollectionSeed QuickTrackAssociatorByHitsImpl::associateSimToReco	(	const edm::Handle< edm::View< TrajectorySeed > > &	,
		const edm::Handle< TrackingParticleCollection > &
	)		const

override

◆ associateSimToReco() [4/4]

reco::SimToRecoCollectionTCandidate QuickTrackAssociatorByHitsImpl::associateSimToReco	(	const edm::Handle< TrackCandidateCollection > &	trackCollectionHandle,
		const edm::RefVector< TrackingParticleCollection > &	trackingParticleCollectionHandle
	)		const

overridevirtual

Reimplemented from reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 793 of file QuickTrackAssociatorByHitsImpl.cc.

References clusterToTPMap_, and hitAssociator_.

                                                                                             {
   // Only pass the one that was successfully created to the templated method.
   if (not clusterToTPMap_)
     return associateSimToRecoImplementation<reco::SimToRecoCollectionTCandidate>(
         trackCollectionHandle, trackingParticleCollectionHandle, nullptr, *hitAssociator_);
   else
     return associateSimToRecoImplementation<reco::SimToRecoCollectionTCandidate>(
         trackCollectionHandle, trackingParticleCollectionHandle, nullptr, *clusterToTPMap_);
 }

◆ associateSimToRecoImplementation()

template<class T_SimToRecoCollection , class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >

T_SimToRecoCollection QuickTrackAssociatorByHitsImpl::associateSimToRecoImplementation	(	const T_TrackCollection &	trackCollection,
		const T_TrackingParticleCollection &	trackingParticleCollection,
		const TrackingParticleRefKeySet *	trackingParticleKeys,
		T_hitOrClusterAssociator	hitOrClusterAssociator
	)		const

private

The method that does the work for both overloads of associateSimToReco.

Parts that actually rely on the type of the collections are delegated out to overloaded functions in the unnamed namespace of the .cc file. Parts that rely on the type of T_hitOrClusterAssociator are delegated out to overloaded methods.

Definition at line 256 of file QuickTrackAssociatorByHitsImpl.cc.

References funct::abs(), absoluteNumberOfHits_, associateTrack(), denomreco, denomsim, getDoubleCount(), mps_fire::i, productGetter_, puritySimToReco_, quality, qualitySimToReco_, simToRecoDenominator_, threeHitTracksAreSpecial_, JetHT_cfg::trackCollection, PrimaryVertexAnalyzer4PUSlimmed_cfi::trackingParticleCollection, and weightedNumberOfTrackClusters().

                                                            {
   T_SimToRecoCollection returnValue(productGetter_);
   if (::collectionSize(trackingParticleCollection) == 0)
     return returnValue;
 
   checkClusterMapProductID(hitOrClusterAssociator, trackingParticleCollection);
 
   size_t collectionSize = ::collectionSize(trackCollection);  // Delegate away type specific part
 
   for (size_t i = 0; i < collectionSize; ++i) {
     // Get a normal pointer for ease of use. This part is type specific so delegate.
     const auto* pTrack = ::getTrackAt(trackCollection, i);
 
     // The return of this function has first as an edm:Ref to the associated TrackingParticle, and second as the number of associated hits
     std::vector<std::pair<edm::Ref<TrackingParticleCollection>, double>> trackingParticleQualityPairs =
         associateTrack(hitOrClusterAssociator,
                        trackingParticleCollection,
                        trackingParticleKeys,
                        pTrack->recHits().begin(),
                        pTrack->recHits().end());
 
     // int nt = 0;
     for (auto iTrackingParticleQualityPair = trackingParticleQualityPairs.begin();
          iTrackingParticleQualityPair != trackingParticleQualityPairs.end();
          ++iTrackingParticleQualityPair) {
       const edm::Ref<TrackingParticleCollection>& trackingParticleRef = iTrackingParticleQualityPair->first;
       double numberOfSharedClusters = iTrackingParticleQualityPair->second;
       double numberOfValidTrackClusters = weightedNumberOfTrackClusters(*pTrack, hitOrClusterAssociator);
       size_t numberOfSimulatedHits = 0;  // Set a few lines below, but only if required.
 
       if (numberOfSharedClusters == 0.0)
         continue;  // No point in continuing if there was no association
 
       if (simToRecoDenominator_ == denomsim ||
           (numberOfSharedClusters < 3.0 &&
            threeHitTracksAreSpecial_))  // the numberOfSimulatedHits is not always required, so can skip counting in some circumstances
       {
         // Note that in the standard TrackAssociatorByHits, all of the hits in associatedTrackingParticleHits are checked for
         // various things.  I'm not sure what these checks are for but they depend on the UseGrouping and UseSplitting settings.
         // This associator works as though both UseGrouping and UseSplitting were set to true, i.e. just counts the number of
         // hits in the tracker.
         numberOfSimulatedHits = trackingParticleRef->numberOfTrackerHits();
       }
 
       //if electron subtract double counting
       if (abs(trackingParticleRef->pdgId()) == 11 &&
           (trackingParticleRef->g4Track_end() - trackingParticleRef->g4Track_begin()) > 1) {
         numberOfSharedClusters -= getDoubleCount(
             hitOrClusterAssociator, pTrack->recHits().begin(), pTrack->recHits().end(), trackingParticleRef);
       }
 
       double purity = numberOfSharedClusters / numberOfValidTrackClusters;
       double quality;
       if (absoluteNumberOfHits_)
         quality = numberOfSharedClusters;
       else if (simToRecoDenominator_ == denomsim && numberOfSimulatedHits != 0)
         quality = numberOfSharedClusters / static_cast<double>(numberOfSimulatedHits);
       else if (simToRecoDenominator_ == denomreco && numberOfValidTrackClusters != 0)
         quality = purity;
       else
         quality = 0;
 
       if (quality > qualitySimToReco_ &&
           !(threeHitTracksAreSpecial_ && numberOfSimulatedHits == 3 && numberOfSharedClusters < 3.0) &&
           (absoluteNumberOfHits_ || (purity > puritySimToReco_))) {
         // Getting the RefToBase is dependent on the type of trackCollection, so delegate that to an overload.
         returnValue.insert(trackingParticleRef, std::make_pair(::getRefToTrackAt(trackCollection, i), quality));
       }
     }
   }
   returnValue.post_insert();
   return returnValue;
 }

◆ associateTrack() [1/2]

template<typename T_TPCollection , typename iter >

std::vector< std::pair< edm::Ref< TrackingParticleCollection >, double > > QuickTrackAssociatorByHitsImpl::associateTrack	(	const TrackerHitAssociator &	hitAssociator,
		const T_TPCollection &	trackingParticles,
		const TrackingParticleRefKeySet *	trackingParticleKeys,
		iter	begin,
		iter	end
	)		const

private

Returns the TrackingParticle that has the most associated hits to the given track.

Return value is a vector of pairs, where first is an edm::Ref to the associated TrackingParticle, and second is the number of associated hits.

Definition at line 335 of file QuickTrackAssociatorByHitsImpl.cc.

References mps_fire::end, getAllSimTrackIdentifiers(), TrackQuality_cff::hitAssociator, mps_fire::i, trackingParticleContainsIdentifier(), and muonClassificationByHits_cfi::trackingParticles.

Referenced by associateRecoToSimImplementation(), and associateSimToRecoImplementation().

                     {
   // The pairs in this vector have a Ref to the associated TrackingParticle as "first" and the weighted number of associated hits as "second"
   std::vector<std::pair<edm::Ref<TrackingParticleCollection>, double>> returnValue;
 
   // The pairs in this vector have first as the sim track identifiers, and second the number of reco hits associated to that sim track.
   // Most reco hits will probably have come from the same sim track, so the number of entries in this vector should be fewer than the
   // number of reco hits.  The pair::second entries should add up to the total number of reco hits though.
   std::vector<std::pair<SimTrackIdentifiers, double>> hitIdentifiers =
       getAllSimTrackIdentifiers(hitAssociator, begin, end);
 
   // Loop over the TrackingParticles
   size_t collectionSize = ::collectionSize(trackingParticles);
 
   for (size_t i = 0; i < collectionSize; ++i) {
     const TrackingParticle* pTrackingParticle = getTrackingParticleAt(trackingParticles, i);
 
     // Historically there was a requirement that pTrackingParticle->numberOfHits() > 0
     // However, in TrackingTruthAccumulator, the numberOfHits is calculated from a subset
     // of the SimHits of the SimTracks of a TrackingParticle (essentially limiting the
     // processType and particleType to those of the "first" hit, and particleType to the pdgId of the SimTrack).
     // But, here the association between tracks and TrackingParticles is done with *all* the hits of
     // TrackingParticle, so we should not rely on the numberOfHits() calculated with a subset of SimHits.
 
     double numberOfAssociatedHits = 0;
     // Loop over all of the sim track identifiers and see if any of them are part of this TrackingParticle. If they are, add
     // the number of reco hits associated to that sim track to the total number of associated hits.
     for (const auto& identifierCountPair : hitIdentifiers) {
       if (trackingParticleContainsIdentifier(pTrackingParticle, identifierCountPair.first))
         numberOfAssociatedHits += identifierCountPair.second;
     }
 
     if (numberOfAssociatedHits > 0) {
       returnValue.push_back(std::make_pair(getRefToTrackingParticleAt(trackingParticles, i), numberOfAssociatedHits));
     }
   }
 
   return returnValue;
 }

◆ associateTrack() [2/2]

template<typename T_TPCollection , typename iter >

std::vector< std::pair< edm::Ref< TrackingParticleCollection >, double > > QuickTrackAssociatorByHitsImpl::associateTrack	(	const ClusterTPAssociation &	clusterToTPMap,
		const T_TPCollection &	trackingParticles,
		const TrackingParticleRefKeySet *	trackingParticleKeys,
		iter	begin,
		iter	end
	)		const

private

Returns the TrackingParticle that has the most associated hits to the given track.

See the notes for the other overload for the return type.

Note that the trackingParticles parameter is not actually required since all the information is in clusterToTPMap, but the method signature has to match the other overload because it is called from a templated method.

Definition at line 380 of file QuickTrackAssociatorByHitsImpl.cc.

References ClusterTPAssociation::empty(), mps_fire::end, ClusterTPAssociation::equal_range(), edm::IndexSet::has(), track_associator::hitsToClusterRefs(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, edm::Ref< C, T, F >::key(), pixelHitWeight_, and FastTimerService_cff::range.

                     {
   // Note that the trackingParticles parameter is not actually required since all the information is in clusterToTPMap,
   // but the method signature has to match the other overload because it is called from a templated method.
 
   // Note further, that we can't completely ignore the
   // trackingParticles parameter, in case it is a subset of those
   // TrackingParticles used to construct clusterToTPMap (via the
   // TrackingParticleRefVector overloads). The trackingParticles
   // parameter is still ignored since looping over it on every call
   // would be expensive, but the keys of the TrackingParticleRefs are
   // cached to an IndexSet (trackingParticleKeys) which is used
   // as a fast search structure.
 
   // The pairs in this vector have a Ref to the associated TrackingParticle as "first" and the weighted number of associated clusters as "second"
   // Note: typedef edm::Ref<TrackingParticleCollection> TrackingParticleRef;
   std::vector<std::pair<edm::Ref<TrackingParticleCollection>, double>> returnValue;
   if (clusterToTPMap.empty())
     return returnValue;
 
   // The pairs in this vector have first as the TP, and second the number of reco clusters associated to that TP.
   // Most reco clusters will probably have come from the same sim track (i.e TP), so the number of entries in this
   // vector should be fewer than the number of clusters. The pair::second entries should add up to the total
   // number of reco clusters though.
   std::vector<OmniClusterRef> oClusters = track_associator::hitsToClusterRefs(begin, end);
 
   std::map<TrackingParticleRef, double> lmap;
   for (std::vector<OmniClusterRef>::const_iterator it = oClusters.begin(); it != oClusters.end(); ++it) {
     auto range = clusterToTPMap.equal_range(*it);
     const double weight = it->isPixel() ? pixelHitWeight_ : 1.0;
     if (range.first != range.second) {
       for (auto ip = range.first; ip != range.second; ++ip) {
         const TrackingParticleRef trackingParticle = (ip->second);
 
         if (trackingParticleKeys && !trackingParticleKeys->has(trackingParticle.key()))
           continue;
 
         // Historically there was a requirement that pTrackingParticle->numberOfHits() > 0
         // However, in TrackingTruthAccumulator, the numberOfHits is calculated from a subset
         // of the SimHits of the SimTracks of a TrackingParticle (essentially limiting the
         // processType and particleType to those of the "first" hit, and particleType to the pdgId of the SimTrack).
         // But, here the association between tracks and TrackingParticles is done with *all* the hits of
         // TrackingParticle, so we should not rely on the numberOfHits() calculated with a subset of SimHits.
 
         /* Alternative implementation to avoid the use of lmap... memory slightly improved but slightly slower...
                  std::pair<edm::Ref<TrackingParticleCollection>,size_t> tpIntPair(trackingParticle, 1);
                  auto tp_range = std::equal_range(returnValue.begin(), returnValue.end(), tpIntPair, tpIntPairGreater);
                  if ((tp_range.second-tp_range.first)>1) {
                  edm::LogError("TrackAssociator") << ">>> Error in counting TPs!" << " file: " << __FILE__ << " line: " << __LINE__;
                  }
                  if(tp_range.first != tp_range.second) {
                  tp_range.first->second++;
                  } else {
                  returnValue.push_back(tpIntPair);
                  std::sort(returnValue.begin(), returnValue.end(), tpIntPairGreater);
                  }
                  */
         auto jpos = lmap.find(trackingParticle);
         if (jpos != lmap.end())
           jpos->second += weight;
         else
           lmap.insert(std::make_pair(trackingParticle, weight));
       }
     }
   }
   // now copy the map to returnValue
   for (auto ip = lmap.begin(); ip != lmap.end(); ++ip) {
     returnValue.push_back(std::make_pair(ip->first, ip->second));
   }
   return returnValue;
 }

◆ getAllSimTrackIdentifiers()

template<typename iter >

std::vector< std::pair< QuickTrackAssociatorByHitsImpl::SimTrackIdentifiers, double > > QuickTrackAssociatorByHitsImpl::getAllSimTrackIdentifiers	(	const TrackerHitAssociator &	hitAssociator,
		iter	begin,
		iter	end
	)		const

private

Returns a vector of pairs where first is a SimTrackIdentifiers (see typedef above) and second is the number of hits that came from that sim track.

This is used so that the TrackingParticle collection only has to be looped over once to search for each sim track, rather than once per hit. E.g. If all the hits in the reco track come from the same sim track, then there will only be one entry with second as the number of hits in the track.

Definition at line 458 of file QuickTrackAssociatorByHitsImpl.cc.

References mps_fire::end, TrackingRecHit::geographicalId(), track_associator::getHitFromIter(), TrackQuality_cff::hitAssociator, sistrip::SpyUtilities::isValid(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixelHitWeight_, DetId::subdetId(), and mps_merge::weight.

Referenced by associateTrack().

                                                                           {
   // The pairs in this vector have first as the sim track identifiers, and second the number of reco hits associated to that sim track.
   std::vector<std::pair<SimTrackIdentifiers, double>> returnValue;
 
   std::vector<SimTrackIdentifiers> simTrackIdentifiers;
   // Loop over all of the rec hits in the track
   //iter tRHIterBeginEnd = getTRHIterBeginEnd( pTrack );
   for (iter iRecHit = begin; iRecHit != end; ++iRecHit) {
     if (track_associator::getHitFromIter(iRecHit)->isValid()) {
       simTrackIdentifiers.clear();
 
       // Get the identifiers for the sim track that this hit came from. There should only be one entry unless clusters
       // have merged (as far as I know).
       hitAssociator.associateHitId(*(track_associator::getHitFromIter(iRecHit)),
                                    simTrackIdentifiers);  // This call fills simTrackIdentifiers
 
       const auto subdetId = track_associator::getHitFromIter(iRecHit)->geographicalId().subdetId();
       const double weight = (subdetId == PixelSubdetector::PixelBarrel || subdetId == PixelSubdetector::PixelEndcap)
                                 ? pixelHitWeight_
                                 : 1.0;
 
       // Loop over each identifier, and add it to the return value only if it's not already in there
       for (std::vector<SimTrackIdentifiers>::const_iterator iIdentifier = simTrackIdentifiers.begin();
            iIdentifier != simTrackIdentifiers.end();
            ++iIdentifier) {
         std::vector<std::pair<SimTrackIdentifiers, double>>::iterator iIdentifierCountPair;
         for (auto iIdentifierCountPair = returnValue.begin(); iIdentifierCountPair != returnValue.end();
              ++iIdentifierCountPair) {
           if (iIdentifierCountPair->first.first == iIdentifier->first &&
               iIdentifierCountPair->first.second == iIdentifier->second) {
             // This sim track identifier is already in the list, so increment the count of how many hits it relates to.
             iIdentifierCountPair->second += weight;
             break;
           }
         }
         if (iIdentifierCountPair == returnValue.end())
           returnValue.push_back(std::make_pair(*iIdentifier, 1.0));
         // This identifier wasn't found, so add it
       }
     }
   }
   return returnValue;
 }

◆ getDoubleCount() [1/2]

template<typename iter >

double QuickTrackAssociatorByHitsImpl::getDoubleCount	(	const TrackerHitAssociator &	hitAssociator,
		iter	begin,
		iter	end,
		TrackingParticleRef	associatedTrackingParticle
	)		const

private

This method was copied almost verbatim from the standard TrackAssociatorByHits.

Modified 01/May/2014 to take the TrackerHitAssociator as a parameter rather than using a member.

Definition at line 522 of file QuickTrackAssociatorByHitsImpl.cc.

References spr::find(), TrackingRecHit::geographicalId(), track_associator::getHitFromIter(), TrackQuality_cff::hitAssociator, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixelHitWeight_, and DetId::subdetId().

Referenced by associateRecoToSimImplementation(), and associateSimToRecoImplementation().

                                                                                                             {
   // This method is largely copied from the standard TrackAssociatorByHits. Once I've tested how much difference
   // it makes I'll go through and comment it properly.
 
   // FIXME: It may be that this piece is not fully correct for
   // counting how many times a single *cluster* is matched to many
   // SimTracks of a single TrackingParticle (see comments in
   // getDoubleCount(ClusterTPAssociation) overload). To be verified
   // some time.
 
   double doubleCount = 0.0;
   std::vector<SimHitIdpr> SimTrackIdsDC;
 
   for (iter iHit = startIterator; iHit != endIterator; iHit++) {
     int idcount = 0;
 
     SimTrackIdsDC.clear();
     hitAssociator.associateHitId(*(track_associator::getHitFromIter(iHit)), SimTrackIdsDC);
     if (SimTrackIdsDC.size() > 1) {
       for (TrackingParticle::g4t_iterator g4T = associatedTrackingParticle->g4Track_begin();
            g4T != associatedTrackingParticle->g4Track_end();
            ++g4T) {
         if (find(SimTrackIdsDC.begin(),
                  SimTrackIdsDC.end(),
                  SimHitIdpr((*g4T).trackId(), SimTrackIdsDC.begin()->second)) != SimTrackIdsDC.end()) {
           idcount++;
         }
       }
     }
     if (idcount > 1) {
       const auto subdetId = track_associator::getHitFromIter(iHit)->geographicalId().subdetId();
       const double weight = (subdetId == PixelSubdetector::PixelBarrel || subdetId == PixelSubdetector::PixelEndcap)
                                 ? pixelHitWeight_
                                 : 1.0;
       doubleCount += weight * (idcount - 1);
     }
   }
 
   return doubleCount;
 }

◆ getDoubleCount() [2/2]

template<typename iter >

double QuickTrackAssociatorByHitsImpl::getDoubleCount	(	const ClusterTPAssociation &	clusterToTPList,
		iter	begin,
		iter	end,
		TrackingParticleRef	associatedTrackingParticle
	)		const

private

Overload for when using cluster to TrackingParticle association list.

Definition at line 567 of file QuickTrackAssociatorByHitsImpl.cc.

References ClusterTPAssociation::equal_range(), TrackingRecHit::geographicalId(), track_associator::getHitFromIter(), track_associator::hitsToClusterRefs(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixelHitWeight_, FastTimerService_cff::range, and DetId::subdetId().

                                                                                                             {
   // This code here was written by Subir Sarkar. I'm just splitting it off into a
   // separate method. - Grimes 01/May/2014
 
   // The point here is that the electron TrackingParticles may contain
   // multiple SimTracks (from the bremsstrahling), and (historically)
   // the each matched hit/cluster has been multiplied by "how many
   // SimTracks from the TrackingParticle" it contains charge from.
   // Here the amount of this double counting is calculated, so that it
   // can be subtracted by the calling code.
   //
   // Note that recently (hence "historically" in the paragraph above)
   // the ClusterTPAssociationProducer was changed to remove the
   // duplicate cluster->TP associations (hence making this function
   // obsolete), but there is more recent proof that there is some
   // duplication left (to be investigated).
 
   double doubleCount = 0;
   std::vector<SimHitIdpr> SimTrackIdsDC;
 
   for (iter iHit = startIterator; iHit != endIterator; iHit++) {
     std::vector<OmniClusterRef> oClusters =
         track_associator::hitsToClusterRefs(iHit, iHit + 1);  //only for the cluster being checked
     for (std::vector<OmniClusterRef>::const_iterator it = oClusters.begin(); it != oClusters.end(); ++it) {
       int idcount = 0;
 
       auto range = clusterToTPList.equal_range(*it);
       if (range.first != range.second) {
         for (auto ip = range.first; ip != range.second; ++ip) {
           const TrackingParticleRef trackingParticle = (ip->second);
           if (associatedTrackingParticle == trackingParticle) {
             idcount++;
           }
         }
       }
 
       if (idcount > 1) {
         const auto subdetId = track_associator::getHitFromIter(iHit)->geographicalId().subdetId();
         const double weight = (subdetId == PixelSubdetector::PixelBarrel || subdetId == PixelSubdetector::PixelEndcap)
                                   ? pixelHitWeight_
                                   : 1.0;
         doubleCount += weight * (idcount - 1);
       }
     }
   }
 
   return doubleCount;
 }

◆ tpIntPairGreater()

static bool QuickTrackAssociatorByHitsImpl::tpIntPairGreater	(	std::pair< edm::Ref< TrackingParticleCollection >, size_t >	i,
		std::pair< edm::Ref< TrackingParticleCollection >, size_t >	j
	)

inlinestaticprivate

Definition at line 116 of file QuickTrackAssociatorByHitsImpl.h.

References mps_fire::i, and dqmiolumiharvest::j.

                                                                                       {
     return (i.first.key() > j.first.key());
   }

◆ trackingParticleContainsIdentifier()

bool QuickTrackAssociatorByHitsImpl::trackingParticleContainsIdentifier	(	const TrackingParticle *	pTrackingParticle,
		const SimTrackIdentifiers &	identifier
	)		const

private

Returns true if the supplied TrackingParticle has the supplied g4 track identifiers.

Definition at line 504 of file QuickTrackAssociatorByHitsImpl.cc.

References TrackingParticle::g4Track_begin(), and TrackingParticle::g4Track_end().

Referenced by associateTrack().

                                                                                                                      {
   // Loop over all of the g4 tracks in the tracking particle
   for (std::vector<SimTrack>::const_iterator iSimTrack = pTrackingParticle->g4Track_begin();
        iSimTrack != pTrackingParticle->g4Track_end();
        ++iSimTrack) {
     // And see if the sim track identifiers match
     if (iSimTrack->eventId() == identifier.second && iSimTrack->trackId() == identifier.first) {
       return true;
     }
   }
 
   // If control has made it this far then none of the identifiers were found in
   // any of the g4 tracks, so return false.
   return false;
 }

◆ weightedNumberOfTrackClusters() [1/5]

template<typename T_Track >

double QuickTrackAssociatorByHitsImpl::weightedNumberOfTrackClusters	(	const T_Track &	track,
		const TrackerHitAssociator &
	)		const

private

Definition at line 816 of file QuickTrackAssociatorByHitsImpl.cc.

References TrackingRecHit::geographicalId(), track_associator::getHitFromIter(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixelHitWeight_, fileCollector::seed, DetId::subdetId(), and mps_merge::weight.

Referenced by associateRecoToSimImplementation(), associateSimToRecoImplementation(), and weightedNumberOfTrackClusters().

                                                                                                         {
   double sum = 0.0;
   for (auto iHit = seed.recHits().begin(); iHit != seed.recHits().end(); ++iHit) {
     const auto subdetId = track_associator::getHitFromIter(iHit)->geographicalId().subdetId();
     const double weight = (subdetId == PixelSubdetector::PixelBarrel || subdetId == PixelSubdetector::PixelEndcap)
                               ? pixelHitWeight_
                               : 1.0;
     sum += weight;
   }
   return sum;
 }

◆ weightedNumberOfTrackClusters() [2/5]

template<typename T_Track >

double QuickTrackAssociatorByHitsImpl::weightedNumberOfTrackClusters	(	const T_Track &	track,
		const ClusterTPAssociation &
	)		const

private

Definition at line 836 of file QuickTrackAssociatorByHitsImpl.cc.

References fileCollector::seed, and weightedNumberOfTrackClusters().

                                                                                                         {
   const auto& hitRange = seed.recHits();
   return weightedNumberOfTrackClusters(hitRange.begin(), hitRange.end());
 }

◆ weightedNumberOfTrackClusters() [3/5]

template<typename iter >

double QuickTrackAssociatorByHitsImpl::weightedNumberOfTrackClusters	(	iter	begin,
		iter	end
	)		const

private

Definition at line 843 of file QuickTrackAssociatorByHitsImpl.cc.

References mps_fire::end, TrackingRecHit::geographicalId(), track_associator::getHitFromIter(), track_associator::hitsToClusterRefs(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LogTrace, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pixelHitWeight_, DetId::rawId(), DetId::subdetId(), and mps_merge::weight.

                                                                                                {
   double weightedClusters = 0.0;
   for (iter iRecHit = begin; iRecHit != end; ++iRecHit) {
     const auto subdetId = track_associator::getHitFromIter(iRecHit)->geographicalId().subdetId();
     const double weight = (subdetId == PixelSubdetector::PixelBarrel || subdetId == PixelSubdetector::PixelEndcap)
                               ? pixelHitWeight_
                               : 1.0;
     LogTrace("QuickTrackAssociatorByHitsImpl")
         << "  detId: " << track_associator::getHitFromIter(iRecHit)->geographicalId().rawId();
     LogTrace("QuickTrackAssociatorByHitsImpl") << "  weight: " << weight;
     std::vector<OmniClusterRef> oClusters =
         track_associator::hitsToClusterRefs(iRecHit, iRecHit + 1);  //only for the cluster being checked
     for (std::vector<OmniClusterRef>::const_iterator it = oClusters.begin(); it != oClusters.end(); ++it) {
       weightedClusters += weight;
     }
   }
   LogTrace("QuickTrackAssociatorByHitsImpl") << "  total weighted clusters: " << weightedClusters;
   return weightedClusters;
 }

◆ weightedNumberOfTrackClusters() [4/5]

template<>

double QuickTrackAssociatorByHitsImpl::weightedNumberOfTrackClusters	(	const reco::Track &	track,
		const TrackerHitAssociator &
	)		const

Definition at line 807 of file QuickTrackAssociatorByHitsImpl.cc.

References AlCaHLTBitMon_ParallelJobs::p, InitialStepPreSplitting_cff::pixelHits, pixelHitWeight_, and HLT_2024v14_cff::track.

                                                                                                         {
   const reco::HitPattern& p = track.hitPattern();
   const auto pixelHits = p.numberOfValidPixelHits();
   const auto otherHits = p.numberOfValidHits() - pixelHits;
   return pixelHits * pixelHitWeight_ + otherHits;
 }

◆ weightedNumberOfTrackClusters() [5/5]

template<>

double QuickTrackAssociatorByHitsImpl::weightedNumberOfTrackClusters	(	const reco::Track &	track,
		const ClusterTPAssociation &
	)		const

Definition at line 831 of file QuickTrackAssociatorByHitsImpl.cc.

References HLT_2024v14_cff::track, and weightedNumberOfTrackClusters().

                                                                                                         {
   return weightedNumberOfTrackClusters(track.recHitsBegin(), track.recHitsEnd());
 }

Member Data Documentation

◆ absoluteNumberOfHits_

bool QuickTrackAssociatorByHitsImpl::absoluteNumberOfHits_

private

Definition at line 244 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSimImplementation(), and associateSimToRecoImplementation().

◆ clusterToTPMap_

const ClusterTPAssociation* QuickTrackAssociatorByHitsImpl::clusterToTPMap_

private

Definition at line 236 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSim(), and associateSimToReco().

◆ cutRecoToSim_

double QuickTrackAssociatorByHitsImpl::cutRecoToSim_

private

Definition at line 241 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSimImplementation().

◆ hitAssociator_

std::unique_ptr<const TrackerHitAssociator> QuickTrackAssociatorByHitsImpl::hitAssociator_

private

Definition at line 235 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSim(), and associateSimToReco().

◆ pixelHitWeight_

double QuickTrackAssociatorByHitsImpl::pixelHitWeight_

private

Definition at line 240 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateTrack(), getAllSimTrackIdentifiers(), getDoubleCount(), and weightedNumberOfTrackClusters().

◆ productGetter_

edm::EDProductGetter const* QuickTrackAssociatorByHitsImpl::productGetter_

private

creates either a ClusterTPAssociation OR a TrackerHitAssociator and stores it in the provided unique_ptr. The other will be null.

A decision is made whether to create a ClusterTPAssociation or a TrackerHitAssociator depending on how this track associator was configured. If the ClusterTPAssociation couldn't be fetched from the event then it falls back to creating a TrackerHitAssociator.

Only one type will be created, never both. The other unique_ptr reference will be null so check for that and decide which to use.

N.B. The value of useClusterTPAssociation_ should not be used to decide which of the two pointers to use. If the cluster to TrackingParticle couldn't be retrieved from the event then pClusterToTPMap will be null but useClusterTPAssociation_ is no longer changed to false.

Definition at line 234 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSimImplementation(), and associateSimToRecoImplementation().