DstCgemTrack.cxx

Go to the documentation of this file.

 
#include "DstEvent/DstCgemTrack.h"
 
  DstCgemTrack::DstCgemTrack()
  : m_trackId(0),      
    m_charge(0),  
    m_poca3D(0,0,0),    
    m_a(5,0),      
    m_Ea(5,0),
    m_stat(0),
    m_chi2(0),
    m_ndof(0),    
    m_firstLayer(0),    
    m_lastLayer(42), 
    m_pxy(0), 
    m_px(0),
    m_py(0),
    m_pz(0),
    m_p(0),
    m_theta(0),      
    m_phi(0),   
    m_x(0),     
    m_y(0),   
    m_z(0),       
    m_r(0)
  {
    for(int i=0;i<5;i++)  m_helix[i] = 0;
    for(int j=0;j<15;j++) m_err[j]   = 999; 
    for(int k=0;k<3;k++)  m_poca[k]  = 0;    
  }
 
 DstCgemTrack::DstCgemTrack( const DstCgemTrack& other)
 : m_trackId(other.m_trackId),
   m_charge(other.m_charge),
   m_poca3D(other.m_poca3D),
   m_a(other.m_a),
   m_Ea(other.m_Ea),
   m_stat(other.m_stat),
   m_chi2(other.m_chi2),
   m_ndof(other.m_ndof),
   m_firstLayer(other.m_firstLayer),
   m_lastLayer(other.m_lastLayer),
   m_pxy(other.m_pxy),
   m_px(other.m_px),
   m_py(other.m_py),
   m_pz(other.m_pz),
   m_p(other.m_p),
   m_theta(other.m_theta),
   m_phi(other.m_phi),
   m_x(other.m_x),
   m_y(other.m_y),
   m_z(other.m_z),
   m_r(other.m_r)
   { 
     for(int i=0;i<5;i++)  m_helix[i] = (other.m_helix)[i];
     for(int j=0;j<15;j++) m_err[j]   = (other.m_err)[j];
     for(int k=0;k<3;k++)  m_poca[k]  = (other.m_poca)[k];
   }
 
 
  DstCgemTrack& DstCgemTrack::operator=(const DstCgemTrack& h){
    if( &h != this ){
      m_trackId    = h.m_trackId;
      m_charge     = h.m_charge;
      m_poca3D     = h.m_poca3D;
      m_a          = h.m_a;
      m_Ea         = h.m_Ea;
      m_stat       = h.m_stat;
      m_chi2       = h.m_chi2;
      m_ndof       = h.m_ndof;
      m_firstLayer = h.m_firstLayer;
      m_lastLayer  = h.m_lastLayer;
      m_pxy        = h.m_pxy;
      m_px         = h.m_px;
      m_py         = h.m_py;
      m_pz         = h.m_pz;
      m_p          = h.m_p;
      m_theta      = h.m_theta;
      m_phi        = h.m_phi;
      m_x          = h.m_x;
      m_y          = h.m_y;
      m_z          = h.m_z;
      m_r          = h.m_r;
    
      for(int i=0;i<5;i++)  m_helix[i] = (h.m_helix)[i];
      for(int j=0;j<15;j++) m_err[j]   = (h.m_err)[j];
      for(int k=0;k<3;k++)  m_poca[k]  = (h.m_poca)[k];
    }
    return *this;
  }
 
  void DstCgemTrack::setHelix(double helix[5]) {   
    for(int i=0; i<5; i++){
      m_helix[i] = helix[i];
      m_a[i] = helix[i]; 
    }
  }
 
  void DstCgemTrack::setPoca(double poca[3]){
    for(int i=0; i<3; i++){
      m_poca[i] = poca[i];
      m_poca3D[i] = poca[i];
    }
  }  
 
  void DstCgemTrack::setError( double err[15]) {           
    for (int i1=0; i1<15; i1++) {
      m_err[i1] = err[i1];
    }
    for(int i=0, k=0; i<5; i++) {
      for(int j=i; j<5; j++) {
    m_Ea[i][j] = m_err[k++];
    m_Ea[j][i] = m_Ea[i][j];
      }
    }
  }      
 
 
  void  DstCgemTrack::setHelix(const HepVector& helix){
    for(int i=0; i<5; i++){        
      m_helix[i] = helix[i];
    }
    m_a = helix; 
  }
  
  void  DstCgemTrack::setPoca(const HepPoint3D& poca){
    for(int i=0; i<3; i++){
      m_poca[i] = poca[i];
    }
    m_poca3D = poca;
  }
  
  void  DstCgemTrack::setError(const HepSymMatrix& err){
    int k = 0;
    for (int i=0; i<5; i++){
      for (int j=i; j<5; j++)
      {
    m_err[k] = err[i][j];
    k++;
      }
    }
    m_Ea = err;
  }
 
  const HepVector DstCgemTrack::helix() const{
    return m_a;
  }
  
  const HepSymMatrix DstCgemTrack::err() const{
    /*
    HepSymMatrix ea(5);
    for(int i = 0, k = 0; i < 5; i++) {
      for(int j = 0; j <= i; j++) {
    ea[i][j] = m_err[k++];
    ea[j][i] = ea[i][j];
      }
    }
    return ea;
   */
    return m_Ea;
  }
    
  const HepLorentzVector DstCgemTrack::p4(const double mass) const{
    return HepLorentzVector(p3(),sqrt(p3()*p3()+mass*mass)); 
  }