CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | Static Protected Attributes
ReferenceTrajectoryBase Class Referenceabstract

#include <ReferenceTrajectoryBase.h>

Inheritance diagram for ReferenceTrajectoryBase:
ReferenceCounted DualReferenceTrajectory ReferenceTrajectory TwoBodyDecayTrajectory DualBzeroReferenceTrajectory BzeroReferenceTrajectory

Public Types

enum  MaterialEffects {
  none, multipleScattering, energyLoss, combined,
  breakPoints, brokenLinesCoarse, brokenLinesFine, localGBL,
  curvlinGBL
}
 
typedef
ReferenceCountingPointer
< ReferenceTrajectoryBase
ReferenceTrajectoryPtr
 

Public Member Functions

virtual ReferenceTrajectoryBaseclone () const =0
 
const AlgebraicMatrixderivatives () const
 
const TMatrixD & gblExtDerivatives () const
 
const TVectorD & gblExtMeasurements () const
 
const TVectorD & gblExtPrecisions () const
 
std::vector< std::pair
< std::vector< gbl::GblPoint >
, TMatrixD > > & 
gblInput ()
 
bool isValid ()
 
const AlgebraicMatrixlocalToTrajectory () const
 
const AlgebraicSymMatrixmeasurementErrors () const
 
const AlgebraicVectormeasurements () const
 
int nominalField () const
 
unsigned int numberOfHitMeas () const
 
unsigned int numberOfHits () const
 
unsigned int numberOfPar () const
 
unsigned int numberOfVirtualMeas () const
 
unsigned int numberOfVirtualPar () const
 
const AlgebraicSymMatrixparameterErrors () const
 
bool parameterErrorsAvailable () const
 
const AlgebraicVectorparameters () const
 
const
TransientTrackingRecHit::ConstRecHitContainer
recHits () const
 
void setParameterErrors (const AlgebraicSymMatrix &error)
 
const AlgebraicSymMatrixtrajectoryPositionErrors () const
 
const AlgebraicVectortrajectoryPositions () const
 
const std::vector
< TrajectoryStateOnSurface > & 
trajectoryStates () const
 
const AlgebraicMatrixtrajectoryToCurv () const
 
virtual ~ReferenceTrajectoryBase ()
 

Protected Member Functions

unsigned int numberOfUsedRecHits (const TransientTrackingRecHit::ConstRecHitContainer &recHits) const
 
 ReferenceTrajectoryBase (unsigned int nPar, unsigned int nHits, unsigned int nVirtualPar, unsigned int nVirtualMeas)
 
bool useRecHit (const TransientTrackingRecHit::ConstRecHitPointer &hitPtr) const
 

Protected Attributes

AlgebraicMatrix theDerivatives
 
TMatrixD theGblExtDerivatives
 
TVectorD theGblExtMeasurements
 
TVectorD theGblExtPrecisions
 
std::vector< std::pair
< std::vector< gbl::GblPoint >
, TMatrixD > > 
theGblInput
 
AlgebraicMatrix theInnerLocalToTrajectory
 
AlgebraicMatrix theInnerTrajectoryToCurvilinear
 
AlgebraicVector theMeasurements
 
AlgebraicSymMatrix theMeasurementsCov
 
int theNomField
 
unsigned int theNumberOfHits
 
unsigned int theNumberOfPars
 
unsigned int theNumberOfVirtualMeas
 
unsigned int theNumberOfVirtualPars
 
bool theParamCovFlag
 
AlgebraicSymMatrix theParameterCov
 
AlgebraicVector theParameters
 
TransientTrackingRecHit::ConstRecHitContainer theRecHits
 
AlgebraicSymMatrix theTrajectoryPositionCov
 
AlgebraicVector theTrajectoryPositions
 
std::vector
< TrajectoryStateOnSurface
theTsosVec
 
bool theValidityFlag
 

Static Protected Attributes

static const unsigned int nMeasPerHit = 2
 

Detailed Description

Author : Gero Flucke (based on code for ORCA by Edmund Widl) date : 2006/09/17 last update:

Date:
2010/09/10 12:11:39

by :

Author:
mussgill

