#include <BzeroReferenceTrajectory.h>

Inheritance diagram for BzeroReferenceTrajectory:

Public Member Functions
	BzeroReferenceTrajectory (const TrajectoryStateOnSurface &tsos, const TransientTrackingRecHit::ConstRecHitContainer &recHits, const MagneticField *magField, const reco::BeamSpot &beamSpot, const ReferenceTrajectoryBase::Config &config)

BzeroReferenceTrajectory *	clone () const override

	~BzeroReferenceTrajectory () override

Public Member Functions inherited from ReferenceTrajectory
ReferenceTrajectory *	clone () const override

	ReferenceTrajectory (const TrajectoryStateOnSurface &referenceTsos, const TransientTrackingRecHit::ConstRecHitContainer &recHits, const MagneticField *magField, const reco::BeamSpot &beamSpot, const ReferenceTrajectoryBase::Config &config)

	~ReferenceTrajectory () override

Public Member Functions inherited from ReferenceTrajectoryBase
const AlgebraicMatrix &	derivatives () 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 AlgebraicMatrix &	localToTrajectory () const

const AlgebraicSymMatrix &	measurementErrors () const

const AlgebraicVector &	measurements () 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 AlgebraicSymMatrix &	parameterErrors () const

bool	parameterErrorsAvailable () const

const AlgebraicVector &	parameters () const

const TransientTrackingRecHit::ConstRecHitContainer &	recHits () const

void	setParameterErrors (const AlgebraicSymMatrix &error)

const AlgebraicSymMatrix &	trajectoryPositionErrors () const

const AlgebraicVector &	trajectoryPositions () const

const std::vector< TrajectoryStateOnSurface > &	trajectoryStates () const

const AlgebraicMatrix &	trajectoryToCurv () const

	~ReferenceTrajectoryBase () override

Private Attributes
double	theMomentumEstimate

Additional Inherited Members
Public Types inherited from ReferenceTrajectory
typedef SurfaceSideDefinition::SurfaceSide	SurfaceSide

Public Types inherited from ReferenceTrajectoryBase
enum	MaterialEffects { none, multipleScattering, energyLoss, combined, breakPoints, brokenLinesCoarse, brokenLinesFine, localGBL, curvlinGBL }

typedef ReferenceCountingPointer< ReferenceTrajectoryBase >	ReferenceTrajectoryPtr

Protected Member Functions inherited from ReferenceTrajectory
virtual bool	addMaterialEffectsBp (const std::vector< AlgebraicMatrix > &allJacobians, const std::vector< AlgebraicMatrix > &allProjections, const std::vector< AlgebraicSymMatrix > &allCurvChanges, const std::vector< AlgebraicSymMatrix > &allDeltaParaCovs, const std::vector< AlgebraicMatrix > &allLocalToCurv)

virtual bool	addMaterialEffectsBrl (const std::vector< AlgebraicMatrix > &allJacobians, const std::vector< AlgebraicMatrix > &allProjections, const std::vector< AlgebraicSymMatrix > &allCurvChanges, const std::vector< AlgebraicSymMatrix > &allDeltaParaCovs, const std::vector< AlgebraicMatrix > &allLocalToCurv, const GlobalTrajectoryParameters &gtp)

virtual bool	addMaterialEffectsBrl (const std::vector< AlgebraicMatrix > &allProjections, const std::vector< AlgebraicSymMatrix > &allDeltaParaCovs, const std::vector< AlgebraicMatrix > &allLocalToCurv, const std::vector< double > &allSteps, const GlobalTrajectoryParameters &gtp, const double minStep=1.0)

virtual bool	addMaterialEffectsCov (const std::vector< AlgebraicMatrix > &allJacobians, const std::vector< AlgebraicMatrix > &allProjections, const std::vector< AlgebraicSymMatrix > &allCurvChanges, const std::vector< AlgebraicSymMatrix > &allDeltaParaCovs)

