KernelSumVisitor.h

Go to the documentation of this file.

// -*- lsst-c++ -*-
#ifndef LSST_IP_DIFFIM_KERNELSUMVISITOR_H
#define LSST_IP_DIFFIM_KERNELSUMVISITOR_H
 
#include <memory>
 
#include "lsst/afw/math.h"
#include "lsst/afw/image.h"
#include "lsst/daf/base/PropertySet.h"
 
namespace lsst {
namespace ip {
namespace diffim {
namespace detail {
 
    template<typename PixelT>
    class KernelSumVisitor : public lsst::afw::math::CandidateVisitor {
    public:
        typedef std::shared_ptr<KernelSumVisitor<PixelT> > Ptr;
 
        enum Mode {AGGREGATE = 0, REJECT = 1};
 
        KernelSumVisitor(lsst::daf::base::PropertySet const& ps);
        virtual ~KernelSumVisitor() {};
 
        void setMode(Mode mode) {_mode = mode;}
        int    getNRejected() {return _nRejected;}
        double getkSumMean()  {return _kSumMean;}
        double getkSumStd()   {return _kSumStd;}
        double getdkSumMax()  {return _dkSumMax;}
        int    getkSumNpts()  {return _kSumNpts;}
 
        void resetKernelSum();
        void processCandidate(lsst::afw::math::SpatialCellCandidate *candidate);
        void processKsumDistribution();
 
    private:
        Mode _mode;                  
        std::vector<double> _kSums;  
        double _kSumMean;            
        double _kSumStd;             
        double _dkSumMax;            
        int    _kSumNpts;            
        int    _nRejected;           
        lsst::daf::base::PropertySet::Ptr _ps;   
    };
 
    template<typename PixelT>
    std::shared_ptr<KernelSumVisitor<PixelT> >
    makeKernelSumVisitor(lsst::daf::base::PropertySet const& ps) {
        return std::shared_ptr<KernelSumVisitor<PixelT>>(new KernelSumVisitor<PixelT>(ps));
    }
 
}}}} // end of namespace lsst::ip::diffim::detail
 
#endif