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, 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