CMS 3D CMS Logo

List of all members | Classes | 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

Classes

struct  Config
 

Public Types

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

Public Member Functions

virtual ReferenceTrajectoryBaseclone () const =0
 
const AlgebraicMatrixderivatives () const
 
const Eigen::MatrixXd & gblExtDerivatives () const
 
const Eigen::VectorXd & gblExtMeasurements () const
 
const Eigen::VectorXd & gblExtPrecisions () const
 
std::vector< std::pair< std::vector< gbl::GblPoint >, Eigen::MatrixXd > > & 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::ConstRecHitContainerrecHits () const
 
void setParameterErrors (const AlgebraicSymMatrix &error)
 
const AlgebraicSymMatrixtrajectoryPositionErrors () const
 
const AlgebraicVectortrajectoryPositions () const
 
const std::vector< TrajectoryStateOnSurface > & trajectoryStates () const
 
const AlgebraicMatrixtrajectoryToCurv () const
 
 ~ReferenceTrajectoryBase () override
 

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
 
Eigen::MatrixXd theGblExtDerivatives
 
Eigen::VectorXd theGblExtMeasurements
 
Eigen::VectorXd theGblExtPrecisions
 
std::vector< std::pair< std::vector< gbl::GblPoint >, Eigen::MatrixXd > > 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< TrajectoryStateOnSurfacetheTsosVec
 
bool theValidityFlag
 

Static Protected Attributes

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

Member Typedef Documentation

Definition at line 109 of file ReferenceTrajectoryBase.h.

Member Enumeration Documentation

Constructor & Destructor Documentation

ReferenceTrajectoryBase::~ReferenceTrajectoryBase ( )
inlineoverride

Definition at line 136 of file ReferenceTrajectoryBase.h.

136 {}
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.

Referenced by nominalField().

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

References theDerivatives.

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

160 { return theDerivatives; }
const Eigen::MatrixXd& ReferenceTrajectoryBase::gblExtDerivatives ( ) const
inline

Returns the GBL external derivatives.

Definition at line 175 of file ReferenceTrajectoryBase.h.

References theGblExtDerivatives.

175 { return theGblExtDerivatives; }
const Eigen::VectorXd& ReferenceTrajectoryBase::gblExtMeasurements ( ) const
inline

Returns the GBL external derivatives.

Definition at line 179 of file ReferenceTrajectoryBase.h.

References theGblExtMeasurements.

179 { return theGblExtMeasurements; }
const Eigen::VectorXd& ReferenceTrajectoryBase::gblExtPrecisions ( ) const
inline

Returns the GBL external derivatives.

Definition at line 183 of file ReferenceTrajectoryBase.h.

References theGblExtPrecisions.

183 { return theGblExtPrecisions; }
std::vector<std::pair<std::vector<gbl::GblPoint>, Eigen::MatrixXd> >& ReferenceTrajectoryBase::gblInput ( )
inline

Returns the GBL input

Definition at line 171 of file ReferenceTrajectoryBase.h.

References theGblInput.

Referenced by TwoBodyDecayTrajectory::construct().

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

Returns the transformation of local to tracjectory parameters

Definition at line 167 of file ReferenceTrajectoryBase.h.

References theInnerLocalToTrajectory.

Referenced by TwoBodyDecayTrajectory::construct().

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

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

Definition at line 146 of file ReferenceTrajectoryBase.h.

References theMeasurementsCov.

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

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

Returns the measurements in local coordinates.

Definition at line 142 of file ReferenceTrajectoryBase.h.

References theMeasurements.

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

142 { return theMeasurements; }
int ReferenceTrajectoryBase::nominalField ( ) const
inline
unsigned int ReferenceTrajectoryBase::numberOfHitMeas ( ) const
inline
unsigned int ReferenceTrajectoryBase::numberOfHits ( ) const
inline

Definition at line 211 of file ReferenceTrajectoryBase.h.

References theNumberOfHits.

Referenced by TwoBodyDecayTrajectory::construct().

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

Definition at line 212 of file ReferenceTrajectoryBase.h.

References theNumberOfPars.

Referenced by TwoBodyDecayTrajectory::construct().

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

Definition at line 25 of file ReferenceTrajectoryBase.cc.

References useRecHit().

Referenced by nominalField().

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

References theNumberOfVirtualPars.

Referenced by TwoBodyDecayTrajectory::construct().

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

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

Definition at line 200 of file ReferenceTrajectoryBase.h.

References theParameterCov.

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

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

Definition at line 192 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 209 of file ReferenceTrajectoryBase.h.

References theRecHits.

Referenced by DualBzeroReferenceTrajectory::clone(), DualReferenceTrajectory::construct(), and nominalField().

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

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

Definition at line 196 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 155 of file ReferenceTrajectoryBase.h.

References theTrajectoryPositionCov.

Referenced by DualReferenceTrajectory::construct().

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

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

Definition at line 151 of file ReferenceTrajectoryBase.h.

References theTrajectoryPositions.

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

151 { 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 204 of file ReferenceTrajectoryBase.h.

References theTsosVec.

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

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

Returns the transformation of tracjectory to curvilinear parameters

Definition at line 164 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(), ReferenceTrajectory::getHitProjectionMatrix(), nominalField(), and numberOfUsedRecHits().

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

Member Data Documentation

unsigned int ReferenceTrajectoryBase::nMeasPerHit {2}
staticprotected
AlgebraicMatrix ReferenceTrajectoryBase::theDerivatives
protected
Eigen::MatrixXd ReferenceTrajectoryBase::theGblExtDerivatives
protected
Eigen::VectorXd ReferenceTrajectoryBase::theGblExtMeasurements
protected
Eigen::VectorXd ReferenceTrajectoryBase::theGblExtPrecisions
protected
std::vector<std::pair<std::vector<gbl::GblPoint>, Eigen::MatrixXd> > 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 229 of file ReferenceTrajectoryBase.h.

Referenced by parameterErrorsAvailable(), and setParameterErrors().

AlgebraicSymMatrix ReferenceTrajectoryBase::theParameterCov
protected

Definition at line 246 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