GsfMultiStateUpdator Class Reference

#include <GsfMultiStateUpdator.h>

Inheritance diagram for GsfMultiStateUpdator:

Public Member Functions
virtual GsfMultiStateUpdator *	clone () const
	GsfMultiStateUpdator ()
TrajectoryStateOnSurface	update (const TrajectoryStateOnSurface &, const TrackingRecHit &) const
Public Member Functions inherited from TrajectoryStateUpdator
	TrajectoryStateUpdator ()
virtual	~TrajectoryStateUpdator ()

Detailed Description

Class which updates a Gaussian mixture trajectory state with the information from a reconstructed hit according to the Gaussian-sum filter (GSF) strategy. The relevant formulas can be derived from those described in R. Fruhwirth, Computer Physics Communications 100 (1997), 1.

Definition at line 16 of file GsfMultiStateUpdator.h.

Constructor & Destructor Documentation

GsfMultiStateUpdator::GsfMultiStateUpdator ( )

inline

Definition at line 20 of file GsfMultiStateUpdator.h.

Referenced by clone().

{}

Member Function Documentation

virtual GsfMultiStateUpdator* GsfMultiStateUpdator::clone ( void  ) const

inlinevirtual

Implements TrajectoryStateUpdator.

Definition at line 25 of file GsfMultiStateUpdator.h.

References GsfMultiStateUpdator().

                                               {
    return new GsfMultiStateUpdator(*this);
  }

TrajectoryStateOnSurface GsfMultiStateUpdator::update ( const TrajectoryStateOnSurface &  tsos, const TrackingRecHit &  aRecHit  ) const

virtual

Implements TrajectoryStateUpdator.

Definition at line 13 of file GsfMultiStateUpdator.cc.

References MultiTrajectoryStateAssembler::addState(), MultiTrajectoryStateAssembler::combinedState(), TrajectoryStateOnSurface::components(), TrajectoryStateOnSurface::globalParameters(), i, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), GlobalTrajectoryParameters::magneticField(), mps_fire::result, TrajectoryStateOnSurface::surface(), KFUpdator::update(), PosteriorWeightsCalculator::weights(), and create_public_pileup_plots::weights.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), GsfConstraintAtVertex::constrainAtPoint(), 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().

                                                                   {
  
  std::vector<TrajectoryStateOnSurface> && predictedComponents = tsos.components();
  if (predictedComponents.empty()) {
    edm::LogError("GsfMultiStateUpdator") << "Trying to update trajectory state with zero components! " ;
    return TrajectoryStateOnSurface();
  }

  std::vector<double> && weights = PosteriorWeightsCalculator(predictedComponents).weights(aRecHit);
  if ( weights.empty() ) {
    edm::LogError("GsfMultiStateUpdator") << " no weights could be retreived. invalid updated state !.";
    return TrajectoryStateOnSurface();
  }

  MultiTrajectoryStateAssembler result;

  int i = 0;
  for (std::vector<TrajectoryStateOnSurface>::const_iterator iter = predictedComponents.begin();
       iter != predictedComponents.end(); iter++) {
    TrajectoryStateOnSurface updatedTSOS = KFUpdator().update(*iter, aRecHit);
    if (updatedTSOS.isValid()){
      result.addState(TrajectoryStateOnSurface(weights[i], 
                                               updatedTSOS.localParameters(),
                           updatedTSOS.localError(), updatedTSOS.surface(), 
                           &(tsos.globalParameters().magneticField()),
                           (*iter).surfaceSide()
                                              ));
    i++;
    }
    else{
      edm::LogError("GsfMultiStateUpdator") << "one of the KF updated state is invalid. skipping.";
    }

  }

  return result.combinedState();
}