#include <DualBzeroReferenceTrajectory.h>

Inheritance diagram for DualBzeroReferenceTrajectory:

Public Types
typedef TransientTrackingRecHit::ConstRecHitContainer	ConstRecHitContainer

Public Types inherited from DualReferenceTrajectory
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
DualBzeroReferenceTrajectory *	clone () const override

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

	~DualBzeroReferenceTrajectory () override

Public Member Functions inherited from DualReferenceTrajectory
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

virtual ReferenceTrajectory *	construct (const TrajectoryStateOnSurface &referenceTsos, const ConstRecHitContainer &recHits, double mass, MaterialEffects materialEffects, const PropagationDirection propDir, const MagneticField *magField, bool useBeamSpot, const reco::BeamSpot &beamSpot) const

AlgebraicVector	extractParameters (const TrajectoryStateOnSurface &referenceTsos) const override

Protected Member Functions inherited from DualReferenceTrajectory
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)

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

Protected Attributes
double	theMomentumEstimate

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

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

Detailed Description

Definition at line 37 of file DualBzeroReferenceTrajectory.h.

Member Typedef Documentation

◆ ConstRecHitContainer

typedef TransientTrackingRecHit::ConstRecHitContainer DualBzeroReferenceTrajectory::ConstRecHitContainer

Definition at line 39 of file DualBzeroReferenceTrajectory.h.

Constructor & Destructor Documentation

◆ DualBzeroReferenceTrajectory()

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

Definition at line 12 of file DualBzeroReferenceTrajectory.cc.

     : DualReferenceTrajectory(tsos.localParameters().mixedFormatVector().kSize - 1,
                               numberOfUsedRecHits(forwardRecHits) + numberOfUsedRecHits(backwardRecHits) - 1,
                               config),
       theMomentumEstimate(config.momentumEstimate) {
   theValidityFlag = construct(tsos, forwardRecHits, backwardRecHits, magField, beamSpot);
 }

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

Referenced by clone().

◆ ~DualBzeroReferenceTrajectory()

DualBzeroReferenceTrajectory::~DualBzeroReferenceTrajectory ( )

inlineoverride

Definition at line 48 of file DualBzeroReferenceTrajectory.h.

48 {}

Member Function Documentation

◆ clone()

DualBzeroReferenceTrajectory* DualBzeroReferenceTrajectory::clone ( void ) const

inlineoverridevirtual

Implements ReferenceTrajectoryBase.

Definition at line 50 of file DualBzeroReferenceTrajectory.h.

50 { return new DualBzeroReferenceTrajectory(*this); }

References DualBzeroReferenceTrajectory().

◆ construct() [1/3]

bool DualReferenceTrajectory::construct

protected

internal method to calculate members

Definition at line 37 of file DualReferenceTrajectory.cc.

                                                                       {
   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/3]

ReferenceTrajectory * DualReferenceTrajectory::construct

protected

Definition at line 122 of file DualReferenceTrajectory.cc.

                                                                                           {
   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);
 }

◆ construct() [3/3]

ReferenceTrajectory * DualBzeroReferenceTrajectory::construct	(	const TrajectoryStateOnSurface &	referenceTsos,
		const ConstRecHitContainer &	recHits,
		double	mass,
		MaterialEffects	materialEffects,
		const PropagationDirection	propDir,
		const MagneticField *	magField,
		bool	useBeamSpot,
		const reco::BeamSpot &	beamSpot
	)		const

protectedvirtual

Definition at line 25 of file DualBzeroReferenceTrajectory.cc.

                                                                                                  {
   if (materialEffects >= breakPoints)
     throw cms::Exception("BadConfig") << "[DualBzeroReferenceTrajectory::construct] Wrong MaterialEffects: "
                                       << materialEffects;
  
   ReferenceTrajectoryBase::Config config(materialEffects, propDir, mass, theMomentumEstimate);
   config.useBeamSpot = useBeamSpot;
   config.hitsAreReverse = false;
   return new BzeroReferenceTrajectory(referenceTsos, recHits, magField, beamSpot, config);
 }

References pwdgSkimBPark_cfi::beamSpot, ReferenceTrajectoryBase::breakPoints, submitPVResolutionJobs::config, Exception, EgHLTOffHistBins_cfi::mass, ReferenceTrajectoryBase::recHits(), theMomentumEstimate, and HLT_FULL_cff::useBeamSpot.

Referenced by DualBzeroReferenceTrajectory().

◆ extractParameters()

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

overrideprotectedvirtual

Reimplemented from DualReferenceTrajectory.

Definition at line 43 of file DualBzeroReferenceTrajectory.cc.

                                                                                                                    {
   AlgebraicVector param = asHepVector<5>(referenceTsos.localParameters().mixedFormatVector());
   return param.sub(2, 5);
 }

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

Member Data Documentation

◆ theMomentumEstimate

double DualBzeroReferenceTrajectory::theMomentumEstimate

protected

Definition at line 66 of file DualBzeroReferenceTrajectory.h.

Referenced by construct().

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ ConstRecHitContainer

Constructor & Destructor Documentation

◆ DualBzeroReferenceTrajectory()

◆ ~DualBzeroReferenceTrajectory()

Member Function Documentation

◆ clone()

◆ construct() [1/3]

◆ construct() [2/3]

◆ construct() [3/3]

◆ extractParameters()

Member Data Documentation

◆ theMomentumEstimate