Base class for reference 'trajectories' of single- or multiparticles stated. Inheriting classes have to calculate all relevant quantities accessed through member functions of this base class:

The local measured x/y coordinates on all crossed detectors as vector:

m = (x1, y1, x2, y2, ..., xN, yN) [transposed vector shown]

their covariance matrix (possibly containing correlations between hits due to material effects which are not taken into account by the ideal trajectory parametrisation, cf. below),

similarly the local x/y cordinates of the reference trajectory with covariance,

the parameters of the (ideal) 'trajectory' (e.g. 5 values for a single track or 9 for a two-track-state with vertex constraint),

the derivatives of the local coordinates of the reference trajectory with respect to the initial 'trajectory'-parameters, e.g. for n parameters 'p' and N hits with 'x/y' coordinates:

D = ( dx1/dp1, dx1/dp2, ..., dx1/dpn,

  dy1/dp1, dy1/dp2, ..., dy1/dpn,

  dx2/dp1, dx2/dp2, ..., dx2/dpn,

  dy2/dp1, dy2/dp2, ..., dy2/dpn,

     .        .             .

     .        .             .

  dxN/dp1, dxN/dp2, ..., dxN/dpn,

  dyN/dp1, dyN/dp2, ..., dyN/dpn )

and finally the TrajectoryStateOnSurface's of the reference trajectory.

Take care to check validity of the calculation (isValid()) before using the results.

.............................................................................

090730 C. Kleinwort: 'Break Points' introduced for better description of multiple scattering (correlations)

For each detector layer (in the hit collection) two ortogonal scattering angles are introduced as new local track parameters ("break points"). To constrain those to the expected mean (=0.) and variance (~1/p^2 X/X0) corresponding measurements are added. Invalid hits may be used to produce break points too (UseInvalidHits = True).

Break points have been implemented for: ReferenceTrajectory, BzeroReferenceTrajectory, DualReferenceTrajectory and DualBzeroReferenceTrajectory For 'Dual' trajectories they make no sense as only the correlations between the hits in a track half are accounted for and not those between the halves!

Break Points are selected by TrajectoryFactory.MaterialEffects = "BreakPoints"

090909 C. Kleinwort: 'Broken Lines' introduced for description of multiple scattering (V. Blobel, Nuclear Instruments and Methods A, 566 (2006), pp. 14-17) Fine Broken Lines are selected by TrajectoryFactory.MaterialEffects = "BrokenLinesFine" (exact derivatives) Coarse Broken Lines are selected by TrajectoryFactory.MaterialEffects = "BrokenLines[Coarse]" (approximate derivatives, closeby hits (ds<1cm) combined)

091127 C. Kleinwort: 1) 'Broken Lines' extended to PCA, required for beamspot constraint and TwoBodyDecayTrajectory, selected with "BrokenLines[Coarse]Pca" or "BrokenLinesFinePca" 2) For coarse Broken Lines linear interpolation is used for combined hits 3) TwoBodyDecayTrajectory implemented for break points and Broken Lines

141103 C. Kleinwort: 'General Broken Lines' introduced for description of multiple scattering (C. Kleinwort, Nuclear Instruments and Methods A, 673 (2012), pp. 107-110) needs GBL version >= V01-13-00 (from svnsrv.desy.de) Selected by TrajectoryFactory.MaterialEffects = "LocalGBL" or = "CurvlinGBL" (for trajectory constructed in local or curvilinear system)

Definition at line 103 of file ReferenceTrajectoryBase.h.

Member Typedef Documentation

Definition at line 108 of file ReferenceTrajectoryBase.h.

Member Enumeration Documentation

Constructor & Destructor Documentation

virtual ReferenceTrajectoryBase::~ReferenceTrajectoryBase ( )
inlinevirtual

Definition at line 113 of file ReferenceTrajectoryBase.h.

113 {}
ReferenceTrajectoryBase::ReferenceTrajectoryBase ( unsigned int  nPar,
unsigned int  nHits,
unsigned int  nVirtualPar,
unsigned int  nVirtualMeas 
)
explicitprotected

Definition at line 3 of file ReferenceTrajectoryBase.cc.

References theRecHits, and theTsosVec.

