PCAxis.h

Go to the documentation of this file.

//
// \file  PCAxis.h
// \brief Definition of a Principal Component Axis object for LArSoft
//
// \author usher@slac.stanford.edu
//
// \see   PCAxis.cxx
//

#ifndef PCAxis_H
#define PCAxis_H

#include <vector>

#include <iosfwd>

namespace recob
{
    //
    // @brief PCAxis is an object containting the results of a Principal Components
    //        Analysis of a group of space points.

    class PCAxis
    {
    public:
    
        typedef std::vector<std::vector<double> > EigenVectors;
    
        PCAxis();
    
    private:
    
        bool         fSvdOK;            
        int          fNumHitsUsed;      
        double       fEigenValues[3];   
        EigenVectors fEigenVectors;     
        double       fAvePosition[3];   
        double       fAveHitDoca;       
        size_t       fID;               
    
    public:
    
        PCAxis(bool ok, int nHits, const double* eigenValues, const EigenVectors& eigenVecs, const double* avePos, const double aveHitDoca = 9999., size_t id=0);
    
        bool                getSvdOK()        const;
        int                 getNumHitsUsed()  const;
        const double*       getEigenValues()  const;
        const EigenVectors& getEigenVectors() const;
        const double*       getAvePosition()  const;
        double              getAveHitDoca()   const;
        size_t              getID()           const;
    
        friend std::ostream&  operator << (std::ostream & o, const PCAxis& a);
        friend bool operator < (const PCAxis& a, const PCAxis& b);
    
    };
    
}


inline bool                               recob::PCAxis::getSvdOK()        const {return fSvdOK;}
inline int                                recob::PCAxis::getNumHitsUsed()  const {return fNumHitsUsed;}
inline const double*                      recob::PCAxis::getEigenValues()  const {return fEigenValues;}
inline const recob::PCAxis::EigenVectors& recob::PCAxis::getEigenVectors() const {return fEigenVectors;}
inline const double*                      recob::PCAxis::getAvePosition()  const {return fAvePosition;}
inline double                             recob::PCAxis::getAveHitDoca()   const {return fAveHitDoca;}
inline size_t                             recob::PCAxis::getID()           const {return fID;}



#endif // PCAxis_H