MuonCalib::TrRelationLookUp Class Reference

#include <TrRelationLookUp.h>

Inheritance diagram for MuonCalib::TrRelationLookUp:

Collaboration diagram for MuonCalib::TrRelationLookUp:

Public Types
using	ParVec = std::vector< double >

Public Member Functions
	TrRelationLookUp (const IRtRelation &rtRelation)
virtual std::string	name () const override final
virtual std::optional< double >	driftTime (const double r) const override final
	Interface method for fetching the drift-time from the radius Returns a nullopt if the time is out of the boundaries. More...
virtual std::optional< double >	driftTimePrime (const double r) const override final
virtual std::optional< double >	driftTime2Prime (const double r) const override final
virtual double	minRadius () const override final
	Returns the minimum drift-radius. More...
virtual double	maxRadius () const override final
	Returns the maximum drift-radius. More...
virtual unsigned	nDoF () const override final
	Returns the number of degrees of freedom of the tr relation. More...
virtual std::string	typeName () const override final
double	getReducedR (const double r) const
	Maps the radius interval [minRadius;maxRadius] to [-1;1] where the minimal radius is on the lower end. More...
double	getReducedRPrime () const
	Returns the derivative of the reduced radisu w.r.t r. More...
unsigned int	nPar () const
const ParVec &	parameters () const
double	par (unsigned int index) const

Private Member Functions
double	getTFromR (const double radius, const IRtRelation &rtRelation) const

Private Attributes
std::vector< double >	m_times {}
std::vector< double >	m_radii {}
double	m_minRadius {}
double	m_maxRadius {}
ParVec	m_parameters {}

Detailed Description

Definition at line 14 of file TrRelationLookUp.h.

Member Typedef Documentation

ParVec

using MuonCalib::CalibFunc::ParVec = std::vector<double>

inherited

Definition at line 35 of file CalibFunc.h.

Constructor & Destructor Documentation

TrRelationLookUp()

MuonCalib::TrRelationLookUp::TrRelationLookUp

(

const IRtRelation &

rtRelation

)

Definition at line 14 of file TrRelationLookUp.cxx.

                                                                    : ITrRelation{ParVec{}} {
 
        constexpr std::size_t numPoints(100);
        double stepSize{rtRelation.radius(rtRelation.tUpper()) / (numPoints - 1)};
        m_times.resize(numPoints);
        m_radii.resize(numPoints);
 
        for(std::size_t i = 0; i < numPoints; ++i){
            m_radii[i] = i * stepSize;
            m_times[i] = getTFromR(m_radii[i], rtRelation);
        }
        m_minRadius = m_radii.front();
        m_maxRadius = m_radii.back();
    }

Member Function Documentation

driftTime()

std::optional< double > MuonCalib::TrRelationLookUp::driftTime ( const double  r ) const

finaloverridevirtual

Interface method for fetching the drift-time from the radius Returns a nullopt if the time is out of the boundaries.

Implements MuonCalib::ITrRelation.

Definition at line 31 of file TrRelationLookUp.cxx.

                                                                        {
        if(r < minRadius()) return m_times.front();
        if(r > maxRadius()) return m_times.back();
        //Find the first element that is not less than r
        auto it = std::lower_bound(m_radii.begin(), m_radii.end(), r);
        //Find the distance to the element right before this one, since r lies between idx and idx+1
        std::size_t idx = std::distance(m_radii.begin(), it) - 1;
        double radius1 = m_radii[idx];
        double radius2 = m_radii[idx + 1];
        if (radius1 == radius2) return m_times[idx];
        double time1 = m_times[idx];
        double time2 = m_times[idx + 1];
        return time1 + (time2 - time1) * (r - radius1) / (radius2 - radius1);
    }

driftTime2Prime()

std::optional< double > MuonCalib::TrRelationLookUp::driftTime2Prime ( const double  r ) const

finaloverridevirtual

Implements MuonCalib::ITrRelation.

Definition at line 59 of file TrRelationLookUp.cxx.

                                                                              {
        if(r < minRadius() || r > maxRadius()) return std::nullopt;
        return 0.0;
    }

driftTimePrime()

std::optional< double > MuonCalib::TrRelationLookUp::driftTimePrime ( const double  r ) const

finaloverridevirtual

Implements MuonCalib::ITrRelation.

