VertexPointEstimator.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef INDETCONVERSIONFINDERTOOLS_VERTEXPOINTESTIMATOR_H
#define INDETCONVERSIONFINDERTOOLS_VERTEXPOINTESTIMATOR_H
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "TrkParameters/TrackParameters.h"
 
namespace InDet {
  
  class VertexPointEstimator : public AthAlgTool {
    
  public:
    VertexPointEstimator (const std::string& type,const std::string& name, const IInterface* parent);
    
    virtual ~VertexPointEstimator();
    
    static const InterfaceID& interfaceID();
 
    virtual StatusCode initialize() override;
    
    virtual StatusCode finalize() override;
 
    typedef std::map<std::string, float> Values_t;
    
    Amg::Vector3D getCirclesIntersectionPoint(const Trk::Perigee* per1,
                                              const Trk::Perigee* per2,
                                              unsigned int flag,
                                              int& errorcode) const;
 
 
    Amg::Vector3D getCirclesIntersectionPoint(const Trk::Perigee* per1,
                                              const Trk::Perigee* per2,
                                              unsigned int flag,
                                              int& errorcode,
                                              Values_t& decors) const;
 
    static std::vector<std::string> decorKeys() ;
 
 
  private:
    Amg::Vector3D intersectionImpl (const Trk::Perigee *per1,
                                    const Trk::Perigee *per2,
                                    unsigned int flag,
                                    int& errorcode,
                                    float& deltaPhi,
                                    float& deltaR) const;
 
 
    static double areaVar(double, double, double, double, double, double, double&) ;
    static double areaVar(double, double, double, double, double, double, double&, double&, double&) ;
    static bool   circleIntersection(double, double, double, 
                  double, double, double, 
                  double&, double&, 
                  double&, double&) ;
    static bool   secondDegree(double, double, double, double&, double&) ;
    static double areaTriangle(double, double, double, double, double, double) ;
    
    static const double s_bmagnt;
    std::vector<double> m_maxDR;         
    std::vector<double> m_maxDZ;         
    std::vector<double> m_maxR;          
    std::vector<double> m_minArcLength;  
    std::vector<double> m_maxArcLength;  
    std::vector<double> m_minDr;         
    std::vector<double> m_maxDr;         
    std::vector<double> m_maxHl;         
    std::vector<double> m_maxPhi;        
    bool m_returnOnError;
  };
  
}
#endif // INDETCONVERSIONFINDERTOOLS_VERTEXPOINTESTIMATOR_H