#include <RecMdcTrack.h>

Inheritance diagram for RecMdcTrack:

Public Member Functions
virtual const CLID &	clID () const

	RecMdcTrack ()

	RecMdcTrack (const RecMdcTrack &)

	RecMdcTrack (const DstMdcTrack &)

RecMdcTrack &	operator= (const RecMdcTrack &)

RecMdcTrack &	operator= (const DstMdcTrack &)

	~RecMdcTrack ()

const double	getVX0 () const

const double	getVY0 () const

const double	getVZ0 () const

const int	getNhits () const

const double	getFiTerm () const

const HepPoint3D &	getPivot () const

const HitRefVec	getVecHits (void) const

const int	getNcluster () const

const ClusterRefVec	getVecClusters () const

const double	getMatchChi2 () const

const map< int, int >	getClusterFitStat ()

void	setVX0 (double x0)

void	setVY0 (double y0)

void	setVZ0 (double z0)

void	setFiTerm (double fiterm)

void	setNhits (int nhits)

void	setPivot (const HepPoint3D &pivot)

void	setVecHits (HitRefVec vechits)

void	setNcluster (int ncluster)

void	setVecClusters (ClusterRefVec vecclusters)

void	setVecClusters (ClusterRefVec vecclusters, map< int, int > clusterFitStat)

void	setMatchChi2 (double MatchChi2)

void	setClusterFitStat (map< int, int > clusterFitStat)

void	updateStatistics ()

Public Member Functions inherited from DstMdcTrack
	DstMdcTrack ()

	DstMdcTrack (const DstMdcTrack &)

DstMdcTrack &	operator= (const DstMdcTrack &)

virtual	~DstMdcTrack ()

const HepVector	helix () const
	......

const HepSymMatrix	err () const

const HepLorentzVector	p4 (const double mass) const

const HepPoint3D	poca () const

const Hep3Vector	p3 () const

const HepPoint3D	x3 () const

const double	helix (const int i) const

const double	err (const int i) const

const int	trackId () const

const int	charge () const

const double	pxy () const

const double	px () const

const double	py () const

const double	pz () const

const double	p () const

const double	theta () const

const double	phi () const

const double	x () const

const double	y () const

const double	z () const

const double	r () const

const int	stat () const

const double	chi2 () const

const int	ndof () const

const int	nster () const

const int	firstLayer () const

const int	lastLayer () const

const int	nlayer () const

const int	nCgemXCluster () const

const int	nCgemVCluster () const

void	setHelix (double helix[5])

void	setPoca (double poca[3])

void	setError (double err[15])

void	setHelix (const HepVector &a)

void	setPoca (const HepPoint3D &poca)

void	setError (const HepSymMatrix &Ea)

void	setTrackId (const int trackId)

void	setCharge (const int charge)

void	setPxy (const double pxy)

void	setPx (const double px)

void	setPy (const double py)

void	setPz (const double pz)

void	setP (const double p)

void	setTheta (const double theta)

void	setPhi (const double phi)

void	setX (const double x)

void	setY (const double y)

void	setZ (const double z)

void	setR (const double r)

void	setStat (const int stat)

void	setChi2 (const double chi)

void	setNdof (const int ndof)

void	setNster (const int ns)

void	setFirstLayer (const int id)

void	setLastLayer (const int id)

void	setNlayer (const int nlayer)

void	setNCgemXCluster (const int n)

void	setNCgemVCluster (const int n)

Static Public Member Functions
static const CLID &	classID ()

Static Public Member Functions inherited from DstMdcTrack
static const CLID &	classID ()

Additional Inherited Members
Protected Attributes inherited from DstMdcTrack
int	m_trackId

int	m_charge

double	m_helix [5]

double	m_err [15]

double	m_poca [3]

HepPoint3D	m_poca3D

HepVector	m_a

HepSymMatrix	m_Ea

int	m_stat

double	m_chi2

int	m_ndof

int	m_nster

int	m_nCgemXClusters

int	m_nCgemVClusters

int	m_firstLayer

int	m_lastLayer

int	m_nlayer

double	m_pxy

double	m_px

