#include <DualReferenceTrajectory.h>

Inheritance diagram for DualReferenceTrajectory:

Public Types
typedef TransientTrackingRecHit::ConstRecHitContainer	ConstRecHitContainer

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

typedef ReferenceCountingPointer< ReferenceTrajectoryBase >	ReferenceTrajectoryPtr

Public Member Functions
DualReferenceTrajectory *	clone () const override

	DualReferenceTrajectory (const TrajectoryStateOnSurface &tsos, const ConstRecHitContainer &forwardRecHits, const ConstRecHitContainer &backwardRecHits, const MagneticField *magField, const reco::BeamSpot &beamSpot, const ReferenceTrajectoryBase::Config &config)

	~DualReferenceTrajectory () 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

Protected Member Functions
virtual bool	construct (const TrajectoryStateOnSurface &referenceTsos, const ConstRecHitContainer &forwardRecHits, const ConstRecHitContainer &backwardRecHits, const MagneticField *magField, const reco::BeamSpot &beamSpot)

virtual ReferenceTrajectory *	construct (const TrajectoryStateOnSurface &referenceTsos, const ConstRecHitContainer &recHits, const MagneticField *magField, const reco::BeamSpot &beamSpot, const bool revertDirection=false) const

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

virtual AlgebraicVector	extractParameters (const TrajectoryStateOnSurface &referenceTsos) const

const PropagationDirection	oppositeDirection (const PropagationDirection propDir) 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

Private Attributes
const double	mass_

const MaterialEffects	materialEffects_

const PropagationDirection	propDir_

const bool	useBeamSpot_

Additional Inherited Members
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

Definition at line 37 of file DualReferenceTrajectory.h.

Member Typedef Documentation

◆ ConstRecHitContainer

typedef TransientTrackingRecHit::ConstRecHitContainer DualReferenceTrajectory::ConstRecHitContainer

Definition at line 39 of file DualReferenceTrajectory.h.

Constructor & Destructor Documentation

◆ DualReferenceTrajectory() [1/2]

DualReferenceTrajectory::DualReferenceTrajectory	(	const TrajectoryStateOnSurface &	tsos,
		const ConstRecHitContainer &	forwardRecHits,
		const ConstRecHitContainer &	backwardRecHits,
		const MagneticField *	magField,
		const reco::BeamSpot &	beamSpot,
		const ReferenceTrajectoryBase::Config &	config
	)

Definition at line 11 of file DualReferenceTrajectory.cc.

References pwdgSkimBPark_cfi::beamSpot, construct(), and ReferenceTrajectoryBase::theValidityFlag.

Referenced by clone().

     : ReferenceTrajectoryBase(tsos.localParameters().mixedFormatVector().kSize,
                               numberOfUsedRecHits(forwardRecHits) + numberOfUsedRecHits(backwardRecHits) - 1,
                               0,
                               0),
       mass_(config.mass),
       materialEffects_(config.materialEffects),
       propDir_(config.propDir),
       useBeamSpot_(config.useBeamSpot) {
   theValidityFlag = this->construct(tsos, forwardRecHits, backwardRecHits, magField, beamSpot);
 }

◆ ~DualReferenceTrajectory()

DualReferenceTrajectory::~DualReferenceTrajectory ( )

inlineoverride

Definition at line 48 of file DualReferenceTrajectory.h.

48 {}

◆ DualReferenceTrajectory() [2/2]

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

protected

Definition at line 28 of file DualReferenceTrajectory.cc.

     : ReferenceTrajectoryBase(nPar, nHits, 0, 0),
       mass_(config.mass),
       materialEffects_(config.materialEffects),
       propDir_(config.propDir),
       useBeamSpot_(config.useBeamSpot) {}

Member Function Documentation

◆ clone()

DualReferenceTrajectory* DualReferenceTrajectory::clone ( void ) const

inlineoverridevirtual

Implements ReferenceTrajectoryBase.

Definition at line 50 of file DualReferenceTrajectory.h.

References DualReferenceTrajectory().

50 { return new DualReferenceTrajectory(*this); }

DualReferenceTrajectory::DualReferenceTrajectory

DualReferenceTrajectory(const TrajectoryStateOnSurface &tsos, const ConstRecHitContainer &forwardRecHits, const ConstRecHitContainer &backwardRecHits, const MagneticField *magField, const reco::BeamSpot &beamSpot, const ReferenceTrajectoryBase::Config &config)

Definition: DualReferenceTrajectory.cc:11

◆ construct() [1/2]

