lsst.jointcal  22.0.1-11-g977d919+31fba84f5c
ProperMotion.h
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 #ifndef LSST_JOINTCAL_PROPER_MOTION_H
26 #define LSST_JOINTCAL_PROPER_MOTION_H
27 
28 #include <iostream>
29 #include <memory>
30 #include <math.h> // atan2
31 
32 namespace lsst {
33 namespace jointcal {
34 
35 class Point;
36 
45 class ProperMotion {
46 public:
47  ProperMotion(double ra, double dec, double raErr, double decErr, double raDecCov = 0)
48  : _ra(ra), _dec(dec), _raErr(raErr), _decErr(decErr), _raDecCov(raDecCov) {
49  _offsetBearing = atan2(dec, ra);
50  };
51 
52  ProperMotion(ProperMotion const &) = default;
53  ProperMotion(ProperMotion &&) = default;
54  ProperMotion &operator=(ProperMotion const &) = default;
56  ~ProperMotion() = default;
57 
68  Point apply(Point star, double timeDeltaYears) const;
69 
70  friend std::ostream &operator<<(std::ostream &stream, ProperMotion const &pm);
71 
72 private:
73  // Proper motion in ra and dec. Units are radians/year
74  double _ra, _dec;
75  double _raErr, _decErr, _raDecCov;
76 
77  // Cache the bearing along-which to offset.
78  double _offsetBearing;
79 };
80 } // namespace jointcal
81 } // namespace lsst
82 
83 #endif // LSST_JOINTCAL_PROPER_MOTION_H
double dec
T atan2(T... args)
A point in a plane.
Definition: Point.h:37
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
ProperMotion(double ra, double dec, double raErr, double decErr, double raDecCov=0)
Definition: ProperMotion.h:47
friend std::ostream & operator<<(std::ostream &stream, ProperMotion const &pm)
Definition: ProperMotion.cc:44
ProperMotion & operator=(ProperMotion &&)=default
ProperMotion & operator=(ProperMotion const &)=default
ProperMotion(ProperMotion &&)=default
ProperMotion(ProperMotion const &)=default
Class for a simple mapping implementing a generic AstrometryTransform.