double	m_py

double	m_pz

double	m_p

double	m_theta

double	m_phi

double	m_x

double	m_y

double	m_z

double	m_r

Detailed Description

Definition at line 30 of file RecMdcTrack.h.

Constructor & Destructor Documentation

◆ RecMdcTrack() [1/3]

RecMdcTrack::RecMdcTrack ( )

Definition at line 20 of file RecMdcTrack.cxx.

  : m_pivot(0,0,0),
    m_fiTerm(0.0),           
    m_nhits(0),
    m_ncluster(0),
    m_matchChi2(9999.0)
   {   
   }

◆ RecMdcTrack() [2/3]

RecMdcTrack::RecMdcTrack ( const RecMdcTrack & rectrk )

Definition at line 29 of file RecMdcTrack.cxx.

  : DstMdcTrack(rectrk),
    m_pivot(rectrk.m_pivot),
    m_fiTerm(rectrk.m_fiTerm),
    m_nhits(rectrk.m_nhits),
    m_vechits(rectrk.m_vechits),
    m_ncluster(rectrk.m_ncluster),
    m_matchChi2(rectrk.m_matchChi2),
    m_vecclusters(rectrk.m_vecclusters)
  {
  }

◆ RecMdcTrack() [3/3]

RecMdcTrack::RecMdcTrack ( const DstMdcTrack & dstrk )

Definition at line 41 of file RecMdcTrack.cxx.

  : DstMdcTrack(dstrk),
    m_pivot(-99.,-99.,-99.),
    m_fiTerm(-99.),
    m_nhits(-99)
  {   
  }

◆ ~RecMdcTrack()

RecMdcTrack::~RecMdcTrack ( )

Definition at line 77 of file RecMdcTrack.cxx.

77 {

78 }

Member Function Documentation

◆ classID()

const CLID & RecMdcTrack::classID ( )

static

Definition at line 16 of file RecMdcTrack.cxx.

                                   {
    return CLID_RecMdcTrack;
  }

Referenced by clID().

◆ clID()

const CLID & RecMdcTrack::clID ( ) const

virtual

Reimplemented from DstMdcTrack.

Definition at line 12 of file RecMdcTrack.cxx.

                                      {
    return RecMdcTrack::classID();
  }

◆ getClusterFitStat()

const map< int, int > RecMdcTrack::getClusterFitStat ( )

inline

Definition at line 72 of file RecMdcTrack.h.

72{return m_clusterFitStat;}

◆ getFiTerm()

const double RecMdcTrack::getFiTerm ( ) const

inline

Definition at line 59 of file RecMdcTrack.h.

59{ return m_fiTerm; }

Referenced by ExtMdcTrack::ReadTrk().

◆ getMatchChi2()

const double RecMdcTrack::getMatchChi2 ( ) const

inline

Definition at line 71 of file RecMdcTrack.h.

71{return m_matchChi2; }

◆ getNcluster()

const int RecMdcTrack::getNcluster ( ) const

inline

Definition at line 69 of file RecMdcTrack.h.

69{ return m_ncluster; }

◆ getNhits()

const int RecMdcTrack::getNhits ( ) const

inline

Definition at line 56 of file RecMdcTrack.h.

56{ return m_nhits; }

Referenced by MdcMergeDups::dumpRecMdcTrack(), MdcTrkRecon::dumpTdsTrack(), TTrackManager::makeTds(), MdcPrintSvc::printRecMdcTrack(), and MdcDedxTrk::set_ExTrk().

◆ getPivot()

const HepPoint3D & RecMdcTrack::getPivot ( ) const

inline

Definition at line 63 of file RecMdcTrack.h.

63{ return m_pivot; }

Referenced by ExtMdcTrack::ReadTrk().

◆ getVecClusters()

const ClusterRefVec RecMdcTrack::getVecClusters ( ) const

inline

Definition at line 70 of file RecMdcTrack.h.

70{ return m_vecclusters;}

Referenced by TestTrack::execute(), and CgemMdcFitAlg::fit().

◆ getVecHits()

const HitRefVec RecMdcTrack::getVecHits ( void ) const

inline

