GblPoint.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file GblPoint.h

/*
 * 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:
        explicit GblPoint(const TMatrixD &aJacobian);
        explicit 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_ */