#include <QuickTrackAssociatorByHitsImpl.h>

Inheritance diagram for QuickTrackAssociatorByHitsImpl:

Public Types
typedef std::vector< std::pair < OmniClusterRef, TrackingParticleRef > >	ClusterTPAssociationList

enum	SimToRecoDenomType { denomnone, denomsim, denomreco }

Public Member Functions
virtual reco::RecoToSimCollection	associateRecoToSim (const edm::Handle< edm::View< reco::Track > > &trackCollectionHandle, const edm::Handle< TrackingParticleCollection > &trackingParticleCollectionHandle) const override
	compare reco to sim the handle of reco::Track and TrackingParticle collections More...

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

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

virtual reco::SimToRecoCollection	associateSimToReco (const edm::Handle< edm::View< reco::Track > > &trackCollectionHandle, const edm::Handle< TrackingParticleCollection > &trackingParticleCollectionHandle) const override
	compare reco to sim the handle of reco::Track and TrackingParticle collections More...

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

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

	QuickTrackAssociatorByHitsImpl (std::shared_ptr< const TrackerHitAssociator > hitAssoc, std::shared_ptr< const ClusterTPAssociationList > clusterToTPMap, bool absoluteNumberOfHits, double qualitySimToReco, double puritySimToReco, double cutRecoToSim, bool threeHitTracksAreSpecial, SimToRecoDenomType simToRecoDenominator)

Public Member Functions inherited from reco::TrackToTrackingParticleAssociatorBaseImpl
virtual reco::RecoToSimCollectionTCandidate	associateRecoToSim (const edm::Handle< TrackCandidateCollection > &, const edm::Handle< TrackingParticleCollection > &) const

