Ext_errmx Class Reference

#include <Ext_errmx.h>

Inheritance diagram for Ext_errmx:

Public Member Functions
	Ext_errmx ()
	Ext_errmx (const Ext_errmx &errmx)
	Ext_errmx (const HepSymMatrix &err)
virtual	~Ext_errmx ()
void	put_err (const double error[])
void	put_err (const HepSymMatrix &err)
const HepSymMatrix &	get_err () const
double	get_plane_err (const Hep3Vector &np, const Hep3Vector &nr) const
const HepVector &	get_plane_errs (const Hep3Vector &np, const Hep3Vector &nr, const Hep3Vector &nt) const
const Hep3Vector *	get_tvs (const int view, const Hep3Vector &pv) const
const Hep3Vector *	get_tvs (const Hep3Vector &pv) const
bool	valid (bool msg) const
Ext_errmx &	operator= (const Ext_errmx &errmx)

Protected Member Functions
void	set_plane_errs (const Hep3Vector &nx, const Hep3Vector &ny, const Hep3Vector &nz) const

Protected Attributes
HepSymMatrix	m_err

Friends
std::ostream &	operator<< (std::ostream &s, const Ext_errmx &errmx)

Detailed Description

Definition at line 36 of file Ext_errmx.h.

Constructor & Destructor Documentation

Ext_errmx() [1/3]

Ext_errmx::Ext_errmx

(

)

Definition at line 29 of file Ext_errmx.cxx.

                     : m_err(Ndim_err,0), m_err3(3,0), m_R(3,3,0), 
    m_err2(2), m_valid(0) {}

Ext_errmx() [2/3]

Ext_errmx::Ext_errmx

(

const Ext_errmx &

errmx

)

Definition at line 34 of file Ext_errmx.cxx.

                                           : m_err(err.m_err),
    m_err3(err.m_err3), m_R(err.m_R), m_err2(err.m_err2), 
    m_valid(err.m_valid) {}

Ext_errmx() [3/3]

Ext_errmx::Ext_errmx

(

const HepSymMatrix &

err

)

Definition at line 38 of file Ext_errmx.cxx.

                                              : m_err(err), 
    m_err3(3,0), m_R(3,3,0), m_err2(2), m_valid(1) {}

~Ext_errmx()

virtual Ext_errmx::~Ext_errmx ( )

inlinevirtual

Definition at line 49 of file Ext_errmx.h.

{};

Member Function Documentation

get_err()

const HepSymMatrix & Ext_errmx::get_err ( ) const

inline

Definition at line 157 of file Ext_errmx.h.

                                                     {
  return m_err;
}

Referenced by ExtSteppingAction::UserSteppingAction().

get_plane_err()

double Ext_errmx::get_plane_err	(	const Hep3Vector &	np,
		const Hep3Vector &	nr
	)		const

Definition at line 85 of file Ext_errmx.cxx.

{
 
// Check the validity of the error matrix.
 
  if( !(*this).valid( 0 ) ){
    //std::cout << "%WARNING at Ext_get_plane_err: You are trying to calculate"
//  << " error \n           using an invalid error matrix." << std::endl;
  }
 
// Construct 3 TCS axes.
 
// x: nx --- unit vector which is perp to np and on the plane of nr and np.
// y: nz x nx
// z: nz( = np ) --- direction of the track.
 
  double nr_np( nr*np );
  double denom( 1.0 - nr_np*nr_np );
  double error( 0.0 );
 
  if( denom > Eps ){    // Track is not parallel to the readout direction.
 
    double fac( 1.0 / sqrt( denom ) );
    Hep3Vector nx( ( nr - nr_np * np ) * fac );
    Hep3Vector ny( np.cross( nx ) );
 
    (*this).set_plane_errs( nx, ny, np );
 
    double sigma2( m_err3( 1, 1 ) );    // Error along nx.
    if( sigma2 > 0 ){
      error = sqrt( sigma2 ) * fac; // Error projection.
    }
 
  } else {  // Track is parallel to the readout direction.
 
    error = Infinite;
 
  }
  return( error );
}

Referenced by ExtSteppingAction::UserSteppingAction().

get_plane_errs()

const HepVector & Ext_errmx::get_plane_errs	(	const Hep3Vector &	np,
		const Hep3Vector &	nr,
		const Hep3Vector &	nt
	)		const

