#include <MTStraightLine.h>

Collaboration diagram for MuonCalib::MTStraightLine:

Public Types
using	MdtHitVec = MuonCalibSegment::MdtHitVec

using	MdtHitPtr = MuonCalibSegment::MdtHitPtr

Public Member Functions
	MTStraightLine ()
	Default constructor: all internal parameters of the straight line are set to 0. More...

	MTStraightLine (const Amg::Vector3D &r_position, const Amg::Vector3D &r_direction, const Amg::Vector3D &r_position_error, const Amg::Vector3D &r_direction_error)
	Constructor: r_position: position vector of the straight line; r_direction: direction vector of the straight line; r_position_error: error on the position vector; r_direction_error: error on the direction vector. More...

	MTStraightLine (const double r_a_x1, const double r_b_x1, const double r_a_x2, const double r_b_x2, const double r_a_x1_err, const double r_b_x1_err, const double r_a_x2_err, const double r_b_x2_err)
	Coordinates: x1, x2, x3. More...

Amg::Vector3D	positionVector () const
	get the position vector of the straight line More...

Amg::Vector3D	directionVector () const
	get the direction vector of the straight line More...

Amg::Vector3D	positionError () const
	get the error on the position vector of the straight line More...

Amg::Vector3D	directionError () const
	get the error on the direction vector of the straight line More...

double	a_x1 () const
	get the slope of the straight line in the x1-x3 plane More...

double	a_x1_error () const
	get the error on the slope of the straight line in the x1-x3 plane More...

double	b_x1 () const
	get the intercept of the straight line in the x1-x3 plane More...

double	b_x1_error () const
	get the error on the intercept of the straight line in the x1-x3 plane More...

double	a_x2 () const
	get the slope of the straight line in the x2-x3 plane More...

double	a_x2_error () const
	get the error on the slope of the straight line in the x2-x3 plane More...

double	b_x2 () const
	get the intercept of the straight line in the x2-x3 plane More...

double	b_x2_error () const
	get the slope of the intercept of the straight line in the x2-x3 plane More...

Amg::Vector3D	pointOnLine (const double lambda) const
	get the point on the line for the given scale factor lambda, point=position_vector+lambda*direction_vector More...

double	signDistFrom (const MTStraightLine &h) const
	get the signed distance of two lines (if both are parallel, dist>0) More...

double	distFromLine (const Amg::Vector3D &point) const
	get the distance of point point from straight line More...

void	setChi2 (double chi2)
	Cache the chi2. More...

double	chi2 () const

void	setNumberOfTrackHits (unsigned int n_hits)
	cache the number of track hits More...

unsigned int	numberOfTrackHits () const

double	chi2PerDegreesOfFreedom () const
	Return chi2 / number of TrackHits - 2. More...

void	setUsedHits (const MdtHitVec &hits)

const MdtHitVec &	trackHits () const

Private Member Functions
void	init ()

void	init (const Amg::Vector3D &r_position, const Amg::Vector3D &r_direction, const Amg::Vector3D &r_position_error, const Amg::Vector3D &r_direction_error)

void	init (const double r_a_x1, const double r_b_x1, const double r_a_x2, const double r_b_x2, const double r_a_x1_err, const double r_b_x1_err, const double r_a_x2_err, const double r_b_x2_err)

Private Attributes
Amg::Vector3D	m_position {0., 0., 0.}

Amg::Vector3D	m_direction {0., 0., 0.}

Amg::Vector3D	m_position_error {0., 0., 0.}

Amg::Vector3D	m_direction_error {0., 0., 0.}

double	m_chi2 {-1}

unsigned int	m_numTrkHits {0}

MdtHitVec	m_used_hits

Detailed Description

Definition at line 16 of file MTStraightLine.h.

Member Typedef Documentation

◆ MdtHitPtr

using MuonCalib::MTStraightLine::MdtHitPtr = MuonCalibSegment::MdtHitPtr

Definition at line 19 of file MTStraightLine.h.

◆ MdtHitVec

using MuonCalib::MTStraightLine::MdtHitVec = MuonCalibSegment::MdtHitVec

Definition at line 18 of file MTStraightLine.h.

Constructor & Destructor Documentation

◆ MTStraightLine() [1/3]

MuonCalib::MTStraightLine::MTStraightLine ( )

inline

Default constructor: all internal parameters of the straight line are set to 0.

Definition at line 22 of file MTStraightLine.h.

