#include <CurvedLine.h>

Collaboration diagram for MuonCalib::CurvedLine:

Public Types
using	MdtHitVec = MuonCalibSegment::MdtHitVec
using	MdtHitPtr = MuonCalibSegment::MdtHitPtr

Public Member Functions
	CurvedLine ()
	Default constructor: a straight line through (0,0,0) pointing in in the local x direction of the chambers is created.
	CurvedLine (std::vector< Amg::Vector3D > &points)
	Constructor.
	CurvedLine (std::vector< Amg::Vector3D > &points, std::vector< Amg::Vector3D > &x_and_y_errors)
	Constructor.
Amg::Vector3D	getPointOnLine (const double loc_z) const
	get the point on the line a the local z coordinate "loc_z"
MTStraightLine	getTangent (const double loc_z) const
	get the tangent to the line a the local z coordinate "loc_z"
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 - 3.
void	setUsedHits (const MdtHitVec &hits)
const MdtHitVec &	trackHits () const

Private Member Functions
void	init (std::vector< Amg::Vector3D > &points, std::vector< Amg::Vector3D > &x_and_y_errors)

Private Attributes
Amg::VectorX	m_coeff_xz
Amg::VectorX	m_coeff_yz
double	m_chi2 {-1}
unsigned int	m_numTrkHits {0}
MdtHitVec	m_used_hits

Detailed Description

Definition at line 30 of file CurvedLine.h.

Member Typedef Documentation

◆ MdtHitPtr

using MuonCalib::CurvedLine::MdtHitPtr = MuonCalibSegment::MdtHitPtr

Definition at line 33 of file CurvedLine.h.

◆ MdtHitVec

using MuonCalib::CurvedLine::MdtHitVec = MuonCalibSegment::MdtHitVec

Definition at line 32 of file CurvedLine.h.

Constructor & Destructor Documentation

◆ CurvedLine() [1/3]

CurvedLine::CurvedLine ( )

Default constructor: a straight line through (0,0,0) pointing in in the local x direction of the chambers is created.

Definition at line 20 of file CurvedLine.cxx.

                       {
    // VARIABLES //
 
    std::vector<Amg::Vector3D> points(3);
    std::vector<Amg::Vector3D> errors(3);
    // FILL VARIABLES //
 
    points[0] = Amg::Vector3D(0.0, 0.0, 0.0);
    errors[0] = Amg::Vector3D(1.0, 1.0, 0.0);
    points[1] = Amg::Vector3D(0.0, 0.0, 1.0);
    errors[1] = Amg::Vector3D(1.0, 1.0, 0.0);
    points[2] = Amg::Vector3D(0.0, 0.0, 2.0);
    errors[2] = Amg::Vector3D(1.0, 1.0, 0.0);
    // INITIALIZATION //
 
    init(points, errors);
}

◆ CurvedLine() [2/3]

CurvedLine::CurvedLine ( std::vector< Amg::Vector3D > & points )

Constructor.

Parameters

points A straight line in the xz and a parabolic line in the yz plane is fitter through the given points. All points get equal weight.

Definition at line 52 of file CurvedLine.cxx.

                                                     {
    // VARIABLES //
 
    std::vector<Amg::Vector3D> errors(points.size(), Amg::Vector3D(1.0, 1.0, 0.0));
    // INITIALIZATION //
 
    init(points, errors);
}

◆ CurvedLine() [3/3]

CurvedLine::CurvedLine	(	std::vector< Amg::Vector3D > &	points,
		std::vector< Amg::Vector3D > &	x_and_y_errors )

Constructor.

Parameters

points	A straight line in the xz and a parabolic line in the yz plane is fitter through the given points.
x_and_y_errors	Errors of the x and y coordinates of the given points

Definition at line 72 of file CurvedLine.cxx.

72{ init(points, x_and_y_errors); }

Member Function Documentation

◆ chi2()

double CurvedLine::chi2 ( ) const

Definition at line 182 of file CurvedLine.cxx.

182{ return m_chi2; }

MuonCalib::CurvedLine::m_chi2

double m_chi2

Definition CurvedLine.h:82

◆ chi2PerDegreesOfFreedom()

double CurvedLine::chi2PerDegreesOfFreedom ( ) const

Return chi2 / number of TrackHits - 3.

Definition at line 186 of file CurvedLine.cxx.

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

MuonCalib::CurvedLine::m_numTrkHits

unsigned int m_numTrkHits

Definition CurvedLine.h:83

◆ getPointOnLine()

Amg::Vector3D CurvedLine::getPointOnLine ( const double loc_z ) const

get the point on the line a the local z coordinate "loc_z"

Definition at line 80 of file CurvedLine.cxx.

                                                                {
    // VARIABLES //
 
    double loc_x{0.0}, loc_y{0.0};
    // CALCULATE THE COORDINATES //
 
    for (int k = 0; k < m_coeff_xz.rows(); k++) { 
        loc_x += m_coeff_xz[k] * std::legendre(k, loc_z); 
    }
 
    for (int k = 0; k < m_coeff_yz.rows(); k++) { 
        loc_y += m_coeff_yz[k] * std::legendre(k, loc_z); 
    }
    // RETURN THE REQUESTED POINT //
 
    return Amg::Vector3D(loc_x, loc_y, loc_z);
}

