CrossDistancesSeedFinder.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRKVERTEXSEEDFINDERTOOLS_CROSSDISTANCESSEEDFINDER_H
#define TRKVERTEXSEEDFINDERTOOLS_CROSSDISTANCESSEEDFINDER_H
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "TrkVertexFitterInterfaces/IVertexSeedFinder.h"
#include "GaudiKernel/ToolHandle.h"
 
#include "TrkVertexSeedFinderUtils/IMode3dFinder.h"
#include "TrkVertexSeedFinderUtils/ITrkDistanceFinder.h"
 
namespace Trk
{
  class Track;
 
  // @author N. Giacinto Piacquadio (Albert-Ludwigs-Universitaet Freiburg - Germany)
  //
  // @ATLAS software
  //
  // This class implements a seed finder for the vertexing algorithm 
  // which is based on the mean points of the distances between all tracks 
  // (taken two by two). There's actually no selection on the tracks.
  // 
  // 
  // A very similar point finder was implemented for the first time by Waltenberger 
  // in CMS. I suggest his PhD thesis (2005) as a reference for the different concepts 
  // of point finders.
  //
  // ---------------------------------------
  // Changes:
  //
  // David Shope <david.richard.shope@cern.ch> (2016-04-19)
  //
  // EDM Migration to xAOD - from Trk::VxCandidate to xAOD::Vertex, 
  //                         from Trk::RecVertex   to xAOD::Vertex,
  //                         from Trk::Vertex      to Amg::Vector3D
 
 
  class CrossDistancesSeedFinder final: public extends<AthAlgTool, IVertexSeedFinder>
  {
  public:
    // Standard Gaudi constructor.
    CrossDistancesSeedFinder (const std::string& t,
                              const std::string& n,
                              const IInterface*  p);
 
    
    virtual ~CrossDistancesSeedFinder();
 
    virtual StatusCode initialize() override;
 
 
    using IVertexSeedFinder::findSeed;

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

    virtual Amg::Vector3D
    findSeed (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;
 

    virtual std::vector<Amg::Vector3D>
    findMultiSeeds (const std::vector<const Trk::TrackParameters*>& perigeeList,
                    const xAOD::Vertex * constraint=0) const override final;
 
 
  private:
 
    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{ this,
                                              "Mode3dFinder",
                                              "Trk::Mode3dTo1dFinder" };
    ToolHandle<ITrkDistanceFinder> m_distancefinder{
      this,
      "TrkDistanceFinder",
      "Trk::SeedNewtonTrkDistanceFinder"
    };
  };
}
#endif