Definition at line 67 of file RecMdcTrack.h.

67{ return m_vechits; }

Referenced by MdcMergeDups::doMergeCurl(), MdcMergeDups::dumpRecMdcTrack(), MdcTrkRecon::dumpTdsTrack(), MdcxCosmicSewer::execute(), CgemMdcFitAlg::fit(), TTrackManager::makeTds(), MdcPrintSvc::printRecMdcTrack(), MdcMergeDups::testByOverlapHit(), RecMdcTrackCnv::TObjectToDataObject(), and CgemMdcFitAlg::updateTracks().

◆ getVX0()

const double RecMdcTrack::getVX0 ( ) const

inline

Definition at line 50 of file RecMdcTrack.h.

50{ return m_pivot[0]; }

Referenced by MdcMergeDups::testByOverlapHit().

◆ getVY0()

const double RecMdcTrack::getVY0 ( ) const

inline

Definition at line 51 of file RecMdcTrack.h.

51{ return m_pivot[1]; }

Referenced by MdcMergeDups::testByOverlapHit().

◆ getVZ0()

const double RecMdcTrack::getVZ0 ( ) const

inline

Definition at line 52 of file RecMdcTrack.h.

52{ return m_pivot[2]; }

Referenced by MdcMergeDups::testByOverlapHit().

◆ operator=() [1/2]

RecMdcTrack & RecMdcTrack::operator= ( const DstMdcTrack & dstrk )

Definition at line 64 of file RecMdcTrack.cxx.

  {
    if( &dstrk != this ){
      DstMdcTrack::operator=(dstrk);  
      m_pivot[0] = -99.;
      m_pivot[1] = -99.;
      m_pivot[2] = -99.;
      m_fiTerm   = -99.;
      m_nhits    = -99;
    }
    return  *this;
  }  

◆ operator=() [2/2]

RecMdcTrack & RecMdcTrack::operator= ( const RecMdcTrack & rectrk )

Definition at line 49 of file RecMdcTrack.cxx.

  {
    if( &rectrk != this ){
      DstMdcTrack::operator=(rectrk);
      m_pivot = rectrk.m_pivot;
      m_fiTerm = rectrk. m_fiTerm;
      m_nhits = rectrk. m_nhits;
      m_vechits = rectrk.m_vechits;
      m_ncluster = rectrk.m_ncluster;
      m_matchChi2 = rectrk.m_matchChi2;
      m_vecclusters = rectrk.m_vecclusters;
    }
    return  *this;
  }

◆ setClusterFitStat()

void RecMdcTrack::setClusterFitStat ( map< int, int > clusterFitStat )

inline

Definition at line 88 of file RecMdcTrack.h.

88{m_clusterFitStat = clusterFitStat;}

◆ setFiTerm()

void RecMdcTrack::setFiTerm ( double fiterm )

inline

Definition at line 77 of file RecMdcTrack.h.

77{ m_fiTerm = fiterm; }

Referenced by TTrackManager::makeTds(), HoughFinder::storeRecTracks(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), RecMdcTrackCnv::TObjectToDataObject(), and CgemMdcFitAlg::updateTracks().

◆ setMatchChi2()

void RecMdcTrack::setMatchChi2 ( double MatchChi2 )

inline

Definition at line 87 of file RecMdcTrack.h.

87{ m_matchChi2 = MatchChi2; }

◆ setNcluster()

void RecMdcTrack::setNcluster ( int ncluster )

inline

Definition at line 83 of file RecMdcTrack.h.

83{ m_ncluster = ncluster;}

Referenced by CgemLineFit::Store_Trk(), MdcTrack::storeTrack(), updateStatistics(), and CgemMdcFitAlg::updateTracks().

◆ setNhits()

void RecMdcTrack::setNhits ( int nhits )

inline

Definition at line 78 of file RecMdcTrack.h.

78{ m_nhits = nhits; }

Referenced by TTrackManager::makeTds(), MdcTrack::storeTrack(), RecMdcTrackCnv::TObjectToDataObject(), updateStatistics(), and CgemMdcFitAlg::updateTracks().

◆ setPivot()

