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_ */