TTClusterAssociationMap< T > Class Template Reference

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

TTClusterAssociationMap()

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.

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

~TTClusterAssociationMap()

template<typename T >

TTClusterAssociationMap< T >::~TTClusterAssociationMap

(

)

Destructor.

Definition at line 107 of file TTClusterAssociationMap.h.

{}

Member Function Documentation

findTrackingParticlePtr()

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.

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

References groupFilesInBlocks::temp.

findTrackingParticlePtrs()

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.

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

findTTClusterRefs()

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.

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

getTrackingParticleToTTClustersMap()

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.

                                             {
    return trackingParticleToClusterVectorMap;
  }

References TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap.

getTTClusterToTrackingParticlesMap()

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.

                                             {
    return clusterToTrackingParticleVectorMap;
  }

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap.

isCombinatoric()

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.

                                                                             {
  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);
}

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

isGenuine()

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 | C NEW 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.

                                                                                                                  {
  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);
}

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

isUnknown()

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.

                                                                                                                  {
  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);
}

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

setTrackingParticleToTTClustersMap()

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.

                {
    trackingParticleToClusterVectorMap = aMap;
  }

References TTClusterAssociationMap< T >::trackingParticleToClusterVectorMap.

setTTClusterToTrackingParticlesMap()

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.

                                                                                                {
    clusterToTrackingParticleVectorMap = aMap;
  }

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap.

Member Data Documentation

clusterToTrackingParticleVectorMap

template<typename T >

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

private

Data members.

Definition at line 80 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 84 of file TTClusterAssociationMap.h.

trackingParticleToClusterVectorMap

template<typename T >

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

private

Definition at line 82 of file TTClusterAssociationMap.h.

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