MuonCalib::TrRelationLookUp Class Reference

#include <TrRelationLookUp.h>

Inheritance diagram for MuonCalib::TrRelationLookUp:

Collaboration diagram for MuonCalib::TrRelationLookUp:

Public Types
using	IRtRelationPtr = GeoModel::TransientConstSharedPtr< IRtRelation >
using	ParVec = std::vector< double >

Public Member Functions
	TrRelationLookUp (const IRtRelationPtr &rtRelation, const ParVec &vec={})
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 std::string	typeName () const override final
unsigned int	nPar () const
const ParVec &	parameters () const
double	par (unsigned int index) const

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

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

Detailed Description

Definition at line 13 of file TrRelationLookUp.h.

Member Typedef Documentation

IRtRelationPtr

using MuonCalib::ITrRelation::IRtRelationPtr = GeoModel::TransientConstSharedPtr<IRtRelation>

inherited

Definition at line 15 of file ITrRelation.h.

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 IRtRelationPtr &	rtRelation,
		const ParVec &	vec = {}
	)

Definition at line 12 of file TrRelationLookUp.cxx.

                                                                                          : ITrRelation{rtRelation, vec} {
 
        std::size_t numPoints(100);
        double stepSize{rtRelation->radius(rtRelation->tUpper()) / (numPoints - 1)};
        m_times = std::vector<double>(numPoints);
        m_radii = std::vector<double>(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 29 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 57 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 44 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);
    }

getTFromR()

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

private

Definition at line 65 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 63 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 62 of file TrRelationLookUp.cxx.

{ return m_minRadius; }

name()

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

finaloverridevirtual

Implements MuonCalib::CalibFunc.

Definition at line 27 of file TrRelationLookUp.cxx.

{ return "TrRelationLookUp"; }

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.

                                                 {
                if (index < nPar())
                    return m_parameters[index];
                else
                    return 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 21 of file ITrRelation.h.

{ return "ITrRelation"; }

Member Data Documentation

m_maxRadius

double MuonCalib::TrRelationLookUp::m_maxRadius {}

private

Definition at line 26 of file TrRelationLookUp.h.

m_minRadius

double MuonCalib::TrRelationLookUp::m_minRadius {}

private

Definition at line 25 of file TrRelationLookUp.h.

m_parameters

ParVec MuonCalib::CalibFunc::m_parameters {}

privateinherited

Definition at line 51 of file CalibFunc.h.

m_radii

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

private

Definition at line 24 of file TrRelationLookUp.h.

m_rt

IRtRelationPtr MuonCalib::ITrRelation::m_rt {}

privateinherited

Definition at line 36 of file ITrRelation.h.

m_times

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

private

Definition at line 23 of file TrRelationLookUp.h.

---

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

• TrRelationLookUp.h
• TrRelationLookUp.cxx