TCgemTrack.cxx

Go to the documentation of this file.

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