#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. More...

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. More...

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. More...

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. More...

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. More...

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. More...

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 29 of file IndexedCrossDistancesSeedFinder.cxx.

   {   
     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 253 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);
 }
  
 std::vector<Amg::Vector3D> IndexedCrossDistancesSeedFinder::findMultiSeeds( 

◆ 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 262 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 80 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 91 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 62 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 71 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 53 of file IndexedCrossDistancesSeedFinder.cxx.

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.

◆ 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.

◆ 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:

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ IndexedCrossDistancesSeedFinder()

◆ ~IndexedCrossDistancesSeedFinder()

Member Function Documentation

◆ findMultiSeeds() [1/2]

◆ findMultiSeeds() [2/2]

◆ findSeed() [1/4]

◆ findSeed() [2/4]

◆ findSeed() [3/4]

◆ findSeed() [4/4]

◆ initialize()

Member Data Documentation

◆ m_constraintcutoff

◆ m_constrainttemp

◆ m_distancefinder

◆ m_maximumDistanceCut

◆ m_maximumTracksNoCut

◆ m_mode3dfinder

◆ m_trackdistcutoff

◆ m_trackdistexppower

◆ m_useweights