/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Mdc/MdcRecEvent/MdcRecEvent-00-05-14/src/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() {

}

---