MuonCalib::RtChebyshev Class Reference

#include <RtChebyshev.h>

Inheritance diagram for MuonCalib::RtChebyshev:

Collaboration diagram for MuonCalib::RtChebyshev:

Public Types
typedef std::vector< double >	ParVec

Public Member Functions
	RtChebyshev (const ParVec &vec)
	initialization constructor, More...
std::string	name (void) const
	get the class name More...
double	radius (double t) const
	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 More...
double	driftvelocity (double t) const
double	tLower (void) const
	< get the lower drift-time bound More...
double	tUpper (void) const
	get the number of parameters used to describe the r(t) relationship More...
unsigned int	numberOfRtParameters (void) const
	get the coefficients of the r(t) polynomial More...
std::vector< double >	rtParameters (void) const
	get the reduced time which is the argument of the Chebyshev polynomial More...
double	get_reduced_time (const double &t) const
virtual std::string	typeName () const
float	GetTmaxDiff () const
	return the difference in total dirft time between the two multilayers (ML1 - ML2) More...
bool	HasTmaxDiff () const
void	SetTmaxDiff (const float &d)
	set the difference in total drift time betwene the two multilayers (ML1 - ML2) More...
unsigned int	nPar () const
const ParVec &	parameters () const
double	par (unsigned int index) const

Public Attributes
float	m_tmax_diff

Private Member Functions
void	_init (void)
	initialization method More...

Private Attributes
const Tschebyscheff_polynomial *	m_Chebyshev
	pointer to the Chebyshev polynomials More...
ParVec	m_parameters

Detailed Description

This class contains the implementation of an r(t) relationship parametrized by a linear combination of Chebyshev polyonomials.

Convention:

\[ r(t) = \sum_{k=0}^{K} p_k*T_k(2*(t-0.5*(tupper+tlower))/(tupper-tlower) \]

where T_k is the Chebyshev polynomial of k-th order, tupper and tlower are upper and lower drift-time bounds.

Units: [t] = ns, [r] = mm, [v] = mm/ns.

Author

Olive.nosp@m.r.Ko.nosp@m.rtner.nosp@m.@cer.nosp@m.n.ch

Date

04.06.2006

Definition at line 52 of file RtChebyshev.h.

Member Typedef Documentation

ParVec

typedef std::vector<double> MuonCalib::CalibFunc::ParVec

inherited

Definition at line 36 of file CalibFunc.h.

Constructor & Destructor Documentation

RtChebyshev()

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

inlineexplicit

initialization constructor,

size of ParVec - 2 = order of the r(t) polynomial,

ParVec[0] = t_low (smallest allowed drift time), ParVec[1] = t_up (largest allowed drift time). ParVec[2...] = parameters of the Chebyshev polynomial

Definition at line 73 of file RtChebyshev.h.

: IRtRelation(vec) { _init(); }

Member Function Documentation

_init()

void RtChebyshev::_init ( void  )

private

initialization method

Definition at line 26 of file RtChebyshev.cxx.

                            {
    // check for consistency //
    if (nPar() < 3) {
        throw std::runtime_error(Form("File: %s, Line: %d\nRtChebyshev::_init() - Not enough parameters!", __FILE__, __LINE__));
    }
    if (parameters()[0] >= parameters()[1]) {
        throw std::runtime_error(
            Form("File: %s, Line: %d\nRtChebyshev::_init() - Lower time boundary >= upper time boundary!", __FILE__, __LINE__));
    }
 
    // pointer to the chebyshev service //
    m_Chebyshev = Tschebyscheff_polynomial::get_Tschebyscheff_polynomial();
 
    return;
}  // end RtChebyshev::_init

driftvelocity()

double RtChebyshev::driftvelocity ( double  t ) const

virtual

Implements MuonCalib::IRtRelation.

Definition at line 85 of file RtChebyshev.cxx.

{ return (radius(t + 1.0) - radius(t)); }

get_reduced_time()

double RtChebyshev::get_reduced_time

(

const double &

t

)

const

Definition at line 125 of file RtChebyshev.cxx.

                                                          {
    return 2 * (t - 0.5 * (parameters()[1] + parameters()[0])) / (parameters()[1] - parameters()[0]);
}

GetTmaxDiff()

float MuonCalib::IRtRelation::GetTmaxDiff ( ) const

inlineinherited

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

Definition at line 27 of file IRtRelation.h.

{ return HasTmaxDiff() ? m_tmax_diff : 0.0; }

HasTmaxDiff()

bool MuonCalib::IRtRelation::HasTmaxDiff ( ) const

inlineinherited

Definition at line 29 of file IRtRelation.h.

{ return (m_tmax_diff < 8e8); }

name()

std::string RtChebyshev::name ( void  ) const

virtual

get the class name

Implements MuonCalib::CalibFunc.

Definition at line 47 of file RtChebyshev.cxx.

{ return std::string("RtChebyshev"); }

nPar()

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

inlineinherited

Definition at line 39 of file CalibFunc.h.

{ return m_parameters.size(); }

numberOfRtParameters()

unsigned int RtChebyshev::numberOfRtParameters

(

void

)

const

get the coefficients of the r(t) polynomial

Definition at line 106 of file RtChebyshev.cxx.

{ return nPar() - 2; }

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

radius()

double RtChebyshev::radius ( double  t ) const

virtual

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

get the drift velocity

Implements MuonCalib::IRtRelation.

Definition at line 54 of file RtChebyshev.cxx.

                                         {
    // INITIAL TIME CHECK //
    if (t < parameters()[0]) return 0.0;
    if (t > parameters()[1]) return 14.6;
    // if x is out of bounds, return 99999 //
    //  if (t<parameters()[0] || t>parameters()[1]) {
    //      return 99999.0;
    //  }
 
    // VARIABLES //
    // argument of the Chebyshev polynomials
    double x(2 * (t - 0.5 * (parameters()[1] + parameters()[0])) / (parameters()[1] - parameters()[0]));
    double rad(0.0);  // auxiliary radius
 
    // CALCULATE r(t) //
    for (unsigned int k = 0; k < nPar() - 2; k++) { rad = rad + parameters()[k + 2] * m_Chebyshev->value(k, x); }
 
    return rad >= 0 ? rad : 0;
}

rtParameters()

std::vector< double > RtChebyshev::rtParameters

(

void

)

const

get the reduced time which is the argument of the Chebyshev polynomial

Definition at line 113 of file RtChebyshev.cxx.

                                                      {
    std::vector<double> alpha(nPar() - 2);
    for (unsigned int k = 0; k < alpha.size(); k++) { alpha[k] = parameters()[k + 2]; }
 
    return alpha;
}

SetTmaxDiff()

void MuonCalib::IRtRelation::SetTmaxDiff ( const float &  d )

inlineinherited

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

Definition at line 32 of file IRtRelation.h.

{ m_tmax_diff = d; }

tLower()

double RtChebyshev::tLower ( void  ) const

virtual

< get the lower drift-time bound

get the upper drift-time bound

Implements MuonCalib::IRtRelation.

Definition at line 92 of file RtChebyshev.cxx.

{ return parameters()[0]; }

tUpper()

double RtChebyshev::tUpper ( void  ) const

virtual

get the number of parameters used to describe the r(t) relationship

Implements MuonCalib::IRtRelation.

Definition at line 99 of file RtChebyshev.cxx.

{ return parameters()[1]; }

typeName()

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

inlinevirtualinherited

Implements MuonCalib::CalibFunc.

Definition at line 18 of file IRtRelation.h.

{ return "IRtRelation"; }

Member Data Documentation

m_Chebyshev

const Tschebyscheff_polynomial* MuonCalib::RtChebyshev::m_Chebyshev

private

pointer to the Chebyshev polynomials

Definition at line 55 of file RtChebyshev.h.

m_parameters

ParVec MuonCalib::CalibFunc::m_parameters

privateinherited

Definition at line 51 of file CalibFunc.h.

m_tmax_diff

float MuonCalib::IRtRelation::m_tmax_diff

inherited

Definition at line 34 of file IRtRelation.h.

---

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

• RtChebyshev.h
• RtChebyshev.cxx