virtual bool	addMaterialEffectsCurvlinGbl (const std::vector< AlgebraicMatrix > &allJacobians, const std::vector< AlgebraicMatrix > &allProjections, const std::vector< AlgebraicSymMatrix > &allCurvChanges, const std::vector< AlgebraicSymMatrix > &allDeltaParaCovs, const std::vector< AlgebraicMatrix > &allLocalToCurv)

virtual bool	addMaterialEffectsLocalGbl (const std::vector< AlgebraicMatrix > &allJacobians, const std::vector< AlgebraicMatrix > &allProjections, const std::vector< AlgebraicSymMatrix > &allCurvatureChanges, const std::vector< AlgebraicSymMatrix > &allDeltaParameterCovs)

virtual bool	construct (const TrajectoryStateOnSurface &referenceTsos, const TransientTrackingRecHit::ConstRecHitContainer &recHits, const MagneticField *magField, const reco::BeamSpot &beamSpot)

MaterialEffectsUpdator *	createUpdator (MaterialEffects materialEffects, double mass) const

virtual void	fillDerivatives (const AlgebraicMatrix &projection, const AlgebraicMatrix &fullJacobian, unsigned int iRow)

virtual void	fillMeasurementAndError (const TransientTrackingRecHit::ConstRecHitPointer &hitPtr, unsigned int iRow, const TrajectoryStateOnSurface &updatedTsos)

virtual void	fillTrajectoryPositions (const AlgebraicMatrix &projection, const AlgebraicVector &mixedLocalParams, unsigned int iRow)

AlgebraicMatrix	getHitProjectionMatrix (const TransientTrackingRecHit::ConstRecHitPointer &recHit) const

template<unsigned int N>
AlgebraicMatrix	getHitProjectionMatrixT (const TransientTrackingRecHit::ConstRecHitPointer &recHit) const

virtual bool	propagate (const Plane &previousSurface, const TrajectoryStateOnSurface &previousTsos, const Plane &newSurface, TrajectoryStateOnSurface &newTsos, AlgebraicMatrix &newJacobian, AlgebraicMatrix &newCurvlinJacobian, double &nextStep, const MagneticField *magField) const

	ReferenceTrajectory (unsigned int nPar, unsigned int nHits, const ReferenceTrajectoryBase::Config &config)

SurfaceSide	surfaceSide (const PropagationDirection dir) const

Protected Member Functions inherited from ReferenceTrajectoryBase
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 inherited from ReferenceTrajectoryBase
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< TrajectoryStateOnSurface >	theTsosVec

bool	theValidityFlag

Static Protected Attributes inherited from ReferenceTrajectoryBase
static constexpr unsigned int	nMeasPerHit {2}

Detailed Description

Class implementing the reference trajectory of a single particle in the absence of a magnetic field (B=0, hence the name), i.e. a straight track with only 4 parameters (no curvature parameter).

Given the TrajectoryStateOnSurface at the first hit and the list of all hits the local measurements, derivatives etc. as described in (and accessed via) ReferenceTrajectoryBase are calculated.

The covariance-matrix of the measurements may include effects of multiple-scattering or energy-loss effects or both. This can be defined in the constructor via the variable 'materialEffects (cf. ReferenceTrajectoryBase):

materialEffects = none/multipleScattering/energyLoss/combined

Since no a-priori momentum estimate exists, but is needed for the calculation of material effects, a fixed value is used (e.g. representing the average momentum of the cosmics spectrum).

By default, the mass is assumed to be the muon-mass, but can be changed via a constructor argument.

LIMITATIONS: So far all input hits have to be valid, but invalid hits would be needed to take into account the material effects in them...

Definition at line 35 of file BzeroReferenceTrajectory.h.

Constructor & Destructor Documentation

◆ BzeroReferenceTrajectory()

BzeroReferenceTrajectory::BzeroReferenceTrajectory	(	const TrajectoryStateOnSurface &	tsos,
		const TransientTrackingRecHit::ConstRecHitContainer &	recHits,
		const MagneticField *	magField,
		const reco::BeamSpot &	beamSpot,
		const ReferenceTrajectoryBase::Config &	config
	)