◆ getTangent()

MTStraightLine CurvedLine::getTangent ( const double loc_z ) const

get the tangent to the line a the local z coordinate "loc_z"

Definition at line 112 of file CurvedLine.cxx.

                                                               {
    // VARIABLES //
 
    Amg::Vector3D null_vec(0.0, 0.0, 0.0);               // auxiliary 0 vector
    Amg::Vector3D point_1(getPointOnLine(loc_z));        // first point on the curved
                                                         // line
    Amg::Vector3D point_2(getPointOnLine(loc_z + 1.0));  // second point on the
                                                         // curved line
    // RETURN THE TANGENT //
 
    return MTStraightLine(point_1, point_2 - point_1, null_vec, null_vec);
}

◆ init()

void CurvedLine::init	(	std::vector< Amg::Vector3D > &	points,
		std::vector< Amg::Vector3D > &	x_and_y_errors )

private

Definition at line 136 of file CurvedLine.cxx.

                                                                                            {
    // CHECK THE NUMBER OF POINTS //
 
    if (points.size() < 3) {
        MsgStream log(Athena::getMessageSvc(), "CurvedLine");
        log << MSG::ERROR << "Class CurvedLine, method init: Not enough points given, must be at least 3 points!" << endmsg;
    }
    // VARIABLES //
 
    BaseFunctionFitter fitter;
    LegendrePolynomial legendre;                            // Legendre polynomial needed by the base
                                                            // function fitter
    std::vector<SamplePoint> sample_points(points.size());  // sample points needed
                                                            // by the base function
                                                            // fitter
                                                            // FILL THE VARIABLES AND PERFORM THE FITS //
 
    // xz plane //
    for (unsigned int k = 0; k < points.size(); k++) {
        sample_points[k].set_x1(points[k].z());
        sample_points[k].set_x2(points[k].x());
        sample_points[k].set_error(x_and_y_errors[k].x());
    }
    fitter.set_number_of_coefficients(2);
    fitter.fit_parameters(sample_points, 1, sample_points.size(), legendre);
    m_coeff_xz = fitter.coefficients();
 
    // yz plane //
    for (unsigned int k = 0; k < points.size(); k++) {
        sample_points[k].set_x1(points[k].z());
        sample_points[k].set_x2(points[k].y());
        sample_points[k].set_error(x_and_y_errors[k].y());
    }
    fitter.set_number_of_coefficients(3);
    fitter.fit_parameters(sample_points, 1, sample_points.size(), legendre);
    m_coeff_yz = fitter.coefficients();
 
}

◆ numberOfTrackHits()

unsigned int CurvedLine::numberOfTrackHits ( ) const

Definition at line 184 of file CurvedLine.cxx.

184{ return m_numTrkHits; }

◆ setChi2()

void CurvedLine::setChi2 ( double chi2 )

Cache the chi2.

Definition at line 181 of file CurvedLine.cxx.

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

MuonCalib::CurvedLine::chi2

double chi2() const

Definition CurvedLine.cxx:182

◆ setNumberOfTrackHits()

void CurvedLine::setNumberOfTrackHits ( unsigned int n_hits )

cache the number of track hits

Definition at line 183 of file CurvedLine.cxx.

183{ m_numTrkHits = n_hits; }

◆ setUsedHits()

void CurvedLine::setUsedHits ( const MdtHitVec & hits )

Definition at line 188 of file CurvedLine.cxx.

188{ m_used_hits = hits; }

MuonCalib::CurvedLine::m_used_hits

MdtHitVec m_used_hits

Definition CurvedLine.h:84

TRTCalib_Extractor.hits

hits

Definition TRTCalib_Extractor.py:35

◆ trackHits()

const CurvedLine::MdtHitVec & CurvedLine::trackHits ( ) const

Definition at line 189 of file CurvedLine.cxx.

189{ return m_used_hits; }

Member Data Documentation

◆ m_chi2

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

private

Definition at line 82 of file CurvedLine.h.

82{-1};

◆ m_coeff_xz

Amg::VectorX MuonCalib::CurvedLine::m_coeff_xz

private

Definition at line 77 of file CurvedLine.h.

◆ m_coeff_yz

Amg::VectorX MuonCalib::CurvedLine::m_coeff_yz

private

Definition at line 79 of file CurvedLine.h.

◆ m_numTrkHits

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

private

Definition at line 83 of file CurvedLine.h.

83{0};

◆ m_used_hits

MdtHitVec MuonCalib::CurvedLine::m_used_hits

private

Definition at line 84 of file CurvedLine.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

◆ CurvedLine() [1/3]

◆ CurvedLine() [2/3]

◆ CurvedLine() [3/3]

Member Function Documentation

◆ chi2()

◆ chi2PerDegreesOfFreedom()

◆ getPointOnLine()

◆ getTangent()

◆ init()

◆ numberOfTrackHits()

◆ setChi2()

◆ setNumberOfTrackHits()

◆ setUsedHits()

◆ trackHits()

Member Data Documentation

◆ m_chi2

◆ m_coeff_xz

◆ m_coeff_yz

◆ m_numTrkHits

◆ m_used_hits