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

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

{}

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 310 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 112 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 99 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 48 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 257 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::findTrackingParticlePtrs(), edm::Ptr< T >::get(), TTClusterAssociationMap< T >::isGenuine(), edm::Ptr< T >::isNull(), TTClusterAssociationMap< T >::isUnknown(), jetUpdater_cfi::sort, 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 148 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::findTrackingParticlePtrs(), edm::Ptr< T >::get(), edm::Ptr< T >::isNull(), TrackingParticle::p4(), jetUpdater_cfi::sort, 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 )
  {
    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 219 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::findTrackingParticlePtrs(), edm::Ptr< T >::get(), edm::Ptr< T >::isNull(), jetUpdater_cfi::sort, 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 53 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 51 of file TTClusterAssociationMap.h.

References TTClusterAssociationMap< T >::clusterToTrackingParticleVectorMap.

      { clusterToTrackingParticleVectorMap = aMap; }

Member Data Documentation

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

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

template<typename T>

int TTClusterAssociationMap< T >::nclus

private

Definition at line 71 of file TTClusterAssociationMap.h.

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

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

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