TTClusterAssociationMap< T > Class Template Reference

Stores association of Truth Particles (TP) to L1 Track-Trigger Clusters. More...

#include <TTClusterAssociationMap.h>

Public Member Functions
const TrackingParticlePtr &	findTrackingParticlePtr (TTClusterRefT< T > aCluster) const
	Get main TP associated to a cluster. (Non-NULL if isGenuine() below is true). More...
const std::vector< TrackingParticlePtr > &	findTrackingParticlePtrs (TTClusterRefT< T > aCluster) const
	Get all TPs associated to a cluster. More...
const std::vector< TTClusterRefT< T > > &	findTTClusterRefs (TrackingParticlePtr aTrackingParticle) const
	Operations. More...
const MapTPToVecClus< T > &	getTrackingParticleToTTClustersMap () const
const MapClusToVecTP< T > &	getTTClusterToTrackingParticlesMap () const
	Get/set cluster <-> truth association maps. More...
bool	isCombinatoric (TTClusterRefT< T > aCluster) const
	Cluster is not "genuine" or "unknown". More...
bool	isGenuine (TTClusterRefT< T > aCluster) const
bool	isUnknown (TTClusterRefT< T > aCluster) const
	Cluster "unknown": i.e. not associated with any TP. More...
void	setTrackingParticleToTTClustersMap (const MapTPToVecClus< T > &aMap)
void	setTTClusterToTrackingParticlesMap (const MapClusToVecTP< T > &aMap)
	TTClusterAssociationMap ()
	Constructors. More...
	~TTClusterAssociationMap ()
	Destructor. More...

Private Attributes
MapClusToVecTP< T >	clusterToTrackingParticleVectorMap_
	Data members. More...
int	nclus
MapTPToVecClus< T >	trackingParticleToClusterVectorMap_

Static Private Attributes
static const TrackingParticlePtr	nullTrackingParticlePtr_
	Close class. More...
static const std::vector< TTClusterRefT< T > >	nullVecClusterRef_
static const std::vector< TrackingParticlePtr >	nullVecTrackingParticlePtr_

Detailed Description

template<typename T>
class TTClusterAssociationMap< T >

Stores association of Truth Particles (TP) to L1 Track-Trigger Clusters.

Contains two maps. One associates each cluster to a vector of all TPs that made its hits. The other associates each TP to a vector of all clusters it contributed to.

(The template structure is used to accomodate types other than PixelDigis, in case they are needed in future).

Author

Nicola Pozzobon

Date

2013, Jul 19 (tidy up: Ian Tomalin, 2020)

Definition at line 44 of file TTClusterAssociationMap.h.

Constructor & Destructor Documentation

TTClusterAssociationMap()

template<typename T >

TTClusterAssociationMap< T >::TTClusterAssociationMap

(

)

Constructors.

Default Constructor NOTE: to be used with setSomething(...) methods

Set default data members

Definition at line 113 of file TTClusterAssociationMap.h.

                                                    {
  nclus = 0;
}

~TTClusterAssociationMap()

template<typename T >

TTClusterAssociationMap< T >::~TTClusterAssociationMap

(

)

Destructor.

Definition at line 120 of file TTClusterAssociationMap.h.

{}

Member Function Documentation

findTrackingParticlePtr()

template<typename T >

const TrackingParticlePtr & TTClusterAssociationMap< T >::findTrackingParticlePtr

(

TTClusterRefT< T >

aCluster

)

const

Get main TP associated to a cluster. (Non-NULL if isGenuine() below is true).

Definition at line 144 of file TTClusterAssociationMap.h.

                                                                                                              {
  if (this->isGenuine(aCluster)) {
    return this->findTrackingParticlePtrs(aCluster).at(0);
  } else {
    return nullTrackingParticlePtr_;
  }
}

findTrackingParticlePtrs()

template<typename T >

const std::vector< TrackingParticlePtr > & TTClusterAssociationMap< T >::findTrackingParticlePtrs

(

TTClusterRefT< T >

aCluster

)

const

Get all TPs associated to a cluster.

Definition at line 134 of file TTClusterAssociationMap.h.

                                     {
  if (clusterToTrackingParticleVectorMap_.find(aCluster) != clusterToTrackingParticleVectorMap_.end()) {
    return clusterToTrackingParticleVectorMap_.find(aCluster)->second;
  } else {
    return nullVecTrackingParticlePtr_;
  }
}

findTTClusterRefs()

template<typename T >

const std::vector< TTClusterRefT< T > > & TTClusterAssociationMap< T >::findTTClusterRefs

(

TrackingParticlePtr

aTrackingParticle

)

const

Operations.

Definition at line 124 of file TTClusterAssociationMap.h.

                                                 {
  if (trackingParticleToClusterVectorMap_.find(aTrackingParticle) != trackingParticleToClusterVectorMap_.end()) {
    return trackingParticleToClusterVectorMap_.find(aTrackingParticle)->second;
  } else {
    return nullVecClusterRef_;
  }
}

