CollinearFitAtTM.h

Go to the documentation of this file.

#ifndef CollinearFitAtTM_h_
#define CollinearFitAtTM_h_

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

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

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

  CollinearFitAtTM ();

  typedef ROOT::Math::SVector<double,6> ResultVector;
  typedef ROOT::Math::SMatrix<double,6,6,ROOT::Math::MatRepSym<double,6> > ResultMatrix;
  bool fit (const TrajectoryMeasurement& tm,
        ResultVector& parameters,
        ResultMatrix& covariance,
        double& chi2);
  bool fit (const AlgebraicVector5& fwdParameters, 
        const AlgebraicSymMatrix55& fwdCovariance,
        const AlgebraicVector5& bwdParameters, 
        const AlgebraicSymMatrix55& bwdCovariance,
        const LocalPoint& hitPosition, const LocalError& hitErrors,
        ResultVector& parameters,
        ResultMatrix& covariance,
        double& chi2);
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;
};

#endif