5  : theValidityFlag(false), theParamCovFlag(false),
6  theNumberOfHits( nHits ), theNumberOfPars( nPar ),
7  theNumberOfVirtualMeas( nVirtualMeas ), theNumberOfVirtualPars( nVirtualPar ),
9  theMeasurements(nMeasPerHit * nHits + nVirtualMeas),
10  theMeasurementsCov(nMeasPerHit * nHits + nVirtualMeas, 0),
13  theParameters(nPar),
14  theParameterCov(nPar, 0),
15  theDerivatives(nMeasPerHit * nHits + nVirtualMeas, nPar + nVirtualPar, 0),
17  theInnerLocalToTrajectory( 5, 5, 0 )
18 {
19  theTsosVec.reserve(nHits);
20  theRecHits.reserve(nHits);
21 }
AlgebraicMatrix theInnerTrajectoryToCurvilinear
AlgebraicMatrix theInnerLocalToTrajectory
static const unsigned int nMeasPerHit
AlgebraicSymMatrix theTrajectoryPositionCov
TransientTrackingRecHit::ConstRecHitContainer theRecHits
AlgebraicSymMatrix theMeasurementsCov
AlgebraicSymMatrix theParameterCov
std::vector< TrajectoryStateOnSurface > theTsosVec

Member Function Documentation

virtual ReferenceTrajectoryBase* ReferenceTrajectoryBase::clone ( ) const
pure virtual
const AlgebraicMatrix& ReferenceTrajectoryBase::derivatives ( void  ) const
inline

Returns the derivatives of the local coordinates of the reference trajectory (i.e. trajectoryPositions) w.r.t. the initial 'trajectory'-parameters.

Definition at line 137 of file ReferenceTrajectoryBase.h.

References theDerivatives.

Referenced by TwoBodyDecayTrajectory::construct(), and DualReferenceTrajectory::construct().

137 { return theDerivatives; }
const TMatrixD& ReferenceTrajectoryBase::gblExtDerivatives ( ) const
inline

Returns the GBL external derivatives.

Definition at line 152 of file ReferenceTrajectoryBase.h.

References theGblExtDerivatives.

const TVectorD& ReferenceTrajectoryBase::gblExtMeasurements ( ) const
inline

Returns the GBL external derivatives.

Definition at line 156 of file ReferenceTrajectoryBase.h.

References theGblExtMeasurements.

const TVectorD& ReferenceTrajectoryBase::gblExtPrecisions ( ) const
inline

Returns the GBL external derivatives.

Definition at line 160 of file ReferenceTrajectoryBase.h.

References theGblExtPrecisions.

std::vector<std::pair<std::vector<gbl::GblPoint>, TMatrixD> >& ReferenceTrajectoryBase::gblInput ( )
inline

Returns the GBL input

Definition at line 148 of file ReferenceTrajectoryBase.h.

References theGblInput.

Referenced by TwoBodyDecayTrajectory::construct().

148 { return theGblInput; }
std::vector< std::pair< std::vector< gbl::GblPoint >, TMatrixD > > theGblInput
bool ReferenceTrajectoryBase::isValid ( void  )
inline
const AlgebraicMatrix& ReferenceTrajectoryBase::localToTrajectory ( ) const
inline

Returns the transformation of local to tracjectory parameters

Definition at line 144 of file ReferenceTrajectoryBase.h.

References theInnerLocalToTrajectory.

Referenced by TwoBodyDecayTrajectory::construct().

144 { return theInnerLocalToTrajectory; }
AlgebraicMatrix theInnerLocalToTrajectory
const AlgebraicSymMatrix& ReferenceTrajectoryBase::measurementErrors ( ) const
inline

Returns the full covariance matrix of the measurements. ORCA-equivalent: covariance()

Definition at line 123 of file ReferenceTrajectoryBase.h.

References theMeasurementsCov.

Referenced by TwoBodyDecayTrajectory::construct(), and DualReferenceTrajectory::construct().

123 { return theMeasurementsCov; }
AlgebraicSymMatrix theMeasurementsCov
const AlgebraicVector& ReferenceTrajectoryBase::measurements ( ) const
inline

Returns the measurements in local coordinates.

