#include <RZLine.h>

Classes
struct	ErrZ2_tag

Public Member Functions
float	chi2 () const

float	cotTheta () const

float	covii () const

float	covsi () const

float	covss () const

float	intercept () const

template<typename P , typename E , typename B >
	RZLine (const P &points, const E &errors, const B &isBarrel)

	RZLine (const std::vector< float > &r, const std::vector< float > &z, const std::vector< float > &errZ)

template<size_t N>
	RZLine (const std::array< float, N > &r, const std::array< float, N > &z, const std::array< float, N > &errZ)

template<typename T >
	RZLine (const T &r, const T &z, const T &errZ2, ErrZ2_tag)

Private Member Functions
template<typename R , typename Z , typename E >
void	calculate (const R &r, const Z &z, const E &errZ2)

template<typename T >
T	sqr (T t)

Private Attributes
float	chi2_

float	cotTheta_

float	covii_

float	covsi_

float	covss_

float	intercept_

Detailed Description

Definition at line 12 of file RZLine.h.

Constructor & Destructor Documentation

template<typename P , typename E , typename B >

RZLine::RZLine	(	const P &	points,
		const E &	errors,
		const B &	isBarrel
	)

inline

Constructor for containers of GlobalPoint, GlobalError, and bool

Template Parameters

P	Container of GlobalPoint
E	Container of GlobalError
B	Container of bool

Container can be e.g. std::vector, std::array, or DynArray.

Although for std::array use this constructor could be specialized to use std::array instead of DynArray for temporary storage.

Definition at line 29 of file RZLine.h.

References calculate(), declareDynArray, mps_fire::i, gen::n, AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), alignCSCRings::r, sqr(), z, and PV3DBase< T, PVType, FrameType >::z().

                                                               {
     const size_t n = points.size();
     declareDynArray(float, n, r);
     declareDynArray(float, n, z);
     declareDynArray(float, n, errZ2);
     for(size_t i=0; i<n; ++i) {
       const GlobalPoint& p = points[i];
       r[i] = p.perp();
       z[i] = p.z();
     }
 
     float simpleCot2 = sqr( (z[n-1]-z[0])/ (r[n-1]-r[0]) );
     for (size_t i=0; i<n; ++i) {
       errZ2[i] = (isBarrel[i]) ? errors[i].czz() :
         errors[i].rerr(points[i])  * simpleCot2;
     }
 
     calculate(r, z, errZ2);
   }

RZLine::RZLine	(	const std::vector< float > &	r,
		const std::vector< float > &	z,
		const std::vector< float > &	errZ
	)

inline

Constructor for std::vector of r, z, and z standard deviation

Definition at line 52 of file RZLine.h.

References calculate(), declareDynArray, mps_fire::i, gen::n, and sqr().

                                         {
     const size_t n = errZ.size();
     declareDynArray(float, n, errZ2);
     for(size_t i=0; i<n; ++i) errZ2[i] = sqr(errZ[i]);
     calculate(r, z, errZ2);
   }

template<size_t N>

RZLine::RZLine	(	const std::array< float, N > &	r,
		const std::array< float, N > &	z,
		const std::array< float, N > &	errZ
	)

inline

Constructor for std::array of r, z, and z standard deviation

Definition at line 65 of file RZLine.h.

References calculate(), mps_fire::i, N, and sqr().

                                          {
     std::array<float, N> errZ2;
     for(size_t i=0; i<N; ++i) errZ2[i] = sqr(errZ[i]);
     calculate(r, z, errZ2);
   }

template<typename T >

RZLine::RZLine	(	const T &	r,
		const T &	z,
		const T &	errZ2,
		ErrZ2_tag
	)

inline

Constructor for container of r, z, and z variance

Template Parameters

T	Container of float

Container can be e.g. std::vector, std::array, or DynArray.

The ErrZ2_tag parameter is used to distinguish this constructor from other 3-parameter constructors.

Passing variance is useful in cases where it is already available to avoid making a square of a square root.

Definition at line 87 of file RZLine.h.

References calculate().

                                                              {
     calculate(r, z, errZ2);
   }

Member Function Documentation

template<typename R , typename Z , typename E >

void RZLine::calculate	(	const R &	r,
		const Z &	z,
		const E &	errZ2
	)

inlineprivate

Definition at line 101 of file RZLine.h.

References chi2_, cotTheta_, covii_, covsi_, covss_, mps_fire::i, intercept_, linearFit(), gen::n, and sqr().

Referenced by RZLine().

                                                          {
     const size_t n = r.size();
     linearFit(r.data(), z.data(), n, errZ2.data(), cotTheta_, intercept_, covss_, covii_, covsi_);
     chi2_ = 0.f;
     for(size_t i=0; i<n; ++i) {
       chi2_ += sqr( ((z[i]-intercept_) - cotTheta_*r[i]) ) / errZ2[i];
     }
   }