◆ MTStraightLine() [2/3]

MuonCalib::MTStraightLine::MTStraightLine	(	const Amg::Vector3D &	r_position,
		const Amg::Vector3D &	r_direction,
		const Amg::Vector3D &	r_position_error,
		const Amg::Vector3D &	r_direction_error
	)

inline

Constructor: r_position: position vector of the straight line; r_direction: direction vector of the straight line; r_position_error: error on the position vector; r_direction_error: error on the direction vector.

Definition at line 26 of file MTStraightLine.h.

◆ MTStraightLine() [3/3]

MuonCalib::MTStraightLine::MTStraightLine	(	const double	r_a_x1,
		const double	r_b_x1,
		const double	r_a_x2,
		const double	r_b_x2,
		const double	r_a_x1_err,
		const double	r_b_x1_err,
		const double	r_a_x2_err,
		const double	r_b_x2_err
	)

inline

Coordinates: x1, x2, x3.

Parametrization of the straight line: x1 = r_a_x1*x3 + r_b_x1; x2 = r_a_x2*x3 + r_b_x2; x3 arbitrary. r_a_x1_err: error on r_a_x1. r_b_x1_err: error on r_b_x1. r_a_x2_err: error on r_a_x2. r_b_x2_err: error on r_b_x2.

Definition at line 36 of file MTStraightLine.h.

Member Function Documentation

◆ a_x1()

double MTStraightLine::a_x1 ( ) const

get the slope of the straight line in the x1-x3 plane

Definition at line 46 of file MTStraightLine.cxx.

46 { return (m_direction.z() != 0.) ? m_direction.x() / m_direction.z() : 0.; }

◆ a_x1_error()

double MTStraightLine::a_x1_error ( ) const

get the error on the slope of the straight line in the x1-x3 plane

Definition at line 48 of file MTStraightLine.cxx.

                                         {
     return m_direction.z() != 0.
                ? std::hypot(m_direction_error.x() / m_direction_error.z(), m_direction_error.y() * a_x1() / m_direction_error.z())
                : 0.;
 }

◆ a_x2()

double MTStraightLine::a_x2 ( ) const

get the slope of the straight line in the x2-x3 plane

Definition at line 65 of file MTStraightLine.cxx.

65 { return (m_direction.z() == 0.0) ? 0. : m_direction.y() / m_direction.z(); }

◆ a_x2_error()

double MTStraightLine::a_x2_error ( ) const

get the error on the slope of the straight line in the x2-x3 plane

Definition at line 66 of file MTStraightLine.cxx.

                                         {
     return (m_direction.z() == 0)
                ? 0.0
                :
  
                std::hypot(m_direction_error.y() / m_direction_error.z(), m_direction_error.y() * a_x2() / m_direction_error.z());
 }

◆ b_x1()

double MTStraightLine::b_x1 ( ) const

get the intercept of the straight line in the x1-x3 plane

Definition at line 53 of file MTStraightLine.cxx.

                                   {
     if (m_direction.z() == 0.0) {
         MsgStream log(Athena::getMessageSvc(), "MTStraightLine");
         log << MSG::WARNING << "Class MTStraightLine, method b_x1: b_x1 not uniquely defined." << endmsg;
         return m_position.x();
     }
     return (m_position.x() - m_position.z() * a_x1());
 }

◆ b_x1_error()

double MTStraightLine::b_x1_error ( ) const

get the error on the intercept of the straight line in the x1-x3 plane

Definition at line 62 of file MTStraightLine.cxx.

                                         {
     return std::hypot(m_position_error.x(), a_x1() * m_position_error.z(), m_position.z() * a_x1_error());
 }

◆ b_x2()

double MTStraightLine::b_x2 ( ) const

get the intercept of the straight line in the x2-x3 plane

Definition at line 74 of file MTStraightLine.cxx.

                                   {
     if (m_direction.z() == 0.0) {
         MsgStream log(Athena::getMessageSvc(), "MTStraightLine");
         log << MSG::WARNING << "Class MTStraightLine, method b_x2: b_x2 not uniquely defined." << endmsg;
         return m_position.x();
     }
     return (m_position.y() - m_position.z() * a_x2());
 }

◆ b_x2_error()

double MTStraightLine::b_x2_error ( ) const

get the slope of the intercept of the straight line in the x2-x3 plane

Definition at line 83 of file MTStraightLine.cxx.

                                         {
     return std::hypot(m_position_error.y(), a_x2() * m_position_error.z(), m_position.z() * a_x2_error());
 }

