#include <BasicMultiTrajectoryState.h>

Inheritance diagram for BasicMultiTrajectoryState:

Public Types
using	Components = BasicTrajectoryState::Components

Public Types inherited from BasicTrajectoryState
typedef BasicTrajectoryState	BTSOS

using	Components = std::vector< TrajectoryStateOnSurface >

typedef Proxy::pointer	pointer

typedef ProxyBase11< BTSOS >	Proxy

typedef SurfaceSideDefinition::SurfaceSide	SurfaceSide

typedef Surface	SurfaceType

Public Member Functions
	BasicMultiTrajectoryState (const std::vector< TSOS > &tsvec)

	BasicMultiTrajectoryState ()

bool	canUpdateLocalParameters () const override

pointer	clone () const override

Components const &	components () const override

void	rescaleError (double factor)

bool	singleState () const override

void	update (const LocalTrajectoryParameters &p, const Surface &aSurface, const MagneticField *field, const SurfaceSide side) override

void	update (double weight, const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const Surface &aSurface, const MagneticField *field, const SurfaceSide side) override

Public Member Functions inherited from BasicTrajectoryState
	BasicTrajectoryState ()

	BasicTrajectoryState (const SurfaceType &aSurface)
	construct invalid trajectory state (without parameters) More...

	BasicTrajectoryState (const FreeTrajectoryState &fts, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)

	BasicTrajectoryState (const FreeTrajectoryState &fts)

	BasicTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)

	BasicTrajectoryState (const LocalTrajectoryParameters &par, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)

	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CartesianTrajectoryError &err, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)

	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CurvilinearTrajectoryError &err, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)

	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)

template<typename... Args>
	BasicTrajectoryState (double iweight, Args &&...args)

const CartesianTrajectoryError	cartesianError () const

TrackCharge	charge () const

const CurvilinearTrajectoryError &	curvilinearError () const

FreeTrajectoryState const *	freeTrajectoryState (bool withErrors=true) const

GlobalVector	globalDirection () const

GlobalVector	globalMomentum () const

const GlobalTrajectoryParameters &	globalParameters () const

GlobalPoint	globalPosition () const

bool	hasError () const

bool	isValid () const

LocalVector	localDirection () const

const LocalTrajectoryError &	localError () const

LocalVector	localMomentum () const

const LocalTrajectoryParameters &	localParameters () const

LocalPoint	localPosition () const

const MagneticField *	magneticField () const

void	rescaleError (double factor)

CurvilinearTrajectoryError &	setCurvilinearError ()

double	signedInverseMomentum () const

const SurfaceType &	surface () const

SurfaceSide	surfaceSide () const
	Position relative to material, defined relative to momentum vector. More...

double	transverseCurvature () const

virtual void	update (const LocalTrajectoryParameters &p, const SurfaceSide side) final

virtual void	update (const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const SurfaceSide side) final

double	weight () const

virtual	~BasicTrajectoryState ()

Private Types
typedef TrajectoryStateOnSurface	TSOS

Private Member Functions
void	combine ()

Private Attributes
Components	theStates

Additional Inherited Members
Static Public Member Functions inherited from BasicTrajectoryState
template<typename T , typename... Args>
static std::shared_ptr< BTSOS >	build (Args &&...args)

template<typename T , typename... Args>
static std::shared_ptr< BTSOS >	churn (Args &&...args)

Detailed Description

Class which combines a set of components of a Gaussian mixture into a single component. Given all the components of a mixture, it calculates the mean and covariance matrix of the entire mixture. This combiner class can also be used in the process of transforming a Gaussian mixture into another Gaussian mixture with a smaller number of components. The relevant formulas can be found in R. Fruhwirth, Computer Physics Communications 100 (1997), 1.

Definition at line 17 of file BasicMultiTrajectoryState.h.

Member Typedef Documentation

using BasicMultiTrajectoryState::Components = BasicTrajectoryState::Components

Definition at line 35 of file BasicMultiTrajectoryState.h.

typedef TrajectoryStateOnSurface BasicMultiTrajectoryState::TSOS

private

Definition at line 18 of file BasicMultiTrajectoryState.h.

Constructor & Destructor Documentation

BasicMultiTrajectoryState::BasicMultiTrajectoryState ( const std::vector< TSOS > & tsvec )

explicit

Definition at line 8 of file BasicMultiTrajectoryState.cc.

References combine(), Exception, mps_fire::i, and UNLIKELY.

     : BasicTrajectoryState(tsvec.front().surface()), theStates(tsvec) {
   // only accept planes!!
   const BoundPlane* bp = dynamic_cast<const BoundPlane*>(&tsvec.begin()->surface());
   if UNLIKELY (bp == nullptr)
     throw cms::Exception("LogicError") << "MultiTrajectoryState constructed on cylinder";
 
   for (auto i = tsvec.begin(); i != tsvec.end(); i++) {
     if UNLIKELY (!i->isValid()) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState constructed with invalid state";
     }
     if UNLIKELY (i->hasError() != tsvec.front().hasError()) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states with and without errors";
     }
     if UNLIKELY (&i->surface() != &tsvec.front().surface()) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states with different surfaces";
     }
     if UNLIKELY (i->surfaceSide() != tsvec.front().surfaceSide()) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states defined before and after material";
     }
     if UNLIKELY (i->localParameters().pzSign() * tsvec.front().localParameters().pzSign() < 0.) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states with different signs of local p_z";
     }
   }
   //
   combine();
 }

BasicMultiTrajectoryState::BasicMultiTrajectoryState ( )

inline

Definition at line 23 of file BasicMultiTrajectoryState.h.

23 {}

Member Function Documentation

