ShapeletPsfCandidate.h

Go to the documentation of this file.

// -*- LSST-C++ -*-
#ifndef LSST_MEAS_ALGORITHMS_SHAPELETPSFCANDIDATE_H
#define LSST_MEAS_ALGORITHMS_SHAPELETPSFCANDIDATE_H
/* 
 * LSST Data Management System
 * Copyright 2008, 2009, 2010 LSST Corporation.
 * 
 * This product includes software developed by the
 * LSST Project (http://www.lsst.org/).
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * 
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the LSST License Statement and 
 * the GNU General Public License along with this program.  If not, 
 * see <http://www.lsstcorp.org/LegalNotices/>.
 */

#include "lsst/afw/math/SpatialCell.h"
#include "lsst/afw/table/Source.h"
#include "boost/shared_ptr.hpp"

namespace lsst {
namespace meas {
namespace algorithms {

    class ShapeletPsfCandidate : 
        public lsst::afw::math::SpatialCellCandidate
    {
    public :
        typedef lsst::afw::math::SpatialCellCandidate base;
        typedef lsst::afw::table::SourceRecord Source;

        typedef boost::shared_ptr<ShapeletPsfCandidate> Ptr;
        typedef boost::shared_ptr<const ShapeletPsfCandidate> ConstPtr;

        inline ShapeletPsfCandidate(
            double x,       
            double y,       
            double size,    
            PTR(Source) source   
        ) : 
            base(x,y), _size(size), _source(source), _rating(1.)
        {}

        inline void setShapelet(Shapelet::ConstPtr shapelet)
        { 
            _shapelet = shapelet; 
            if (_shapelet->hasCovariance())
                _rating = _shapelet->getValues()(0) / sqrt((*_shapelet->getCovariance())(0,0));
        }

        inline double getX() const { return base::getXCenter(); }
        inline double getY() const { return base::getYCenter(); }

        inline double getSize() const { return _size; }

        inline PTR(Source) getSource() const 
        { 
            assert(_source);
            return _source; 
        }

        inline Shapelet::ConstPtr getShapelet() const 
        { 
            assert(_shapelet);
            return _shapelet; 
        }

        inline bool hasShapelet() const { return static_cast<bool>(_shapelet); }

        inline double getCandidateRating() const { return _rating; }

        inline void setBad() 
        { 
            lsst::afw::math::SpatialCellCandidate::setStatus(
                lsst::afw::math::SpatialCellCandidate::BAD);
        }

    private :

        double _size;
        PTR(Source) _source;
        Shapelet::ConstPtr _shapelet;
        double _rating;
    };

}}}

#endif