CMS 3D CMS Logo

Public Member Functions

KFUpdator Class Reference

#include <KFUpdator.h>

Inheritance diagram for KFUpdator:
TrajectoryStateUpdator

List of all members.

Public Member Functions

virtual KFUpdatorclone () const
 KFUpdator ()
template<unsigned int D>
TrajectoryStateOnSurface update (const TrajectoryStateOnSurface &, const TransientTrackingRecHit &) const
TrajectoryStateOnSurface update (const TrajectoryStateOnSurface &, const TransientTrackingRecHit &) const

Detailed Description

Update trajectory state by combining predicted state and measurement as prescribed in the Kalman Filter algorithm (see R. Fruhwirth, NIM A262 (1987) 444).

x_filtered = x_predicted + K * (measurement - H * x_predicted)

x_filtered, x_predicted filtered and predicted state vectors
measurement measurement vector
H "measurement matrix" projects state vector onto measurement space
K Kalman gain matrix
(formulae for K and error matrix of filtered state not shown)

This implementation works for measurements of all dimensions. It relies on CLHEP double precision vectors and matrices for matrix calculations.

Arguments: TrajectoryState & predicted state
RecHit & reconstructed hit

Initial author: P.Vanlaer 25.02.1999 Ported from ORCA.

Date:
2007/05/09 13:50:25
Revision:
1.3
Author:
vanlaer, cerati

Definition at line 34 of file KFUpdator.h.


Constructor & Destructor Documentation

KFUpdator::KFUpdator ( ) [inline]

Definition at line 40 of file KFUpdator.h.

Referenced by clone().

{}

Member Function Documentation

virtual KFUpdator* KFUpdator::clone ( void  ) const [inline, virtual]

Implements TrajectoryStateUpdator.

Definition at line 48 of file KFUpdator.h.

References KFUpdator().

                                    {
    return new KFUpdator(*this);
  }
TrajectoryStateOnSurface KFUpdator::update ( const TrajectoryStateOnSurface tsos,
const TransientTrackingRecHit aRecHit 
) const [virtual]

Implements TrajectoryStateUpdator.

Definition at line 10 of file KFUpdator.cc.

References TrackingRecHit::dimension(), Exception, TransientTrackingRecHit::hit(), and mergeVDriftHistosByStation::name.

Referenced by CRackTrajectoryBuilder::AddHit(), SiStripElectronSeedGenerator::altCheckHitsAndTSOS(), SeedForPhotonConversion1Leg::buildSeed(), SeedFromConsecutiveHitsCreator::buildSeed(), SeedForPhotonConversionFromQuadruplets::buildSeed(), SiStripElectronSeedGenerator::checkHitsAndTSOS(), SeedFromNuclearInteraction::construct(), OutInConversionSeedFinder::createSeed(), InOutConversionSeedFinder::createSeed(), CosmicMuonSmoother::fit(), ConvBremSeedProducer::produce(), KFBasedPixelFitter::run(), SeedGeneratorForCRack::seeds(), SimpleCosmicBONSeeder::seeds(), SeedGeneratorForCosmics::seeds(), CosmicMuonSmoother::smooth(), KFSwitching1DUpdator::update(), GsfMultiStateUpdator::update(), and KFSwitchingUpdator::update().

                                                                                         {
    switch (aRecHit.dimension()) {
        case 1: return update<1>(tsos,aRecHit);
        case 2: return update<2>(tsos,aRecHit);
        case 3: return update<3>(tsos,aRecHit);
        case 4: return update<4>(tsos,aRecHit);
        case 5: return update<5>(tsos,aRecHit);
    }
    throw cms::Exception("Rec hit of invalid dimension (not 1,2,3,4,5)") <<
         "The value was " << aRecHit.dimension() << 
        ", type is " << typeid(aRecHit).name() << 
        ", persistent rechit type " << (aRecHit.hit() ? typeid(*aRecHit.hit()).name() : "NULL") << "\n";
}
template<unsigned int D>
TrajectoryStateOnSurface KFUpdator::update ( const TrajectoryStateOnSurface ,
const TransientTrackingRecHit  
) const [virtual]