bool BasicMultiTrajectoryState::canUpdateLocalParameters ( ) const

inlineoverridevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 39 of file BasicMultiTrajectoryState.h.

39 { return false; }

pointer BasicMultiTrajectoryState::clone ( void ) const

inlineoverridevirtual

Implements BasicTrajectoryState.

Definition at line 33 of file BasicMultiTrajectoryState.h.

33 { return build<BasicMultiTrajectoryState>(*this); }

void BasicMultiTrajectoryState::combine ( )

private

Definition at line 47 of file BasicMultiTrajectoryState.cc.

References change_name::diff, SiStripPI::mean, theStates, UNLIKELY, BasicTrajectoryState::update(), and BasicTrajectoryState::weight().

Referenced by BasicMultiTrajectoryState(), and rescaleError().

                                         {
   const std::vector<TrajectoryStateOnSurface>& tsos = theStates;
 
   if UNLIKELY (tsos.empty()) {
     edm::LogError("MultiTrajectoryStateCombiner") << "Trying to collapse empty set of trajectory states!";
     return;
   }
 
   double pzSign = tsos.front().localParameters().pzSign();
   for (std::vector<TrajectoryStateOnSurface>::const_iterator it = tsos.begin(); it != tsos.end(); it++) {
     if UNLIKELY (it->localParameters().pzSign() != pzSign) {
       edm::LogError("MultiTrajectoryStateCombiner")
           << "Trying to collapse trajectory states with different signs on p_z!";
       return;
     }
   }
 
   if UNLIKELY (tsos.size() == 1) {
     BasicTrajectoryState::update(tsos.front().weight(),
                                  tsos.front().localParameters(),
                                  tsos.front().localError(),
                                  tsos.front().surface(),
                                  tsos.front().magneticField(),
                                  tsos.front().surfaceSide());
     return;
   }
 
   double sumw = 0.;
   //int dim = tsos.front().localParameters().vector().num_row();
   AlgebraicVector5 mean;
   AlgebraicSymMatrix55 covarPart1, covarPart2, covtmp;
   for (auto it1 = tsos.begin(); it1 != tsos.end(); it1++) {
     auto weight = it1->weight();
     auto const& param = it1->localParameters().vector();
     sumw += weight;
     mean += weight * param;
     covarPart1 += weight * it1->localError().matrix();
     for (auto it2 = it1 + 1; it2 != tsos.end(); it2++) {
       AlgebraicVector5 diff = param - it2->localParameters().vector();
       ROOT::Math::AssignSym::Evaluate(covtmp, ROOT::Math::TensorProd(diff, diff));
       covarPart2 += (weight * it2->weight()) * covtmp;
     }
   }
   double sumwI = 1.0 / sumw;
   mean *= sumwI;
   covarPart1 *= sumwI;
   covarPart2 *= (sumwI * sumwI);
   AlgebraicSymMatrix55 covar = covarPart1 + covarPart2;
 
   BasicTrajectoryState::update(sumw,
                                LocalTrajectoryParameters(mean, pzSign),
                                LocalTrajectoryError(covar),
                                tsos.front().surface(),
                                tsos.front().magneticField(),
                                tsos.front().surfaceSide());
 }

Components const& BasicMultiTrajectoryState::components ( ) const

inlineoverridevirtual

Implements BasicTrajectoryState.

Definition at line 36 of file BasicMultiTrajectoryState.h.

References theStates.

36 { return theStates; }

BasicMultiTrajectoryState::theStates

Components theStates

Definition: BasicMultiTrajectoryState.h:53

void BasicMultiTrajectoryState::rescaleError ( double factor )

Rescaling the error of the mixture with a given factor. Please note that this rescaling is imposed on each of the components of the mixture and does therefore not exactly correspond to rescaling theCombinedState with the same factor.

Definition at line 36 of file BasicMultiTrajectoryState.cc.

References combine(), theStates, and UNLIKELY.

                                                           {
   if UNLIKELY (theStates.empty()) {
     edm::LogError("BasicMultiTrajectoryState") << "Trying to rescale errors of empty MultiTrajectoryState!";
     return;
   }
 
   for (auto& is : theStates)
     is.rescaleError(factor);
   combine();
 }

bool BasicMultiTrajectoryState::singleState ( ) const

inlineoverridevirtual

Implements BasicTrajectoryState.

Definition at line 37 of file BasicMultiTrajectoryState.h.

37 { return false; }

void BasicMultiTrajectoryState::update	(	const LocalTrajectoryParameters &	p,
		const Surface &	aSurface,
		const MagneticField *	field,
		const SurfaceSide	side
	)

overridevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 104 of file BasicMultiTrajectoryState.cc.

References Exception.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                                {
   throw cms::Exception("LogicError",
                        "BasicMultiTrajectoryState::update(LocalTrajectoryParameters, Surface, ...) called even if "
                        "canUpdateLocalParameters() is false");
 }

void BasicMultiTrajectoryState::update	(	double	weight,
		const LocalTrajectoryParameters &	p,
		const LocalTrajectoryError &	err,
		const Surface &	aSurface,
		const MagneticField *	field,
		const SurfaceSide	side
	)

overridevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 113 of file BasicMultiTrajectoryState.cc.

References Exception.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                                {
   throw cms::Exception("LogicError",
                        "BasicMultiTrajectoryState::update(LocalTrajectoryParameters, LocalTrajectoryError, ...) called "
                        "even if canUpdateLocalParameters() is false");
 }

Member Data Documentation

Components BasicMultiTrajectoryState::theStates

private

Definition at line 53 of file BasicMultiTrajectoryState.h.

Referenced by combine(), components(), and rescaleError().

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation