MuonCalib::RtSimplePolynomial Class Reference

#include <RtSimplePolynomial.h>

Inheritance diagram for MuonCalib::RtSimplePolynomial:

Collaboration diagram for MuonCalib::RtSimplePolynomial:

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

Public Member Functions
	RtSimplePolynomial (const ParVec &vec)
virtual std::string	name () const override final
	get the class name More...
virtual double	radius (double t) const override final
	get the drift velocity More...
virtual double	driftVelocity (double t) const override final
	get the drift acceleration More...
virtual double	driftAcceleration (double t) const override final
	Returns the acceleration of the r-t relation. More...
virtual double	tLower () const override final
	< get the lower drift-time bound More...
virtual double	tUpper () const override final
	Returns the upper time covered by the r-t. More...
virtual double	tBinWidth () const override final
	Returns the step-size for the sampling. More...
virtual unsigned	nDoF () const override final
	get the coefficients of the r(t) polynomial More...
std::vector< double >	rtParameters () const
virtual std::string	typeName () const override final
double	GetTmaxDiff () const
	return the difference in total dirft time between the two multilayers (ML1 - ML2) More...
bool	hasTmaxDiff () const
void	SetTmaxDiff (const double d)
	set the difference in total drift time betwene the two multilayers (ML1 - ML2) More...
double	getReducedTime (const double t) const
	map the in the interval [tLower;tUpper] onto the interval [-1. More...
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 Attributes
std::optional< double >	m_tmax_diff {std::nullopt}
ParVec	m_parameters {}

Detailed Description

Definition at line 15 of file RtSimplePolynomial.h.

Member Typedef Documentation

ParVec

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

inherited

Definition at line 35 of file CalibFunc.h.

Constructor & Destructor Documentation

RtSimplePolynomial()

RtSimplePolynomial::RtSimplePolynomial ( const ParVec &  vec )

explicit

Definition at line 9 of file RtSimplePolynomial.cxx.

                                                        : 
    IRtRelation(vec) { 
    // check for consistency //
    if (nPar() < 3) {
        THROW_EXCEPTION("RtSimplePolynomial::_init() - Not enough parameters!");
    }
    if (tLower() >= tUpper()) {
        THROW_EXCEPTION("Lower time boundary ("<<tLower()<<")>= upper time ("<<tUpper()<<") boundary!");
    }
}  // end RtSimplePolynomial::_init

Member Function Documentation

CalibFunc()

MuonCalib::CalibFunc::CalibFunc

inlineexplicitinherited

Definition at line 36 of file CalibFunc.h.

                                                  : 
                    m_parameters{vec} {}

dReducedTimeDt()

double MuonCalib::IRtRelation::dReducedTimeDt ( ) const

inlineinherited

Definition at line 53 of file IRtRelation.h.

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

driftAcceleration()

double RtSimplePolynomial::driftAcceleration ( double  t ) const

finaloverridevirtual

Returns the acceleration of the r-t relation.

Implements MuonCalib::IRtRelation.

Definition at line 56 of file RtSimplePolynomial.cxx.

                                                           {
    double acc{0.};
    const double x = getReducedTime(t);
     // Chain rule
    const double dx_dt = dReducedTimeDt();
    for (unsigned int k = 2; k < nDoF(); ++k) {
        acc += par(k+2) *k * (k-1) * std::pow(x,k-2)*std::pow(dx_dt,2);
    }
    return acc * t;
}

driftVelocity()

double RtSimplePolynomial::driftVelocity ( double  t ) const

finaloverridevirtual

get the drift acceleration

Implements MuonCalib::IRtRelation.

Definition at line 42 of file RtSimplePolynomial.cxx.

                                                       { 
    // Set derivative to 0 outside of the bounds
    if (t < tLower() || t > tUpper()) return 0.0;
 
    // Chain rule
    double drdt{0.};
    const double x = getReducedTime(t);
     // Chain rule
    const double dx_dt = dReducedTimeDt();
    for (unsigned int k = 1; k < nDoF(); ++k) {
        drdt += par(k+2) * k * std::pow(x,k-1)*dx_dt;
    }
    return drdt; 
}

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.

{ return m_tmax_diff.value_or(0.); }

hasTmaxDiff()

bool MuonCalib::IRtRelation::hasTmaxDiff ( ) const

inlineinherited

Definition at line 42 of file IRtRelation.h.

{ return m_tmax_diff.has_value(); }

name()

std::string RtSimplePolynomial::name ( ) const

finaloverridevirtual

get the class name

get the radius corresponding to the drift time t; if t is not within [t_low, t_up] an unphysical radius of 99999 is returned

Implements MuonCalib::CalibFunc.

Definition at line 20 of file RtSimplePolynomial.cxx.

{ return "RtSimplePolynomial"; }

nDoF()

unsigned RtSimplePolynomial::nDoF ( ) const

finaloverridevirtual

get the coefficients of the r(t) polynomial

Implements MuonCalib::IRtRelation.

Definition at line 68 of file RtSimplePolynomial.cxx.

{ return nPar() -2; }

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

radius()

double RtSimplePolynomial::radius ( double  t ) const

finaloverridevirtual

get the drift velocity

Implements MuonCalib::IRtRelation.

Definition at line 23 of file RtSimplePolynomial.cxx.

                                                {
    // INITIAL TIME CHECK //
    if (t < tLower()) return 0.0;
    if (t > tUpper()) return 14.6;
 
    double rad{0.0};  // auxiliary radius
    const double x = getReducedTime(t);
    // CALCULATE r(t) //
    for (unsigned int k = 0; k < nDoF(); k++) { 
        rad += par(k+2) * std::pow(x,k); 
    }
    return std::max(rad, 0.);
}

rtParameters()

std::vector< double > RtSimplePolynomial::rtParameters

(

)

const

Definition at line 70 of file RtSimplePolynomial.cxx.

                                                         {
    return std::vector<double>{parameters().begin() +2, parameters().end()};
}

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.

{ m_tmax_diff = d; }

tBinWidth()

double RtSimplePolynomial::tBinWidth ( ) const

finaloverridevirtual

Returns the step-size for the sampling.

Implements MuonCalib::IRtRelation.

Definition at line 21 of file RtSimplePolynomial.cxx.

{ return s_tBinWidth; }

tLower()

double RtSimplePolynomial::tLower ( ) const

finaloverridevirtual

< get the lower drift-time bound

get the upper drift-time bound

Implements MuonCalib::IRtRelation.

Definition at line 66 of file RtSimplePolynomial.cxx.

{ return par(0); }

tUpper()

double RtSimplePolynomial::tUpper ( ) const

finaloverridevirtual

Returns the upper time covered by the r-t.

Implements MuonCalib::IRtRelation.

Definition at line 67 of file RtSimplePolynomial.cxx.

{ return par(1); }

typeName()

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

inlinefinaloverridevirtualinherited

Implements MuonCalib::CalibFunc.

Definition at line 23 of file IRtRelation.h.

{ return "IRtRelation"; }

Member Data Documentation

m_parameters

ParVec MuonCalib::CalibFunc::m_parameters {}

privateinherited

Definition at line 48 of file CalibFunc.h.

m_tmax_diff

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

privateinherited

Definition at line 58 of file IRtRelation.h.

s_tBinWidth

constexpr 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:

• RtSimplePolynomial.h
• RtSimplePolynomial.cxx