TrackParameters.h

Go to the documentation of this file.

#ifndef ATLFAST_TRACKPARAMETERS_H
#define ATLFAST_TRACKPARAMETERS_H
//

#include "AtlfastEvent/ITrackParameters.h"
#include "AtlfastEvent/Phi.h"
#include <cmath>

namespace Atlfast
{
  class TrackParameters: public ITrackParameters
  {
    public:
      
      TrackParameters() {}
      TrackParameters( double eta,
                       Phi phi,
                       double pT,
                       double impactParameter,
                       double zPerigee,
                       double cotTheta,
                       double invPtCharge
                     ) :
        m_eta(eta),
        m_phi(phi),
        m_pT(pT),
        m_impactParameter(impactParameter),
        m_zPerigee(zPerigee),
        m_cotTheta(cotTheta),
        m_invPtCharge(invPtCharge) {}
      
      virtual double  eta()             const { return m_eta; }
      virtual Phi     phi()             const { return m_phi; } 
      virtual double  pT()              const { return m_pT; }
      virtual double  impactParameter() const { return m_impactParameter; }
      virtual double  zPerigee()        const { return m_zPerigee; }
      virtual double  cotTheta()        const { return m_cotTheta; }
      virtual double  invPtCharge()     const { return m_invPtCharge; }  
      
      // derived quantities
      double theta() const;
      double invPCharge() const;
 
    private:
    
      double m_eta;
      Phi    m_phi;
      double m_pT;
      double m_impactParameter;
      double m_zPerigee;
      double m_cotTheta;
      double m_invPtCharge;
    
  };

  //               inline (no virtual) methods           //
  inline double TrackParameters::theta() const
  {
    double theta = ( m_cotTheta  ?  std::atan( 1.0 / m_cotTheta )  :  M_PI/2 );
    if ( theta < 0.0 )  theta += M_PI;
    return theta;
  }
  
  inline double TrackParameters::invPCharge() const
  {
    double sinTheta = ( m_cotTheta  ?  std::sqrt( 1.0 / ( 1.0  + m_cotTheta * m_cotTheta ) )  :  1.0 );
    return m_invPtCharge * sinTheta;
  }      


}
#endif

---