Trk::IndexedCrossDistancesSeedFinder Class Reference

#include <IndexedCrossDistancesSeedFinder.h>

Inheritance diagram for Trk::IndexedCrossDistancesSeedFinder:

Collaboration diagram for Trk::IndexedCrossDistancesSeedFinder:

Public Member Functions
	IndexedCrossDistancesSeedFinder (const std::string &t, const std::string &n, const IInterface *p)
virtual	~IndexedCrossDistancesSeedFinder ()
virtual StatusCode	initialize () 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 Amg::Vector3D	findSeed (const double vx, const double vy, 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 Amg::Vector3D	findSeed (const double vx, const double vy, std::unique_ptr< Trk::IMode3dInfo > &info, 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::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.

Private Attributes
bool	m_useweights
float	m_trackdistcutoff
int	m_trackdistexppower
float	m_constraintcutoff
float	m_constrainttemp
unsigned int	m_maximumTracksNoCut
double	m_maximumDistanceCut
ToolHandle< IMode3dFinder >	m_mode3dfinder
ToolHandle< ITrkDistanceFinder >	m_distancefinder

Detailed Description

Definition at line 35 of file IndexedCrossDistancesSeedFinder.h.

Constructor & Destructor Documentation

IndexedCrossDistancesSeedFinder()

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

Definition at line 27 of file IndexedCrossDistancesSeedFinder.cxx.

                                                                                                                               : 
    base_class(t,n,p),
    m_useweights(true),
    m_trackdistcutoff(0.020),
    m_trackdistexppower(2),
    m_constraintcutoff(9.),
    m_constrainttemp(9.),
    m_maximumTracksNoCut(20),
    m_maximumDistanceCut(3.)
  {   
    declareProperty("useweights",m_useweights);
    declareProperty("trackdistcutoff",m_trackdistcutoff);
    declareProperty("trackdistexppower",m_trackdistexppower);
    declareProperty("constrainttemp",m_constrainttemp);
    declareProperty("constraintcutoff",m_constraintcutoff);
    declareProperty("maximumTracksNoCut",m_maximumTracksNoCut);
    declareProperty("maximumDistanceCut",m_maximumDistanceCut);
  }

~IndexedCrossDistancesSeedFinder()

Trk::IndexedCrossDistancesSeedFinder::~IndexedCrossDistancesSeedFinder ( )

virtualdefault

Member Function Documentation

findMultiSeeds() [1/2]

std::vector< Amg::Vector3D > Trk::IndexedCrossDistancesSeedFinder::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 251 of file IndexedCrossDistancesSeedFinder.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::IndexedCrossDistancesSeedFinder::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 260 of file IndexedCrossDistancesSeedFinder.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]

Amg::Vector3D Trk::IndexedCrossDistancesSeedFinder::findSeed ( const double vx, const double vy, 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. Must specify the primary vertex position.

Definition at line 78 of file IndexedCrossDistancesSeedFinder.cxx.

  {
    std::unique_ptr<Trk::IMode3dInfo> info;
    return findSeed (vx, vy, info, perigeeList, constraint);
  }

findSeed() [2/4]

Amg::Vector3D Trk::IndexedCrossDistancesSeedFinder::findSeed ( const double vx, const double vy, std::unique_ptr< Trk::IMode3dInfo > & info, 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. Must specify the primary vertex position. Allows returning additional diagnostic information.

Definition at line 89 of file IndexedCrossDistancesSeedFinder.cxx.

  {
    ATH_MSG_DEBUG( " Enter IndexedCrossDistancesSeedFinder  " );
 
    bool useCutOnDistance=false;
    if (perigeeList.size()>m_maximumTracksNoCut)
    {
      useCutOnDistance=true;
    }
    
    //now implement the code you already had in the standalone code...
 
    //Calculate and cache the covariance matrix for the constraint
    AmgSymMatrix(3) weightMatrixPositionConstraint;
    weightMatrixPositionConstraint.setIdentity(); //very arbitrary
    if (constraint != nullptr) {
      weightMatrixPositionConstraint = constraint->covariancePosition().inverse();
    }
    
    //Prepare the vector of points, on which the 3d mode has later to be calculated
    std::vector<PositionAndWeight> CrossingPointsAndWeights;
    std::vector<Amg::Vector3D> CrossingPoints;
    
    //Implement here SeedPointFinder algorithm, acting with the track vector.... 
    const std::vector<const Trk::TrackParameters*>::const_iterator begin=perigeeList.begin();
    const std::vector<const Trk::TrackParameters*>::const_iterator end=perigeeList.end();
 
    ATH_MSG_DEBUG( " Loop pairs of TrackParameters for modes " );
 
    std::vector< std::pair <int, int> > trkidx ;
    int idx_i = 0 ;
    for (std::vector<const Trk::TrackParameters*>::const_iterator i=begin;i!=end-1;++i) 
    {
      idx_i ++ ;
      const Trk::Perigee* MyI=dynamic_cast<const Trk::Perigee*>(*i);
      if (MyI==nullptr) {
    ATH_MSG_WARNING( "Neutrals not supported for seeding. Rejecting this track..." );
    continue;
      } 
      
      int idx_j = idx_i ;   //  1 has been added to idx_i
      for (std::vector<const Trk::TrackParameters*>::const_iterator j=i+1;j!=end;++j++) {
 
        idx_j ++ ;
        const Trk::Perigee* MyJ=dynamic_cast<const Trk::Perigee*>(*j);
 
        if (MyJ==nullptr) {
      ATH_MSG_WARNING( "Neutrals not supported for seeding. Rejecting this track..." );
      continue;
        }       
 
    
    try {
      
          std::optional<ITrkDistanceFinder::TwoPoints> result
            = m_distancefinder->CalculateMinimumDistance(*MyI,*MyJ);
      if (!result) {ATH_MSG_DEBUG("Problem with distance finder: THIS POINT WILL BE SKIPPED!"); 
          } 
          else 
          {
        //Get the points which connect the minimum distance between the two tracks
            double distance = Amg::distance (result->first, result->second);
#ifdef CROSSDISTANCESSEEDFINDER_DEBUG
        ATH_MSG_DEBUG("Point 1 x: " << result->first.x() << 
              " y: " << result->first.y() << 
              " z: " << result->first.z() );
        ATH_MSG_DEBUG("Point 2 x: " << result->second.x() << 
              " y: " << result->second.y() << 
              " z: " << result->second.z() );
        ATH_MSG_DEBUG("distance is: " << distance );
#endif
        
        Amg::Vector3D thepoint((result->first+result->second)/2.);
 
            if (m_useweights)
            {
              PositionAndWeight thispoint(thepoint,
                                          1./pow(m_trackdistcutoff+distance,m_trackdistexppower));
              
              ATH_MSG_VERBOSE("distance weight: " << 1./pow(m_trackdistcutoff+distance,m_trackdistexppower) );
              
              if (constraint!=nullptr) {
                
        Amg::Vector3D DeltaP(thepoint-constraint->position());
        Amg::Vector3D DeltaPConv;
                ATH_MSG_DEBUG("position x: " << DeltaP.x() << "position y: " << DeltaP.y() << "position z: " << DeltaP.z() );
                DeltaPConv[0]=DeltaP.x();
                DeltaPConv[1]=DeltaP.y();
                DeltaPConv[2]=DeltaP.z();
 
 
                double chi2=DeltaPConv.transpose()*weightMatrixPositionConstraint*DeltaPConv;
 
                ATH_MSG_DEBUG(" chi: " << chi2  << 
                  " beam weight " << 1./(1.+exp((chi2-m_constraintcutoff)/m_constrainttemp)) );
                
                thispoint.second=thispoint.second*1./(1.+exp((chi2-m_constraintcutoff)/m_constrainttemp));
                
              }
 
              if ((!useCutOnDistance || distance<m_maximumDistanceCut) && thispoint.second > 1e-10)
              {
                CrossingPointsAndWeights.push_back( thispoint );
 
                trkidx.emplace_back( idx_i - 1 , idx_j - 1 );
 
                ATH_MSG_VERBOSE( " crossing with track pair : " << idx_i - 1 <<" " 
                 << MyI->parameters()[Trk::d0] <<" "<< idx_j - 1 <<" " 
                 << MyJ->parameters()[Trk::d0] );
 
              }
            }
            else
            {
          const Amg::Vector3D& thispoint(thepoint);
              if ( !useCutOnDistance || distance<m_maximumDistanceCut )
              {
                CrossingPoints.push_back( thispoint );
              }
            }
          }
        } catch (...) {
          ATH_MSG_ERROR("Something wrong in distance calculation: please report..." );
        }
      }
      //to be understood...
    }
 
    //Now all points have been collected (N*(N-1)/2) and 
    //the mode has to be calculated
 
    if (    ( CrossingPoints.empty() && ! m_useweights  ) 
         || ( m_useweights && CrossingPointsAndWeights.empty() ) )
    {
      return Amg::Vector3D(0,0,0);
    }
   
    ATH_MSG_DEBUG(" crossing points prepared : " << CrossingPointsAndWeights.size() ); 
 
    Amg::Vector3D myresult;
 
    if (m_useweights)
    {
      myresult=m_mode3dfinder->getMode(vx, vy, CrossingPointsAndWeights, info);
    }
    else
    {
      myresult=m_mode3dfinder->getMode(vx, vy, CrossingPoints, info);
    }
 
    info->setTrkidx (std::move (trkidx));
    ATH_MSG_DEBUG(" 3D modes found ! " ); 
 
    return myresult;
    
  }

findSeed() [3/4]

Amg::Vector3D Trk::IndexedCrossDistancesSeedFinder::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 60 of file IndexedCrossDistancesSeedFinder.cxx.

  {
    ATH_MSG_ERROR ("Need to supply a primary vertex.");
    return Amg::Vector3D(0.,0.,0.);
  }

findSeed() [4/4]

Amg::Vector3D Trk::IndexedCrossDistancesSeedFinder::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 69 of file IndexedCrossDistancesSeedFinder.cxx.

  {
    ATH_MSG_ERROR ("Need to supply a primary vertex.");
    return Amg::Vector3D(0.,0.,0.);
  }

initialize()

StatusCode Trk::IndexedCrossDistancesSeedFinder::initialize ( )

overridevirtual

Definition at line 51 of file IndexedCrossDistancesSeedFinder.cxx.

  { 
    ATH_CHECK( m_mode3dfinder.retrieve() );
    ATH_CHECK( m_distancefinder.retrieve() );
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_constraintcutoff

float Trk::IndexedCrossDistancesSeedFinder::m_constraintcutoff

private

Definition at line 124 of file IndexedCrossDistancesSeedFinder.h.

m_constrainttemp

float Trk::IndexedCrossDistancesSeedFinder::m_constrainttemp

private

Definition at line 125 of file IndexedCrossDistancesSeedFinder.h.

m_distancefinder

ToolHandle<ITrkDistanceFinder> Trk::IndexedCrossDistancesSeedFinder::m_distancefinder

private

Initial value:

{
      this,
      "TrkDistanceFinder",
      "Trk::SeedNewtonTrkDistanceFinder"
    }

Definition at line 132 of file IndexedCrossDistancesSeedFinder.h.

                                                   {
      this,
      "TrkDistanceFinder",
      "Trk::SeedNewtonTrkDistanceFinder"
    };

m_maximumDistanceCut

double Trk::IndexedCrossDistancesSeedFinder::m_maximumDistanceCut

private

Definition at line 127 of file IndexedCrossDistancesSeedFinder.h.

m_maximumTracksNoCut

unsigned int Trk::IndexedCrossDistancesSeedFinder::m_maximumTracksNoCut

private

Definition at line 126 of file IndexedCrossDistancesSeedFinder.h.

m_mode3dfinder

ToolHandle<IMode3dFinder> Trk::IndexedCrossDistancesSeedFinder::m_mode3dfinder

private

Initial value:

{ this,
                                              "Mode3dFinder",
                                              "Trk::Mode3dFromFsmw1dFinder" }

Definition at line 129 of file IndexedCrossDistancesSeedFinder.h.

                                            { this,
                                              "Mode3dFinder",
                                              "Trk::Mode3dFromFsmw1dFinder" };

m_trackdistcutoff

float Trk::IndexedCrossDistancesSeedFinder::m_trackdistcutoff

private

Definition at line 122 of file IndexedCrossDistancesSeedFinder.h.

m_trackdistexppower

int Trk::IndexedCrossDistancesSeedFinder::m_trackdistexppower

private

Definition at line 123 of file IndexedCrossDistancesSeedFinder.h.

m_useweights

bool Trk::IndexedCrossDistancesSeedFinder::m_useweights

private

Definition at line 121 of file IndexedCrossDistancesSeedFinder.h.

---

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

• IndexedCrossDistancesSeedFinder.h
• IndexedCrossDistancesSeedFinder.cxx