bool DualReferenceTrajectory::construct	(	const TrajectoryStateOnSurface &	referenceTsos,
		const ConstRecHitContainer &	forwardRecHits,
		const ConstRecHitContainer &	backwardRecHits,
		const MagneticField *	magField,
		const reco::BeamSpot &	beamSpot
	)

protectedvirtual

internal method to calculate members

Definition at line 37 of file DualReferenceTrajectory.cc.

Referenced by DualReferenceTrajectory().

                                                                       {
   if (materialEffects_ >= breakPoints)
     throw cms::Exception("BadConfig") << "[DualReferenceTrajectory::construct] Wrong MaterialEffects: "
                                       << materialEffects_;
 
   ReferenceTrajectoryBase* fwdTraj = construct(refTsos, forwardRecHits, magField, beamSpot);
 
   // set flag for opposite direction to true
   ReferenceTrajectoryBase* bwdTraj = construct(refTsos, backwardRecHits, magField, beamSpot, true);
 
   if (!(fwdTraj->isValid() && bwdTraj->isValid())) {
     delete fwdTraj;
     delete bwdTraj;
     return false;
   }
 
   //
   // Combine both reference trajactories to a dual reference trajectory
   //
 
   const std::vector<TrajectoryStateOnSurface>& fwdTsosVec = fwdTraj->trajectoryStates();
   const std::vector<TrajectoryStateOnSurface>& bwdTsosVec = bwdTraj->trajectoryStates();
   theTsosVec.insert(theTsosVec.end(), fwdTsosVec.begin(), fwdTsosVec.end());
   theTsosVec.insert(theTsosVec.end(), ++bwdTsosVec.begin(), bwdTsosVec.end());
 
   const ConstRecHitContainer& fwdRecHits = fwdTraj->recHits();
   const ConstRecHitContainer& bwdRecHits = bwdTraj->recHits();
   theRecHits.insert(theRecHits.end(), fwdRecHits.begin(), fwdRecHits.end());
   theRecHits.insert(theRecHits.end(), ++bwdRecHits.begin(), bwdRecHits.end());
 
   theParameters = extractParameters(refTsos);
 
   unsigned int nParam = theNumberOfPars;
   unsigned int nFwdMeas = fwdTraj->numberOfHitMeas();
   unsigned int nBwdMeas = bwdTraj->numberOfHitMeas();
   unsigned int nFwdBP = fwdTraj->numberOfVirtualMeas();
   unsigned int nBwdBP = bwdTraj->numberOfVirtualMeas();
   unsigned int nMeas = nFwdMeas + nBwdMeas - nMeasPerHit;
 
   theMeasurements.sub(1, fwdTraj->measurements().sub(1, nFwdMeas));
   theMeasurements.sub(nFwdMeas + 1, bwdTraj->measurements().sub(nMeasPerHit + 1, nBwdMeas));
 
   theMeasurementsCov.sub(1, fwdTraj->measurementErrors().sub(1, nFwdMeas));
   theMeasurementsCov.sub(nFwdMeas + 1, bwdTraj->measurementErrors().sub(nMeasPerHit + 1, nBwdMeas));
 
   theTrajectoryPositions.sub(1, fwdTraj->trajectoryPositions());
   theTrajectoryPositions.sub(nFwdMeas + 1, bwdTraj->trajectoryPositions().sub(nMeasPerHit + 1, nBwdMeas));
 
   theTrajectoryPositionCov.sub(1, fwdTraj->trajectoryPositionErrors());
   theTrajectoryPositionCov.sub(nFwdMeas + 1, bwdTraj->trajectoryPositionErrors().sub(nMeasPerHit + 1, nBwdMeas));
 
   theDerivatives.sub(1, 1, fwdTraj->derivatives().sub(1, nFwdMeas, 1, nParam));
   theDerivatives.sub(nFwdMeas + 1, 1, bwdTraj->derivatives().sub(nMeasPerHit + 1, nBwdMeas, 1, nParam));
 
   // for the break points
   // DUAL with break points makes no sense: (MS) correlations between the two parts are lost !
   if (nFwdBP > 0) {
     theMeasurements.sub(nMeas + 1, fwdTraj->measurements().sub(nFwdMeas + 1, nFwdMeas + nFwdBP));
     theMeasurementsCov.sub(nMeas + 1, fwdTraj->measurementErrors().sub(nFwdMeas + 1, nFwdMeas + nFwdBP));
     theDerivatives.sub(1, nParam + 1, fwdTraj->derivatives().sub(1, nFwdMeas, nParam + 1, nParam + nFwdBP));
     theDerivatives.sub(nMeas + 1,
                        nParam + 1,
                        fwdTraj->derivatives().sub(nFwdMeas + 1, nFwdMeas + nFwdBP, nParam + 1, nParam + nFwdBP));
   }
   if (nBwdBP > 0) {
     theMeasurements.sub(nMeas + nFwdBP + 1, bwdTraj->measurements().sub(nBwdMeas + 1, nBwdMeas + nBwdBP));
     theMeasurementsCov.sub(nMeas + nFwdBP + 1, bwdTraj->measurementErrors().sub(nBwdMeas + 1, nBwdMeas + nBwdBP));
     theDerivatives.sub(nFwdMeas + 1,
                        nParam + nFwdBP + 1,
                        bwdTraj->derivatives().sub(nMeasPerHit + 1, nBwdMeas, nParam + 1, nParam + nBwdBP));
     theDerivatives.sub(nMeas + nFwdBP + 1,
                        nParam + nFwdBP + 1,
                        bwdTraj->derivatives().sub(nBwdMeas + 1, nBwdMeas + nBwdBP, nParam + 1, nParam + nBwdBP));
   }
 
   delete fwdTraj;
   delete bwdTraj;
 
   return true;
 }

