ReferencePosition.cc

Go to the documentation of this file.

// -*- lsst-c++ -*-

/* 
 * 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/>.
 */

#ifndef LSST_AP_MATCH_REFERENCEPOSITION_CC
#define LSST_AP_MATCH_REFERENCEPOSITION_CC

#include "ReferencePosition.h"


namespace lsst { namespace ap { namespace match {

ReferencePosition::ReferencePosition(
    int64_t id,   
    lsst::afw::geom::Angle const ra,    
    lsst::afw::geom::Angle const decl,  
    double epoch  
   ) :
    _sc(ra, decl),
    _id(id),
    _epoch(epoch),
    _p(_sc.getVector().asEigen()),
    _v(Eigen::Vector3d::Zero()),
    _parallax(0.0),
    _minDecl(decl),
    _maxDecl(decl),
    _minRa(ra),
    _maxRa(ra),
    _flags(0)
{ }

inline int64_t ReferencePosition::getId() const {
    return _id;
}

inline double ReferencePosition::getEpoch() const {
    return _epoch;
}

inline int ReferencePosition::getFlags() const {
   return _flags;
}

inline lsst::afw::coord::IcrsCoord const & ReferencePosition::getSphericalCoords() const {
    return _sc;
}

inline Eigen::Vector3d const & ReferencePosition::getPosition() const {
    return _p; 
}

inline Eigen::Vector3d const & ReferencePosition::getVelocity() const {
    return _v;
}

inline Eigen::Vector3d const ReferencePosition::getPosition(double epoch) const {
    if ((_flags & MOVING) == 0) {
        return _p;
    }
    Eigen::Vector3d p = _p + _v*(epoch - _epoch);
    if ((_flags & PARALLAX_COR) != 0) {
        p -= lsst::ap::utils::earthPosition(epoch);
    }
    return p.normalized();
}

inline Eigen::Vector3d const ReferencePosition::getPosition(
    double epoch,
    Eigen::Vector3d const &origin // barycentric coordinate of origin, AU
) const {
    if ((_flags & MOVING) == 0) {
        return _p;
    }
    Eigen::Vector3d p = _p + _v*(epoch - _epoch);
    if ((_flags & PARALLAX_COR) != 0) {
        p -= origin;
    }
    return p.normalized();
}

}}} // namespace lsst::ap::match

#endif // LSST_AP_MATCH_REFERENCEPOSITION_CC