MeasurementExtractor.h

Go to the documentation of this file.

#ifndef MeasurementExtractor_H
#define MeasurementExtractor_H
 
#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHit.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
 
class MeasurementExtractor {
public:
  // construct
  MeasurementExtractor(const TrajectoryStateOnSurface &aTSoS) : theTSoS(aTSoS) {}
 
  // access
 
  // Following methods can be overloaded against their argument
  // thus allowing one to have different behaviour for different RecHit types
 
  AlgebraicVector measuredParameters(const TrackingRecHit &);
  AlgebraicSymMatrix measuredError(const TrackingRecHit &);
 
  template <unsigned int D>
  typename AlgebraicROOTObject<D>::Vector measuredParameters(const TrackingRecHit &hit) {
    typedef typename AlgebraicROOTObject<D, 5>::Matrix Mat;
    AlgebraicVector5 par5(theTSoS.localParameters().vector());
    Mat H = asSMatrix<D, 5>(hit.projectionMatrix());
    return H * par5;
  }
 
  template <unsigned int D>
  typename AlgebraicROOTObject<D>::SymMatrix measuredError(const TrackingRecHit &hit) {
    typedef typename AlgebraicROOTObject<D, 5>::Matrix Mat;
    const AlgebraicSymMatrix55 &err5 = theTSoS.localError().matrix();
    Mat H = asSMatrix<D, 5>(hit.projectionMatrix());
    return ROOT::Math::Similarity(H, err5);
  }
 
private:
  const TrajectoryStateOnSurface &theTSoS;
};
 
#endif