CentroidUtilities.h

Go to the documentation of this file.

// -*- lsst-c++ -*-
/*
 * LSST Data Management System
 * Copyright 2008-2014 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/>.
 */

#ifndef LSST_MEAS_BASE_CentroidUtilities_h_INCLUDED
#define LSST_MEAS_BASE_CentroidUtilities_h_INCLUDED

#include "lsst/meas/base/constants.h"
#include "lsst/afw/table/aggregates.h"

namespace lsst { namespace meas { namespace base {

struct CentroidResult {
    CentroidElement x; 
    CentroidElement y; 
    ErrElement xSigma; 
    ErrElement ySigma; 
    ErrElement x_y_Cov; 

    CentroidResult();

    Centroid const getCentroid() const;

    void setCentroid(Centroid const & centroid);

    CentroidCov const getCentroidErr() const;

    void setCentroidErr(CentroidCov const & matrix);

};

class CentroidResultKey : public afw::table::FunctorKey<CentroidResult> {
public:

    static CentroidResultKey addFields(
        afw::table::Schema & schema,
        std::string const & name,
        std::string const & doc,
        UncertaintyEnum uncertainty
    );

    CentroidResultKey() : _centroid(), _centroidErr() {}

    CentroidResultKey(
        afw::table::PointKey<CentroidElement> const & centroid,
        afw::table::CovarianceMatrixKey<ErrElement,2> const & centroidErr
    ) :
        _centroid(centroid), _centroidErr(centroidErr)
    {}

    CentroidResultKey(afw::table::SubSchema const & s);

    virtual CentroidResult get(afw::table::BaseRecord const & record) const;

    virtual void set(afw::table::BaseRecord & record, CentroidResult const & value) const;

    bool operator==(CentroidResultKey const & other) const {
        return _centroid == other._centroid && _centroidErr == other._centroidErr;
    }
    bool operator!=(CentroidResultKey const & other) const { return !(*this == other); }

    bool isValid() const { return _centroid.isValid() && _centroidErr.isValid(); }

    afw::table::PointKey<CentroidElement> getCentroid() const { return _centroid; }

    afw::table::CovarianceMatrixKey<ErrElement,2> getCentroidErr() const { return _centroidErr; }

    afw::table::Key<CentroidElement> getX() const { return _centroid.getX(); }

    afw::table::Key<CentroidElement> getY() const { return _centroid.getY(); }

private:
    afw::table::PointKey<CentroidElement> _centroid;
    afw::table::CovarianceMatrixKey<ErrElement,2> _centroidErr;
};

}}} // lsst::meas::base

#endif // !LSST_MEAS_BASE_CentroidUtilities_h_INCLUDED