CMS 3D CMS Logo

Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | Static Protected Attributes

ReferenceTrajectoryBase Class Reference

#include <ReferenceTrajectoryBase.h>

Inheritance diagram for ReferenceTrajectoryBase:
ReferenceCounted DualReferenceTrajectory ReferenceTrajectory TwoBodyDecayTrajectory DualBzeroReferenceTrajectory BzeroReferenceTrajectory

List of all members.

Public Types

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

Public Member Functions

virtual ReferenceTrajectoryBaseclone () const =0
const AlgebraicMatrixderivatives () const
bool isValid ()
const AlgebraicMatrixlocalToTrajectory () const
const AlgebraicSymMatrixmeasurementErrors () const
const AlgebraicVectormeasurements () 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
AlgebraicMatrix theInnerLocalToTrajectory
AlgebraicMatrix theInnerTrajectoryToCurvilinear
AlgebraicVector theMeasurements
AlgebraicSymMatrix theMeasurementsCov
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:
2011/09/06 13:50:18

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

Definition at line 94 of file ReferenceTrajectoryBase.h.


Member Typedef Documentation

Definition at line 99 of file ReferenceTrajectoryBase.h.


Member Enumeration Documentation

Enumerator:
none 
multipleScattering 
energyLoss 
combined 
breakPoints 
brokenLinesCoarse 
brokenLinesFine 

Definition at line 101 of file ReferenceTrajectoryBase.h.


Constructor & Destructor Documentation

virtual ReferenceTrajectoryBase::~ReferenceTrajectoryBase ( ) [inline, virtual]

Definition at line 104 of file ReferenceTrajectoryBase.h.

{}
ReferenceTrajectoryBase::ReferenceTrajectoryBase ( unsigned int  nPar,
unsigned int  nHits,
unsigned int  nVirtualPar,
unsigned int  nVirtualMeas 
) [explicit, protected]

Definition at line 3 of file ReferenceTrajectoryBase.cc.

References theRecHits, and theTsosVec.

  : theValidityFlag(false), theParamCovFlag(false),
    theNumberOfHits( nHits ), theNumberOfPars( nPar ), 
    theNumberOfVirtualMeas( nVirtualMeas ), theNumberOfVirtualPars( nVirtualPar ),
    theTsosVec(), theRecHits(),
    theMeasurements(nMeasPerHit * nHits + nVirtualMeas), 
    theMeasurementsCov(nMeasPerHit * nHits + nVirtualMeas, 0),
    theTrajectoryPositions(nMeasPerHit * nHits), 
    theTrajectoryPositionCov(nMeasPerHit * nHits, 0),
    theParameters(nPar), 
    theParameterCov(nPar, 0),
    theDerivatives(nMeasPerHit * nHits + nVirtualMeas, nPar + nVirtualPar, 0),
    theInnerTrajectoryToCurvilinear( 5, 5, 0 ),
    theInnerLocalToTrajectory( 5, 5, 0 )    
{
  theTsosVec.reserve(nHits);
  theRecHits.reserve(nHits);
}

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 128 of file ReferenceTrajectoryBase.h.

References theDerivatives.

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

{ return theDerivatives; }
bool ReferenceTrajectoryBase::isValid ( void  ) [inline]
const AlgebraicMatrix& ReferenceTrajectoryBase::localToTrajectory ( ) const [inline]

Returns the transformation of local to tracjectory parameters

Definition at line 135 of file ReferenceTrajectoryBase.h.

References theInnerLocalToTrajectory.

Referenced by TwoBodyDecayTrajectory::construct().

const AlgebraicSymMatrix& ReferenceTrajectoryBase::measurementErrors ( ) const [inline]

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

Definition at line 114 of file ReferenceTrajectoryBase.h.

References theMeasurementsCov.

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

{ return theMeasurementsCov; }
const AlgebraicVector& ReferenceTrajectoryBase::measurements ( ) const [inline]

Returns the measurements in local coordinates.

Definition at line 110 of file ReferenceTrajectoryBase.h.

References theMeasurements.

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

