TrajectoryStateCombiner Class Reference

#include <TrajectoryStateCombiner.h>

Public Types
typedef TrajectoryStateOnSurface	TSOS

Public Member Functions
TSOS	combine (const TSOS &pTsos1, const TSOS &pTsos2) const
TSOS	operator() (const TSOS &pTsos1, const TSOS &pTsos2) const

Detailed Description

Combines the information from two trajectory states via a weighted mean. The input states should not be correlated. Ported from ORCA

Author

todorov, cerati

Definition at line 13 of file TrajectoryStateCombiner.h.

Member Typedef Documentation

TSOS

typedef TrajectoryStateOnSurface TrajectoryStateCombiner::TSOS

Definition at line 15 of file TrajectoryStateCombiner.h.

Member Function Documentation

combine()

TrajectoryStateOnSurface TrajectoryStateCombiner::combine	(	const TSOS &	pTsos1,
		const TSOS &	pTsos2
	)		const

Definition at line 6 of file TrajectoryStateCombiner.cc.

                                                                                                    {
  auto pzSign = Tsos1.localParameters().pzSign();
  AlgebraicVector5&& x1 = Tsos1.localParameters().vector();
  AlgebraicVector5&& x2 = Tsos2.localParameters().vector();
  const AlgebraicSymMatrix55& C1 = (Tsos1.localError().matrix());
  const AlgebraicSymMatrix55& C2 = (Tsos2.localError().matrix());
 
  AlgebraicSymMatrix55&& Csum = C1 + C2;
  bool ok = invertPosDefMatrix(Csum);
  AlgebraicMatrix55&& K = C1 * Csum;
 
  if (!ok) {
    if (!(C1(0, 0) == 0.0 && C2(0, 0) == 0.0))  //do not make noise about obviously bad input
      edm::LogWarning("MatrixInversionFailure")
          << "the inversion of the combined error matrix failed. Impossible to get a combined state."
          << "\nmatrix 1:" << C1 << "\nmatrix 2:" << C2;
    return TSOS();
  }
 
  AlgebraicVector5&& xcomb = x1 + K * (x2 - x1);
  //AlgebraicSymMatrix55 Ccomb; Ccomb.assign(K*C2);
  AlgebraicSymMatrix55&& Ccomb = (AlgebraicMatrix55(K * C2)).LowerBlock();
 
  return TSOS(LocalTrajectoryParameters(xcomb, pzSign),
              LocalTrajectoryError(Ccomb),
              Tsos1.surface(),
              &(Tsos1.globalParameters().magneticField()));
}

References TrajectoryStateOnSurface::globalParameters(), invertPosDefMatrix(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), GlobalTrajectoryParameters::magneticField(), LocalTrajectoryError::matrix(), convertSQLiteXML::ok, LocalTrajectoryParameters::pzSign(), TrajectoryStateOnSurface::surface(), LocalTrajectoryParameters::vector(), testProducerWithPsetDescEmpty_cfi::x1, and testProducerWithPsetDescEmpty_cfi::x2.

Referenced by OverlapValidation::analyze(), CSCPairResidualsConstraint::dphidzFromTrack(), AlignmentMonitorTemplate::event(), AlignmentMonitorMuonResiduals::event(), operator()(), AlignmentMuonHIPTrajectorySelector::produce(), TrackInfoProducerAlgorithm::run(), HIPAlignmentAlgorithm::run(), and TrajectoryAtInvalidHit::TrajectoryAtInvalidHit().

operator()()

TSOS TrajectoryStateCombiner::operator() ( const TSOS &  pTsos1, const TSOS &  pTsos2  ) const

inline

Definition at line 19 of file TrajectoryStateCombiner.h.

{ return combine(pTsos1, pTsos2); }

References combine().