◆ chi2()

double MTStraightLine::chi2 ( ) const

Definition at line 135 of file MTStraightLine.cxx.

135 { return m_chi2; }

◆ chi2PerDegreesOfFreedom()

double MTStraightLine::chi2PerDegreesOfFreedom ( ) const

Return chi2 / number of TrackHits - 2.

Definition at line 139 of file MTStraightLine.cxx.

139 { return m_chi2 / (m_numTrkHits > 2 ? m_numTrkHits - 2 : 0.01); }

◆ directionError()

Amg::Vector3D MTStraightLine::directionError ( ) const

get the error on the direction vector of the straight line

Definition at line 45 of file MTStraightLine.cxx.

45 { return m_direction_error; }

◆ directionVector()

Amg::Vector3D MTStraightLine::directionVector ( ) const

get the direction vector of the straight line

Definition at line 43 of file MTStraightLine.cxx.

43 { return m_direction; }

◆ distFromLine()

double MTStraightLine::distFromLine ( const Amg::Vector3D & point ) const

get the distance of point point from straight line

Definition at line 120 of file MTStraightLine.cxx.

                                                                   {
     //:::::::::::::::::::::::::
     //:: AUXILIARY VARIABLES ::
     //:::::::::::::::::::::::::
  
     Amg::Vector3D u = m_direction;
  
     //::::::::::::::::::::::::
     //:: CALCULATE DISTANCE ::
     //::::::::::::::::::::::::
  
     return std::sqrt((point - m_position).dot(point - m_position) -
                      ((point - m_position).dot(u.unit())) * ((point - m_position).dot(u.unit())));
 }

◆ init() [1/3]

void MTStraightLine::init ( )

private

Definition at line 14 of file MTStraightLine.cxx.

                           {
     m_position = Amg::Vector3D(0.0, 0.0, 0.0);
     m_direction = Amg::Vector3D(0.0, 0.0, 0.0);
     m_position_error = Amg::Vector3D(0.0, 0.0, 0.0);
     m_direction_error = Amg::Vector3D(0.0, 0.0, 0.0);
  
     return;
 }

◆ init() [2/3]

void MTStraightLine::init	(	const Amg::Vector3D &	r_position,
		const Amg::Vector3D &	r_direction,
		const Amg::Vector3D &	r_position_error,
		const Amg::Vector3D &	r_direction_error
	)

private

Definition at line 23 of file MTStraightLine.cxx.

                                                                 {
     m_position = r_position;
     m_direction = r_direction;
     m_position_error = r_position_error;
     m_direction_error = r_direction_error;
  
     return;
 }

◆ init() [3/3]

void MTStraightLine::init	(	const double	r_a_x1,
		const double	r_b_x1,
		const double	r_a_x2,
		const double	r_b_x2,
		const double	r_a_x1_err,
		const double	r_b_x1_err,
		const double	r_a_x2_err,
		const double	r_b_x2_err
	)

private

Definition at line 33 of file MTStraightLine.cxx.

                                                                                                         {
     m_position = Amg::Vector3D(r_b_x1, r_b_x2, 0.0);
     m_direction = Amg::Vector3D(r_a_x1, r_a_x2, 1.0);
     m_position_error = Amg::Vector3D(r_b_x1_err, r_b_x2_err, 0.0);
     m_direction_error = Amg::Vector3D(r_a_x1_err, r_a_x2_err, 1.0);
  
     return;
 }

◆ numberOfTrackHits()

unsigned int MTStraightLine::numberOfTrackHits ( ) const

Definition at line 137 of file MTStraightLine.cxx.

137 { return m_numTrkHits; }

◆ pointOnLine()

Amg::Vector3D MTStraightLine::pointOnLine ( const double lambda ) const

get the point on the line for the given scale factor lambda, point=position_vector+lambda*direction_vector

Definition at line 87 of file MTStraightLine.cxx.

87 { return m_position + lambda * m_direction; }

◆ positionError()

Amg::Vector3D MTStraightLine::positionError ( ) const

get the error on the position vector of the straight line

Definition at line 44 of file MTStraightLine.cxx.

44 { return m_position_error; }

◆ positionVector()

Amg::Vector3D MTStraightLine::positionVector ( ) const

get the position vector of the straight line

Definition at line 42 of file MTStraightLine.cxx.

42 { return m_position; }

