Trk::CrossDistancesSeedFinder Class Reference

#include <CrossDistancesSeedFinder.h>

Inheritance diagram for Trk::CrossDistancesSeedFinder:

Collaboration diagram for Trk::CrossDistancesSeedFinder:

Public Member Functions
	CrossDistancesSeedFinder (const std::string &t, const std::string &n, const IInterface *p)
virtual	~CrossDistancesSeedFinder ()
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 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
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 42 of file CrossDistancesSeedFinder.h.

Constructor & Destructor Documentation

CrossDistancesSeedFinder()

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

Definition at line 42 of file CrossDistancesSeedFinder.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);
  }

~CrossDistancesSeedFinder()

Trk::CrossDistancesSeedFinder::~CrossDistancesSeedFinder ( )

virtualdefault

Member Function Documentation

findMultiSeeds() [1/2]

std::vector< Amg::Vector3D > Trk::CrossDistancesSeedFinder::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 244 of file CrossDistancesSeedFinder.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::CrossDistancesSeedFinder::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 252 of file CrossDistancesSeedFinder.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::CrossDistancesSeedFinder::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 73 of file CrossDistancesSeedFinder.cxx.

                                                                                                                                 {
    
 
    //create perigees from track list
    std::vector<const TrackParameters*> perigeeList;
    for (const auto *iter : VectorTrk) {
      if (std::isnan(iter->perigeeParameters()->parameters()[Trk::d0])) {
    continue;
      }  
      perigeeList.push_back(iter->perigeeParameters());
    }
 
    if (perigeeList.size()<2)
    {
      return Amg::Vector3D(0.,0.,0.);
    }
   
    //create seed from perigee list
    return findSeed(perigeeList,constraint);
    
  }

findSeed() [4/4]

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

                                                                                                                                             {
 
    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();
 
    for (std::vector<const Trk::TrackParameters*>::const_iterator i=begin;i!=end-1;++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;
      } 
      
      for (std::vector<const Trk::TrackParameters*>::const_iterator j=i+1;j!=end;++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;
      }     
 
#ifdef CROSSDISTANCESSEEDFINDER_DEBUG
    
      ATH_MSG_DEBUG( "Track 1 d0: " << MyI->parameters()[Trk::d0] << " Track2 d0: " << MyJ->parameters()[Trk::d0]  );
      ATH_MSG_DEBUG( "Track 1 z0: " << MyI->parameters()[Trk::z0] << "Track2 z0: " << MyJ->parameters()[Trk::z0]  );
#endif
    
    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 = dist (result.value());
#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_DEBUG( "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);
              }
            }
            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() && CrossingPointsAndWeights.empty())
    {
      return Amg::Vector3D(0,0,0);
    }
   
    Amg::Vector3D myresult;
 
    if (m_useweights)
    {
      myresult=m_mode3dfinder->getMode(0, 0, CrossingPointsAndWeights);
    }
    else
    {
      myresult=m_mode3dfinder->getMode(0, 0, CrossingPoints);
    }
    
#ifdef CROSSDISTANCESSEEDFINDER_DEBUG
    ATH_MSG_INFO( "Resulting mean POINT FOUND:  x: " << myresult.x() << 
                  " y: " << myresult.y() << 
                  " z: " << myresult.z()  );
#endif
    
    
    return myresult;
    
    
    
    
  }

initialize()

StatusCode Trk::CrossDistancesSeedFinder::initialize ( )

overridevirtual

Definition at line 66 of file CrossDistancesSeedFinder.cxx.

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

Member Data Documentation

m_constraintcutoff

float Trk::CrossDistancesSeedFinder::m_constraintcutoff

private

Definition at line 104 of file CrossDistancesSeedFinder.h.

m_constrainttemp

float Trk::CrossDistancesSeedFinder::m_constrainttemp

private

Definition at line 105 of file CrossDistancesSeedFinder.h.

m_distancefinder

ToolHandle<ITrkDistanceFinder> Trk::CrossDistancesSeedFinder::m_distancefinder

private

Initial value:

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

Definition at line 112 of file CrossDistancesSeedFinder.h.

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

m_maximumDistanceCut

double Trk::CrossDistancesSeedFinder::m_maximumDistanceCut

private

Definition at line 107 of file CrossDistancesSeedFinder.h.

m_maximumTracksNoCut

unsigned int Trk::CrossDistancesSeedFinder::m_maximumTracksNoCut

private

Definition at line 106 of file CrossDistancesSeedFinder.h.

m_mode3dfinder

ToolHandle<IMode3dFinder> Trk::CrossDistancesSeedFinder::m_mode3dfinder

private

Initial value:

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

Definition at line 109 of file CrossDistancesSeedFinder.h.

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

m_trackdistcutoff

float Trk::CrossDistancesSeedFinder::m_trackdistcutoff

private

Definition at line 102 of file CrossDistancesSeedFinder.h.

m_trackdistexppower

int Trk::CrossDistancesSeedFinder::m_trackdistexppower

private

Definition at line 103 of file CrossDistancesSeedFinder.h.

m_useweights

bool Trk::CrossDistancesSeedFinder::m_useweights

private

Definition at line 101 of file CrossDistancesSeedFinder.h.

---

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

• CrossDistancesSeedFinder.h
• CrossDistancesSeedFinder.cxx