BasicMultiTrajectoryState Class Reference

#include <BasicMultiTrajectoryState.h>

Inheritance diagram for BasicMultiTrajectoryState:

Public Member Functions
	BasicMultiTrajectoryState (const std::vector< TSOS > &tsvec)
	BasicMultiTrajectoryState ()
virtual bool	canUpdateLocalParameters () const
pointer	clone () const override
virtual std::vector < TrajectoryStateOnSurface >	components () const
void	rescaleError (double factor)
virtual void	update (const LocalTrajectoryParameters &p, const Surface &aSurface, const MagneticField *field, const SurfaceSide side) override
virtual 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...
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)
virtual void	update (const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const SurfaceSide side)
double	weight () const
virtual	~BasicTrajectoryState ()

Private Types
typedef TrajectoryStateOnSurface	TSOS

Private Member Functions
void	combine ()

Private Attributes
std::vector< TSOS >	theStates

Additional Inherited Members
Public Types inherited from BasicTrajectoryState
typedef BasicTrajectoryState	BTSOS
typedef Proxy::pointer	pointer
typedef ProxyBase11< BTSOS >	Proxy
typedef SurfaceSideDefinition::SurfaceSide	SurfaceSide
typedef Surface	SurfaceType

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

typedef TrajectoryStateOnSurface BasicMultiTrajectoryState::TSOS

private

Definition at line 19 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(), edm::hlt::Exception, 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==0 )
           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 25 of file BasicMultiTrajectoryState.h.

{}

Member Function Documentation

virtual bool BasicMultiTrajectoryState::canUpdateLocalParameters ( ) const

inlinevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 45 of file BasicMultiTrajectoryState.h.

{ return false; }

pointer BasicMultiTrajectoryState::clone ( void  ) const

inlineoverridevirtual

Implements BasicTrajectoryState.

Definition at line 36 of file BasicMultiTrajectoryState.h.

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

void BasicMultiTrajectoryState::combine ( )

private

Definition at line 56 of file BasicMultiTrajectoryState.cc.

References diffTreeTool::diff, timingPdfMaker::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++) {
    double weight = it1->weight();
    AlgebraicVector5 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()
                              );
}

virtual std::vector<TrajectoryStateOnSurface> BasicMultiTrajectoryState::components ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 40 of file BasicMultiTrajectoryState.h.

References theStates.

                                                                 {
    return theStates;
  }

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 42 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 (std::vector<TSOS>::iterator it = theStates.begin(); it != theStates.end(); it++) {
    it->rescaleError(factor);
  }
  combine();
}

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

overridevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 118 of file BasicMultiTrajectoryState.cc.

References edm::hlt::Exception.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
  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 129 of file BasicMultiTrajectoryState.cc.

References edm::hlt::Exception.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

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

Member Data Documentation

std::vector<TSOS> BasicMultiTrajectoryState::theStates

private

Definition at line 59 of file BasicMultiTrajectoryState.h.

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