#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.
	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.
	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.
Amg::Vector3D	positionVector () const
	get the position vector of the straight line
Amg::Vector3D	directionVector () const
	get the direction vector of the straight line
Amg::Vector3D	positionError () const
	get the error on the position vector of the straight line
Amg::Vector3D	directionError () const
	get the error on the direction vector of the straight line
double	a_x1 () const
	get the slope of the straight line in the x1-x3 plane
double	a_x1_error () const
	get the error on the slope of the straight line in the x1-x3 plane
double	b_x1 () const
	get the intercept of the straight line in the x1-x3 plane
double	b_x1_error () const
	get the error on the intercept of the straight line in the x1-x3 plane
double	a_x2 () const
	get the slope of the straight line in the x2-x3 plane
double	a_x2_error () const
	get the error on the slope of the straight line in the x2-x3 plane
double	b_x2 () const
	get the intercept of the straight line in the x2-x3 plane
double	b_x2_error () const
	get the slope of the intercept of the straight line in the x2-x3 plane
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
double	signDistFrom (const MTStraightLine &h) const
	get the signed distance of two lines (if both are parallel, dist>0)
double	distFromLine (const Amg::Vector3D &point) const
	get the distance of point point from straight line
void	setChi2 (double chi2)
	Cache the chi2.
double	chi2 () const
void	setNumberOfTrackHits (unsigned int n_hits)
	cache the number of track hits
unsigned int	numberOfTrackHits () const
double	chi2PerDegreesOfFreedom () const
	Return chi2 / number of TrackHits - 2.
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.; }

MuonCalib::MTStraightLine::m_direction

Amg::Vector3D m_direction

Definition MTStraightLine.h:117

◆ 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; }

MuonCalib::MTStraightLine::m_chi2

double m_chi2

Definition MTStraightLine.h:123

◆ 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); }

MuonCalib::MTStraightLine::m_numTrkHits

unsigned int m_numTrkHits

Definition MTStraightLine.h:124

◆ 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; }

MuonCalib::MTStraightLine::chi2

double chi2() const

Definition MTStraightLine.cxx:135

◆ 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; }

MuonCalib::MTStraightLine::m_used_hits

MdtHitVec m_used_hits

Definition MTStraightLine.h:125

TRTCalib_Extractor.hits

hits

Definition TRTCalib_Extractor.py:35

◆ 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.

123{-1};

◆ m_direction

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

private

Definition at line 117 of file MTStraightLine.h.

117{0., 0., 0.}; // direction vector of the straight line

◆ m_direction_error

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

private

Definition at line 121 of file MTStraightLine.h.

121{0., 0., 0.}; // error on the direction vector

◆ m_numTrkHits

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

private

Definition at line 124 of file MTStraightLine.h.

124{0};

◆ m_position

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

private

Definition at line 116 of file MTStraightLine.h.

116{0., 0., 0.}; // position vector of the straight line

◆ m_position_error

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

private

Definition at line 120 of file MTStraightLine.h.

120{0., 0., 0.}; // error on the position vector

◆ 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