---

TrTrack.cpp

Go to the documentation of this file.

// $Id: $ // Include files

// local
#include "Event/TrTrack.h"

//-----------------------------------------------------------------------------
// Implementation file for class : TrTrack
//
// 2004-10-07 : Eduardo Rodrigues
//-----------------------------------------------------------------------------


//=============================================================================
// Retrieve the position and momentum vectors and the corresponding
// 6D covariance matrix (pos:1->3,mom:4-6) of a track closest to the beam-line
//=============================================================================
StatusCode TrTrack::positionAndMomentum( HepPoint3D &pos,
                                         HepVector3D &mom,
                                         HepSymMatrix &cov6D )
{
  m_physicsState.positionAndMomentum( pos, mom, cov6D );

  return StatusCode::SUCCESS;
}

//=============================================================================
// Retrieve the 3D-position (+ errors) of a track closest to the beam-line
//=============================================================================
StatusCode TrTrack::position( HepPoint3D &pos,
                              HepSymMatrix &errPos )
{
  pos    = m_physicsState.position();
  errPos = m_physicsState.errPosition();

  return StatusCode::SUCCESS;
};

//=============================================================================
// Retrieve the slopes (dx/dz,dy/dz,1) of a track closest to the beam-line
//=============================================================================
StatusCode TrTrack::slopes( HepVector3D &slopes,
                            HepSymMatrix &errSlopes )
{
  slopes    = m_physicsState.slopes();
  errSlopes = m_physicsState.errSlopes();

  return StatusCode::SUCCESS;
};

//=============================================================================
// Retrieve the momentum of a track closest to the beam-line
//=============================================================================
double TrTrack::p()
{
  return m_physicsState.p();
};

//=============================================================================
// Retrieve the transverse momentum of a track closest to the beam-line
//=============================================================================
double TrTrack::pt()
{
  return m_physicsState.pt();
};

//=============================================================================
// Retrieve the momentum vector (+ errors) of a track closest to the beam-line
//=============================================================================
StatusCode TrTrack::momentum( HepVector3D &mom,
                              HepSymMatrix &errMom )
{
  mom    = m_physicsState.momentum();
  errMom = m_physicsState.errMomentum();

  return StatusCode::SUCCESS;
};

//=============================================================================
// Retrieve the 6D covariance matrix (x,y,z,px,py,pz)
// of a track closest to the beam-line
//=============================================================================
StatusCode TrTrack::posMomCovariance( HepSymMatrix &cov6D )
{
  cov6D = m_physicsState.posMomCovariance();

  return StatusCode::SUCCESS;
};

//=============================================================================
// Retrieve the pointer to the state closest to the given z-position
//=============================================================================
TrState* TrTrack::closestState( double z ) {
  double minDist = 999999999.;
  TrState* best = 0;
  for ( std::vector<TrState*>::iterator it = m_states.begin() ;
        m_states.end() != it; it++ ) {
    if ( minDist > fabs( z - (*it)->z() ) ) {
      minDist = fabs( z-(*it)->z() );
      best    = *it;
    }
  }
  return best;
}

//=============================================================================
// Retrieve the pointer to the state closest to the given plane
//=============================================================================
TrState* TrTrack::closestState( HepPlane3D &plane ) {
  double minDist = 999999999.;
  double dist;
  TrState* best = 0;
  for ( std::vector<TrState*>::iterator it = m_states.begin() ;
        m_states.end() != it; it++ ) {
    dist = plane.distance( ((*it) -> position()) );
    if ( minDist > dist ) {
      minDist = dist;
      best    = *it;
    }
  }
  return best;
}

//=============================================================================
// Clone the track
//=============================================================================
TrTrack* TrTrack::clone() const
{
  return new TrTrack(*this);
};

//=============================================================================

---