GFTools.h

Go to the documentation of this file.

/* Copyright 2011, Technische Universitaet Muenchen,
   Author: Karl Bicker

   This file is part of GENFIT.

   GENFIT is free software: you can redistribute it and/or modify
   it under the terms of the GNU Lesser General Public License as published
   by the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.

   GENFIT is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public License
   along with GENFIT.  If not, see <http://www.gnu.org/licenses/>.
*/
/* @{ */

#ifndef GFTOOLS_H
#define GFTOOLS_H

#include <cmath>
#include <GFAbsTrackRep.h>
#include <GFDetPlane.h>
#include <GFException.h>
#include <GFTrack.h>
#include <TMath.h>
#include <TMatrixT.h>
#include <TDecompSVD.h>
//#include <TDecompBK.h>
//#include <TDecompChol.h>

namespace GFTools {
        
        TMatrixT<double> getSmoothedPos(GFTrack* trk, unsigned int irep, unsigned int ihit);

        TMatrixT<double> getSmoothedCov(GFTrack* trk, unsigned int irep, unsigned int ihit);

        bool getSmoothedData(GFTrack* trk, unsigned int irep, unsigned int ihit, TMatrixT<double>& smoothed_state, TMatrixT<double>& smoothed_cov);

        bool getSmoothedData(GFTrack* trk, unsigned int irep, unsigned int ihit, TMatrixT<double>& smoothed_state, TMatrixT<double>& smoothed_cov, GFDetPlane& smoothing_plane);

        bool getSmoothedData(GFTrack* trk, unsigned int irep, unsigned int ihit, TMatrixT<double>& smoothed_state, TMatrixT<double>& smoothed_cov, GFDetPlane& smoothing_plane, TMatrixT<double>& auxInfo);

        GFDetPlane getSmoothingPlane(GFTrack* trk, unsigned int irep, unsigned int ihit);

        bool getBiasedSmoothedData(GFTrack* trk, unsigned int irep, unsigned int ihit, TMatrixT<double>& smoothed_state, TMatrixT<double>& smoothed_cov);

        bool getBiasedSmoothedData(GFTrack* trk, unsigned int irep, unsigned int ihit, TMatrixT<double>& smoothed_state, TMatrixT<double>& smoothed_cov, GFDetPlane& smoothing_plane);

        bool getBiasedSmoothedData(GFTrack* trk, unsigned int irep, unsigned int ihit, TMatrixT<double>& smoothed_state, TMatrixT<double>& smoothed_cov, GFDetPlane& smoothing_plane, TMatrixT<double>& auxInfo);

        TMatrixT<double> getBiasedSmoothedPos(GFTrack* trk, unsigned int irep, unsigned int ihit);

        TMatrixT<double> getBiasedSmoothedCov(GFTrack* trk, unsigned int irep, unsigned int ihit);


        void invertMatrix(const TMatrixT<double>& mat, TMatrixT<double>& inv);
        double getSmoothedChiSqu(GFTrack* const trk, unsigned int irep, unsigned int ihit);

}

#endif