GFTools.h

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

/* 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(const GFTrack* trk, unsigned int irep, unsigned int ihit);

  TVector3 getSmoothedPosXYZ(const GFTrack* trk, unsigned int irep, unsigned int ihit, bool *ret = 0);
  TVector3 getBiasedSmoothedPosXYZ(const GFTrack* trk, unsigned int irep, unsigned int ihit, bool *ret = 0);

  //TVector3 getSmoothedMomXYZ(const GFTrack* trk, unsigned int irep, unsigned int ihit);
  TVector3 getBiasedSmoothedMomXYZ(const GFTrack* trk, unsigned int irep, unsigned int ihit, bool *ret = 0);

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

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

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

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

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

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

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

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

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

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

  double getTrackLength(const GFTrack* trk, unsigned int irep, unsigned int startHit=0, unsigned int endHit=0);

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

  void updateRepSmoothed(GFTrack* trk, unsigned int irep, unsigned int ihit);

  unsigned int getClosestHit(const GFTrack* trk, unsigned int irep, const TVector3& pos, double& distance, bool checkEveryHit = false);

}

#endif