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 33 of file Ext_errmx.h.

Constructor & Destructor Documentation

Ext_errmx() [1/3]

Ext_errmx::Ext_errmx

(

)

Definition at line 28 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 ) {}

Referenced by Ext_errmx(), Ext_xp_err::Ext_xp_err(), operator<<, Ext_xp_err::operator<<, operator=(), Ext_xp_err::operator=(), and Ext_xp_err::set_err().

Ext_errmx() [2/3]

Ext_errmx::Ext_errmx

(

const Ext_errmx &

errmx

)

Definition at line 33 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 40 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 46 of file Ext_errmx.h.

{};

Member Function Documentation

get_err()

const HepSymMatrix & Ext_errmx::get_err ( ) const

inline

Definition at line 149 of file Ext_errmx.h.

{ return m_err; }

get_plane_err()

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

Definition at line 83 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 );
}

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 294 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 222 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 330 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 47 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().

put_err() [2/2]

void Ext_errmx::put_err ( const HepSymMatrix & err )

inline

Definition at line 141 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 315 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 Symbol Documentation

operator<<

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

friend

Definition at line 347 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 128 of file Ext_errmx.h.

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

---

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

• Ext_errmx.h
• Ext_errmx.cxx