getTrackingParticleToTTClustersMap()

template<typename T >

const MapTPToVecClus<T>& TTClusterAssociationMap< T >::getTrackingParticleToTTClustersMap ( ) const

inline

Definition at line 55 of file TTClusterAssociationMap.h.

{ return trackingParticleToClusterVectorMap_; }

References TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap_.

getTTClusterToTrackingParticlesMap()

template<typename T >

const MapClusToVecTP<T>& TTClusterAssociationMap< T >::getTTClusterToTrackingParticlesMap ( ) const

inline

Get/set cluster <-> truth association maps.

Definition at line 54 of file TTClusterAssociationMap.h.

{ return clusterToTrackingParticleVectorMap_; }

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap_.

isCombinatoric()

template<typename T >

bool TTClusterAssociationMap< T >::isCombinatoric

(

TTClusterRefT< T >

aCluster

)

const

Cluster is not "genuine" or "unknown".

Get the TrackingParticles

If the vector is empty, then the cluster is UNKNOWN

If we are here, it means there are some TrackingParticles

Loop over the TrackingParticles

Get the TrackingParticle

Count the NULL TrackingParticles

Store the pointers (addresses) of the TrackingParticle to be able to count how many different there are

Count how many different TrackingParticle there are

COMBINATORIC means no NULLs and more than one good TP OR, in alternative, only one good TP but non-zero NULLS

Definition at line 270 of file TTClusterAssociationMap.h.

                                                                               {
  const std::vector<TrackingParticlePtr>& theseTrackingParticles = this->findTrackingParticlePtrs(aCluster);
 
  if (theseTrackingParticles.empty())
    return false;
 
  bool genuineClu = this->isGenuine(aCluster);
  bool unknownClu = this->isUnknown(aCluster);
 
  if (genuineClu || unknownClu)
    return false;
 
  return true;
 
  unsigned int nullTPs = 0;
  unsigned int goodDifferentTPs = 0;
  std::vector<const TrackingParticle*> tpAddressVector;
 
  for (unsigned int itp = 0; itp < theseTrackingParticles.size(); itp++) {
    TrackingParticlePtr curTP = theseTrackingParticles.at(itp);
 
    if (curTP.isNull()) {
      nullTPs++;
    } else {
      tpAddressVector.push_back(curTP.get());
    }
  }
 
  std::sort(tpAddressVector.begin(), tpAddressVector.end());
  tpAddressVector.erase(std::unique(tpAddressVector.begin(), tpAddressVector.end()), tpAddressVector.end());
  goodDifferentTPs = tpAddressVector.size();
 
  //return ( ( nullTPs == 0 && goodDifferentTPs > 1 ) || ( nullTPs > 0 && goodDifferentTPs > 0 ) );
  return (goodDifferentTPs > 1);
}

References edm::Ptr< T >::get(), edm::Ptr< T >::isNull(), jetUpdater_cfi::sort, and tier0::unique().

isGenuine()

template<typename T >

bool TTClusterAssociationMap< T >::isGenuine

(

TTClusterRefT< T >

aCluster

)

const

— Get quality of L1 cluster based on truth info. (exactly 1 of following 3 functions is always true) Cluster "genuine": i.e. cluster associated to exactly 1 TP. (If other TPs are associated, but have in total < 1% of Pt of main TP, or if they are null, then they are neglected here).

MC truth Table to define Genuine, Combinatoric and Unknown

N = number of NULL TP pointers D = number of GOOD TP pointers different from each other

OLD DEFINITION

N / D--> | 0 | 1 | >1

0 | U | G | C

>0 | U | C | C NEW DEFINITION SV 060617

N / D--> | 0 | 1 | >1 (with 1 TP getting >99% of the total pT) | >1

0 | U | G | G | C

>0 | U | G | G | C

Get the TrackingParticles

If the vector is empty, then the cluster is UNKNOWN

If we are here, it means there are some TrackingParticles

Loop over the TrackingParticles

Get the TrackingParticle

Count the NULL TrackingParticles

Get the TrackingParticle

Count the NULL TrackingParticles

Store the pointers (addresses) of the TrackingParticle to be able to count how many different there are

Count how many different TrackingParticle there are

Definition at line 177 of file TTClusterAssociationMap.h.

                                                                          {
  const std::vector<TrackingParticlePtr>& theseTrackingParticles = this->findTrackingParticlePtrs(aCluster);
 
  if (theseTrackingParticles.empty())
    return false;
 
  unsigned int nullTPs = 0;
  unsigned int goodDifferentTPs = 0;
  std::vector<const TrackingParticle*> tpAddressVector;
 
  std::vector<float> tp_mom;
 
  float tp_tot = 0;
 
  for (const auto& tp : theseTrackingParticles) {
    const TrackingParticlePtr& curTP = tp;
 
    if (curTP.isNull()) {
      //      nullTPs++;
      tp_mom.push_back(0);
    } else {
      tp_mom.push_back(curTP.get()->p4().pt());
      tp_tot += curTP.get()->p4().pt();
    }
  }
 
  if (tp_tot == 0)
    return false;
 
  for (unsigned int itp = 0; itp < theseTrackingParticles.size(); itp++) {
    TrackingParticlePtr curTP = theseTrackingParticles.at(itp);
 
    if (tp_mom.at(itp) <= 0.01 * tp_tot) {
      nullTPs++;
    } else {
      tpAddressVector.push_back(curTP.get());
    }
  }
 
  std::sort(tpAddressVector.begin(), tpAddressVector.end());
  tpAddressVector.erase(std::unique(tpAddressVector.begin(), tpAddressVector.end()), tpAddressVector.end());
  goodDifferentTPs = tpAddressVector.size();
 
  return (goodDifferentTPs == 1);
}

