GsfMultiStateUpdator Class Reference

#include <GsfMultiStateUpdator.h>

Inheritance diagram for GsfMultiStateUpdator:

Public Member Functions
GsfMultiStateUpdator *	clone () const override
	GsfMultiStateUpdator ()
TrajectoryStateOnSurface	update (const TrajectoryStateOnSurface &, const TrackingRecHit &) const override
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::GsfMultiStateUpdator ( )

inline

Definition at line 18 of file GsfMultiStateUpdator.h.

Referenced by clone().

{}

Member Function Documentation

clone()

GsfMultiStateUpdator* GsfMultiStateUpdator::clone ( void  ) const

inlineoverridevirtual

Implements TrajectoryStateUpdator.

Definition at line 22 of file GsfMultiStateUpdator.h.

References GsfMultiStateUpdator().

{ return new GsfMultiStateUpdator(*this); }

update()

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

overridevirtual

Implements TrajectoryStateUpdator.

Definition at line 13 of file GsfMultiStateUpdator.cc.

References TrajectoryStateOnSurface::curvilinearError(), TrajectoryStateOnSurface::globalParameters(), mps_fire::i, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), GlobalTrajectoryParameters::magneticField(), CurvilinearTrajectoryError::matrix(), LocalTrajectoryError::posDef(), mps_fire::result, TrajectoryStateOnSurface::surface(), KFUpdator::update(), LocalTrajectoryError::valid(), PosteriorWeightsCalculator::weights(), and HLT_2022v12_cff::weights.

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

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

  auto&& 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 (auto const& tsosI : predictedComponents) {
    TrajectoryStateOnSurface updatedTSOS = KFUpdator().update(tsosI, aRecHit);

    if (double det;
        updatedTSOS.isValid() && updatedTSOS.localError().valid() && updatedTSOS.localError().posDef() &&
        (det = 0., updatedTSOS.curvilinearError().matrix().Sub<AlgebraicSymMatrix22>(0, 0).Det(det) && det > 0) &&
        (det = 0., updatedTSOS.curvilinearError().matrix().Sub<AlgebraicSymMatrix33>(0, 0).Det(det) && det > 0) &&
        (det = 0., updatedTSOS.curvilinearError().matrix().Sub<AlgebraicSymMatrix44>(0, 0).Det(det) && det > 0) &&
        (det = 0., updatedTSOS.curvilinearError().matrix().Det2(det) && det > 0)) {
      result.addState(TrajectoryStateOnSurface(weights[i],
                                               updatedTSOS.localParameters(),
                                               updatedTSOS.localError(),
                                               updatedTSOS.surface(),
                                               &(tsos.globalParameters().magneticField()),
                                               tsosI.surfaceSide()));
    } else {
      edm::LogError("GsfMultiStateUpdator") << "KF updated state " << i << " is invalid. skipping.";
    }
    ++i;
  }

  return result.combinedState();
}