KHitMulti.h

Go to the documentation of this file.

#ifndef KHITMULTI_H
#define KHITMULTI_H

#include "lardata/RecoObjects/KHit.h"

namespace trkf {

  class KHitMulti : public KHitBase
  {
  public:

    KHitMulti();

    KHitMulti(const std::shared_ptr<const Surface>& psurf);

    virtual ~KHitMulti();

    // Accessors.

    int getMeasDim() const {return fMeasDim;}

    const std::vector<std::shared_ptr<const KHit<1> > >& getMeasVec() const {return fMeasVec;}

    const ublas::vector<double>& getMeasVector() const {return fMvec;}

    const ublas::symmetric_matrix<double>& getMeasError() const {return fMerr;}

    const ublas::vector<double>& getPredVector() const {return fPvec;}

    const ublas::symmetric_matrix<double>& getPredError() const {return fPerr;}

    const ublas::vector<double>& getResVector() const {return fRvec;}

    const ublas::symmetric_matrix<double>& getResError() const {return fRerr;}

    const ublas::symmetric_matrix<double>& getResInvError() const {return fRinv;}

    const ublas::matrix<double>& getH() const {return fH;}

    double getChisq() const {return fChisq;}

    // Modifiers.

    void addMeas(const std::shared_ptr<const KHitBase>& pmeas);

    void addMeas(const std::shared_ptr<const KHit<1> >& pmeas);

    // Overrides.

    bool predict(const KETrack& tre, const Propagator* prop = 0,
                 const KTrack* ref = 0) const;

    void update(KETrack& tre) const;

    virtual std::ostream& Print(std::ostream& out, bool doTitle = true) const;

  private:

    // Attributes.

    int fMeasDim;                                   
    std::vector<std::shared_ptr<const KHit<1> > > fMeasVec;   
    mutable ublas::vector<double> fMvec;                    
    mutable ublas::symmetric_matrix<double> fMerr;          
    mutable ublas::vector<double> fPvec;            
    mutable ublas::symmetric_matrix<double> fPerr;  
    mutable ublas::vector<double> fRvec;            
    mutable ublas::symmetric_matrix<double> fRerr;  
    mutable ublas::symmetric_matrix<double> fRinv;  
    mutable ublas::matrix<double> fH;               
    mutable double fChisq;                          
  };
}

#endif