LSSTApplications  17.0+11,17.0+34,17.0+56,17.0+57,17.0+59,17.0+7,17.0-1-g377950a+33,17.0.1-1-g114240f+2,17.0.1-1-g4d4fbc4+28,17.0.1-1-g55520dc+49,17.0.1-1-g5f4ed7e+52,17.0.1-1-g6dd7d69+17,17.0.1-1-g8de6c91+11,17.0.1-1-gb9095d2+7,17.0.1-1-ge9fec5e+5,17.0.1-1-gf4e0155+55,17.0.1-1-gfc65f5f+50,17.0.1-1-gfc6fb1f+20,17.0.1-10-g87f9f3f+1,17.0.1-11-ge9de802+16,17.0.1-16-ga14f7d5c+4,17.0.1-17-gc79d625+1,17.0.1-17-gdae4c4a+8,17.0.1-2-g26618f5+29,17.0.1-2-g54f2ebc+9,17.0.1-2-gf403422+1,17.0.1-20-g2ca2f74+6,17.0.1-23-gf3eadeb7+1,17.0.1-3-g7e86b59+39,17.0.1-3-gb5ca14a,17.0.1-3-gd08d533+40,17.0.1-30-g596af8797,17.0.1-4-g59d126d+4,17.0.1-4-gc69c472+5,17.0.1-6-g5afd9b9+4,17.0.1-7-g35889ee+1,17.0.1-7-gc7c8782+18,17.0.1-9-gc4bbfb2+3,w.2019.22
LSSTDataManagementBasePackage
FittedStar.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_FITTED_STAR_H
26 #define LSST_JOINTCAL_FITTED_STAR_H
27 
28 #include <iostream>
29 #include <fstream>
30 
31 #include "Eigen/Core"
32 #include "lsst/jointcal/BaseStar.h"
33 #include "lsst/jointcal/StarList.h"
34 
35 namespace lsst {
36 namespace jointcal {
37 
38 class MeasuredStar;
39 class RefStar;
40 class AstrometryTransform;
41 
44 
46 struct PmBlock {
47  // proper motion in x and y. Units depend on how you fit them
48  double pmx, pmy;
49  double epmx, epmy, epmxy;
50  double color; // OK it is unrelated, but associated in practice
51  bool mightMove;
52 
53  PmBlock() : pmx(0), pmy(0), epmx(0), epmy(0), epmxy(0), color(0), mightMove(false){};
54 };
55 
61 class FittedStar : public BaseStar, public PmBlock {
62 public:
63  FittedStar() : BaseStar(), _indexInMatrix(-1), _measurementCount(0), _refStar(nullptr) {}
64 
65  FittedStar(const BaseStar& baseStar)
66  : BaseStar(baseStar), _indexInMatrix(0), _measurementCount(0), _refStar(nullptr) {}
67 
69  FittedStar(const MeasuredStar& measuredStar);
70 
73  FittedStar(FittedStar const&) = default;
74  FittedStar(FittedStar&&) = delete;
75  FittedStar& operator=(FittedStar const&) = delete;
76  FittedStar& operator=(FittedStar&&) = delete;
77 
80  _indexInMatrix = -1;
81  _measurementCount = 0;
82  _refStar = nullptr;
83  _flux = 0;
84  _fluxErr = 0;
85  _mag = 0;
86  _magErr = 0;
87  }
88 
90  void dump(std::ostream& stream = std::cout) const {
92  stream << " mcount: " << _measurementCount;
93  }
94 
96  int getMeasurementCount() const { return _measurementCount; }
97  int& getMeasurementCount() { return _measurementCount; }
98 
100  void addMagMeasurement(double magValue, double magWeight);
101 
103  void setIndexInMatrix(Eigen::Index const index) { _indexInMatrix = index; };
104 
106  Eigen::Index getIndexInMatrix() const { return _indexInMatrix; }
107 
109  void setRefStar(const RefStar* _refStar);
110 
112  const RefStar* getRefStar() const { return _refStar; };
113 
114 private:
115  Eigen::Index _indexInMatrix;
116  int _measurementCount;
117  const RefStar* _refStar;
118 };
119 
120 /****** FittedStarList */
121 
123 class FittedStarList : public StarList<FittedStar> {
124 public:
126 
128  FittedStarList() { inTangentPlaneCoordinates = true; }
129 };
130 
131 typedef FittedStarList::const_iterator FittedStarCIterator;
132 typedef FittedStarList::iterator FittedStarIterator;
133 
136 const BaseStarList& Fitted2Base(const FittedStarList& This);
137 const BaseStarList* Fitted2Base(const FittedStarList* This);
138 } // namespace jointcal
139 } // namespace lsst
140 
141 #endif // LSST_JOINTCAL_FITTED_STAR_H
Objects used as position anchors, typically USNO stars.
Definition: RefStar.h:39
virtual void dump(std::ostream &stream=std::cout) const
utility
Definition: BaseStar.h:82
objects whose position is going to be fitted. Coordinates in Common Tangent Plane.
Definition: FittedStar.h:46
int getMeasurementCount() const
Definition: FittedStar.h:96
void setIndexInMatrix(Eigen::Index const index)
index is a value that a fit can set and reread....
Definition: FittedStar.h:103
FittedStarList::const_iterator FittedStarCIterator
Definition: FittedStar.h:131
The base class for handling stars. Used by all matching routines.
Definition: BaseStar.h:50
Eigen::Index getIndexInMatrix() const
Definition: FittedStar.h:106
std::lists of Stars.
Definition: StarList.h:58
A base class for image defects.
FittedStarList::iterator FittedStarIterator
Definition: FittedStar.h:132
A list of FittedStar s. Such a list is typically constructed by Associations.
Definition: FittedStar.h:123
objects measured on actual images.
Definition: MeasuredStar.h:46
void dump(std::ostream &stream=std::cout) const
utility
Definition: FittedStar.h:90
BaseStarList & Fitted2Base(FittedStarList &This)
Definition: FittedStar.cc:64
const RefStar * getRefStar() const
Get the astrometric reference star associated with this star.
Definition: FittedStar.h:112
STL class.
The objects which have been measured several times.
Definition: FittedStar.h:61
FittedStar(const BaseStar &baseStar)
Definition: FittedStar.h:65