TwoBodyDecayTrajectoryState Class Reference

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

References construct().

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

~TwoBodyDecayTrajectoryState()

TwoBodyDecayTrajectoryState::~TwoBodyDecayTrajectoryState ( void  )

inline

Definition at line 26 of file TwoBodyDecayTrajectoryState.h.

{}

Member Function Documentation

construct()

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

private

Definition at line 42 of file TwoBodyDecayTrajectoryState.cc.

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

Referenced by TwoBodyDecayTrajectoryState().

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

decayParameters()

const TwoBodyDecayParameters& TwoBodyDecayTrajectoryState::decayParameters ( void  ) const

inline

Definition at line 31 of file TwoBodyDecayTrajectoryState.h.

References theParameters.

{ return theParameters; }

derivatives()

const Derivatives& TwoBodyDecayTrajectoryState::derivatives ( void  ) const

inline

Definition at line 35 of file TwoBodyDecayTrajectoryState.h.

References theDerivatives.

Referenced by construct().

{ return theDerivatives; }

isValid()

bool TwoBodyDecayTrajectoryState::isValid ( void  ) const

inline

Definition at line 28 of file TwoBodyDecayTrajectoryState.h.

References theValidityFlag.

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

{ return theValidityFlag; }

particleMass()

double TwoBodyDecayTrajectoryState::particleMass ( void  ) const

inline

Definition at line 30 of file TwoBodyDecayTrajectoryState.h.

References theParticleMass.

{ return theParticleMass; }

primaryMass()

double TwoBodyDecayTrajectoryState::primaryMass ( void  ) const

inline

Definition at line 39 of file TwoBodyDecayTrajectoryState.h.

References thePrimaryMass.

{ return thePrimaryMass; }

primaryWidth()

double TwoBodyDecayTrajectoryState::primaryWidth ( void  ) const

inline

Definition at line 40 of file TwoBodyDecayTrajectoryState.h.

References thePrimaryWidth.

{ return 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.

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

Referenced by construct().

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

rescaleError()

void TwoBodyDecayTrajectoryState::rescaleError

(

double

scale

)

Definition at line 35 of file TwoBodyDecayTrajectoryState.cc.

References pfClustersFromCombinedCaloHF_cfi::scale, theOriginalTsos, and theRefittedTsos.

Referenced by TwoBodyDecayConstraintProducer::produce().

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

setError()

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

private

Definition at line 133 of file TwoBodyDecayTrajectoryState.cc.

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

Referenced by construct().

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

trajectoryStates()

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

inline

Definition at line 32 of file TwoBodyDecayTrajectoryState.h.

References theOriginalTsos, and theRefittedTsos.

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

                                                                                   {
    return useRefittedState ? theRefittedTsos : theOriginalTsos;
  }

Member Data Documentation

nDecayParam

const unsigned int TwoBodyDecayTrajectoryState::nDecayParam = TwoBodyDecayParameters::dimension

staticprivate

Definition at line 68 of file TwoBodyDecayTrajectoryState.h.

nLocalParam

const unsigned int TwoBodyDecayTrajectoryState::nLocalParam = 5

staticprivate

Definition at line 67 of file TwoBodyDecayTrajectoryState.h.

theDerivatives

Derivatives TwoBodyDecayTrajectoryState::theDerivatives

private

Definition at line 60 of file TwoBodyDecayTrajectoryState.h.

Referenced by construct(), and derivatives().

theOriginalTsos

TsosContainer TwoBodyDecayTrajectoryState::theOriginalTsos

private

Definition at line 61 of file TwoBodyDecayTrajectoryState.h.

Referenced by construct(), rescaleError(), and trajectoryStates().

theParameters

TwoBodyDecayParameters TwoBodyDecayTrajectoryState::theParameters

private

Definition at line 59 of file TwoBodyDecayTrajectoryState.h.

Referenced by construct(), decayParameters(), and setError().

theParticleMass

double TwoBodyDecayTrajectoryState::theParticleMass

private

Definition at line 57 of file TwoBodyDecayTrajectoryState.h.

Referenced by construct(), and particleMass().

thePrimaryMass

double TwoBodyDecayTrajectoryState::thePrimaryMass

private

Definition at line 64 of file TwoBodyDecayTrajectoryState.h.

Referenced by primaryMass().

thePrimaryWidth

double TwoBodyDecayTrajectoryState::thePrimaryWidth

private

Definition at line 65 of file TwoBodyDecayTrajectoryState.h.

Referenced by primaryWidth().

theRefittedTsos

TsosContainer TwoBodyDecayTrajectoryState::theRefittedTsos

private

Definition at line 62 of file TwoBodyDecayTrajectoryState.h.

Referenced by construct(), rescaleError(), and trajectoryStates().

theValidityFlag

bool TwoBodyDecayTrajectoryState::theValidityFlag

private

Definition at line 55 of file TwoBodyDecayTrajectoryState.h.

Referenced by construct(), and isValid().