Constructor with Tsos at first hit (in physical order) and list of hits [if (hitsAreReverse) ==> order of hits is in opposite direction compared to the flight of particle, but note that BzeroReferenceTrajectory::recHits() returns the hits always in order of flight], the material effects to be considered and a particle mass and a momentum extimate, the magnetic field of the event is needed for propagations etc.

Definition at line 10 of file BzeroReferenceTrajectory.cc.

References asHepVector(), pwdgSkimBPark_cfi::beamSpot, LocalTrajectoryParameters::charge(), ReferenceTrajectory::construct(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LocalTrajectoryParameters::mixedFormatVector(), LocalTrajectoryParameters::pzSign(), ReferenceTrajectoryBase::recHits(), TrajectoryStateOnSurface::surface(), ReferenceTrajectory::surfaceSide(), ReferenceTrajectoryBase::theDerivatives, theMomentumEstimate, ReferenceTrajectoryBase::theParameters, and ReferenceTrajectoryBase::theValidityFlag.

Referenced by clone().

     : ReferenceTrajectory(tsos.localParameters().mixedFormatVector().kSize, recHits.size(), config),
       theMomentumEstimate(config.momentumEstimate) {
   // no check against magField == 0
 
   // No estimate for momentum of cosmics available -> set to default value.
   theParameters = asHepVector(tsos.localParameters().mixedFormatVector());
   theParameters[0] = 1. / theMomentumEstimate;
 
   LocalTrajectoryParameters locParamWithFixedMomentum(
       asSVector<5>(theParameters), tsos.localParameters().pzSign(), tsos.localParameters().charge());
 
   const TrajectoryStateOnSurface refTsosWithFixedMomentum(
       locParamWithFixedMomentum, tsos.localError(), tsos.surface(), magField, surfaceSide(config.propDir));
 
   if (config.hitsAreReverse) {
     TransientTrackingRecHit::ConstRecHitContainer fwdRecHits;
     fwdRecHits.reserve(recHits.size());
 
     for (TransientTrackingRecHit::ConstRecHitContainer::const_reverse_iterator it = recHits.rbegin();
          it != recHits.rend();
          ++it)
       fwdRecHits.push_back(*it);
 
     theValidityFlag = this->construct(refTsosWithFixedMomentum, fwdRecHits, magField, beamSpot);
   } else {
     theValidityFlag = this->construct(refTsosWithFixedMomentum, recHits, magField, beamSpot);
   }
 
   // Exclude momentum from the parameters and also the derivatives of the measurements w.r.t. the momentum.
   theParameters = theParameters.sub(2, 5);
   theDerivatives = theDerivatives.sub(1, theDerivatives.num_row(), 2, theDerivatives.num_col());
 }

◆ ~BzeroReferenceTrajectory()

BzeroReferenceTrajectory::~BzeroReferenceTrajectory ( )

inlineoverride

Definition at line 50 of file BzeroReferenceTrajectory.h.

50 {}

Member Function Documentation

◆ clone()

BzeroReferenceTrajectory* BzeroReferenceTrajectory::clone ( void ) const

inlineoverridevirtual

Implements ReferenceTrajectoryBase.

Definition at line 52 of file BzeroReferenceTrajectory.h.

References BzeroReferenceTrajectory().

52 { return new BzeroReferenceTrajectory(*this); }

BzeroReferenceTrajectory::BzeroReferenceTrajectory

BzeroReferenceTrajectory(const TrajectoryStateOnSurface &tsos, const TransientTrackingRecHit::ConstRecHitContainer &recHits, const MagneticField *magField, const reco::BeamSpot &beamSpot, const ReferenceTrajectoryBase::Config &config)

Definition: BzeroReferenceTrajectory.cc:10

Member Data Documentation

◆ theMomentumEstimate

double BzeroReferenceTrajectory::theMomentumEstimate

private

Definition at line 55 of file BzeroReferenceTrajectory.h.

Referenced by BzeroReferenceTrajectory().

Public Member Functions

Private Attributes

Additional Inherited Members