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>
 
///////////////////////change//////////
// #include "Helix.h"
#include "TrackUtil/Helix.h"
///////////////////////////////////////
#include "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_ */