References edm::Ptr< T >::get(), edm::Ptr< T >::isNull(), TrackingParticle::p4(), jetUpdater_cfi::sort, cmsswSequenceInfo::tp, and tier0::unique().

isUnknown()

template<typename T >

bool TTClusterAssociationMap< T >::isUnknown

(

TTClusterRefT< T >

aCluster

)

const

Cluster "unknown": i.e. not associated with any TP.

Get the TrackingParticles

If the vector is empty, then the cluster is UNKNOWN

If we are here, it means there are some TrackingParticles

Loop over the TrackingParticles

Get the TrackingParticle

Count the non-NULL TrackingParticles

Store the pointers (addresses) of the TrackingParticle to be able to count how many different there are

Count how many different TrackingParticle there are

UNKNOWN means no good TP is found

Definition at line 235 of file TTClusterAssociationMap.h.

                                                                          {
  const std::vector<TrackingParticlePtr>& theseTrackingParticles = this->findTrackingParticlePtrs(aCluster);
 
  if (theseTrackingParticles.empty())
    return true;
 
  unsigned int goodDifferentTPs = 0;
  std::vector<const TrackingParticle*> tpAddressVector;
 
  for (unsigned int itp = 0; itp < theseTrackingParticles.size(); itp++) {
    TrackingParticlePtr curTP = theseTrackingParticles.at(itp);
 
    if (!curTP.isNull()) {
      tpAddressVector.push_back(curTP.get());
    }
  }
 
  std::sort(tpAddressVector.begin(), tpAddressVector.end());
  tpAddressVector.erase(std::unique(tpAddressVector.begin(), tpAddressVector.end()), tpAddressVector.end());
  goodDifferentTPs = tpAddressVector.size();
 
  return (goodDifferentTPs == 0);
}

References edm::Ptr< T >::get(), edm::Ptr< T >::isNull(), jetUpdater_cfi::sort, and tier0::unique().

setTrackingParticleToTTClustersMap()

template<typename T >

void TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap ( const MapTPToVecClus< T > &  aMap )

inline

Definition at line 58 of file TTClusterAssociationMap.h.

{ trackingParticleToClusterVectorMap_ = aMap; }

References TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap_.

setTTClusterToTrackingParticlesMap()

template<typename T >

void TTClusterAssociationMap< T >::setTTClusterToTrackingParticlesMap ( const MapClusToVecTP< T > &  aMap )

inline

Definition at line 57 of file TTClusterAssociationMap.h.

{ clusterToTrackingParticleVectorMap_ = aMap; }

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap_.

Member Data Documentation

clusterToTrackingParticleVectorMap_

template<typename T >

MapClusToVecTP<T> TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap_

private

Data members.

Definition at line 83 of file TTClusterAssociationMap.h.

Referenced by TTClusterAssociationMap< T >::getTTClusterToTrackingParticlesMap(), and TTClusterAssociationMap< T >::setTTClusterToTrackingParticlesMap().

nclus

template<typename T >

int TTClusterAssociationMap< T >::nclus

private

Definition at line 86 of file TTClusterAssociationMap.h.

nullTrackingParticlePtr_

template<typename T >

const TrackingParticlePtr TTClusterAssociationMap< T >::nullTrackingParticlePtr_

staticprivate

Close class.

Implementation of methods

Here, in the header file, the methods which do not depend on the specific type <T> that can fit the template. Other methods, with type-specific features, are implemented in the source file.

Definition at line 89 of file TTClusterAssociationMap.h.

nullVecClusterRef_

template<typename T >

const std::vector< TTClusterRefT< T > > TTClusterAssociationMap< T >::nullVecClusterRef_

staticprivate

Definition at line 91 of file TTClusterAssociationMap.h.

nullVecTrackingParticlePtr_

template<typename T >

const std::vector< TrackingParticlePtr > TTClusterAssociationMap< T >::nullVecTrackingParticlePtr_

staticprivate

Definition at line 90 of file TTClusterAssociationMap.h.

trackingParticleToClusterVectorMap_

template<typename T >

MapTPToVecClus<T> TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap_

private

Definition at line 84 of file TTClusterAssociationMap.h.

Referenced by TTClusterAssociationMap< T >::getTrackingParticleToTTClustersMap(), and TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap().