ImpactPoint3dEstimator.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRKVERTEXFITTERUTILS_IMPACTPOINT3DESTIMATOR_H
#define TRKVERTEXFITTERUTILS_IMPACTPOINT3DESTIMATOR_H
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ToolHandle.h"
#include "TrkVertexFitterInterfaces/IImpactPoint3dEstimator.h"
#include "MagFieldConditions/AtlasFieldCacheCondObj.h"
#include "TrkSurfaces/PlaneSurface.h"
#include <memory>
 

 
namespace Trk
{
  class IExtrapolator;
 
 
  class ImpactPoint3dEstimator : public extends<AthAlgTool, IImpactPoint3dEstimator>
  {
  public:
    virtual StatusCode initialize() override;
    virtual StatusCode finalize() override;
    
    ImpactPoint3dEstimator(const std::string& t, const std::string& n, const IInterface*  p);
 
    virtual ~ImpactPoint3dEstimator();
    
    virtual std::unique_ptr<PlaneSurface>
    Estimate3dIP(const EventContext& ctx,
                 const Trk::TrackParameters* trackPerigee,
                 const Amg::Vector3D* theVertex,
                 double& distance) const override;
 
    virtual std::unique_ptr<PlaneSurface>
    Estimate3dIP(const EventContext& ctx,
                 const Trk::NeutralParameters* neutralPerigee,
                 const Amg::Vector3D* theVertex,
                 double& distance) const override;
   
    virtual bool addIP3dAtaPlane(const EventContext& ctx, VxTrackAtVertex &,const Amg::Vector3D & vertex) const override;

    virtual const Trk::AtaPlane * IP3dAtaPlane(const EventContext& ctx,VxTrackAtVertex & vtxTrack,const Amg::Vector3D & vertex) const override;
 
    //Same for neutrals
    virtual const Trk::NeutralAtaPlane * IP3dNeutralAtaPlane(const EventContext& ctx,const NeutralParameters * initNeutPerigee,const Amg::Vector3D & vertex) const override;
 
  private:
    template<typename T>
    std::unique_ptr<PlaneSurface>
    Estimate3dIPNoCurvature(const T*,
                            const Amg::Vector3D* theVertex,
                            double& distance) const;
 
    
    
    ToolHandle< Trk::IExtrapolator > m_extrapolator;
    SG::ReadCondHandleKey<AtlasFieldCacheCondObj> m_fieldCacheCondObjInputKey 
          {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"};
 
    int m_maxiterations;
    double m_precision;
    
  };
 
}
 
#endif