lsst.jointcal  22.0.1-25-g5bf6245+c12e52a65d
ProperMotion.cc
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 /*
3  * This file is part of jointcal.
4  *
5  * Developed for the LSST Data Management System.
6  * This product includes software developed by the LSST Project
7  * (https://www.lsst.org).
8  * See the COPYRIGHT file at the top-level directory of this distribution
9  * for details of code ownership.
10  *
11  * This program is free software: you can redistribute it and/or modify
12  * it under the terms of the GNU General Public License as published by
13  * the Free Software Foundation, either version 3 of the License, or
14  * (at your option) any later version.
15  *
16  * This program is distributed in the hope that it will be useful,
17  * but WITHOUT ANY WARRANTY; without even the implied warranty of
18  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19  * GNU General Public License for more details.
20  *
21  * You should have received a copy of the GNU General Public License
22  * along with this program. If not, see <https://www.gnu.org/licenses/>.
23  */
24 
25 #include <cmath>
26 #include <iostream>
27 #include "lsst/geom/SpherePoint.h"
29 #include "lsst/jointcal/Point.h"
30 
31 namespace lsst {
32 namespace jointcal {
33 
34 Point ProperMotion::apply(Point star, double timeDeltaYears) const {
35  geom::SpherePoint spherePoint(star.x, star.y, geom::degrees);
36  double amount = std::hypot(_ra * timeDeltaYears, _dec * timeDeltaYears);
37  auto result = spherePoint.offset(_offsetBearing * geom::radians, amount * geom::radians);
38  Point newStar(star);
39  newStar.x = result.getRa().asDegrees();
40  newStar.y = result.getDec().asDegrees();
41  return newStar;
42 }
43 
45  stream << "pm_ra*cos(dec)=" << pm._ra << "rad/yr, pm_dec=" << pm._dec << "rad/yr, pm_raErr=" << pm._raErr
46  << "rad/yr, pm_decErr=" << pm._decErr << "rad/yr, pm_raDecCov=" << pm._raDecCov;
47  return stream;
48 }
49 
50 } // namespace jointcal
51 } // namespace lsst
SpherePoint offset(Angle const &bearing, Angle const &amount) const
A point in a plane.
Definition: Point.h:37
double x
coordinate
Definition: Point.h:42
Proper motion data for a reference star or fitted star.
Definition: ProperMotion.h:45
Point apply(Point star, double timeDeltaYears) const
Apply proper motion correction to the input star, returning a star with PM-corrected coordinates.
Definition: ProperMotion.cc:34
T hypot(T... args)
constexpr AngleUnit degrees
constexpr AngleUnit radians
std::ostream & operator<<(std::ostream &stream, AstrometryMapping const &mapping)
Class for a simple mapping implementing a generic AstrometryTransform.