SKalman.hh

Go to the documentation of this file.

#ifndef __SKALMAN_HH
#define __SKALMAN_HH
// ********************************************************************
//
// source:
//
// type:      source code
//
// created:   10. Apr 2001
//
// author:    Thorsten Glebe
//            HERA-B Collaboration
//            Max-Planck-Institut fuer Kernphysik
//            Saupfercheckweg 1
//            69117 Heidelberg
//            Germany
//            E-mail: T.Glebe@mpi-hd.mpg.de
//
// Description: A fixed size class for Kalman filter Vertex fit
//
// changes:
// 10 Apr 2001 (TG) creation
// 12 Apr 2001 (TG) p(), tvec(), evec(), pvec(), F(), ES(), DS() added
// 18 Apr 2001 (TG) doc++ comments added, removed use_mom template parameter
// 24 Apr 2001 (TG) tx(), ty() added
// 04 Mai 2001 (TG) added Track inheritance, some functions are dummy
// 11 Jun 2001 (TG) added COV(), cov_tx(), cov_ty(), cov_p()
// 29 Jun 2001 (TG) added charge(), pz(), energy(), xf(), rap(), xfabs() removed
// 05 Jul 2001 (TG) collect() added
// 09 Jul 2001 (TG) added isValid(), valid(), invalid() from Track
// 10 Jul 2001 (TG) operator==() added
// 09 Jan 2002 (TG) changed return type of propagate() to bool
//
// ********************************************************************
#include <cmath>
#include <iosfwd>
#include "vt++/VtTrack.hh"
 
template <unsigned int NTR> class SVertex;
 
//==============================================================================
// SKalman
//==============================================================================
template <unsigned int NTR>
class SKalman : public Track {
public:
  SKalman();
  SKalman(const Track& t, const SVertex<NTR>& v);
 
  float x()   const;
  float y()   const;
  float x(float z)   const;
  float y(float z)   const;
  float z()   const;
  float tx() const;
  float ty() const;
  float p() const;
  float pt()  const;
  float pz()    const;
  float chi2() const;
  float phi() const;
  float theta() const;
  float eta() const;
  int   charge() const;
 
  float energy(double mass = 0.) const;
  float xf(double mass = 0.) const;
  float rap(double mass = 0.) const;
 
  SVector<double,3> xvec() const;
  SVector<double,3> tvec() const;
  SVector<double,3> evec() const;
  SVector<double,3> pvec() const;
 
  void collect(vector<Track*>& c) const;
 
  bool isValid() const;
  void valid();
  void invalid();
  const SMatrix<double,5>& COV()  const;
  const SMatrix<double,5>& CINV() const;
  bool propagate(const double z);
  unsigned short int ndf()  const;
  float cov_x(double dz=0.)        const;
  float cov_y(double dz=0.)        const;
  float cov_tx()                   const;
  float cov_ty()                   const;
  float cov_p()                    const;
 
  bool operator==(const Track& rhs) const;
 
  const SMatrix<double,3>& KCOV()  const;
  const SMatrix<double,3>& KCINV() const;
  const SVector<double,3>& xv()    const;
  const SMatrix<double,3>& F()     const;
  const SMatrix<double,3>& ES()    const;
  const SMatrix<double,3>& DS()    const;
 
  bool filter(const unsigned int I);
  bool smooth();
  bool smoothC();
  std::ostream& print( std::ostream& ) const;
  
private:
  const Track*      track_;
  const SVertex<NTR>* vtx_;
  double k_tx;
  double k_ty;
  double k_chi2;
  SVector<double,3>   k_qvs;
  SVector<double,3>   k_xv;
  SVector<double,3>   k_qv;
  SMatrix<double,3>   k_C;
  SMatrix<double,3>   k_CINV;
  SVector<double,5>   k_pc;
  SMatrix<double,3,5> k_WBG;
  SMatrix<double,3>   k_W;
  SMatrix<double,3>   k_F;
  SMatrix<double,3>   k_ES;
  SMatrix<double,3>   k_DS;
};
 
 
//==============================================================================
// operator <<
//==============================================================================
template <unsigned int NTR>
inline std::ostream& operator<< (std::ostream& os, const SKalman<NTR>& vt) {
  return vt.print( os );
}
 
#include "SKalman.icc"
#endif