#include <TwoBodyDecayTrajectoryState.h>

Public Types
typedef std::pair< AlgebraicMatrix, AlgebraicMatrix >	Derivatives

typedef std::pair< TrajectoryStateOnSurface, TrajectoryStateOnSurface >	TsosContainer

Public Member Functions
const TwoBodyDecayParameters &	decayParameters (void) const

const Derivatives &	derivatives (void) const

bool	isValid (void) const

double	particleMass (void) const

double	primaryMass (void) const

double	primaryWidth (void) const

void	rescaleError (double scale)

const TsosContainer &	trajectoryStates (bool useRefittedState=true) const

	TwoBodyDecayTrajectoryState (const TsosContainer &tsos, const TwoBodyDecay &tbd, double particleMass, const MagneticField *magField, bool propagateErrors=false)

	~TwoBodyDecayTrajectoryState (void)

Private Member Functions
void	construct (const MagneticField *magField, bool propagateErrors)

bool	propagateSingleState (const FreeTrajectoryState &fts, const GlobalTrajectoryParameters &gtp, const AlgebraicMatrix &startDeriv, const Surface &surface, const MagneticField *magField, TrajectoryStateOnSurface &tsos, AlgebraicMatrix &endDeriv) const

void	setError (FreeTrajectoryState &fts, AlgebraicMatrix &derivative) const

Private Attributes
Derivatives	theDerivatives

TsosContainer	theOriginalTsos

TwoBodyDecayParameters	theParameters

double	theParticleMass

double	thePrimaryMass

double	thePrimaryWidth

TsosContainer	theRefittedTsos

bool	theValidityFlag

Static Private Attributes
static const unsigned int	nDecayParam = TwoBodyDecayParameters::dimension

static const unsigned int	nLocalParam = 5

Detailed Description

Definition at line 12 of file TwoBodyDecayTrajectoryState.h.

Member Typedef Documentation

◆ Derivatives

typedef std::pair<AlgebraicMatrix, AlgebraicMatrix> TwoBodyDecayTrajectoryState::Derivatives

Definition at line 15 of file TwoBodyDecayTrajectoryState.h.

◆ TsosContainer

typedef std::pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface> TwoBodyDecayTrajectoryState::TsosContainer

Definition at line 14 of file TwoBodyDecayTrajectoryState.h.

Constructor & Destructor Documentation

◆ TwoBodyDecayTrajectoryState()

TwoBodyDecayTrajectoryState::TwoBodyDecayTrajectoryState	(	const TsosContainer &	tsos,
		const TwoBodyDecay &	tbd,
		double	particleMass,
		const MagneticField *	magField,
		bool	propagateErrors = `false`
	)

The constructor takes the two trajectory states that are to be updated (typically the innermost trajectory states of two tracks) and the decay parameters.

Definition at line 20 of file TwoBodyDecayTrajectoryState.cc.

     : theValidityFlag(false),
       theParticleMass(particleMass),
       theParameters(tbd.decayParameters()),
       theDerivatives(AlgebraicMatrix(nLocalParam, nDecayParam), AlgebraicMatrix(nLocalParam, nDecayParam)),
       theOriginalTsos(tsos),
       thePrimaryMass(tbd.primaryMass()),
       thePrimaryWidth(tbd.primaryWidth()) {
   construct(magField, propagateErrors);
 }

References construct().

◆ ~TwoBodyDecayTrajectoryState()

TwoBodyDecayTrajectoryState::~TwoBodyDecayTrajectoryState ( void )

inline

Definition at line 26 of file TwoBodyDecayTrajectoryState.h.

26 {}

Member Function Documentation

◆ construct()

void TwoBodyDecayTrajectoryState::construct	(	const MagneticField *	magField,
		bool	propagateErrors
	)

private