void RecMdcTrack::setPivot ( const HepPoint3D & pivot )

inline

Definition at line 79 of file RecMdcTrack.h.

79{ m_pivot = pivot; }

Referenced by TTrackManager::makeTds(), HoughFinder::storeRecTracks(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

◆ setVecClusters() [1/2]

void RecMdcTrack::setVecClusters ( ClusterRefVec vecclusters )

Definition at line 86 of file RecMdcTrack.cxx.

{
    m_vecclusters = vecclusters;
    updateStatistics();
}

Referenced by CgemLineFit::Store_Trk(), HoughFinder::storeRecTracks(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), RecMdcTrackCnv::TObjectToDataObject(), and CgemMdcFitAlg::updateTracks().

◆ setVecClusters() [2/2]

void RecMdcTrack::setVecClusters	(	ClusterRefVec	vecclusters,
		map< int, int >	clusterFitStat )

Definition at line 92 of file RecMdcTrack.cxx.

{
    m_vecclusters = vecclusters;
    m_clusterFitStat = clusterFitStat;
    updateStatistics();
}

◆ setVecHits()

void RecMdcTrack::setVecHits ( HitRefVec vechits )

Definition at line 80 of file RecMdcTrack.cxx.

{
    m_vechits = vechits;
    updateStatistics();
}

Referenced by TTrackManager::makeTds(), HoughFinder::storeRecTracks(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), RecMdcTrackCnv::TObjectToDataObject(), and CgemMdcFitAlg::updateTracks().

◆ setVX0()

void RecMdcTrack::setVX0 ( double x0 )

inline

Definition at line 74 of file RecMdcTrack.h.

74{ m_pivot[0] = x0; }

Referenced by TTrackManager::makeTds(), HoughFinder::storeRecTracks(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), RecMdcTrackCnv::TObjectToDataObject(), and CgemMdcFitAlg::updateTracks().

◆ setVY0()

void RecMdcTrack::setVY0 ( double y0 )

inline

Definition at line 75 of file RecMdcTrack.h.

75{ m_pivot[1] = y0; }

Referenced by TTrackManager::makeTds(), HoughFinder::storeRecTracks(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), RecMdcTrackCnv::TObjectToDataObject(), and CgemMdcFitAlg::updateTracks().

◆ setVZ0()

void RecMdcTrack::setVZ0 ( double z0 )

inline

Definition at line 76 of file RecMdcTrack.h.

76{ m_pivot[2] = z0; }

Referenced by TTrackManager::makeTds(), HoughFinder::storeRecTracks(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), RecMdcTrackCnv::TObjectToDataObject(), and CgemMdcFitAlg::updateTracks().

◆ updateStatistics()

void RecMdcTrack::updateStatistics ( )

Definition at line 99 of file RecMdcTrack.cxx.

{
    set<int> mdc_layers;
    set<int> cgem_layers;
 
    int nDof = 0;
    int nSter = 0;
    int maxLayerId = -1;
    int minLayerId = 43;
    int nLayer = 0;
    double fiTerm = 999.;
 
    MdcID* mdcID = new MdcID();
    int nHits = m_vechits.size();
    HitRefVec::iterator iterHit = m_vechits.begin();
    for( ; iterHit != m_vechits.end(); iterHit++){
        Identifier mdcid = (*iterHit)->getMdcId();
        int layerId = MdcID::layer(mdcid);
        //count nDof
        //cout<< (*iterHit)->getStat()<<endl;
        if ((*iterHit)->getStat()==1){
            nDof++;
        }
        //count nSter
        //cout<<layerId<<"  "<<mdcID->is_axial(mdcid)<<endl;
        if (!(mdcID->is_axial(mdcid))){
        //if (layerId >=20 &&layerId <36){
            nSter++;
        }
        // find maxLayerId and minLayerId
        if (layerId >= maxLayerId){
            maxLayerId = layerId;
            //fiTerm=(1./sqrt(1.+m_helix[4]*m_helix[4]))*(*iterHit)->getFltLen()*m_helix[3]/333.567;
        }
        if (layerId < minLayerId){
            minLayerId = layerId;
        }
        //count nlayer
        mdc_layers.insert(layerId);
        //count nHits
        //nHits++;
    }
 
    //cout<<__FILE__ <<"  "<<__LINE__ << endl;
    //cout<<"nDof "<<nDof-5<<endl;
    //cout<<"nSter "<<nSter<<endl;
    //cout<<"maxLayerId "<<maxLayerId<<endl;
    //cout<<"minLayerId "<<minLayerId<<endl;
    //cout<<"nHits "<<nHits<<endl;
    //cout<<endl;
 
    int nCgemXClusters=0;
    int nCgemVClusters=0; 
    int nCluster=m_vecclusters.size();
    ClusterRefVec::iterator iterCluster = m_vecclusters.begin();
    for( ; iterCluster != m_vecclusters.end(); iterCluster++){
        (*iterCluster)->setTrkId(trackId());
        int flag = (*iterCluster)->getflag();
        int clusterid = (*iterCluster)->getclusterid();
        map<int,int>::iterator iter = m_clusterFitStat.find(clusterid);
        int stat =1;
        if(iter != m_clusterFitStat.end())stat = iter->second;
        switch(flag){
            case 0:
                nCgemXClusters++;
                if(stat)nDof++;
                break;
            case 1:
                nCgemVClusters++;
                if(stat)nDof++;
                nSter++;
                break;
            case 2:
                nCgemXClusters++;
                nCgemVClusters++;
                if(stat)nDof+=2;
                nSter++;
                break;
            case 3:
                nCgemXClusters++;
                nCgemVClusters++;
                if(stat)nDof+=2;
                nSter++;
                break;
            default: break;
        }
        int layerId = (*iterCluster)->getlayerid();
        if (layerId >= maxLayerId){
            maxLayerId = layerId;
            //fiTerm=(1./sqrt(1.+m_helix[4]*m_helix[4]))*(*iterHit)->getFltLen()*m_helix[3]/333.567;
        }
        if (layerId < minLayerId){
            minLayerId = layerId;
        }
        cgem_layers.insert(layerId);
        nCluster++;
    }
 
    nLayer = mdc_layers.size() + cgem_layers.size();
 
    //cout<<__FILE__ <<"  "<<__LINE__ << endl;
    //cout<<"nDof "<<nDof-5<<endl;
    //cout<<"nSter "<<nSter<<endl;
    //cout<<"maxLayerId "<<maxLayerId<<endl;
    //cout<<"minLayerId "<<minLayerId<<endl;
    //cout<<"nCluster "<<nCluster<<endl;
    //cout<<"nLayer "<<nLayer<<endl;
    //cout<<endl;
 
    setNdof(nDof-5);
    setNster(nSter);
    setFirstLayer(minLayerId);
    setLastLayer(maxLayerId);
    setNlayer(nLayer);
    setNCgemXCluster(nCgemXClusters);
    setNCgemVCluster(nCgemVClusters);
    setNhits(nHits);
    setNcluster(nCluster);
 
    //cout<<"mdc nlayer:"<<mdc_layers.size()<<endl;
    //cout<<"cgem nlayer:"<<cgem_layers.size()<<endl;
}

Referenced by setVecClusters(), setVecClusters(), and setVecHits().

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RecMdcTrack() [1/3]

◆ RecMdcTrack() [2/3]

◆ RecMdcTrack() [3/3]

◆ ~RecMdcTrack()

Member Function Documentation

◆ classID()

◆ clID()

◆ getClusterFitStat()

◆ getFiTerm()

◆ getMatchChi2()

◆ getNcluster()

◆ getNhits()

◆ getPivot()

◆ getVecClusters()

◆ getVecHits()

◆ getVX0()

◆ getVY0()

◆ getVZ0()

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ setClusterFitStat()

◆ setFiTerm()

◆ setMatchChi2()

◆ setNcluster()

◆ setNhits()

◆ setPivot()

◆ setVecClusters() [1/2]

◆ setVecClusters() [2/2]

◆ setVecHits()

◆ setVX0()

◆ setVY0()

◆ setVZ0()

◆ updateStatistics()