NOTE: this is needed even if it seems not. More...

#include <TTClusterAssociationMap.h>

Public Member Functions
edm::Ptr< TrackingParticle >	findTrackingParticlePtr (edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster) const

std::vector< edm::Ptr< TrackingParticle > >	findTrackingParticlePtrs (edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster) const
	Operations. More...

std::vector< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > >	findTTClusterRefs (edm::Ptr< TrackingParticle > aTrackingParticle) const
	Operations. More...

std::map< edm::Ptr< TrackingParticle >, std::vector< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > > >	getTrackingParticleToTTClustersMap () const

std::map< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > >, std::vector< edm::Ptr< TrackingParticle > > >	getTTClusterToTrackingParticlesMap () const
	Maps. More...

bool	isCombinatoric (edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster) const

bool	isGenuine (edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster) const
	MC Truth methods. More...

bool	isUnknown (edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster) const

void	setTrackingParticleToTTClustersMap (std::map< edm::Ptr< TrackingParticle >, std::vector< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > > > aMap)

void	setTTClusterToTrackingParticlesMap (std::map< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > >, std::vector< edm::Ptr< TrackingParticle > > > aMap)

	TTClusterAssociationMap ()
	Constructors. More...

	~TTClusterAssociationMap ()
	Destructor. More...

Private Attributes
std::map< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > >, std::vector< edm::Ptr< TrackingParticle > > >	clusterToTrackingParticleVectorMap
	Data members. More...

int	nclus

std::map< edm::Ptr< TrackingParticle >, std::vector< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > > >	trackingParticleToClusterVectorMap

Detailed Description

template<typename T>
class TTClusterAssociationMap< T >

NOTE: this is needed even if it seems not.

Class to store the MC truth of L1 Track Trigger clusters.

After moving from SimDataFormats to DataFormats, the template structure of the class was maintained in order to accomodate any types other than PixelDigis in case there is such a need in the future.

Author: Nicola Pozzobon

Date: 2013, Jul 19

Definition at line 33 of file TTClusterAssociationMap.h.

Constructor & Destructor Documentation

template<typename T >

TTClusterAssociationMap< T >::TTClusterAssociationMap ( )

Constructors.

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.Default Constructor NOTE: to be used with setSomething(...) methods

Set default data members

Definition at line 98 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap, TTClusterAssociationMap< T >::nclus, and TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap.

                                                     {
   clusterToTrackingParticleVectorMap.clear();
   trackingParticleToClusterVectorMap.clear();
   nclus = 0;
 }

template<typename T >

TTClusterAssociationMap< T >::~TTClusterAssociationMap ( )

Destructor.

Definition at line 107 of file TTClusterAssociationMap.h.

107 {}

Member Function Documentation

template<typename T >

edm::Ptr< TrackingParticle > TTClusterAssociationMap< T >::findTrackingParticlePtr ( edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster ) const

Definition at line 299 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::findTrackingParticlePtrs(), TTClusterAssociationMap< T >::isGenuine(), and groupFilesInBlocks::temp.

Referenced by TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap().

                                                                              {
   if (this->isGenuine(aCluster)) {
     return this->findTrackingParticlePtrs(aCluster).at(0);
   }
 
   edm::Ptr<TrackingParticle>* temp = new edm::Ptr<TrackingParticle>();
   return *temp;
 }

template<typename T >

std::vector< edm::Ptr< TrackingParticle > > TTClusterAssociationMap< T >::findTrackingParticlePtrs ( edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster ) const

Operations.

Definition at line 123 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap.

Referenced by TTClusterAssociationMap< T >::findTrackingParticlePtr(), TTClusterAssociationMap< T >::isCombinatoric(), TTClusterAssociationMap< T >::isGenuine(), TTClusterAssociationMap< T >::isUnknown(), and TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap().

                                                                              {
   if (clusterToTrackingParticleVectorMap.find(aCluster) != clusterToTrackingParticleVectorMap.end()) {
     return clusterToTrackingParticleVectorMap.find(aCluster)->second;
   }
 
   std::vector<edm::Ptr<TrackingParticle> > tempVector;
   tempVector.clear();
   return tempVector;
 }

template<typename T >