Definition at line 42 of file TwoBodyDecayTrajectoryState.cc.

                                                                                                {
   // construct global trajectory parameters at the starting point
   TwoBodyDecayModel tbdDecayModel(theParameters[TwoBodyDecayParameters::mass], theParticleMass);
   pair<AlgebraicVector, AlgebraicVector> secondaryMomenta = tbdDecayModel.cartesianSecondaryMomenta(theParameters);
  
   GlobalPoint vtx(theParameters[TwoBodyDecayParameters::x],
                   theParameters[TwoBodyDecayParameters::y],
                   theParameters[TwoBodyDecayParameters::z]);
  
   GlobalVector p1(secondaryMomenta.first[0], secondaryMomenta.first[1], secondaryMomenta.first[2]);
  
   GlobalVector p2(secondaryMomenta.second[0], secondaryMomenta.second[1], secondaryMomenta.second[2]);
  
   GlobalTrajectoryParameters gtp1(vtx, p1, theOriginalTsos.first.charge(), magField);
   FreeTrajectoryState fts1(gtp1);
  
   GlobalTrajectoryParameters gtp2(vtx, p2, theOriginalTsos.second.charge(), magField);
   FreeTrajectoryState fts2(gtp2);
  
   // contruct derivatives at the starting point
   TwoBodyDecayDerivatives tbdDerivatives(theParameters[TwoBodyDecayParameters::mass], theParticleMass);
   pair<AlgebraicMatrix, AlgebraicMatrix> derivatives = tbdDerivatives.derivatives(theParameters);
  
   AlgebraicMatrix deriv1(6, 9, 0);
   deriv1.sub(1, 1, AlgebraicMatrix(3, 3, 1));
   deriv1.sub(4, 4, derivatives.first);
  
   AlgebraicMatrix deriv2(6, 9, 0);
   deriv2.sub(1, 1, AlgebraicMatrix(3, 3, 1));
   deriv2.sub(4, 4, derivatives.second);
  
   // compute errors of initial states
   if (propagateErrors) {
     setError(fts1, deriv1);
     setError(fts2, deriv2);
   }
  
   // propgate states and derivatives from the starting points to the end points
   bool valid1 = propagateSingleState(
       fts1, gtp1, deriv1, theOriginalTsos.first.surface(), magField, theRefittedTsos.first, theDerivatives.first);
  
   bool valid2 = propagateSingleState(
       fts2, gtp2, deriv2, theOriginalTsos.second.surface(), magField, theRefittedTsos.second, theDerivatives.second);
  
   theValidityFlag = valid1 && valid2;
  
   return;
 }

References TwoBodyDecayModel::cartesianSecondaryMomenta(), TwoBodyDecayDerivatives::derivatives(), derivatives(), TwoBodyDecayParameters::mass, p1, p2, propagateSingleState(), setError(), theDerivatives, theOriginalTsos, theParameters, theParticleMass, theRefittedTsos, theValidityFlag, extraflags_cff::vtx, TwoBodyDecayParameters::x, TwoBodyDecayParameters::y, and TwoBodyDecayParameters::z.

Referenced by TwoBodyDecayTrajectoryState().

◆ decayParameters()

const TwoBodyDecayParameters& TwoBodyDecayTrajectoryState::decayParameters ( void ) const

inline

Definition at line 31 of file TwoBodyDecayTrajectoryState.h.

31 { return theParameters; }

References theParameters.

◆ derivatives()

const Derivatives& TwoBodyDecayTrajectoryState::derivatives ( void ) const

inline

Definition at line 35 of file TwoBodyDecayTrajectoryState.h.

35 { return theDerivatives; }

References theDerivatives.

Referenced by construct().

◆ isValid()

bool TwoBodyDecayTrajectoryState::isValid ( void ) const

inline

Definition at line 28 of file TwoBodyDecayTrajectoryState.h.

28 { return theValidityFlag; }

References theValidityFlag.

Referenced by ntupleDataFormat._Object::_checkIsValid(), TwoBodyDecayTrajectoryFactory::constructTrajectories(), TwoBodyDecayMomConstraintProducer::momentaAtInnermostSurface(), and core.AutoHandle.AutoHandle::ReallyLoad().

◆ particleMass()

double TwoBodyDecayTrajectoryState::particleMass ( void ) const

inline

Definition at line 30 of file TwoBodyDecayTrajectoryState.h.

30 { return theParticleMass; }

References theParticleMass.

◆ primaryMass()

double TwoBodyDecayTrajectoryState::primaryMass ( void ) const

inline

Definition at line 39 of file TwoBodyDecayTrajectoryState.h.

39 { return thePrimaryMass; }

References thePrimaryMass.

◆ primaryWidth()

double TwoBodyDecayTrajectoryState::primaryWidth ( void ) const

inline

Definition at line 40 of file TwoBodyDecayTrajectoryState.h.

