Trk::TrackDensitySeedFinder Class Reference

#include <TrackDensitySeedFinder.h>

Inheritance diagram for Trk::TrackDensitySeedFinder:

Collaboration diagram for Trk::TrackDensitySeedFinder:

Public Member Functions
	TrackDensitySeedFinder (const std::string &t, const std::string &n, const IInterface *p)
virtual	~TrackDensitySeedFinder ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual Amg::Vector3D	findSeed (const std::vector< const Trk::Track * > &vectorTrk, const xAOD::Vertex *constraint=0) const override final
	Finds a linearization point out of a vector of tracks and returns it as an Amg::Vector3D object.
virtual Amg::Vector3D	findSeed (const std::vector< const Trk::TrackParameters * > &perigeeList, const xAOD::Vertex *constraint=0) const override final
	Finds a linearization point out of a vector of TrackParameters and returns it as an Amg::Vector3D object.
virtual std::pair< Amg::Vector3D, Amg::MatrixX >	findAnalyticSeed (const std::vector< const Trk::TrackParameters * > &perigeeList, const xAOD::Vertex *constraint=0) const override final
virtual std::vector< Amg::Vector3D >	findMultiSeeds (const std::vector< const Trk::Track * > &vectorTrk, const xAOD::Vertex *constraint=0) const override final
	Finds full vector of linearization points from a vector of tracks and returns it as an Amg::Vector3D object.
virtual std::vector< Amg::Vector3D >	findMultiSeeds (const std::vector< const Trk::TrackParameters * > &perigeeList, const xAOD::Vertex *constraint=0) const override final
	Finds full vector of linearization points from a vector of TrackParameters and returns it as an Amg::Vector3D object.
virtual Amg::Vector3D	findSeed (const double, const double, const std::vector< const Trk::TrackParameters * > &perigeeList, const xAOD::Vertex *constraint=0) const
	Finds a linearization point out of a vector of TrackParameters and returns it as an Amg::Vector3D object.
virtual Amg::Vector3D	findSeed (const double, const double, std::unique_ptr< Trk::IMode3dInfo > &, const std::vector< const Trk::TrackParameters * > &perigeeList, const xAOD::Vertex *constraint=0) const
	Finds a linearization point out of a vector of TrackParameters and returns it as an Amg::Vector3D object.

Private Attributes
ToolHandle< IVertexTrackDensityEstimator >	m_densityEstimator

Detailed Description

Definition at line 29 of file TrackDensitySeedFinder.h.

Constructor & Destructor Documentation

TrackDensitySeedFinder()

Trk::TrackDensitySeedFinder::TrackDensitySeedFinder	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Definition at line 25 of file TrackDensitySeedFinder.cxx.

                                                                                                             : 
    base_class(t,n,p)
  {   
  }

~TrackDensitySeedFinder()

Trk::TrackDensitySeedFinder::~TrackDensitySeedFinder ( )

virtualdefault

Member Function Documentation

finalize()

StatusCode Trk::TrackDensitySeedFinder::finalize ( )

overridevirtual

Definition at line 42 of file TrackDensitySeedFinder.cxx.

  {
    ATH_MSG_DEBUG("Finalize successful");
    return StatusCode::SUCCESS;
  }

findAnalyticSeed()

std::pair< Amg::Vector3D, Amg::MatrixX > Trk::TrackDensitySeedFinder::findAnalyticSeed ( const std::vector< const Trk::TrackParameters * > & perigeeList, const xAOD::Vertex * constraint = 0 ) const

finaloverridevirtual

Definition at line 93 of file TrackDensitySeedFinder.cxx.

                                                                                                                                                                       {
      std::pair<double,double> zResult {0.,0.};
    if ( !perigeeList.empty() ) 
    {
      zResult = m_densityEstimator->globalMaximumWithWidth (perigeeList);
      ATH_MSG_DEBUG("Highest density Z position found: " << zResult.first << "with width: " << zResult.second);
    }
    else
    {
      ATH_MSG_DEBUG("No tracks with sufficient weight; return z position = 0");
    }
 
    if (constraint)
    {
      Amg::Vector3D positionVector = Amg::Vector3D(constraint->position().x(), constraint->position().y(), zResult.first + constraint->position().z());
      Amg::MatrixX covarianceMatrix = constraint->covariancePosition() ;
      
      ATH_MSG_DEBUG("The vertex seed width is " << zResult.second);
     
     // if no seed is found, we get a NaN here. 
     // This should not propagate through the vertex fit.
     // The upstream code aborts if position[z] == 0, which we need to ensure.
     // In this case, the cov matrix will not be used. 
      if(std::isnan(zResult.second)){
        positionVector[2] = 0.; 
      }
      else covarianceMatrix.fillSymmetric(2,2,std::pow(zResult.second,2.));
 
      return std::make_pair(positionVector,covarianceMatrix) ;
    }
    else
    {
      Amg::Vector3D positionVector = Amg::Vector3D(0.,0.,zResult.first);
      Amg::MatrixX covarianceMatrix;
      for(const auto i : {0,1,2} ){
        for(const auto j : {0,1,2}){
          covarianceMatrix(i,j)=0;
        }
      }
      return {positionVector,covarianceMatrix};
    }  
}

findMultiSeeds() [1/2]

std::vector< Amg::Vector3D > Trk::TrackDensitySeedFinder::findMultiSeeds ( const std::vector< const Trk::Track * > & vectorTrk, const xAOD::Vertex * constraint = 0 ) const

