KernelSumVisitor.h

Go to the documentation of this file.

// -*- lsst-c++ -*-
#ifndef LSST_IP_DIFFIM_KERNELSUMVISITOR_H
#define LSST_IP_DIFFIM_KERNELSUMVISITOR_H

#include "boost/shared_ptr.hpp"
#include "lsst/afw/math.h"
#include "lsst/afw/image.h"
#include "lsst/pex/policy/Policy.h"

namespace lsst { 
namespace ip { 
namespace diffim { 
namespace detail {

    template<typename PixelT>
    class KernelSumVisitor : public lsst::afw::math::CandidateVisitor {
    public:
        typedef boost::shared_ptr<KernelSumVisitor<PixelT> > Ptr;

        enum Mode {AGGREGATE = 0, REJECT = 1};
        
        KernelSumVisitor(lsst::pex::policy::Policy const& policy);
        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::pex::policy::Policy _policy;   
    };    
    
    template<typename PixelT>
    boost::shared_ptr<KernelSumVisitor<PixelT> >
    makeKernelSumVisitor(lsst::pex::policy::Policy const& policy) {
        return typename KernelSumVisitor<PixelT>::Ptr(new KernelSumVisitor<PixelT>(policy));
    }

}}}} // end of namespace lsst::ip::diffim::detail

#endif