{ return theMeasurements; }
unsigned int ReferenceTrajectoryBase::numberOfHitMeas ( ) const [inline]
unsigned int ReferenceTrajectoryBase::numberOfHits ( void  ) const [inline]

Definition at line 162 of file ReferenceTrajectoryBase.h.

References theNumberOfHits.

Referenced by TwoBodyDecayTrajectory::construct().

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

Definition at line 163 of file ReferenceTrajectoryBase.h.

References theNumberOfPars.

Referenced by TwoBodyDecayTrajectory::construct().

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

Definition at line 25 of file ReferenceTrajectoryBase.cc.

References useRecHit().

{
  unsigned int nUsedHits = 0;
  TransientTrackingRecHit::ConstRecHitContainer::const_iterator itHit;
  for ( itHit = recHits.begin(); itHit != recHits.end(); ++itHit ) if ( useRecHit( *itHit ) ) ++nUsedHits;
  return nUsedHits;
}
unsigned int ReferenceTrajectoryBase::numberOfVirtualMeas ( ) const [inline]
unsigned int ReferenceTrajectoryBase::numberOfVirtualPar ( ) const [inline]

Definition at line 165 of file ReferenceTrajectoryBase.h.

References theNumberOfVirtualPars.

Referenced by TwoBodyDecayTrajectory::construct().

const AlgebraicSymMatrix& ReferenceTrajectoryBase::parameterErrors ( ) const [inline]

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

Definition at line 151 of file ReferenceTrajectoryBase.h.

References theParameterCov.

{ return theParameterCov; }
bool ReferenceTrajectoryBase::parameterErrorsAvailable ( ) const [inline]

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

Definition at line 143 of file ReferenceTrajectoryBase.h.

References theParamCovFlag.

{ return theParamCovFlag; }
const AlgebraicVector& ReferenceTrajectoryBase::parameters ( void  ) const [inline]

Returns the set of 'track'-parameters.

Definition at line 139 of file ReferenceTrajectoryBase.h.

References theParameters.

Referenced by ReferenceTrajectory::fillDerivatives().

{ return theParameters; }
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 160 of file ReferenceTrajectoryBase.h.

References theRecHits.

Referenced by DualReferenceTrajectory::construct().

{ return theRecHits; }
void ReferenceTrajectoryBase::setParameterErrors ( const AlgebraicSymMatrix error) [inline]

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

Definition at line 147 of file ReferenceTrajectoryBase.h.

References 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 123 of file ReferenceTrajectoryBase.h.

References theTrajectoryPositionCov.

Referenced by DualReferenceTrajectory::construct().

const AlgebraicVector& ReferenceTrajectoryBase::trajectoryPositions ( ) const [inline]

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

Definition at line 119 of file ReferenceTrajectoryBase.h.

References theTrajectoryPositions.

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

const std::vector<TrajectoryStateOnSurface>& ReferenceTrajectoryBase::trajectoryStates ( ) const [inline]

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

Definition at line 155 of file ReferenceTrajectoryBase.h.

References theTsosVec.

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

{ return theTsosVec; }
const AlgebraicMatrix& ReferenceTrajectoryBase::trajectoryToCurv ( ) const [inline]

Returns the transformation of tracjectory to curvilinear parameters

Definition at line 132 of file ReferenceTrajectoryBase.h.

References theInnerTrajectoryToCurvilinear.

bool ReferenceTrajectoryBase::useRecHit ( const TransientTrackingRecHit::ConstRecHitPointer hitPtr) const [protected]

Definition at line 35 of file ReferenceTrajectoryBase.cc.

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

{
  return hitPtr->isValid();
}

Member Data Documentation

const unsigned int ReferenceTrajectoryBase::nMeasPerHit = 2 [static, protected]
unsigned int ReferenceTrajectoryBase::theNumberOfHits [protected]
unsigned int ReferenceTrajectoryBase::theNumberOfPars [protected]

Definition at line 179 of file ReferenceTrajectoryBase.h.

Referenced by parameterErrorsAvailable(), and setParameterErrors().

Definition at line 196 of file ReferenceTrajectoryBase.h.

Referenced by parameterErrors(), and setParameterErrors().