LSSTApplications  17.0+124,17.0+14,17.0+73,18.0.0+37,18.0.0+80,18.0.0-4-g68ffd23+4,18.1.0-1-g0001055+12,18.1.0-1-g03d53ef+5,18.1.0-1-g1349e88+55,18.1.0-1-g2505f39+44,18.1.0-1-g5315e5e+4,18.1.0-1-g5e4b7ea+14,18.1.0-1-g7e8fceb+4,18.1.0-1-g85f8cd4+48,18.1.0-1-g8ff0b9f+4,18.1.0-1-ga2c679d+1,18.1.0-1-gd55f500+35,18.1.0-10-gb58edde+2,18.1.0-11-g0997b02+4,18.1.0-13-gfe4edf0b+12,18.1.0-14-g259bd21+21,18.1.0-19-gdb69f3f+2,18.1.0-2-g5f9922c+24,18.1.0-2-gd3b74e5+11,18.1.0-2-gfbf3545+32,18.1.0-26-g728bddb4+5,18.1.0-27-g6ff7ca9+2,18.1.0-3-g52aa583+25,18.1.0-3-g8ea57af+9,18.1.0-3-gb69f684+42,18.1.0-3-gfcaddf3+6,18.1.0-32-gd8786685a,18.1.0-4-gf3f9b77+6,18.1.0-5-g1dd662b+2,18.1.0-5-g6dbcb01+41,18.1.0-6-gae77429+3,18.1.0-7-g9d75d83+9,18.1.0-7-gae09a6d+30,18.1.0-9-gc381ef5+4,w.2019.45
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