TMdcTrack.cxx

Go to the documentation of this file.

#include "RootEventData/TMdcTrack.h"
#include <cmath>
 
using namespace std;
 
ClassImp( TMdcTrack );
 
//************************************************
//
TMdcTrack::TMdcTrack() { Clear(); }
//************************************************
//
TMdcTrack::~TMdcTrack() { Clear(); }
Int_t TMdcTrack::charge() const {
  Int_t    charge;
  Double_t temp = m_helix[2];
  if ( temp > 0.0000000001 ) charge = 1;
  else if ( temp < -0.0000000001 ) charge = -1;
  else charge = 0;
  return charge;
}
 
Double_t TMdcTrack::x() const { return m_helix[0] * cos( m_helix[1] ); }
 
Double_t TMdcTrack::y() const { return m_helix[0] * sin( m_helix[1] ); }
 
Double_t TMdcTrack::z() const { return m_helix[3]; }
 
Double_t TMdcTrack::r() const {
  //    return  sqrt(x()*x()+y()*y());
  return fabs( m_helix[0] );
}
 
Double_t TMdcTrack::pxy() const {
  if ( m_helix[2] != 0 ) return 1. / fabs( m_helix[2] );
  else return 0.;
}
 
Double_t TMdcTrack::px() const { return pxy() * ( -sin( m_helix[1] ) ); }
 
Double_t TMdcTrack::py() const { return pxy() * cos( m_helix[1] ); }
 
Double_t TMdcTrack::pz() const { return pxy() * m_helix[4]; }
 
Double_t TMdcTrack::p() const {
  //    return  sqrt(px()*px() + py()*py() + pz()*pz());
  return pxy() * sqrt( 1. + m_helix[4] * m_helix[4] );
}
 
Double_t TMdcTrack::theta() const { return acos( pz() / p() ); }
 
Double_t TMdcTrack::phi() const { return atan2( py(), px() ); }
 
void TMdcTrack::setHelix( const Double_t helix[5] ) {
  for ( Int_t i = 0; i < 5; i++ ) { m_helix[i] = helix[i]; }
}
 
void TMdcTrack::setErr( const Double_t er[15] ) {
  for ( Int_t i = 0; i < 15; i++ ) { m_err[i] = er[i]; }
}