Ext_xp_err.h

Go to the documentation of this file.

//
// File: Ext_xp_err.h
//
// Descrition: Handle extrapolation error matrix( x, y, z, px, py, pz ).
//             The used coordinate system is the cartesian BESIII coordinate system.
//
//             Modified from BELLE:  " File: Ext_xp_err.h
//                                    Creation: 29-Dec-1997
//                                    Version: 01-Mar-1999
//                                    $Id: Ext_xp_err.h,v 1.4.4.1 2013/02/21 02:30:58 lich Exp $"
//
// Data: 2005.3.30 by    Wang Liangliang
//
 
 
#ifndef _Ext_XP_Err_Flag_
#define _Ext_XP_Err_Flag_
#include <iosfwd>
 
///////////////////////change//////////
//#include  "TrkExtAlg/Helix.h"
#include "TrackUtil/Helix.h"
///////////////////////////////////////
#include    "TrkExtAlg/Ext_errmx.h"
 
#include        <vector>
 
 
class Ext_xp_err : public Ext_errmx {
 
public:
  Ext_xp_err();// Constructor
  Ext_xp_err( const Ext_xp_err &xp_err );// Copy constructor
  ~Ext_xp_err(){};// Destructor
   
  // Set symmetry matrix parameters. 
  void set_err(const HepSymMatrix &err,const Hep3Vector &xv,const Hep3Vector &pv, const double &q, 
           const double &mass );
                    
  inline const Hep3Vector & get_x() const;// Get coordinate at the current point.
  inline const Hep3Vector & get_p() const;// Get momentum at the current point.
  inline const double & get_q() const; // Get the charge of the particle.
  inline double get_mass() const;// Get the mass of the particle.
  bool move( const Hep3Vector &xv1, const Hep3Vector &pv1, const Hep3Vector &B,
         const int ms_on, const double chi_cc ); // Move x -> x'.
  Ext_xp_err &operator=( const Ext_xp_err &xp_err );
  friend std::ostream &operator<<( std::ostream &s, const Ext_xp_err &xp_err );
  inline  void set_pos(const Hep3Vector &pos);
  inline  void set_mom(const Hep3Vector &mom);   
 
private:
  Hep3Vector    m_xv;       // Position(x,y,z) in cm unit.
  Hep3Vector    m_pv;       // Momentum(px,py,pz) in GeV/c.
  HepMatrix m_xp_jcb;   // Jacobian for x -> x'.
  HepMatrix m_h2xp_jcb; // Jacobian for helix -> xp.
  double    m_q;        // Charge of the particle.
  double    m_mass2;    // Mass**2 of the particle.
 
};
 
 
//  Get the current coordinate.
inline const Hep3Vector & Ext_xp_err::get_x() const {
  return m_xv;
}
 
 
//  Get the momentum at the current position.
inline const Hep3Vector & Ext_xp_err::get_p() const {
  return m_pv;
}
 
 
//  Get the charge of the particle.
inline const double & Ext_xp_err::get_q() const {
  return m_q;
}
 
 
//  Get the mass of the particle.
inline double Ext_xp_err::get_mass() const {
  return sqrt(m_mass2);
}
 
inline  void Ext_xp_err::set_pos(const Hep3Vector &pos)
{
m_xv = pos;
}
 
inline  void Ext_xp_err::set_mom(const Hep3Vector &mom)
{
m_pv = mom;
}
 
 
#endif /* _Ext_XP_Err_Flag_ */