Definition at line 119 of file ReferenceTrajectoryBase.h.

References theMeasurements.

Referenced by TwoBodyDecayTrajectory::construct(), and DualReferenceTrajectory::construct().

119 { return theMeasurements; }
int ReferenceTrajectoryBase::nominalField ( ) const
inline

Definition at line 193 of file ReferenceTrajectoryBase.h.

References theNomField.

Referenced by TwoBodyDecayTrajectory::construct().

unsigned int ReferenceTrajectoryBase::numberOfHitMeas ( ) const
inline
unsigned int ReferenceTrajectoryBase::numberOfHits ( ) const
inline

Definition at line 188 of file ReferenceTrajectoryBase.h.

References theNumberOfHits.

Referenced by TwoBodyDecayTrajectory::construct().

188 { return theNumberOfHits; }
unsigned int ReferenceTrajectoryBase::numberOfPar ( ) const
inline

Definition at line 189 of file ReferenceTrajectoryBase.h.

References theNumberOfPars.

Referenced by TwoBodyDecayTrajectory::construct().

189 { return theNumberOfPars; }
unsigned int ReferenceTrajectoryBase::numberOfUsedRecHits ( const TransientTrackingRecHit::ConstRecHitContainer recHits) const
protected

Definition at line 25 of file ReferenceTrajectoryBase.cc.

References useRecHit().

26 {
27  unsigned int nUsedHits = 0;
28  TransientTrackingRecHit::ConstRecHitContainer::const_iterator itHit;
29  for ( itHit = recHits.begin(); itHit != recHits.end(); ++itHit ) if ( useRecHit( *itHit ) ) ++nUsedHits;
30  return nUsedHits;
31 }
const TransientTrackingRecHit::ConstRecHitContainer & recHits() const
bool useRecHit(const TransientTrackingRecHit::ConstRecHitPointer &hitPtr) const
unsigned int ReferenceTrajectoryBase::numberOfVirtualMeas ( ) const
inline
unsigned int ReferenceTrajectoryBase::numberOfVirtualPar ( ) const
inline

Definition at line 191 of file ReferenceTrajectoryBase.h.

References theNumberOfVirtualPars.

Referenced by TwoBodyDecayTrajectory::construct().

191 { return theNumberOfVirtualPars; }
const AlgebraicSymMatrix& ReferenceTrajectoryBase::parameterErrors ( ) const
inline

Returns the covariance matrix of the 'track'-parameters.

Definition at line 177 of file ReferenceTrajectoryBase.h.

References theParameterCov.

177 { return theParameterCov; }
AlgebraicSymMatrix theParameterCov
bool ReferenceTrajectoryBase::parameterErrorsAvailable ( ) const
inline

Returns true if the covariance matrix of the 'track'-parameters is set.

Definition at line 169 of file ReferenceTrajectoryBase.h.

References theParamCovFlag.

const AlgebraicVector& ReferenceTrajectoryBase::parameters ( void  ) const
inline
const TransientTrackingRecHit::ConstRecHitContainer& ReferenceTrajectoryBase::recHits ( ) const
inline