std::vector< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > > TTClusterAssociationMap< T >::findTTClusterRefs ( edm::Ptr< TrackingParticle > aTrackingParticle ) const

Operations.

Definition at line 112 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap.

Referenced by TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap().

                                                                                               {
   if (trackingParticleToClusterVectorMap.find(aTrackingParticle) != trackingParticleToClusterVectorMap.end()) {
     return trackingParticleToClusterVectorMap.find(aTrackingParticle)->second;
   }
 
   std::vector<edm::Ref<edmNew::DetSetVector<TTCluster<T> >, TTCluster<T> > > tempVector;
   tempVector.clear();
   return tempVector;
 }

template<typename T>

std::map<edm::Ptr<TrackingParticle>, std::vector<edm::Ref<edmNew::DetSetVector<TTCluster<T> >, TTCluster<T> > > > TTClusterAssociationMap< T >::getTrackingParticleToTTClustersMap ( ) const

inline

Definition at line 50 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap.

                                              {
     return trackingParticleToClusterVectorMap;
   }

template<typename T>

std::map<edm::Ref<edmNew::DetSetVector<TTCluster<T> >, TTCluster<T> >, std::vector<edm::Ptr<TrackingParticle> > > TTClusterAssociationMap< T >::getTTClusterToTrackingParticlesMap ( ) const

inline

Maps.

Data members: getABC( ... ) Helper methods: findABC( ... )

Definition at line 46 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap.

                                              {
     return clusterToTrackingParticleVectorMap;
   }

template<typename T >

bool TTClusterAssociationMap< T >::isCombinatoric ( edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster ) const

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 250 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::findTrackingParticlePtrs(), edm::Ptr< T >::get(), TTClusterAssociationMap< T >::isGenuine(), edm::Ptr< T >::isNull(), TTClusterAssociationMap< T >::isUnknown(), and tier0::unique().

Referenced by TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap().

                                                                              {
   std::vector<edm::Ptr<TrackingParticle> > 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++) {
     edm::Ptr<TrackingParticle> 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);
 }

template<typename T >

bool TTClusterAssociationMap< T >::isGenuine ( edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster ) const

MC Truth methods.

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

N / D–> | 0 | 1 | >1

0 | U | G | C

>0 | U | C | CNEW 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 157 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::findTrackingParticlePtrs(), edm::Ptr< T >::get(), edm::Ptr< T >::isNull(), TrackingParticle::p4(), and tier0::unique().

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

                                                                                                                   {
   std::vector<edm::Ptr<TrackingParticle> > 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 edm::Ptr<TrackingParticle>& 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++) {
     edm::Ptr<TrackingParticle> 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);
 }

template<typename T >

bool TTClusterAssociationMap< T >::isUnknown ( edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > aCluster ) const

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 215 of file TTClusterAssociationMap.h.

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

Referenced by TTClusterAssociationMap< T >::isCombinatoric(), TTStubAssociationMap< T >::isUnknown(), and TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap().

                                                                                                                   {
   std::vector<edm::Ptr<TrackingParticle> > 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++) {
     edm::Ptr<TrackingParticle> 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);
 }

template<typename T>

void TTClusterAssociationMap< T >::setTrackingParticleToTTClustersMap ( std::map< edm::Ptr< TrackingParticle >, std::vector< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > > > > aMap )

inline

Definition at line 58 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::findTrackingParticlePtr(), TTClusterAssociationMap< T >::findTrackingParticlePtrs(), TTClusterAssociationMap< T >::findTTClusterRefs(), TTClusterAssociationMap< T >::isCombinatoric(), TTClusterAssociationMap< T >::isGenuine(), TTClusterAssociationMap< T >::isUnknown(), and TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap.

                 {
     trackingParticleToClusterVectorMap = aMap;
   }

template<typename T>

void TTClusterAssociationMap< T >::setTTClusterToTrackingParticlesMap ( std::map< edm::Ref< edmNew::DetSetVector< TTCluster< T > >, TTCluster< T > >, std::vector< edm::Ptr< TrackingParticle > > > aMap )

inline

Definition at line 54 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap.

                                                                                                 {
     clusterToTrackingParticleVectorMap = aMap;
   }