/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/TrkExtAlg/TrkExtAlg-00-00-64/TrkExtAlg/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.5 2015/07/29 00:28:14 sunss Exp $"
//
// Data: 2005.3.30 by    Wang Liangliang
//


#ifndef _Ext_XP_Err_Flag_
#define _Ext_XP_Err_Flag_
#include <iosfwd>

//#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_ */

---