GblMeasurement.h

Go to the documentation of this file.

/*
 * GblMeasurement.h
 *
 *  Created on: 31 Mar 2021
 *      Author: kleinwrt
 */
 
#ifndef GBLMEASUREMENT_H_
#define GBLMEASUREMENT_H_
 
#include<iostream>
#include<vector>
#include <array>
#include<math.h>
#include "VMatrix.h"
 
#ifdef GBL_EIGEN_SUPPORT_ROOT
#include "TVectorD.h"
#include "TMatrixD.h"
#include "TMatrixDSym.h"
#include "TMatrixDSymEigen.h"
#endif
 
#include "Eigen/Dense"
 
namespace gbl {
 
typedef Eigen::Matrix<double, 5, 1> Vector5d;
typedef Eigen::Matrix<double, 2, 3> Matrix23d;
typedef Eigen::Matrix<double, 2, 5> Matrix25d;
typedef Eigen::Matrix<double, 3, 2> Matrix32d;
typedef Eigen::Matrix<double, 5, 5> Matrix5d;
 
 
class GblMeasurement {
public:
        GblMeasurement(const Eigen::MatrixXd &aProjection,
                        const Eigen::VectorXd &aResiduals,
                        const Eigen::MatrixXd &aPrecision, double minPrecision = 0.);
        GblMeasurement(const Eigen::VectorXd &aResiduals,
                        const Eigen::MatrixXd &aPrecision, double minPrecision = 0.);
        GblMeasurement(const GblMeasurement&) = default;
        GblMeasurement& operator=(const GblMeasurement&) = default;
        GblMeasurement(GblMeasurement&&) = default;
        GblMeasurement& operator=(GblMeasurement&&) = default;
        virtual ~GblMeasurement();
#ifdef GBL_EIGEN_SUPPORT_ROOT
        // input via ROOT
        GblMeasurement(const TMatrixD &aProjection, const TVectorD &aResiduals,
                        const TVectorD &aPrecision, double minPrecision = 0.);
        GblMeasurement(const TMatrixD &aProjection, const TVectorD &aResiduals,
                        const TMatrixDSym &aPrecision, double minPrecision = 0.);
        GblMeasurement(const TVectorD &aResiduals, const TVectorD &aPrecision,
                        double minPrecision = 0.);
        GblMeasurement(const TVectorD &aResiduals, const TMatrixDSym &aPrecision,
                        double minPrecision = 0.);
        void addLocals(const TMatrixD &aDerivatives);
        void addGlobals(const std::vector<int> &aLabels,
                        const TMatrixD &aDerivatives);
#endif
 
        // input via Eigen
        void addLocals(const Eigen::MatrixXd &aDerivatives);
        void addGlobals(const std::vector<int> &aLabels,
                        const Eigen::MatrixXd &aDerivatives);
        void setEnabled(bool flag);
        bool isEnabled() const;
        unsigned int getMeasDim() const;
        double getMeasPrecMin() const;
        void getMeasurement(Matrix5d &aProjection, Vector5d &aResiduals,
                        Vector5d &aPrecision) const;
        void getMeasTransformation(Eigen::MatrixXd &aTransformation) const;
        unsigned int getNumLocals() const;
        const Eigen::MatrixXd& getLocalDerivatives() const;
        unsigned int getNumGlobals() const;
        void printMeasurement(unsigned int level = 0) const;
        void getGlobalLabels(std::vector<int> &aLabels) const;
        void getGlobalDerivatives(Eigen::MatrixXd &aDerivatives) const;
        void getGlobalLabelsAndDerivatives(unsigned int aRow,
                        std::vector<int> &aLabels, std::vector<double> &aDerivatives) const;
 
private:
        bool enabled; 
        unsigned int measDim; 
        double measPrecMin; 
        Matrix5d measProjection; 
 
        Vector5d measResiduals; 
        Vector5d measPrecision; 
        bool transFlag; 
        Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic,
                        Eigen::ColMajor /* default */, 5, 5> measTransformation; 
        Eigen::MatrixXd localDerivatives; 
        std::vector<int> globalLabels; 
        Eigen::MatrixXd globalDerivatives; 
};
}
 
#endif /* GBLMEASUREMENT_H_ */