Get the 2D projected error vector on the detector plane. \begin{itemize} \item np: track direction unit vector. \item nr: readout direction unit vector. \item nt: ( np x nr ) unit vector. \item HepVector(1)= error(sigma) along the nr. \item HepVector(2)= error(sigma) along the nt. \end{itemize}

Definition at line 140 of file Ext_errmx.cxx.

{
 
// Check the validity of the error matrix.
 
  if( !(*this).valid( 1 ) ){
//    std::cout << "%WARNING at Ext_get_plane_errs: You are trying to calculate"
//  << " error \n           using an invalid error matrix." << std::endl;
  }
 
  double nr_np( nr*np );
  double denom_r( 1.0 - nr_np*nr_np );
 
  if( denom_r > Eps ){  // nr is not parallel to the track direction: np.
 
    double nt_np( nt*np );
    double denom_t( 1.0 - nt_np*nt_np );
    double fac_r( 1.0 / sqrt( denom_r ) );
    Hep3Vector nx( ( nr - nr_np * np ) * fac_r );
    Hep3Vector ny( np.cross( nx ) );
 
    (*this).set_plane_errs( nx, ny, np );
 
    double sigma2( m_err3( 1, 1 ) );    // Error along nx.
    if( sigma2 > 0 ){
      m_err2( 1 ) = sqrt( sigma2 ) * fac_r; // Error projection.
    } else {
      m_err2( 1 ) = 0.0;
    }
 
    if( denom_t > Eps ){   // nt is not parallel to the track direction: np.
      double fac_t( 1.0 / sqrt( denom_t ) );
      sigma2 = m_err3( 2, 2 );
      if( sigma2 > 0 ){
    m_err2( 2 ) = sqrt( sigma2 ) * fac_t;   // Error projection.
      } else {
    m_err2( 2 ) = 0.0;
      }
    } else {    // nt is parallel to the track direction: np.
      m_err2( 2 ) = (*this).get_plane_err( np, nt );
    }
  } else {  // nr is parallel to the track direction: np.
    m_err2( 1 ) = (*this).get_plane_err( np, nr );
    m_err2( 2 ) = (*this).get_plane_err( np, nt );
  }
  return( m_err2 );
}

get_tvs() [1/2]

const Hep3Vector * Ext_errmx::get_tvs

(

const Hep3Vector &

pv

)

const

Get the 2D projected track error perpendicular to a given vector at the current point. pv: a given vector, momentum vector, for example. \begin{itemize} \item Hep3Vector(1)= error(1 sigma) along the direction which is perpendicular to pv on the plane formed by pv and z-axis. \item Hep3Vector(2)= error(1 sigma) along the direction which is (pv) x (HepVector(1) direction). \end{itemize}

Definition at line 278 of file Ext_errmx.cxx.

{
 
  Hep3Vector    np( pv.unit() );
  Hep3Vector    nz( 0.0, 0.0, 1.0 );
  Hep3Vector    nt( (nz.cross(np)).unit() );
  Hep3Vector    nr( nt.cross(np) );
 
  const HepVector & err_v = (*this).get_plane_errs( np, nr, nt );
  *m_nv     = err_v[0]*nr;
  *(m_nv+1) = err_v[1]*nt;
  return( m_nv );
 
}

get_tvs() [2/2]

const Hep3Vector * Ext_errmx::get_tvs	(	const int	view,
		const Hep3Vector &	pv
	)		const

Get 2D projected track error perpendicular to a given vector at the current point. pv: a given vector, momentum vector for example. \begin{itemize} \item view=1. \begin{itemize} \item Hep3Vector(1)= error(1 sigma) along the direction which is perpendicular to pv on the xy plane. \item Hep3Vector(2)= error(1 sigma) along the direction which is (pv) x (HepVector(1) direction). \end{itemize} \item view=2. \begin{itemize} \item Hep3Vector(1)= error(1 sigma) along the direction which is perpendicular to pv on the zy plane. \item Hep3Vector(2)= error(1 sigma) along the direction which is (pv) x (HepVector(1) direction). \end{itemize} \item view=3. \begin{itemize} \item Hep3Vector(1)= error(1 sigma) along the direction which is perpendicular to pv on the zx plane. \item Hep3Vector(2)= error(1 sigma) along the direction which is (pv) x (HepVector(1) direction). \end{itemize} \end{itemize}

Definition at line 215 of file Ext_errmx.cxx.