◆ construct() [2/2]

ReferenceTrajectory * DualReferenceTrajectory::construct	(	const TrajectoryStateOnSurface &	referenceTsos,
		const ConstRecHitContainer &	recHits,
		const MagneticField *	magField,
		const reco::BeamSpot &	beamSpot,
		const bool	revertDirection = `false`
	)		const

protectedvirtual

Definition at line 122 of file DualReferenceTrajectory.cc.

References pwdgSkimBPark_cfi::beamSpot, ReferenceTrajectoryBase::breakPoints, submitPVResolutionJobs::config, Exception, mass_, materialEffects_, oppositeDirection(), propDir_, ReferenceTrajectoryBase::recHits(), and useBeamSpot_.

                                                                                           {
   if (materialEffects_ >= breakPoints)
     throw cms::Exception("BadConfig") << "[DualReferenceTrajectory::construct] Wrong MaterialEffects: "
                                       << materialEffects_;
 
   ReferenceTrajectoryBase::Config config(
       materialEffects_, (revertDirection ? oppositeDirection(propDir_) : propDir_), mass_);
   config.useBeamSpot = useBeamSpot_;
   config.hitsAreReverse = false;
   return new ReferenceTrajectory(referenceTsos, recHits, magField, beamSpot, config);
 }

◆ extractParameters()

AlgebraicVector DualReferenceTrajectory::extractParameters ( const TrajectoryStateOnSurface & referenceTsos ) const

protectedvirtual

Reimplemented in DualBzeroReferenceTrajectory.

Definition at line 138 of file DualReferenceTrajectory.cc.

References TrajectoryStateOnSurface::localParameters(), and LocalTrajectoryParameters::mixedFormatVector().

Referenced by construct().

                                                                                                               {
   return asHepVector<5>(referenceTsos.localParameters().mixedFormatVector());
 }

◆ oppositeDirection()

const PropagationDirection DualReferenceTrajectory::oppositeDirection ( const PropagationDirection propDir ) const

inlineprotected

Definition at line 71 of file DualReferenceTrajectory.h.

References alongMomentum, anyDirection, and oppositeToMomentum.

Referenced by construct().

                                                                                                 {
     return (propDir == anyDirection) ? anyDirection : ((propDir == alongMomentum) ? oppositeToMomentum : alongMomentum);
   }

Member Data Documentation

◆ mass_

const double DualReferenceTrajectory::mass_

private

Definition at line 76 of file DualReferenceTrajectory.h.

Referenced by construct().

◆ materialEffects_

const MaterialEffects DualReferenceTrajectory::materialEffects_

private

Definition at line 77 of file DualReferenceTrajectory.h.

Referenced by construct().

◆ propDir_

const PropagationDirection DualReferenceTrajectory::propDir_

private

Definition at line 78 of file DualReferenceTrajectory.h.

Referenced by construct().

◆ useBeamSpot_

const bool DualReferenceTrajectory::useBeamSpot_

private

Definition at line 79 of file DualReferenceTrajectory.h.

Referenced by construct().

Public Types

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ ConstRecHitContainer

Constructor & Destructor Documentation

◆ DualReferenceTrajectory() [1/2]

◆ ~DualReferenceTrajectory()

◆ DualReferenceTrajectory() [2/2]

Member Function Documentation

◆ clone()

◆ construct() [1/2]

◆ construct() [2/2]

◆ extractParameters()

◆ oppositeDirection()

Member Data Documentation

◆ mass_

◆ materialEffects_

◆ propDir_

◆ useBeamSpot_