GblPoint.h

Go to the documentation of this file.

/*
 * GblPoint.h
 *
 *  Created on: Aug 18, 2011
 *      Author: kleinwrt
 */
 
#ifndef GBLPOINT_H_
#define GBLPOINT_H_
 
#include<iostream>
#include<vector>
#include<math.h>
#include <stdexcept>
#include "TVectorD.h"
#include "TMatrixD.h"
#include "TMatrixDSym.h"
#include "TMatrixDSymEigen.h"
 
#include "Math/SMatrix.h"
#include "Math/SVector.h"
typedef ROOT::Math::SMatrix<double, 2> SMatrix22;
typedef ROOT::Math::SMatrix<double, 2, 3> SMatrix23;
typedef ROOT::Math::SMatrix<double, 2, 5> SMatrix25;
typedef ROOT::Math::SMatrix<double, 2, 7> SMatrix27;
typedef ROOT::Math::SMatrix<double, 3, 2> SMatrix32;
typedef ROOT::Math::SMatrix<double, 3> SMatrix33;
typedef ROOT::Math::SMatrix<double, 5> SMatrix55;
typedef ROOT::Math::SVector<double, 2> SVector2;
typedef ROOT::Math::SVector<double, 5> SVector5;
 
namespace gbl {
 
 
class GblPoint {
public:
  GblPoint(const TMatrixD &aJacobian);
  GblPoint(const SMatrix55 &aJacobian);
  virtual ~GblPoint();
  void addMeasurement(const TMatrixD &aProjection, const TVectorD &aResiduals,
      const TVectorD &aPrecision, double minPrecision = 0.);
  void addMeasurement(const TMatrixD &aProjection, const TVectorD &aResiduals,
      const TMatrixDSym &aPrecision, double minPrecision = 0.);
  void addMeasurement(const TVectorD &aResiduals, const TVectorD &aPrecision,
      double minPrecision = 0.);
  void addMeasurement(const TVectorD &aResiduals,
      const TMatrixDSym &aPrecision, double minPrecision = 0.);
  unsigned int hasMeasurement() const;
  void getMeasurement(SMatrix55 &aProjection, SVector5 &aResiduals,
      SVector5 &aPrecision) const;
  void getMeasTransformation(TMatrixD &aTransformation) const;
  void addScatterer(const TVectorD &aResiduals, const TVectorD &aPrecision);
  void addScatterer(const TVectorD &aResiduals,
      const TMatrixDSym &aPrecision);
  bool hasScatterer() const;
  void getScatterer(SMatrix22 &aTransformation, SVector2 &aResiduals,
      SVector2 &aPrecision) const;
  void getScatTransformation(TMatrixD &aTransformation) const;
  void addLocals(const TMatrixD &aDerivatives);
  unsigned int getNumLocals() const;
  const TMatrixD& getLocalDerivatives() const;
  void addGlobals(const std::vector<int> &aLabels,
      const TMatrixD &aDerivatives);
  unsigned int getNumGlobals() const;
  std::vector<int> getGlobalLabels() const;
  const TMatrixD& getGlobalDerivatives() const;
  void setLabel(unsigned int aLabel);
  unsigned int getLabel() const;
  void setOffset(int anOffset);
  int getOffset() const;
  const SMatrix55& getP2pJacobian() const;
  void addPrevJacobian(const SMatrix55 &aJac);
  void addNextJacobian(const SMatrix55 &aJac);
  void getDerivatives(int aDirection, SMatrix22 &matW, SMatrix22 &matWJ,
      SVector2 &vecWd) const;
  void printPoint(unsigned int level = 0) const;
 
private:
  unsigned int theLabel; 
  int theOffset; 
  SMatrix55 p2pJacobian; 
  SMatrix55 prevJacobian; 
  SMatrix55 nextJacobian; 
  unsigned int measDim; 
  SMatrix55 measProjection; 
  SVector5 measResiduals; 
  SVector5 measPrecision; 
  bool transFlag; 
  TMatrixD measTransformation; 
  bool scatFlag; 
  SMatrix22 scatTransformation; 
  SVector2 scatResiduals; 
  SVector2 scatPrecision; 
  TMatrixD localDerivatives; 
  std::vector<int> globalLabels; 
  TMatrixD globalDerivatives; 
};
}
#endif /* GBLPOINT_H_ */