{
 
  Hep3Vector    np( pv.unit() );
  Hep3Vector    nr;
 
  switch( view ){
 
  case 1:   // xy view
 
    if( np.x() != 0 || np.y() != 0 ){
      nr.setX(  np.y() );
      nr.setY( -np.x() );
      nr = nr.unit();
    } else {        // Pointing to z-direction.
      nr.setX( 1 );
    }
    break;
 
  case 2:   // zy view
 
    if( np.y() != 0 || np.z() != 0 ){
      nr.setY( -np.z() );
      nr.setZ(  np.y() );
      nr = nr.unit();
    } else {        // Pointing to x-direction.
      nr.setZ( 1 );
    }
    break;
 
  case 3:   // zx view
 
    if( np.z() != 0 || np.x() != 0 ){
      nr.setX( -np.z() );
      nr.setZ(  np.x() );
      nr = nr.unit();
    } else {        // Pointing to z-direction.
      nr.setZ( 1 );
    }
    break;
  } /* End of switch */
 
  Hep3Vector nt( np.cross( nr ) );
  const HepVector & err_v = (*this).get_plane_errs( np, nr, nt );
  *m_nv     = err_v[0]*nr;
  *(m_nv+1) = err_v[1]*nt;
  return( m_nv );
}

operator=()

Ext_errmx & Ext_errmx::operator=

(

const Ext_errmx &

errmx

)

Definition at line 316 of file Ext_errmx.cxx.

{
  if( this != &err ){
    m_err   = err.m_err;
    m_err3  = err.m_err3;
    m_err2  = err.m_err2;
    m_R     = err.m_R;
    m_valid = err.m_valid;
    *m_nv   = *err.m_nv;
    *(m_nv+1)   = *(err.m_nv+1);
  }
  return *this;
}

put_err() [1/2]

void Ext_errmx::put_err

(

const double

error[]

)

Definition at line 45 of file Ext_errmx.cxx.

{
  int ne = 0;
 
  for( int i = 1; i <= m_err.num_col(); i++ ){
    for( int j = 1; j <= i; j++ ){
      m_err.fast( i, j ) = error[ ne++ ];
    }
  }
  m_valid = 1;
}

Referenced by Ext_xp_err::set_err(), and ExtSteppingAction::UserSteppingAction().

put_err() [2/2]

void Ext_errmx::put_err ( const HepSymMatrix &  err )

inline

Definition at line 148 of file Ext_errmx.h.

{
  m_err = err;
  m_valid = 1;
}

set_plane_errs()

void Ext_errmx::set_plane_errs ( const Hep3Vector &  nx, const Hep3Vector &  ny, const Hep3Vector &  nz  ) const

protected

Definition at line 60 of file Ext_errmx.cxx.

{
// Setup the rotation matrix.
 
  m_R( 1, 1 ) = nx.x();
  m_R( 1, 2 ) = nx.y();
  m_R( 1, 3 ) = nx.z();
  m_R( 2, 1 ) = ny.x();
  m_R( 2, 2 ) = ny.y();
  m_R( 2, 3 ) = ny.z();
  m_R( 3, 1 ) = nz.x();
  m_R( 3, 2 ) = nz.y();
  m_R( 3, 3 ) = nz.z();
 
// Get the 3x3 sub matrix and Rotate the error matrix.
 
  m_err3 = (m_err.sub( 1, 3 )).similarity( m_R );
 
}

valid()

bool Ext_errmx::valid

(

bool

msg

)

const

Definition at line 300 of file Ext_errmx.cxx.

{
  if( m_valid ){
    if( Ext_err_valid( msg, m_err, 6 ) ){
      return 1;
    } else {
      m_valid = 0;
      return 0;
    }
  } else {
    return 0;
  }
}

Friends And Related Function Documentation

operator<<

std::ostream & operator<< ( std::ostream &  s, const Ext_errmx &  errmx  )

friend

Definition at line 333 of file Ext_errmx.cxx.

{
  s << " m_valid: " << err.m_valid << '\n'
    << "m_err: " << err.m_err << " m_err3: " << err.m_err3 
    << " m_R: " << err.m_R << " m_err2: " << err.m_err2
    << " *m_nv: " << *err.m_nv 
    << " *(m_nv+1): " << *(err.m_nv+1)
    << std::endl;
  return s;
}

Member Data Documentation

m_err

HepSymMatrix Ext_errmx::m_err

protected

Definition at line 133 of file Ext_errmx.h.

Referenced by get_err(), Ext_xp_err::move(), operator=(), put_err(), set_plane_errs(), and valid().

---

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

• Ext_errmx.h
• Ext_errmx.cxx