CollinearFitAtTM.h

Go to the documentation of this file.

#ifndef CollinearFitAtTM2_h_
#define CollinearFitAtTM2_h_

#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/Measurement1D.h"
#include "DataFormats/Math/interface/Vector.h"
#include "DataFormats/Math/interface/Error.h"

// #include "Workspace/TrajectoryMeasurementFits/interface/RandomVector.h"

class CollinearFitAtTM2 {
public:
  enum { ParQpIn=0, ParQpOut, ParDxDz, ParDyDz, ParX, ParY };

  CollinearFitAtTM2 (const TrajectoryMeasurement& tm);
  CollinearFitAtTM2 (const AlgebraicVector5& fwdParameters, 
            const AlgebraicSymMatrix55& fwdCovariance,
            const AlgebraicVector5& bwdParameters, 
            const AlgebraicSymMatrix55& bwdCovariance,
            const LocalPoint& hitPosition, 
            const LocalError& hitErrors);

  typedef ROOT::Math::SVector<double,6> ResultVector;
  typedef ROOT::Math::SMatrix<double,6,6,ROOT::Math::MatRepSym<double,6> > ResultMatrix;
  bool valid () const {return valid_;}
  double chi2 () const {return chi2_;}
  int ndof () const {return ndof_;}
  const ResultVector& parameters () const {return parameters_;}
  const ResultMatrix& covariance () const {return covariance_;}
  Measurement1D deltaP () const;

private:
  void initJacobian ();
  bool fit (const AlgebraicVector5& fwdParameters, 
        const AlgebraicSymMatrix55& fwdCovariance,
        const AlgebraicVector5& bwdParameters, 
        const AlgebraicSymMatrix55& bwdCovariance,
        const LocalPoint& hitPosition, const LocalError& hitErrors);
  
private:
  ROOT::Math::SMatrix<double,12,6> jacobian_;
  ROOT::Math::SVector<double,12> measurements_;  
  ROOT::Math::SMatrix<double,12,12,ROOT::Math::MatRepSym<double,12> > weightMatrix_;
  ROOT::Math::SVector<double,6> projectedMeasurements_;
//   RandomVector randomGenerator;

  bool valid_;
  ResultVector parameters_;
  ResultMatrix covariance_;
  double chi2_;
  int ndof_;
};

#endif