Helix.h

Go to the documentation of this file.

 
#ifndef PANDORA_HELIX_H
#define PANDORA_HELIX_H 1
 
#include "Objects/CartesianVector.h"
 
#include "Pandora/StatusCodes.h"
 
namespace pandora
{
 
class Helix
{
public:
    Helix(const float phi0, const float d0, const float z0, const float omega, const float tanlambda, const float bField);
 
    Helix(const CartesianVector &position, const CartesianVector &momentum, const float charge, const float bField);
 
    StatusCode GetPointInXY(const float x0, const float y0, const float ax, const float ay, const CartesianVector &referencePoint,
        CartesianVector &intersectionPoint) const;
 
    StatusCode GetPointInXY(const float x0, const float y0, const float ax, const float ay, const CartesianVector &referencePoint,
        CartesianVector &intersectionPoint, float &genericTime) const;
 
    StatusCode GetPointInZ(const float zPlane, const CartesianVector &referencePoint, CartesianVector &intersectionPoint) const;
 
    StatusCode GetPointInZ(const float zPlane, const CartesianVector &referencePoint, CartesianVector &intersectionPoint,
        float &genericTime) const;
 
    StatusCode GetPointOnCircle(const float radius, const CartesianVector &referencePoint, CartesianVector &intersectionPoint) const;
 
    StatusCode GetPointOnCircle(const float radius, const CartesianVector &referencePoint, CartesianVector &intersectionPoint,
        float &genericTime) const;
 
    StatusCode GetDistanceToPoint(const CartesianVector &point, CartesianVector &distance) const;
 
    StatusCode GetDistanceToPoint(const CartesianVector &point, CartesianVector &distance, float &genericTime) const;
 
    StatusCode GetDistanceToHelix(const Helix *const pHelix, CartesianVector &positionOfClosestApproach, CartesianVector &v0momentum,
        float &helixDistance) const;
 
    CartesianVector GetExtrapolatedMomentum(const CartesianVector &position) const;
 
    const CartesianVector &GetMomentum() const;
 
    const CartesianVector &GetReferencePoint() const;
 
    float GetPhi0() const;
 
    float GetD0() const;
 
    float GetZ0() const;
 
    float GetOmega() const;
 
    float GetTanLambda() const;
 
    float GetPxy() const;
 
    float GetCharge() const;
 
    float GetXCentre() const;
 
    float GetYCentre() const;
 
    float GetRadius() const;
 
private:
    static const float FCT;
    static const float TWO_PI;
    static const float HALF_PI;
 
    CartesianVector     m_referencePoint;       
    CartesianVector     m_momentum;             
 
    float               m_phi0;                 
    float               m_d0;                   
    float               m_z0;                   
    float               m_omega;                
    float               m_tanLambda;            
    float               m_pxy;                  
    float               m_charge;               
    float               m_xCentre;              
    float               m_yCentre;              
    float               m_radius;               
    float               m_phiRefPoint;          
    float               m_phiAtPCA;             
    float               m_xAtPCA;               
    float               m_yAtPCA;               
    float               m_pxAtPCA;              
    float               m_pyAtPCA;              
    float               m_phiMomRefPoint;       
};
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline StatusCode Helix::GetPointInXY(const float x0, const float y0, const float ax, const float ay, const CartesianVector &referencePoint,
    CartesianVector &intersectionPoint) const
{
    float genericTime;
    return this->GetPointInXY(x0, y0, ax, ay, referencePoint, intersectionPoint, genericTime);
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline StatusCode Helix::GetPointInZ(const float zPlane, const CartesianVector &referencePoint, CartesianVector &intersectionPoint) const
{
    float genericTime;
    return this->GetPointInZ(zPlane, referencePoint, intersectionPoint, genericTime);
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline StatusCode Helix::GetPointOnCircle(const float radius, const CartesianVector &referencePoint, CartesianVector &intersectionPoint) const
{
    float genericTime;
    return this->GetPointOnCircle(radius, referencePoint, intersectionPoint, genericTime);
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline StatusCode Helix::GetDistanceToPoint(const CartesianVector &point, CartesianVector &distance) const
{
    float genericTime;
    return this->GetDistanceToPoint(point, distance, genericTime);
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline const CartesianVector &Helix::GetMomentum() const
{
    return m_momentum;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline const CartesianVector &Helix::GetReferencePoint() const
{
    return m_referencePoint;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetPhi0() const
{
    if (m_phi0 < 0.)
        return m_phi0 + TWO_PI;
 
    return m_phi0;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetD0() const
{
    return m_d0;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetZ0() const
{
    return m_z0;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetOmega() const
{
    return m_omega;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetTanLambda() const
{
    return m_tanLambda;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetPxy() const
{
    return m_pxy;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetCharge() const
{
    return m_charge;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetXCentre() const
{
    return m_xCentre;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetYCentre() const
{
    return m_yCentre;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
inline float Helix::GetRadius() const
{
    return m_radius;
}
 
} // namespace pandora
 
#endif // #ifndef PANDORA_HELIX_H