RecMdcKalTrack.cxx

Go to the documentation of this file.

/*
 *
 *  this class models "Kalman Track" used in Mdc Reconstruction
 *
 * ********************************************************/
 
#include "MdcRecEvent/RecMdcKalTrack.h"
// #include "MdcRecEvent/MdcKalHelixSeg.h"
 
#include <cstdio>
#include <fstream>
 
RecMdcKalTrack::RecMdcKalTrack() : DstMdcKalTrack() {
  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;
 
  pidSeg = 4;
 
  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_nster[0][pid]      = -1;
    m_nster[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.; }
}
 
RecMdcKalTrack::RecMdcKalTrack( const RecMdcKalTrack& trk ) : DstMdcKalTrack( trk ) {
  for ( int i = 0; i < 5; i++ ) m_vechelixsegs[i] = trk.m_vechelixsegs[i];
}
 
RecMdcKalTrack::RecMdcKalTrack( const DstMdcKalTrack& dstrk ) : DstMdcKalTrack( 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_nster[0][pid]      = -1;
    m_nster[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.; }
}
 
RecMdcKalTrack& RecMdcKalTrack::operator=( const DstMdcKalTrack& dstrk ) {
  if ( &dstrk != this )
  {
    DstMdcKalTrack::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;
}
 
// 2011-05-06
HelixSegRefVec RecMdcKalTrack::getVecHelixSegs( int pid ) const {
  if ( pid >= 0 && pid < 5 ) return m_vechelixsegs[pid];
  else { return m_vechelixsegs[pidSeg]; }
}
 
// 2011-05-06
void RecMdcKalTrack::setVecHelixSegs( const HelixSegRefVec& vechelixsegs, int pid ) {
  if ( pid >= 0 && pid < 5 )
  {
    m_vechelixsegs[pid] = vechelixsegs;
    pidSeg              = pid;
  }
  else std::cout << "RecMdcKalTrack::setVecHelixSegs error, wrong particle id" << std::endl;
}
 
RecMdcKalTrack::~RecMdcKalTrack() {}