Definition at line 46 of file TrRelationLookUp.cxx.

                                                                             {
        if(r < minRadius() || r > maxRadius()) return 0.;
        //do the same thing as above but return dt/dr
        auto it = std::lower_bound(m_radii.begin(), m_radii.end(), r);
        std::size_t idx = std::distance(m_radii.begin(), it) - 1;
        double radius1 = m_radii[idx];
        double radius2 = m_radii[idx + 1];
        if (radius1 == radius2) return std::nullopt;
        double time1 = m_times[idx];
        double time2 = m_times[idx + 1];
        return (time2 - time1) / (radius2 - radius1);
    }

getReducedR()

double MuonCalib::ITrRelation::getReducedR ( const double  r ) const

inlineinherited

Maps the radius interval [minRadius;maxRadius] to [-1;1] where the minimal radius is on the lower end.

Definition at line 42 of file ITrRelation.h.

                                                     {
                return  mapToUnitInterval(r, minRadius(), maxRadius());
            }

getReducedRPrime()

double MuonCalib::ITrRelation::getReducedRPrime ( ) const

inlineinherited

Returns the derivative of the reduced radisu w.r.t r.

Definition at line 46 of file ITrRelation.h.

                                            {
                return unitIntervalPrime(minRadius(), maxRadius());
            }

getTFromR()

double MuonCalib::TrRelationLookUp::getTFromR ( const double  radius, const IRtRelation &  rtRelation  ) const

private

Definition at line 67 of file TrRelationLookUp.cxx.

                                                                                               {
        double precision{0.001};
        double tMax{rtRelation.tUpper()};
        double tMin{rtRelation.tLower()};
 
        //Search for the drift time
        while (tMax - tMin > 0.1 and std::abs(rtRelation.radius(0.5 * (tMin + tMax)) - radius) > precision) {
            double midPoint = 0.5 * (tMin + tMax);
            if (rtRelation.radius(midPoint) > radius) {
                tMax = midPoint;
            } else {
                tMin = midPoint;
            }
        }
        return 0.5 * (tMin + tMax);
    }

maxRadius()

double MuonCalib::TrRelationLookUp::maxRadius ( ) const

finaloverridevirtual

Returns the maximum drift-radius.

Implements MuonCalib::ITrRelation.

Definition at line 65 of file TrRelationLookUp.cxx.

{ return m_maxRadius; }

minRadius()

double MuonCalib::TrRelationLookUp::minRadius ( ) const

finaloverridevirtual

Returns the minimum drift-radius.

Implements MuonCalib::ITrRelation.

Definition at line 64 of file TrRelationLookUp.cxx.

{ return m_minRadius; }

name()

std::string MuonCalib::TrRelationLookUp::name ( ) const

finaloverridevirtual

Implements MuonCalib::CalibFunc.

Definition at line 29 of file TrRelationLookUp.cxx.

{ return "TrRelationLookUp"; }

nDoF()

unsigned MuonCalib::TrRelationLookUp::nDoF ( ) const

finaloverridevirtual

Returns the number of degrees of freedom of the tr relation.

Implements MuonCalib::ITrRelation.

Definition at line 11 of file TrRelationLookUp.cxx.

                                          {
        return m_times.size();
    }

nPar()

unsigned int MuonCalib::CalibFunc::nPar ( ) const

inlineinherited

Definition at line 39 of file CalibFunc.h.

{ return m_parameters.size(); }

par()

double MuonCalib::CalibFunc::par ( unsigned int  index ) const

inlineinherited

Definition at line 41 of file CalibFunc.h.

                                                 {
                return index < nPar() ? m_parameters[index] : 0.;
            }

parameters()

const ParVec& MuonCalib::CalibFunc::parameters ( ) const

inlineinherited

Definition at line 40 of file CalibFunc.h.

{ return m_parameters; }

typeName()

virtual std::string MuonCalib::ITrRelation::typeName ( ) const

inlinefinaloverridevirtualinherited

Implements MuonCalib::CalibFunc.

Definition at line 24 of file ITrRelation.h.

{ return "ITrRelation"; }

Member Data Documentation

m_maxRadius

double MuonCalib::TrRelationLookUp::m_maxRadius {}

private

Definition at line 28 of file TrRelationLookUp.h.

m_minRadius

double MuonCalib::TrRelationLookUp::m_minRadius {}

private

Definition at line 27 of file TrRelationLookUp.h.

m_parameters

ParVec MuonCalib::CalibFunc::m_parameters {}

privateinherited

Definition at line 48 of file CalibFunc.h.

m_radii

std::vector<double> MuonCalib::TrRelationLookUp::m_radii {}

private

Definition at line 26 of file TrRelationLookUp.h.

m_times

std::vector<double> MuonCalib::TrRelationLookUp::m_times {}

private

Definition at line 25 of file TrRelationLookUp.h.

---

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

• TrRelationLookUp.h
• TrRelationLookUp.cxx