◆ setChi2()

void MTStraightLine::setChi2 ( double chi2 )

Cache the chi2.

Definition at line 134 of file MTStraightLine.cxx.

134 { m_chi2 = std::isnan(chi2) ? -1 : chi2; }

◆ setNumberOfTrackHits()

void MTStraightLine::setNumberOfTrackHits ( unsigned int n_hits )

cache the number of track hits

Definition at line 136 of file MTStraightLine.cxx.

136 { m_numTrkHits = n_hits; }

◆ setUsedHits()

void MTStraightLine::setUsedHits ( const MdtHitVec & hits )

Definition at line 141 of file MTStraightLine.cxx.

141 { m_used_hits = hits; }

◆ signDistFrom()

double MTStraightLine::signDistFrom ( const MTStraightLine & h ) const

get the signed distance of two lines (if both are parallel, dist>0)

Definition at line 89 of file MTStraightLine.cxx.

                                                                  {
     //:::::::::::::::::::::::::
     //:: AUXILIARY VARIABLES ::
     //:::::::::::::::::::::::::
  
     Amg::Vector3D a = m_position, u = m_direction;
     // position and direction vectors of g
     Amg::Vector3D b = h.positionVector(), v = h.directionVector();
     // position and direction vectors of h
     Amg::Vector3D n;  // normal vector of plane spanned by g and h
     Amg::Vector3D d;  // distance vector
  
     //:::::::::::::
     //:: PROGRAM ::
     //:::::::::::::
  
     //::::::::::::::::::::::::
     //:: collinearity check ::
     //::::::::::::::::::::::::
  
     n = u.cross(v);
     d = a - b;
     if (n.dot(n) == 0.0) {
         // straight lines are parallel
         // no sign given
         return std::sqrt(d.dot(d) - (u.unit().dot(d)) * (u.unit().dot(d)));
     }
  
     return (d.dot(n.unit()));
 }

◆ trackHits()

const MTStraightLine::MdtHitVec & MTStraightLine::trackHits ( ) const

Definition at line 142 of file MTStraightLine.cxx.

142 { return m_used_hits; }

Member Data Documentation

◆ m_chi2

double MuonCalib::MTStraightLine::m_chi2 {-1}

private

Definition at line 123 of file MTStraightLine.h.

◆ m_direction

Amg::Vector3D MuonCalib::MTStraightLine::m_direction {0., 0., 0.}

private

Definition at line 117 of file MTStraightLine.h.

◆ m_direction_error

Amg::Vector3D MuonCalib::MTStraightLine::m_direction_error {0., 0., 0.}

private

Definition at line 121 of file MTStraightLine.h.

◆ m_numTrkHits

unsigned int MuonCalib::MTStraightLine::m_numTrkHits {0}

private

Definition at line 124 of file MTStraightLine.h.

◆ m_position

Amg::Vector3D MuonCalib::MTStraightLine::m_position {0., 0., 0.}

private

Definition at line 116 of file MTStraightLine.h.

◆ m_position_error

Amg::Vector3D MuonCalib::MTStraightLine::m_position_error {0., 0., 0.}

private

Definition at line 120 of file MTStraightLine.h.

◆ m_used_hits

MdtHitVec MuonCalib::MTStraightLine::m_used_hits

private

Definition at line 125 of file MTStraightLine.h.

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

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ MdtHitPtr

◆ MdtHitVec

Constructor & Destructor Documentation

◆ MTStraightLine() [1/3]

◆ MTStraightLine() [2/3]

◆ MTStraightLine() [3/3]

Member Function Documentation

◆ a_x1()

◆ a_x1_error()

◆ a_x2()

◆ a_x2_error()

◆ b_x1()

◆ b_x1_error()

◆ b_x2()

◆ b_x2_error()

◆ chi2()

◆ chi2PerDegreesOfFreedom()

◆ directionError()

◆ directionVector()

◆ distFromLine()

◆ init() [1/3]

◆ init() [2/3]

◆ init() [3/3]

◆ numberOfTrackHits()

◆ pointOnLine()

◆ positionError()

◆ positionVector()

◆ setChi2()

◆ setNumberOfTrackHits()

◆ setUsedHits()

◆ signDistFrom()

◆ trackHits()

Member Data Documentation

◆ m_chi2

◆ m_direction

◆ m_direction_error

◆ m_numTrkHits

◆ m_position

◆ m_position_error

◆ m_used_hits