DualReferenceTrajectory Class Reference

#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.

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

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

Referenced by clone().

~DualReferenceTrajectory()

DualReferenceTrajectory::~DualReferenceTrajectory ( )

inlineoverride

Definition at line 48 of file DualReferenceTrajectory.h.

{}

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.

{ return new DualReferenceTrajectory(*this); }

References DualReferenceTrajectory().

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.

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

References pwdgSkimBPark_cfi::beamSpot, ReferenceTrajectoryBase::breakPoints, ReferenceTrajectoryBase::derivatives(), Exception, extractParameters(), ReferenceTrajectoryBase::isValid(), materialEffects_, ReferenceTrajectoryBase::measurementErrors(), ReferenceTrajectoryBase::measurements(), ReferenceTrajectoryBase::nMeasPerHit, pfTracksFromL1Tracks_cfi::nParam, ReferenceTrajectoryBase::numberOfHitMeas(), ReferenceTrajectoryBase::numberOfVirtualMeas(), ReferenceTrajectoryBase::recHits(), ReferenceTrajectoryBase::theDerivatives, ReferenceTrajectoryBase::theMeasurements, ReferenceTrajectoryBase::theMeasurementsCov, ReferenceTrajectoryBase::theNumberOfPars, ReferenceTrajectoryBase::theParameters, ReferenceTrajectoryBase::theRecHits, ReferenceTrajectoryBase::theTrajectoryPositionCov, ReferenceTrajectoryBase::theTrajectoryPositions, ReferenceTrajectoryBase::theTsosVec, ReferenceTrajectoryBase::trajectoryPositionErrors(), ReferenceTrajectoryBase::trajectoryPositions(), and ReferenceTrajectoryBase::trajectoryStates().

Referenced by DualReferenceTrajectory().

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.

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

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

extractParameters()

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

protectedvirtual

Reimplemented in DualBzeroReferenceTrajectory.

Definition at line 138 of file DualReferenceTrajectory.cc.

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

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

Referenced by construct().

oppositeDirection()

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

inlineprotected

Definition at line 71 of file DualReferenceTrajectory.h.

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

References alongMomentum, anyDirection, and oppositeToMomentum.

Referenced by construct().

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().