virtual reco::SimToRecoCollectionTCandidate	associateSimToReco (const edm::Handle< TrackCandidateCollection > &, 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...

Private Member Functions
template<class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >
reco::RecoToSimCollection	associateRecoToSimImplementation (T_TrackCollection trackCollection, T_TrackingParticleCollection trackingParticleCollection, T_hitOrClusterAssociator hitOrClusterAssociator) const
	The method that does the work for both overloads of associateRecoToSim. More...

template<class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >
reco::SimToRecoCollection	associateSimToRecoImplementation (T_TrackCollection trackCollection, T_TrackingParticleCollection trackingParticleCollection, 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 > , size_t > >	associateTrack (const TrackerHitAssociator &hitAssociator, T_TPCollection trackingParticles, 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 > , size_t > >	associateTrack (const ClusterTPAssociationList &clusterToTPMap, T_TPCollection trackingParticles, 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, size_t > >	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 >
int	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 >
int	getDoubleCount (const ClusterTPAssociationList &clusterToTPList, iter begin, iter end, TrackingParticleRef associatedTrackingParticle) const
	Overload for when using cluster to TrackingParticle association list. More...

const TrackingRecHit *	getHitFromIter (trackingRecHit_iterator iter) const

const TrackingRecHit *	getHitFromIter (TrackingRecHitCollection::const_iterator iter) const

template<typename iter >
std::vector< OmniClusterRef >	getMatchedClusters (iter begin, iter end) const

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

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_

std::shared_ptr< const ClusterTPAssociationList >	clusterToTPMap_

double	cutRecoToSim_

std::shared_ptr< const TrackerHitAssociator >	hitAssociator_
	creates either a ClusterTPAssociationList 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

Definition at line 61 of file QuickTrackAssociatorByHitsImpl.h.

Member Typedef Documentation

typedef std::vector<std::pair<OmniClusterRef, TrackingParticleRef> > QuickTrackAssociatorByHitsImpl::ClusterTPAssociationList

Definition at line 64 of file QuickTrackAssociatorByHitsImpl.h.

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

private

This is enough information to uniquely identify a sim track.

Definition at line 101 of file QuickTrackAssociatorByHitsImpl.h.

Member Enumeration Documentation

enum QuickTrackAssociatorByHitsImpl::SimToRecoDenomType

Enumerator
denomnone
denomsim
denomreco

Definition at line 65 of file QuickTrackAssociatorByHitsImpl.h.

65 {denomnone,denomsim,denomreco};

QuickTrackAssociatorByHitsImpl::denomsim

Definition: QuickTrackAssociatorByHitsImpl.h:65

QuickTrackAssociatorByHitsImpl::denomnone

Definition: QuickTrackAssociatorByHitsImpl.h:65

QuickTrackAssociatorByHitsImpl::denomreco

Definition: QuickTrackAssociatorByHitsImpl.h:65

Constructor & Destructor Documentation

QuickTrackAssociatorByHitsImpl::QuickTrackAssociatorByHitsImpl	(	std::shared_ptr< const TrackerHitAssociator >	hitAssoc,
		std::shared_ptr< const ClusterTPAssociationList >	clusterToTPMap,
		bool	absoluteNumberOfHits,
		double	qualitySimToReco,
		double	puritySimToReco,
		double	cutRecoToSim,
		bool	threeHitTracksAreSpecial,
		SimToRecoDenomType	simToRecoDenominator
	)

Definition at line 88 of file QuickTrackAssociatorByHitsImpl.cc.

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

Member Function Documentation

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

overridevirtual

compare reco to sim the handle of reco::Track and TrackingParticle collections

Reimplemented from reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 108 of file QuickTrackAssociatorByHitsImpl.cc.

References associateRecoToSimImplementation(), clusterToTPMap_, and hitAssociator_.

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

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 124 of file QuickTrackAssociatorByHitsImpl.cc.

References associateRecoToSimImplementation(), clusterToTPMap_, and hitAssociator_.

 {
     // Only pass the one that was successfully created to the templated method.
     if( not clusterToTPMap_ ) return associateRecoToSimImplementation( trackCollection, trackingParticleCollection, *hitAssociator_ );
     else return associateRecoToSimImplementation( trackCollection, trackingParticleCollection, *clusterToTPMap_ );
 }

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

overridevirtual

Reimplemented from reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 516 of file QuickTrackAssociatorByHitsImpl.cc.

References funct::abs(), absoluteNumberOfHits_, associateTrack(), clusterToTPMap_, cutRecoToSim_, getDoubleCount(), hitAssociator_, i, edm::AssociationMap< Tag >::insert(), LogTrace, edm::AssociationMap< Tag >::post_insert(), HLT_25ns14e33_v1_cff::quality, TrajectorySeed::recHits(), edm::AssociationMap< Tag >::size(), and threeHitTracksAreSpecial_.

 {
 
     edm::LogVerbatim( "TrackAssociator" ) << "Starting TrackAssociatorByHitsImpl::associateRecoToSim - #seeds=" << pSeedCollectionHandle_->size()
             << " #TPs=" << trackingParticleCollectionHandle->size();
 
     reco::RecoToSimCollectionSeed returnValue;
 
     size_t collectionSize=pSeedCollectionHandle_->size();
 
     for( size_t i=0; i < collectionSize; ++i )
     {
         const TrajectorySeed* pSeed= &( *pSeedCollectionHandle_)[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>,size_t> > trackingParticleQualityPairs=
                 (clusterToTPMap_) ? associateTrack( *clusterToTPMap_, trackingParticleCollectionHandle, pSeed->recHits().first, pSeed->recHits().second ) : associateTrack( *hitAssociator_, trackingParticleCollectionHandle, pSeed->recHits().first, pSeed->recHits().second );
         for( std::vector<std::pair<edm::Ref<TrackingParticleCollection>,size_t> >::const_iterator iTrackingParticleQualityPair=
                 trackingParticleQualityPairs.begin(); iTrackingParticleQualityPair != trackingParticleQualityPairs.end();
                 ++iTrackingParticleQualityPair )
         {
             const edm::Ref<TrackingParticleCollection>& trackingParticleRef=iTrackingParticleQualityPair->first;
             size_t numberOfSharedHits=iTrackingParticleQualityPair->second;
             size_t numberOfValidTrackHits=pSeed->recHits().second - pSeed->recHits().first;
 
             if( numberOfSharedHits == 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 )
             {
                 if( clusterToTPMap_ ) numberOfSharedHits-=getDoubleCount( *clusterToTPMap_, pSeed->recHits().first, pSeed->recHits().second, trackingParticleRef );
                 else numberOfSharedHits-=getDoubleCount( *hitAssociator_, pSeed->recHits().first, pSeed->recHits().second, trackingParticleRef );
             }
 
             double quality;
             if( absoluteNumberOfHits_ ) quality=static_cast<double>( numberOfSharedHits );
             else if( numberOfValidTrackHits != 0 ) quality=
                     (static_cast<double>( numberOfSharedHits ) / static_cast<double>( numberOfValidTrackHits ));
             else quality=0;
 
             if( quality > cutRecoToSim_ && !(threeHitTracksAreSpecial_ && numberOfValidTrackHits == 3 && numberOfSharedHits < 3) )
             {
                 returnValue.insert( edm::RefToBase < TrajectorySeed > (pSeedCollectionHandle_, i), std::make_pair( trackingParticleRef, quality ) );
             }
         }
     }
 
     LogTrace( "TrackAssociator" ) << "% of Assoc Seeds=" << ((double)returnValue.size()) / ((double)pSeedCollectionHandle_->size());
     returnValue.post_insert();
     return returnValue;
 
 }

template<class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >

reco::RecoToSimCollection QuickTrackAssociatorByHitsImpl::associateRecoToSimImplementation	(	T_TrackCollection	trackCollection,
		T_TrackingParticleCollection	trackingParticleCollection,
		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 142 of file QuickTrackAssociatorByHitsImpl.cc.

References funct::abs(), absoluteNumberOfHits_, associateTrack(), cutRecoToSim_, reco::Track::found(), getDoubleCount(), i, edm::AssociationMap< Tag >::insert(), HLT_25ns14e33_v1_cff::quality, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), and threeHitTracksAreSpecial_.

Referenced by associateRecoToSim().

 {
     reco::RecoToSimCollection returnValue;
 
     size_t collectionSize=::collectionSize(trackCollection); // Delegate away type specific part
 
     for( size_t i=0; i < collectionSize; ++i )
     {
         const reco::Track* pTrack=::getTrackAt(trackCollection,i); // Get a normal pointer for ease of use. This part is type specific so delegate.
 
         // 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>,size_t> > trackingParticleQualityPairs=associateTrack( hitOrClusterAssociator, trackingParticleCollection, pTrack->recHitsBegin(), pTrack->recHitsEnd() );
 
         // int nt = 0;
         for( std::vector<std::pair<edm::Ref<TrackingParticleCollection>,size_t> >::const_iterator iTrackingParticleQualityPair=
                 trackingParticleQualityPairs.begin(); iTrackingParticleQualityPair != trackingParticleQualityPairs.end();
                 ++iTrackingParticleQualityPair )
         {
             const edm::Ref<TrackingParticleCollection>& trackingParticleRef=iTrackingParticleQualityPair->first;
             size_t numberOfSharedHits=iTrackingParticleQualityPair->second;
             size_t numberOfValidTrackHits=pTrack->found();
 
             if( numberOfSharedHits == 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 )
             {
                 numberOfSharedHits-=getDoubleCount( hitOrClusterAssociator, pTrack->recHitsBegin(), pTrack->recHitsEnd(), trackingParticleRef );
             }
 
             double quality;
             if( absoluteNumberOfHits_ ) quality=static_cast<double>( numberOfSharedHits );
             else if( numberOfValidTrackHits != 0 ) quality=
                     (static_cast<double>( numberOfSharedHits ) / static_cast<double>( numberOfValidTrackHits ));
             else quality=0;
             if( quality > cutRecoToSim_ && !(threeHitTracksAreSpecial_ && numberOfValidTrackHits == 3 && numberOfSharedHits < 3) )
             {
                 // 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 ) );
             }
         }
     }
     return returnValue;
 }

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

overridevirtual

compare reco to sim the handle of reco::Track and TrackingParticle collections

Reimplemented from reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 116 of file QuickTrackAssociatorByHitsImpl.cc.

References associateSimToRecoImplementation(), clusterToTPMap_, and hitAssociator_.

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

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 132 of file QuickTrackAssociatorByHitsImpl.cc.

References associateSimToRecoImplementation(), clusterToTPMap_, and hitAssociator_.

 {
     // Only pass the one that was successfully created to the templated method.
     if( not clusterToTPMap_ ) return associateSimToRecoImplementation( trackCollection, trackingParticleCollection, *hitAssociator_ );
     else return associateSimToRecoImplementation( trackCollection, trackingParticleCollection, *clusterToTPMap_ );
 }

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

overridevirtual

Reimplemented from reco::TrackToTrackingParticleAssociatorBaseImpl.

Definition at line 570 of file QuickTrackAssociatorByHitsImpl.cc.

References funct::abs(), absoluteNumberOfHits_, associateTrack(), clusterToTPMap_, denomreco, denomsim, getDoubleCount(), hitAssociator_, i, edm::AssociationMap< Tag >::insert(), LogTrace, edm::AssociationMap< Tag >::post_insert(), puritySimToReco_, HLT_25ns14e33_v1_cff::quality, qualitySimToReco_, TrajectorySeed::recHits(), simToRecoDenominator_, edm::AssociationMap< Tag >::size(), and threeHitTracksAreSpecial_.

 {
 
     edm::LogVerbatim( "TrackAssociator" ) << "Starting TrackAssociatorByHitsImpl::associateSimToReco - #seeds=" << pSeedCollectionHandle_->size()
             << " #TPs=" << trackingParticleCollectionHandle->size();
 
     reco::SimToRecoCollectionSeed returnValue;
 
     size_t collectionSize=pSeedCollectionHandle_->size();
 
     for( size_t i=0; i < collectionSize; ++i )
     {
         const TrajectorySeed* pSeed= &( *pSeedCollectionHandle_)[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>,size_t> > trackingParticleQualityPairs=
                 (clusterToTPMap_) ? associateTrack( *clusterToTPMap_, trackingParticleCollectionHandle, pSeed->recHits().first, pSeed->recHits().second ) : associateTrack( *hitAssociator_, trackingParticleCollectionHandle, pSeed->recHits().first, pSeed->recHits().second );
         for( std::vector<std::pair<edm::Ref<TrackingParticleCollection>,size_t> >::const_iterator iTrackingParticleQualityPair=
                 trackingParticleQualityPairs.begin(); iTrackingParticleQualityPair != trackingParticleQualityPairs.end();
                 ++iTrackingParticleQualityPair )
         {
             const edm::Ref<TrackingParticleCollection>& trackingParticleRef=iTrackingParticleQualityPair->first;
             size_t numberOfSharedHits=iTrackingParticleQualityPair->second;
             size_t numberOfValidTrackHits=pSeed->recHits().second - pSeed->recHits().first;
             size_t numberOfSimulatedHits=0; // Set a few lines below, but only if required.
 
             if( numberOfSharedHits == 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 )
             {
                 if( clusterToTPMap_ ) numberOfSharedHits-=getDoubleCount( *clusterToTPMap_, pSeed->recHits().first, pSeed->recHits().second, trackingParticleRef );
                 else numberOfSharedHits-=getDoubleCount( *hitAssociator_, pSeed->recHits().first, pSeed->recHits().second, trackingParticleRef );
             }
 
             if( simToRecoDenominator_ == denomsim || (numberOfSharedHits < 3 && 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();
             }
 
             double purity=static_cast<double>( numberOfSharedHits ) / static_cast<double>( numberOfValidTrackHits );
             double quality;
             if( absoluteNumberOfHits_ ) quality=static_cast<double>( numberOfSharedHits );
             else if( simToRecoDenominator_ == denomsim && numberOfSimulatedHits != 0 ) quality=static_cast<double>( numberOfSharedHits )
                     / static_cast<double>( numberOfSimulatedHits );
             else if( simToRecoDenominator_ == denomreco && numberOfValidTrackHits != 0 ) quality=purity;
             else quality=0;
 
             if( quality > qualitySimToReco_ && !(threeHitTracksAreSpecial_ && numberOfSimulatedHits == 3 && numberOfSharedHits < 3)
                     && (absoluteNumberOfHits_ || (purity > puritySimToReco_)) )
             {
                 returnValue.insert( trackingParticleRef, std::make_pair( edm::RefToBase < TrajectorySeed > (pSeedCollectionHandle_, i), quality ) );
             }
         }
     }
     return returnValue;
 
     LogTrace("TrackAssociator") << "% of Assoc TPs=" << ((double)returnValue.size())/((double)trackingParticleCollectionHandle->size());
     returnValue.post_insert();
     return returnValue;
 }

template<class T_TrackCollection , class T_TrackingParticleCollection , class T_hitOrClusterAssociator >

reco::SimToRecoCollection QuickTrackAssociatorByHitsImpl::associateSimToRecoImplementation	(	T_TrackCollection	trackCollection,
		T_TrackingParticleCollection	trackingParticleCollection,
		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 188 of file QuickTrackAssociatorByHitsImpl.cc.

References funct::abs(), absoluteNumberOfHits_, associateTrack(), denomreco, denomsim, reco::Track::found(), getDoubleCount(), i, edm::AssociationMap< Tag >::insert(), puritySimToReco_, HLT_25ns14e33_v1_cff::quality, qualitySimToReco_, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), simToRecoDenominator_, and threeHitTracksAreSpecial_.

Referenced by associateSimToReco().

 {
     reco::SimToRecoCollection returnValue;
 
     size_t collectionSize=::collectionSize(trackCollection); // Delegate away type specific part
 
     for( size_t i=0; i<collectionSize; ++i )
     {
         const reco::Track* pTrack=::getTrackAt(trackCollection,i); // Get a normal pointer for ease of use. This part is type specific so delegate.
 
         // 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>,size_t> > trackingParticleQualityPairs=associateTrack( hitOrClusterAssociator, trackingParticleCollection, pTrack->recHitsBegin(), pTrack->recHitsEnd() );
 
         // int nt = 0;
         for( std::vector< std::pair<edm::Ref<TrackingParticleCollection>,size_t> >::const_iterator iTrackingParticleQualityPair=trackingParticleQualityPairs.begin();
                 iTrackingParticleQualityPair!=trackingParticleQualityPairs.end(); ++iTrackingParticleQualityPair )
         {
             const edm::Ref<TrackingParticleCollection>& trackingParticleRef=iTrackingParticleQualityPair->first;
             size_t numberOfSharedHits=iTrackingParticleQualityPair->second;
             size_t numberOfValidTrackHits=pTrack->found();
             size_t numberOfSimulatedHits=0; // Set a few lines below, but only if required.
 
             if( numberOfSharedHits==0 ) continue; // No point in continuing if there was no association
 
             if( simToRecoDenominator_==denomsim || (numberOfSharedHits<3 && 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 )
             {
                 numberOfSharedHits -= getDoubleCount( hitOrClusterAssociator, pTrack->recHitsBegin(), pTrack->recHitsEnd(), trackingParticleRef );
             }
 
             double purity=static_cast<double>(numberOfSharedHits)/static_cast<double>(numberOfValidTrackHits);
             double quality;
             if( absoluteNumberOfHits_ ) quality=static_cast<double>(numberOfSharedHits);
             else if( simToRecoDenominator_==denomsim && numberOfSimulatedHits != 0 ) quality=static_cast<double>(numberOfSharedHits)/static_cast<double>(numberOfSimulatedHits);
             else if( simToRecoDenominator_==denomreco && numberOfValidTrackHits != 0 ) quality=purity;
             else quality=0;
 
             if( quality>qualitySimToReco_ && !( threeHitTracksAreSpecial_ && numberOfSimulatedHits==3 && numberOfSharedHits<3 ) && ( 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 ) );
             }
         }
     }
     return returnValue;
 
 }

template<typename T_TPCollection , typename iter >

std::vector< std::pair< edm::Ref< TrackingParticleCollection >, size_t > > QuickTrackAssociatorByHitsImpl::associateTrack	(	const TrackerHitAssociator &	hitAssociator,
		T_TPCollection	trackingParticles,
		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 245 of file QuickTrackAssociatorByHitsImpl.cc.

References getAllSimTrackIdentifiers(), i, TrackingParticle::numberOfHits(), and trackingParticleContainsIdentifier().

Referenced by associateRecoToSim(), associateRecoToSimImplementation(), associateSimToReco(), and associateSimToRecoImplementation().

 {
     // The pairs in this vector have a Ref to the associated TrackingParticle as "first" and the number of associated hits as "second"
     std::vector< std::pair<edm::Ref<TrackingParticleCollection>,size_t> > 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,size_t> > 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 );
 
         // Ignore TrackingParticles with no hits
         if( pTrackingParticle->numberOfHits()==0 ) continue;
 
         size_t 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( std::vector< std::pair<SimTrackIdentifiers,size_t> >::const_iterator iIdentifierCountPair=hitIdentifiers.begin(); iIdentifierCountPair!=hitIdentifiers.end(); ++iIdentifierCountPair )
         {
             if( trackingParticleContainsIdentifier( pTrackingParticle, iIdentifierCountPair->first ) ) numberOfAssociatedHits+=iIdentifierCountPair->second;
         }
 
         if( numberOfAssociatedHits>0 )
         {
             returnValue.push_back( std::make_pair( getRefToTrackingParticleAt(trackingParticles,i), numberOfAssociatedHits ) );
         }
     }
 
     return returnValue;
 }

template<typename T_TPCollection , typename iter >

std::vector< std::pair< edm::Ref< TrackingParticleCollection >, size_t > > QuickTrackAssociatorByHitsImpl::associateTrack	(	const ClusterTPAssociationList &	clusterToTPMap,
		T_TPCollection	trackingParticles,
		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 282 of file QuickTrackAssociatorByHitsImpl.cc.

References clusterTPAssociationListGreater(), and getMatchedClusters().

 {
     // 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.
 
     // The pairs in this vector have a Ref to the associated TrackingParticle as "first" and the number of associated clusters as "second"
     // Note: typedef edm::Ref<TrackingParticleCollection> TrackingParticleRef;
     std::vector < std::pair<edm::Ref<TrackingParticleCollection>,size_t> > 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=getMatchedClusters( begin, end );
 
     std::map < TrackingParticleRef, size_t > lmap;
     for( std::vector<OmniClusterRef>::const_iterator it=oClusters.begin(); it != oClusters.end(); ++it )
     {
 
         std::pair < OmniClusterRef, TrackingParticleRef > clusterTPpairWithDummyTP( *it, TrackingParticleRef() ); //TP is dummy: for clusterTPAssociationListGreater sorting only the cluster is needed
         auto range=std::equal_range( clusterToTPMap.begin(), clusterToTPMap.end(), clusterTPpairWithDummyTP, clusterTPAssociationListGreater );
         if( range.first != range.second )
         {
             for( auto ip=range.first; ip != range.second; ++ip )
             {
 
                 const TrackingParticleRef trackingParticle=(ip->second);
 
                 // Ignore TrackingParticles with no hits
                 if( trackingParticle->numberOfHits() == 0 ) continue;
 
                 /* 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;
                 else lmap.insert( std::make_pair( trackingParticle, 1 ) );
             }
         }
     }
     // 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;
 }

template<typename iter >

std::vector< std::pair< QuickTrackAssociatorByHitsImpl::SimTrackIdentifiers, size_t > > 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 391 of file QuickTrackAssociatorByHitsImpl.cc.

References TrackerHitAssociator::associateHitId(), end, getHitFromIter(), and getDQMSummary::iter.

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,size_t> > 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( 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( *(getHitFromIter( iRecHit )), simTrackIdentifiers ); // This call fills simTrackIdentifiers
             // 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,size_t> >::iterator iIdentifierCountPair;
                 for( 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;
                         break;
                     }
                 }
                 if( iIdentifierCountPair == returnValue.end() ) returnValue.push_back( std::make_pair( *iIdentifier, 1 ) );
                 // This identifier wasn't found, so add it
             }
         }
     }
     return returnValue;
 }

template<typename iter >

int 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 448 of file QuickTrackAssociatorByHitsImpl.cc.

References TrackerHitAssociator::associateHitId(), spr::find(), getHitFromIter(), and getDQMSummary::iter.

Referenced by associateRecoToSim(), associateRecoToSimImplementation(), associateSimToReco(), 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.
 
     int doubleCount=0;
     std::vector < SimHitIdpr > SimTrackIdsDC;
 
     for( iter iHit=startIterator; iHit != endIterator; iHit++ )
     {
         int idcount=0;
 
         SimTrackIdsDC.clear();
         hitAssociator.associateHitId( *(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 ) doubleCount+=(idcount - 1);
     }
 
     return doubleCount;
 }

template<typename iter >

int QuickTrackAssociatorByHitsImpl::getDoubleCount	(	const ClusterTPAssociationList &	clusterToTPList,
		iter	begin,
		iter	end,
		TrackingParticleRef	associatedTrackingParticle
	)		const

private

Overload for when using cluster to TrackingParticle association list.

Definition at line 480 of file QuickTrackAssociatorByHitsImpl.cc.

References clusterTPAssociationListGreater(), getMatchedClusters(), and getDQMSummary::iter.

 {
     // This code here was written by Subir Sarkar. I'm just splitting it off into a
     // separate method. - Grimes 01/May/2014
 
     int doubleCount=0;
     std::vector < SimHitIdpr > SimTrackIdsDC;
 
     for( iter iHit=startIterator; iHit != endIterator; iHit++ )
     {
         int idcount=0;
 
         std::vector < OmniClusterRef > oClusters=getMatchedClusters( iHit, iHit + 1 );  //only for the cluster being checked
         for( std::vector<OmniClusterRef>::const_iterator it=oClusters.begin(); it != oClusters.end(); ++it )
         {
             std::pair<OmniClusterRef,TrackingParticleRef> clusterTPpairWithDummyTP( *it, TrackingParticleRef() ); //TP is dummy: for clusterTPAssociationListGreater sorting only the cluster is needed
             auto range=std::equal_range( clusterToTPList.begin(), clusterToTPList.end(), clusterTPpairWithDummyTP, clusterTPAssociationListGreater );
             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 ) doubleCount+=(idcount - 1);
     }
 
     return doubleCount;
 }

const TrackingRecHit* QuickTrackAssociatorByHitsImpl::getHitFromIter ( trackingRecHit_iterator iter ) const

inlineprivate

Definition at line 164 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by getAllSimTrackIdentifiers(), getDoubleCount(), and getMatchedClusters().

                                                                            {
     return &(**iter);
   }

const TrackingRecHit* QuickTrackAssociatorByHitsImpl::getHitFromIter ( TrackingRecHitCollection::const_iterator iter ) const

inlineprivate

Definition at line 168 of file QuickTrackAssociatorByHitsImpl.h.

                                                                                           {
     return &(*iter);
   }

template<typename iter >

std::vector< OmniClusterRef > QuickTrackAssociatorByHitsImpl::getMatchedClusters	(	iter	begin,
		iter	end
	)		const

private

Definition at line 341 of file QuickTrackAssociatorByHitsImpl.cc.

References SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), SiPixelRecHit::cluster(), compareJSON::const, end, TrackingRecHit::geographicalId(), getHitFromIter(), edm::Ref< C, T, F >::isNonnull(), TrackingRecHit::isValid(), getDQMSummary::iter, SiStripMatchedRecHit2D::monoClusterRef(), SiStripMatchedRecHit2D::monoHit(), TrackerSingleRecHit::omniClusterRef(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, SiStripMatchedRecHit2D::stereoClusterRef(), SiStripMatchedRecHit2D::stereoHit(), DetId::subdetId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, and SiStripDetId::TOB.

Referenced by associateTrack(), and getDoubleCount().

 {
   std::vector<OmniClusterRef> returnValue;
   for (iter iRecHit = begin; iRecHit != end; ++iRecHit) {
     const TrackingRecHit* rhit = getHitFromIter(iRecHit);
     if (rhit->isValid()) {
       int subdetid = rhit->geographicalId().subdetId();
       if (subdetid==PixelSubdetector::PixelBarrel||subdetid==PixelSubdetector::PixelEndcap) {
     const SiPixelRecHit* pRHit = dynamic_cast<const SiPixelRecHit*>(rhit);
     if (!pRHit->cluster().isNonnull())
       edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!" << " file: " << __FILE__ << " line: " << __LINE__;
     returnValue.push_back(pRHit->omniClusterRef());
       }
       else if (subdetid==SiStripDetId::TIB||subdetid==SiStripDetId::TOB||subdetid==SiStripDetId::TID||subdetid==SiStripDetId::TEC) {
     const std::type_info &tid = typeid(*rhit);
     if (tid == typeid(SiStripMatchedRecHit2D)) {
       const SiStripMatchedRecHit2D* sMatchedRHit = dynamic_cast<const SiStripMatchedRecHit2D*>(rhit);
       if (!sMatchedRHit->monoHit().cluster().isNonnull() || !sMatchedRHit->stereoHit().cluster().isNonnull())
         edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!" << " file: " << __FILE__ << " line: " << __LINE__;
       returnValue.push_back(sMatchedRHit->monoClusterRef());
       returnValue.push_back(sMatchedRHit->stereoClusterRef());
     }
     else if (tid == typeid(SiStripRecHit2D)) {
       const SiStripRecHit2D* sRHit = dynamic_cast<const SiStripRecHit2D*>(rhit);
       if (!sRHit->cluster().isNonnull())
         edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!" << " file: " << __FILE__ << " line: " << __LINE__;
       returnValue.push_back(sRHit->omniClusterRef());
     }
     else if (tid == typeid(SiStripRecHit1D)) {
       const SiStripRecHit1D* sRHit = dynamic_cast<const SiStripRecHit1D*>(rhit);
       if (!sRHit->cluster().isNonnull())
         edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!" << " file: " << __FILE__ << " line: " << __LINE__;
       returnValue.push_back(sRHit->omniClusterRef());
     }
     else {
       auto const & thit = static_cast<BaseTrackerRecHit const&>(*rhit);
       if ( thit.isProjected() ) {
       } else {
         edm::LogError("TrackAssociator") << ">>> getMatchedClusters: TrackingRecHit not associated to any SiStripCluster! subdetid = " << subdetid;
       }
     }
       }
       else {
     edm::LogError("TrackAssociator") << ">>> getMatchedClusters: TrackingRecHit not associated to any cluster! subdetid = " << subdetid;
       }
     }
   }
   return returnValue;
 }

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

inlinestaticprivate

Definition at line 104 of file QuickTrackAssociatorByHitsImpl.h.

References i, and j.

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

i

int i

Definition: DBlmapReader.cc:9

j

int j

Definition: DBlmapReader.cc:9

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 430 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;
 }

Member Data Documentation

bool QuickTrackAssociatorByHitsImpl::absoluteNumberOfHits_

private

Definition at line 195 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSim(), associateRecoToSimImplementation(), associateSimToReco(), and associateSimToRecoImplementation().

std::shared_ptr<const ClusterTPAssociationList> QuickTrackAssociatorByHitsImpl::clusterToTPMap_

private

Definition at line 188 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSim(), and associateSimToReco().

double QuickTrackAssociatorByHitsImpl::cutRecoToSim_

private

Definition at line 192 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSim(), and associateRecoToSimImplementation().

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

private

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

A decision is made whether to create a ClusterTPAssociationList or a TrackerHitAssociator depending on how this track associator was configured. If the ClusterTPAssociationList 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 187 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateRecoToSim(), and associateSimToReco().

double QuickTrackAssociatorByHitsImpl::puritySimToReco_

private

Definition at line 191 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateSimToReco(), and associateSimToRecoImplementation().

double QuickTrackAssociatorByHitsImpl::qualitySimToReco_

private

Definition at line 190 of file QuickTrackAssociatorByHitsImpl.h.

Referenced by associateSimToReco(), and associateSimToRecoImplementation().

SimToRecoDenomType QuickTrackAssociatorByHitsImpl::simToRecoDenominator_

private