RecCgemKalTrack.cxx

Go to the documentation of this file.

/*
*                                                                              
*  this class models "Kalman Track" used in Cgem Reconstruction        
*
* ********************************************************/
 
#include "CgemRecEvent/RecCgemKalTrack.h"
//#include "CgemRecEvent/CgemKalHelixSeg.h"
 
RecCgemKalTrack::RecCgemKalTrack() 
:DstCgemKalTrack() 
{
  HepVector a(5,0);
  HepSymMatrix Ea(5,0);
  HepPoint3D poca(0,0,0);
  HepPoint3D point(0.,0.,0);
  HepPoint3D pivot(0.,0.,0);
  
  m_ta(5,0);
  m_tEa(5,0);
  m_trackId = -1;
  for(int pid=0; pid<5; pid++){
    m_lhelixs.push_back(a);
    m_lerrors.push_back(Ea);
//    m_fhelixs.push_back(a);
//    m_ferrors.push_back(Ea);
    m_lpoints.push_back(point);
    m_lpivots.push_back(pivot);
    
    m_length[pid] = 0.0;
    m_tof[pid] = 0.0;
    m_nhits[pid] = 0;
    m_pathSM[pid] = 0.0;
    m_fiTerm[pid] = 0.0;
    m_stat[0][pid] = -1;
    m_stat[1][pid] = -1;
    m_firstLayer[0][pid] = 0;
    m_firstLayer[1][pid] = 0;
    m_lastLayer[0][pid] = 43; 
    m_lastLayer[1][pid] = 43;
    m_chisq[0][pid] = 999.0;    
    m_chisq[1][pid] = 999.0;   
    m_ndf[0][pid] = 0;
    m_ndf[1][pid] = 0;
  }
  for( int i=0; i<43; i++) {   
      m_pathl[i] = 0.0;
  }
  for(int j=0; j<5; j++) {
      m_thelix[j] = 0.0;
  }
  for(int k=0; k<15; k++){
      m_terror[k]= 0.;
  }
}
 
 
    RecCgemKalTrack::RecCgemKalTrack( const RecCgemKalTrack& trk )
:DstCgemKalTrack(trk),
    m_vechelixsegs(trk.m_vechelixsegs)
{   
}
 
    RecCgemKalTrack::RecCgemKalTrack( const DstCgemKalTrack& dstrk )
:DstCgemKalTrack(dstrk)
{
    HepVector a(5,-99);
    HepSymMatrix Ea(5,-99);
    HepPoint3D poca(-99.,-99.,-99.);
    HepPoint3D point(0.,0.,0);
    HepPoint3D pivot(0.,0.,0);
 
    m_ta(5,0);
    m_tEa(5,0);
 
    //m_trackId = -1;
    for(int pid=0; pid<5; pid++){
        //m_pocas.push_back(poca);
        //m_zhelixs.push_back(a);
    //m_zerrors.push_back(Ea);
    m_lhelixs.push_back(a);
    m_lerrors.push_back(Ea);
 //   m_fhelixs.push_back(a);
 //   m_ferrors.push_back(Ea);
    m_lpoints.push_back(point);
    m_lpivots.push_back(pivot);
    
    m_length[pid] = 0.0;
    m_tof[pid] = 0.0;
    m_nhits[pid] = 0;
    m_pathSM[pid] = 0.0;
    m_fiTerm[pid] = 0.0;
     
    m_stat[0][pid] = -1;
    m_stat[1][pid] = -1;
    m_firstLayer[0][pid] = 0;
    m_firstLayer[1][pid] = 0;
    m_lastLayer[0][pid] = 43;
    m_lastLayer[1][pid] = 43;
    m_chisq[0][pid] = 999.0;
    m_chisq[1][pid] = 999.0;
    m_ndf[0][pid] = 0;
    m_ndf[1][pid] = 0;
  }
  for(int i=0; i<43; i++) {
    m_pathl[i] = 0.0;
  }
  for(int j=0; j<5; j++) {
      m_thelix[j] = 0.0;
  }
  for(int k=0; k<15; k++){
      m_terror[k]= 0.;
  }
}
 
RecCgemKalTrack& RecCgemKalTrack::operator=(const DstCgemKalTrack& dstrk)
{
  if( &dstrk != this ){
    DstCgemKalTrack::operator=(dstrk);
    HepVector a(5,0);
    HepSymMatrix Ea(5,0);
    HepPoint3D poca(-99.,-99.,-99.);
    HepPoint3D point(0.,0.,0.); 
    HepPoint3D pivot(0.,0.,0.);
    
    //m_trackId = -1;
    //m_pocas.clear();
    //m_zhelixs.clear();
    //m_zerrors.clear();
    m_lhelixs.clear();
    m_lerrors.clear();
 //   m_fhelixs.clear();
 //   m_ferrors.clear();
    m_lpoints.clear();
    m_lpivots.clear();
    m_ta(5,0) ;  
    m_tEa(5,0);
    for(int pid=0; pid<5; pid++){
      //m_pocas.push_back(poca);
      //m_zhelixs.push_back(a);
      //m_zerrors.push_back(Ea);
      m_lhelixs.push_back(a);
      m_lerrors.push_back(Ea);
 //     m_fhelixs.push_back(a);
 //     m_ferrors.push_back(Ea);
      m_lpoints.push_back(point);
  
      
      m_length[pid] = 0.0;
      m_tof[pid] = 0.0;
      m_nhits[pid] = 0;
      m_pathSM[pid] = 0.;
      m_fiTerm[pid] = 0.;
    }
    
    for(int i=0; i<43; i++) {
      m_pathl[i] = 0.0;
    }
    for(int j=0; j<5; j++) {
      m_thelix[j] = 0.0;
    }
    for(int k=0; k<15; k++){
        m_terror[k]= 0.;
    }
  }
  return  *this;
}
 
RecCgemHelixSegRefVec RecCgemKalTrack::getVecHelixSegs(void) const {
  return m_vechelixsegs;
}
                    
void RecCgemKalTrack::setVecHelixSegs(const RecCgemHelixSegRefVec& vechelixsegs){
  m_vechelixsegs = vechelixsegs;
}
 
RecCgemKalTrack::~RecCgemKalTrack() {
 
}