Returns the TransientTrackingRecHits (as ConstRecHitPointer's), order might depend on concrete implementation of inheriting class

Definition at line 186 of file ReferenceTrajectoryBase.h.

References theRecHits.

Referenced by DualReferenceTrajectory::construct().

186 { return theRecHits; }
TransientTrackingRecHit::ConstRecHitContainer theRecHits
void ReferenceTrajectoryBase::setParameterErrors ( const AlgebraicSymMatrix error)
inline

Set the covariance matrix of the 'track'-parameters.

Definition at line 173 of file ReferenceTrajectoryBase.h.

References relativeConstraints::error, theParamCovFlag, and theParameterCov.

Referenced by TwoBodyDecayTrajectoryFactory::constructTrajectories().

const AlgebraicSymMatrix& ReferenceTrajectoryBase::trajectoryPositionErrors ( ) const
inline

Returns the covariance matrix of the reference trajectory.

Definition at line 132 of file ReferenceTrajectoryBase.h.

References theTrajectoryPositionCov.

Referenced by DualReferenceTrajectory::construct().

132 { return theTrajectoryPositionCov; }
AlgebraicSymMatrix theTrajectoryPositionCov
const AlgebraicVector& ReferenceTrajectoryBase::trajectoryPositions ( ) const
inline

Returns the local coordinates of the reference trajectory. ORCA-equivalent: referenceTrack()

Definition at line 128 of file ReferenceTrajectoryBase.h.

References theTrajectoryPositions.

Referenced by TwoBodyDecayTrajectory::construct(), and DualReferenceTrajectory::construct().

128 { return theTrajectoryPositions; }
const std::vector<TrajectoryStateOnSurface>& ReferenceTrajectoryBase::trajectoryStates ( ) const
inline

Returns the Tsos at the surfaces of the hits, parallel to recHits()

Definition at line 181 of file ReferenceTrajectoryBase.h.

References theTsosVec.

Referenced by TwoBodyDecayTrajectory::construct(), DualReferenceTrajectory::construct(), and TwoBodyDecayTrajectory::constructTsosVecWithErrors().

181 { return theTsosVec; }
std::vector< TrajectoryStateOnSurface > theTsosVec
const AlgebraicMatrix& ReferenceTrajectoryBase::trajectoryToCurv ( ) const
inline

Returns the transformation of tracjectory to curvilinear parameters

Definition at line 141 of file ReferenceTrajectoryBase.h.

References theInnerTrajectoryToCurvilinear.

AlgebraicMatrix theInnerTrajectoryToCurvilinear
bool ReferenceTrajectoryBase::useRecHit ( const TransientTrackingRecHit::ConstRecHitPointer hitPtr) const
protected

Definition at line 35 of file ReferenceTrajectoryBase.cc.

Referenced by ReferenceTrajectory::construct(), and numberOfUsedRecHits().

36 {
37  return hitPtr->isValid();
38 }

Member Data Documentation

const unsigned int ReferenceTrajectoryBase::nMeasPerHit = 2
staticprotected
AlgebraicMatrix ReferenceTrajectoryBase::theDerivatives
protected
TMatrixD ReferenceTrajectoryBase::theGblExtDerivatives
protected
TVectorD ReferenceTrajectoryBase::theGblExtMeasurements
protected
TVectorD ReferenceTrajectoryBase::theGblExtPrecisions
protected
std::vector<std::pair<std::vector<gbl::GblPoint>, TMatrixD> > ReferenceTrajectoryBase::theGblInput
protected
AlgebraicMatrix ReferenceTrajectoryBase::theInnerLocalToTrajectory
protected
AlgebraicMatrix ReferenceTrajectoryBase::theInnerTrajectoryToCurvilinear
protected
AlgebraicVector ReferenceTrajectoryBase::theMeasurements
protected
AlgebraicSymMatrix ReferenceTrajectoryBase::theMeasurementsCov
protected
int ReferenceTrajectoryBase::theNomField
protected
unsigned int ReferenceTrajectoryBase::theNumberOfHits
protected
unsigned int ReferenceTrajectoryBase::theNumberOfPars
protected
unsigned int ReferenceTrajectoryBase::theNumberOfVirtualMeas
protected
unsigned int ReferenceTrajectoryBase::theNumberOfVirtualPars
protected
bool ReferenceTrajectoryBase::theParamCovFlag
protected

Definition at line 206 of file ReferenceTrajectoryBase.h.

Referenced by parameterErrorsAvailable(), and setParameterErrors().

AlgebraicSymMatrix ReferenceTrajectoryBase::theParameterCov
protected

Definition at line 223 of file ReferenceTrajectoryBase.h.

Referenced by parameterErrors(), and setParameterErrors().

AlgebraicVector ReferenceTrajectoryBase::theParameters
protected
TransientTrackingRecHit::ConstRecHitContainer ReferenceTrajectoryBase::theRecHits
protected
AlgebraicSymMatrix ReferenceTrajectoryBase::theTrajectoryPositionCov
protected
AlgebraicVector ReferenceTrajectoryBase::theTrajectoryPositions
protected
std::vector<TrajectoryStateOnSurface> ReferenceTrajectoryBase::theTsosVec
protected
bool ReferenceTrajectoryBase::theValidityFlag
protected