40 { return thePrimaryWidth; }

References thePrimaryWidth.

◆ propagateSingleState()

bool TwoBodyDecayTrajectoryState::propagateSingleState	(	const FreeTrajectoryState &	fts,
		const GlobalTrajectoryParameters &	gtp,
		const AlgebraicMatrix &	startDeriv,
		const Surface &	surface,
		const MagneticField *	magField,
		TrajectoryStateOnSurface &	tsos,
		AlgebraicMatrix &	endDeriv
	)		const

private

Definition at line 91 of file TwoBodyDecayTrajectoryState.cc.

                                                                                         {
   AnalyticalPropagator propagator(magField);
  
   // propagate state
   pair<TrajectoryStateOnSurface, double> tsosWithPath = propagator.propagateWithPath(fts, surface);
  
   // check if propagation was successful
   if (!tsosWithPath.first.isValid())
     return false;
  
   // jacobian for transformation from cartesian to curvilinear frame at the starting point
   JacobianCartesianToCurvilinear cartToCurv(gtp);
   const AlgebraicMatrix56& matCartToCurv = cartToCurv.jacobian();
  
   // jacobian in curvilinear frame for propagation from the starting point to the end point
   AnalyticalCurvilinearJacobian curvJac(
       gtp, tsosWithPath.first.globalPosition(), tsosWithPath.first.globalMomentum(), tsosWithPath.second);
   const AlgebraicMatrix55& matCurvJac = curvJac.jacobian();
  
   // jacobian for transformation from curvilinear to local frame at the end point
   JacobianCurvilinearToLocal curvToLoc(surface, tsosWithPath.first.localParameters(), *magField);
   const AlgebraicMatrix55& matCurvToLoc = curvToLoc.jacobian();
  
   AlgebraicMatrix56 tmpDeriv = matCurvToLoc * matCurvJac * matCartToCurv;
   AlgebraicMatrix hepMatDeriv(asHepMatrix(tmpDeriv));
   //AlgebraicMatrix hepMatDeriv = asHepMatrix< 5, 6 >( tmpDeriv );
  
   // replace original state with new state
   tsos = tsosWithPath.first;
  
   // propagate derivative matrix
   endDeriv = hepMatDeriv * startDeriv;
  
   return true;
 }

References AnalyticalPropagator_cfi::AnalyticalPropagator, asHepMatrix(), JacobianCartesianToCurvilinear::jacobian(), JacobianCurvilinearToLocal::jacobian(), AnalyticalCurvilinearJacobian::jacobian(), and TrackCandidateProducer_cfi::propagator.

Referenced by construct().

◆ rescaleError()

void TwoBodyDecayTrajectoryState::rescaleError ( double scale )

Definition at line 35 of file TwoBodyDecayTrajectoryState.cc.

                                                            {
   theOriginalTsos.first.rescaleError(scale);
   theOriginalTsos.second.rescaleError(scale);
   theRefittedTsos.first.rescaleError(scale);
   theRefittedTsos.second.rescaleError(scale);
 }

References L1EGammaCrystalsEmulatorProducer_cfi::scale, theOriginalTsos, and theRefittedTsos.

◆ setError()

void TwoBodyDecayTrajectoryState::setError	(	FreeTrajectoryState &	fts,
		AlgebraicMatrix &	derivative
	)		const

private

Definition at line 133 of file TwoBodyDecayTrajectoryState.cc.

                                                                                                       {
   AlgebraicSymMatrix ftsCartesianError(theParameters.covariance().similarity(derivative));
   fts.setCartesianError(asSMatrix<6>(ftsCartesianError));
 }

References TwoBodyDecayParameters::covariance(), funct::derivative(), FreeTrajectoryState::setCartesianError(), and theParameters.

Referenced by construct().

◆ trajectoryStates()

const TsosContainer& TwoBodyDecayTrajectoryState::trajectoryStates ( bool useRefittedState = true ) const

inline

Definition at line 32 of file TwoBodyDecayTrajectoryState.h.

                                                                                    {
     return useRefittedState ? theRefittedTsos : theOriginalTsos;
   }

References theOriginalTsos, and theRefittedTsos.

Referenced by TwoBodyDecayTrajectoryFactory::constructTrajectories(), and TwoBodyDecayMomConstraintProducer::momentaAtInnermostSurface().