Equidistant look up table for rt-relations with the time as key. More...

#include <RtRelationLookUp.h>

Inheritance diagram for MuonCalib::RtRelationLookUp:

Collaboration diagram for MuonCalib::RtRelationLookUp:

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

Public Member Functions
	RtRelationLookUp (const ParVec &vec)
virtual std::string	name () const override final
virtual double	radius (double t) const override final
	returns drift radius for a given time
virtual double	driftVelocity (double t) const override final
	returns drift velocity for a given time
virtual double	driftAcceleration (double t) const override final
	returns the acceleration for a given time
virtual double	tLower () const override final
	return rt range
virtual double	tUpper () const override final
	Returns the upper time covered by the r-t.
virtual double	tBinWidth () const override final
	Returns the step-size for the sampling.
virtual unsigned	nDoF () const override final
	Returns the number of degrees of freedom of the relation function.
virtual std::string	typeName () const override final
double	GetTmaxDiff () const
	return the difference in total dirft time between the two multilayers (ML1 - ML2)
bool	hasTmaxDiff () const
void	SetTmaxDiff (const double d)
	set the difference in total drift time betwene the two multilayers (ML1 - ML2)
double	getReducedTime (const double t) const
	map the in the interval [tLower;tUpper] onto the interval [-1.;1.
double	dReducedTimeDt () const
	CalibFunc (const ParVec &vec)
unsigned int	nPar () const
const ParVec &	parameters () const
double	par (unsigned int index) const

Static Protected Attributes
static constexpr double	s_tBinWidth = 1.e-3

Private Member Functions
int	getBin (double t) const
double	getRadius (int bin) const
int	binInRtRange (double t) const

Private Attributes
double	m_t_min {0.}
double	m_bin_size {0.}
std::optional< double >	m_tmax_diff {std::nullopt}
ParVec	m_parameters {}

Detailed Description

Equidistant look up table for rt-relations with the time as key.

The first parameter should be the time corresponding to the first bin. The second parameter should be the binsize.

The r value is calculated by linear interpolation.

Definition at line 23 of file RtRelationLookUp.h.

Member Typedef Documentation

◆ ParVec

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

inherited

Definition at line 35 of file CalibFunc.h.

Constructor & Destructor Documentation

◆ RtRelationLookUp()

MuonCalib::RtRelationLookUp::RtRelationLookUp ( const ParVec & vec )

explicit

Definition at line 10 of file RtRelationLookUp.cxx.

                                                        : 
        IRtRelation(vec) {
        if (vec.size() < 4) {
            m_t_min = 9e9;
            m_bin_size = 1.0;  // will be always out of range
            MsgStream log(Athena::getMessageSvc(), "RtRelationLookUp");
            log << MSG::WARNING << "<to few parameters>" << endmsg;
        } else {
            m_t_min = par(0);
            m_bin_size = par(1);
            if (m_bin_size == 0) {
                MsgStream log(Athena::getMessageSvc(), "RtRelationLookUp");
                log << MSG::WARNING << "<bin size=0>" << endmsg;
            }
        }
    }

Member Function Documentation

◆ binInRtRange()

int MuonCalib::RtRelationLookUp::binInRtRange ( double t ) const

inlineprivate

Definition at line 76 of file RtRelationLookUp.cxx.

                                                            {
        // get bin
        int bin = getBin(t);
 
        // if t corresponds to a negativ bin return first
        if (bin < 0) bin = 0;
 
        // if t corresponds to a bin outside the range of the lookup table return the last bin
        if (static_cast<unsigned>(bin) >= nDoF()) bin = nDoF() - 1;
 
        return bin;
    }

◆ CalibFunc()

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

inlineexplicitinherited

Definition at line 36 of file CalibFunc.h.

36 :

37 m_parameters{vec} {}

MuonCalib::CalibFunc::m_parameters

ParVec m_parameters

Definition CalibFunc.h:48

◆ dReducedTimeDt()

double MuonCalib::IRtRelation::dReducedTimeDt ( ) const

inlineinherited

Definition at line 53 of file IRtRelation.h.

                                      {
            return unitIntervalPrime(tLower(), tUpper());
        }

◆ driftAcceleration()

double MuonCalib::RtRelationLookUp::driftAcceleration ( double t ) const

finaloverridevirtual

returns the acceleration for a given time

Implements MuonCalib::IRtRelation.

Definition at line 72 of file RtRelationLookUp.cxx.

                                                             {
        return (driftVelocity(t+1) - driftVelocity(t-1)) / 2.;
    }

◆ driftVelocity()

double MuonCalib::RtRelationLookUp::driftVelocity ( double t ) const

finaloverridevirtual

returns drift velocity for a given time

Implements MuonCalib::IRtRelation.

Definition at line 57 of file RtRelationLookUp.cxx.

                                                         {
        // get best matching bin in rt range
        int bin = binInRtRange(t);
 
        // shift bin so we are using the last two bins for extrapolation
        if (static_cast<unsigned>(bin) >= nDoF() - 1) bin = nDoF() - 2;
 
        double r1 = getRadius(bin);      // get bin value
        double r2 = getRadius(bin + 1);  // get value of next bin
        double dr = r2 - r1;
 
        double v = dr / m_bin_size;
 
        return v;
    }

◆ getBin()

int MuonCalib::RtRelationLookUp::getBin ( double t ) const

inlineprivate

Definition at line 30 of file RtRelationLookUp.cxx.

                                                      {
        double t_minus_tmin{t - m_t_min};
        double rel = t_minus_tmin / m_bin_size;
        if (rel < static_cast<double>(INT_MIN)) return INT_MIN;
        if (rel > static_cast<double>(INT_MAX)) return INT_MAX;
        return static_cast<int>(rel);
    }

◆ getRadius()

double MuonCalib::RtRelationLookUp::getRadius ( int bin ) const

inlineprivate

Definition at line 47 of file RtRelationLookUp.h.

47{ return par(bin + 2); }

◆ getReducedTime()

double MuonCalib::IRtRelation::getReducedTime ( const double t ) const

inlineinherited

map the in the interval [tLower;tUpper] onto the interval [-1.;1.

] where tLower is mapped to -1. & tUpper to 1.;

Definition at line 49 of file IRtRelation.h.

                                                     {
           return mapToUnitInterval(t, tLower(), tUpper());
        }

◆ GetTmaxDiff()

double MuonCalib::IRtRelation::GetTmaxDiff ( ) const

inlineinherited

return the difference in total dirft time between the two multilayers (ML1 - ML2)

Definition at line 40 of file IRtRelation.h.

40{ return m_tmax_diff.value_or(0.); }

MuonCalib::IRtRelation::m_tmax_diff

std::optional< double > m_tmax_diff

Definition IRtRelation.h:58

◆ hasTmaxDiff()

bool MuonCalib::IRtRelation::hasTmaxDiff ( ) const

inlineinherited

Definition at line 42 of file IRtRelation.h.

42{ return m_tmax_diff.has_value(); }

◆ name()

virtual std::string MuonCalib::RtRelationLookUp::name ( ) const

inlinefinaloverridevirtual

Implements MuonCalib::CalibFunc.

Definition at line 27 of file RtRelationLookUp.h.

27{ return "RtRelationLookUp"; }

◆ nDoF()

unsigned MuonCalib::RtRelationLookUp::nDoF ( ) const

finaloverridevirtual

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

Implements MuonCalib::IRtRelation.

Definition at line 26 of file RtRelationLookUp.cxx.

                                         {
        return nPar() -2;
    }

◆ nPar()

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

inlineinherited

Definition at line 39 of file CalibFunc.h.

39{ 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.

40{ return m_parameters; }

◆ radius()

double MuonCalib::RtRelationLookUp::radius ( double t ) const

finaloverridevirtual

returns drift radius for a given time

Implements MuonCalib::IRtRelation.

Definition at line 38 of file RtRelationLookUp.cxx.

                                                  {
        // get best matching bin in rt range
        int bin = binInRtRange(t);
 
        // shift bin so we are using the last two bins for extrapolation
        if (static_cast<unsigned>(bin) >= nDoF() - 1) bin = nDoF() - 2;
 
        double r1 = getRadius(bin);      // get bin value
        double r2 = getRadius(bin + 1);  // get value of next bin
        double dr = r2 - r1;
 
        // scale factor for interpolation
        double scale = (t - m_t_min) / m_bin_size - (double)bin;
 
        double r = r1 + dr * scale;
 
        return r >= 0 ? r : 0;
    }

◆ SetTmaxDiff()

void MuonCalib::IRtRelation::SetTmaxDiff ( const double d )

inlineinherited

set the difference in total drift time betwene the two multilayers (ML1 - ML2)

Definition at line 45 of file IRtRelation.h.

45{ m_tmax_diff = d; }

hist_file_dump.d

d

Definition hist_file_dump.py:142

◆ tBinWidth()

double MuonCalib::RtRelationLookUp::tBinWidth ( ) const

finaloverridevirtual

Returns the step-size for the sampling.

Implements MuonCalib::IRtRelation.

Definition at line 91 of file RtRelationLookUp.cxx.

91{return m_bin_size; }

◆ tLower()

double MuonCalib::RtRelationLookUp::tLower ( ) const

finaloverridevirtual

return rt range

Implements MuonCalib::IRtRelation.

Definition at line 89 of file RtRelationLookUp.cxx.

89{ return m_t_min; }

◆ tUpper()

double MuonCalib::RtRelationLookUp::tUpper ( ) const

finaloverridevirtual

Returns the upper time covered by the r-t.

Implements MuonCalib::IRtRelation.

Definition at line 90 of file RtRelationLookUp.cxx.

90{ return m_t_min + m_bin_size * nDoF(); }

◆ typeName()

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

inlinefinaloverridevirtualinherited

Implements MuonCalib::CalibFunc.

Definition at line 23 of file IRtRelation.h.

23{ return "IRtRelation"; }

Member Data Documentation

◆ m_bin_size

double MuonCalib::RtRelationLookUp::m_bin_size {0.}

private

Definition at line 52 of file RtRelationLookUp.h.

52{0.};

◆ m_parameters

ParVec MuonCalib::CalibFunc::m_parameters {}

privateinherited

Definition at line 48 of file CalibFunc.h.

48{};

◆ m_t_min

double MuonCalib::RtRelationLookUp::m_t_min {0.}

private

Definition at line 51 of file RtRelationLookUp.h.

51{0.};

◆ m_tmax_diff

std::optional<double> MuonCalib::IRtRelation::m_tmax_diff {std::nullopt}

privateinherited

Definition at line 58 of file IRtRelation.h.

58{std::nullopt};

◆ s_tBinWidth

double MuonCalib::IRtRelation::s_tBinWidth = 1.e-3

staticconstexprprotectedinherited

Definition at line 60 of file IRtRelation.h.

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

Public Types

Public Member Functions

Static Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ ParVec

Constructor & Destructor Documentation

◆ RtRelationLookUp()

Member Function Documentation

◆ binInRtRange()

◆ CalibFunc()

◆ dReducedTimeDt()

◆ driftAcceleration()

◆ driftVelocity()

◆ getBin()

◆ getRadius()

◆ getReducedTime()

◆ GetTmaxDiff()

◆ hasTmaxDiff()

◆ name()

◆ nDoF()

◆ nPar()

◆ par()

◆ parameters()

◆ radius()

◆ SetTmaxDiff()

◆ tBinWidth()

◆ tLower()

◆ tUpper()

◆ typeName()

Member Data Documentation

◆ m_bin_size

◆ m_parameters

◆ m_t_min

◆ m_tmax_diff

◆ s_tBinWidth