finaloverridevirtual

Finds full vector of linearization points from a vector of tracks and returns it as an Amg::Vector3D object.

Intended for seed finders that produce all at once. If you want an additional constraint can be taken into account.

Definition at line 136 of file TrackDensitySeedFinder.cxx.

  {
    //implemented to satisfy inheritance but this algorithm only supports one seed at a time
    ATH_MSG_WARNING("Multi-seeding requested but seed finder not able to operate in that mode, returning no seeds");
    return std::vector<Amg::Vector3D>(0);
  }

findMultiSeeds() [2/2]

std::vector< Amg::Vector3D > Trk::TrackDensitySeedFinder::findMultiSeeds ( const std::vector< const Trk::TrackParameters * > & perigeeList, const xAOD::Vertex * constraint = 0 ) const

finaloverridevirtual

Finds full vector of linearization points from a vector of TrackParameters and returns it as an Amg::Vector3D object.

Intended for seed finders that produce all at once. If you want an additional constraint can be taken into account.

Definition at line 144 of file TrackDensitySeedFinder.cxx.

  {
     //implemented to satisfy inheritance but this algorithm only supports one seed at a time
    ATH_MSG_WARNING("Multi-seeding requested but seed finder not able to operate in that mode, returning no seeds");
    return std::vector<Amg::Vector3D>(0);
  }

findSeed() [1/4]

virtual Amg::Vector3D Trk::IVertexSeedFinder::findSeed ( const double , const double , const std::vector< const Trk::TrackParameters * > & perigeeList, const xAOD::Vertex * constraint = 0 ) const

inline

Finds a linearization point out of a vector of TrackParameters and returns it as an Amg::Vector3D object.

If you want an additional constraint can be taken into account. Must specify the primary vertex position.

Definition at line 79 of file IVertexSeedFinder.h.

    {
      return findSeed (perigeeList, constraint);
    }

findSeed() [2/4]

virtual Amg::Vector3D Trk::IVertexSeedFinder::findSeed ( const double , const double , std::unique_ptr< Trk::IMode3dInfo > & , const std::vector< const Trk::TrackParameters * > & perigeeList, const xAOD::Vertex * constraint = 0 ) const

inline

Finds a linearization point out of a vector of TrackParameters and returns it as an Amg::Vector3D object.

If you want an additional constraint can be taken into account. Must specify the primary vertex position. Allows returning additional diagnostic information.

Definition at line 96 of file IVertexSeedFinder.h.

    {
      return findSeed (perigeeList, constraint);
    }

findSeed() [3/4]

Amg::Vector3D Trk::TrackDensitySeedFinder::findSeed ( const std::vector< const Trk::Track * > & vectorTrk, const xAOD::Vertex * constraint = 0 ) const

finaloverridevirtual

Finds a linearization point out of a vector of tracks and returns it as an Amg::Vector3D object.

If you want an additional constraint can be taken into account.

Definition at line 48 of file TrackDensitySeedFinder.cxx.

  {
    
    //create perigees from track list
    std::vector<const TrackParameters*> perigeeList;
    
    std::vector<const Trk::Track*>::const_iterator begin=VectorTrk.begin();
    std::vector<const Trk::Track*>::const_iterator end=VectorTrk.end();
    
    for (std::vector<const Trk::Track*>::const_iterator iter=begin;
     iter!=end;++iter) 
    {
      perigeeList.push_back((*iter)->perigeeParameters());
    }
   
    //create seed from perigee list
    return findSeed(perigeeList,constraint);  
  }

findSeed() [4/4]

Amg::Vector3D Trk::TrackDensitySeedFinder::findSeed ( const std::vector< const Trk::TrackParameters * > & perigeeList, const xAOD::Vertex * constraint = 0 ) const

finaloverridevirtual

Finds a linearization point out of a vector of TrackParameters and returns it as an Amg::Vector3D object.

If you want an additional constraint can be taken into account.

Definition at line 68 of file TrackDensitySeedFinder.cxx.

  {
    double zResult {0.};
    if ( !perigeeList.empty() ) 
    {
      zResult = m_densityEstimator->globalMaximum (perigeeList);
      ATH_MSG_DEBUG("Highest density Z position found: " << zResult);
    }
    else
    {
      ATH_MSG_DEBUG("No tracks with sufficient weight; return z position = 0");
    }
 
    if (constraint)
    {
      return Amg::Vector3D(constraint->position().x(), constraint->position().y(), zResult + constraint->position().z());
    }
    else
    {
      return Amg::Vector3D(0.,0.,zResult);
    }  
  }

initialize()

StatusCode Trk::TrackDensitySeedFinder::initialize ( )

overridevirtual

Definition at line 35 of file TrackDensitySeedFinder.cxx.

  { 
    ATH_CHECK( m_densityEstimator.retrieve() );
    ATH_MSG_DEBUG("Initialize successful");
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_densityEstimator

ToolHandle<IVertexTrackDensityEstimator> Trk::TrackDensitySeedFinder::m_densityEstimator

private

Initial value:

{
      this,
      "DensityEstimator",
      "Trk::GaussianTrackDensity",
      "Track density tool"
    }

Definition at line 92 of file TrackDensitySeedFinder.h.

                                                               {
      this,
      "DensityEstimator",
      "Trk::GaussianTrackDensity",
      "Track density tool"
    };

---

The documentation for this class was generated from the following files:

• TrackDensitySeedFinder.h
• TrackDensitySeedFinder.cxx