LSST Applications  21.0.0-172-gfb10e10a+18fedfabac,22.0.0+297cba6710,22.0.0+80564b0ff1,22.0.0+8d77f4f51a,22.0.0+a28f4c53b1,22.0.0+dcf3732eb2,22.0.1-1-g7d6de66+2a20fdde0d,22.0.1-1-g8e32f31+297cba6710,22.0.1-1-geca5380+7fa3b7d9b6,22.0.1-12-g44dc1dc+2a20fdde0d,22.0.1-15-g6a90155+515f58c32b,22.0.1-16-g9282f48+790f5f2caa,22.0.1-2-g92698f7+dcf3732eb2,22.0.1-2-ga9b0f51+7fa3b7d9b6,22.0.1-2-gd1925c9+bf4f0e694f,22.0.1-24-g1ad7a390+a9625a72a8,22.0.1-25-g5bf6245+3ad8ecd50b,22.0.1-25-gb120d7b+8b5510f75f,22.0.1-27-g97737f7+2a20fdde0d,22.0.1-32-gf62ce7b1+aa4237961e,22.0.1-4-g0b3f228+2a20fdde0d,22.0.1-4-g243d05b+871c1b8305,22.0.1-4-g3a563be+32dcf1063f,22.0.1-4-g44f2e3d+9e4ab0f4fa,22.0.1-42-gca6935d93+ba5e5ca3eb,22.0.1-5-g15c806e+85460ae5f3,22.0.1-5-g58711c4+611d128589,22.0.1-5-g75bb458+99c117b92f,22.0.1-6-g1c63a23+7fa3b7d9b6,22.0.1-6-g50866e6+84ff5a128b,22.0.1-6-g8d3140d+720564cf76,22.0.1-6-gd805d02+cc5644f571,22.0.1-8-ge5750ce+85460ae5f3,master-g6e05de7fdc+babf819c66,master-g99da0e417a+8d77f4f51a,w.2021.48
LSST Data Management Base Package
AstrometryFit.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_ASTROMETRY_FIT_H
26 #define LSST_JOINTCAL_ASTROMETRY_FIT_H
27 
28 #include <string>
29 #include <iostream>
30 #include <sstream>
31 
33 #include "lsst/jointcal/CcdImage.h"
34 #include "lsst/jointcal/Chi2.h"
39 
40 namespace lsst {
41 namespace jointcal {
42 
78 class AstrometryFit : public FitterBase {
79 public:
82  std::shared_ptr<AstrometryModel> astrometryModel, double posError);
83 
85  AstrometryFit(AstrometryFit const &) = delete;
89 
108  void assignIndices(std::string const &whatToFit) override;
109 
117  void freezeErrorTransform() { _astrometryModel->freezeErrorTransform(); }
118 
119  void offsetParams(Eigen::VectorXd const &delta) override;
120 
122  std::shared_ptr<AstrometryModel> getModel() const { return _astrometryModel; }
123 
127  void checkStuff();
128 
129 protected:
131  void saveChi2MeasContributions(std::string const &filename) const override;
132 
134  void saveChi2RefContributions(std::string const &filename) const override;
135 
136 private:
137  bool _fittingDistortions, _fittingPos, _fittingPM;
138  std::shared_ptr<AstrometryModel> _astrometryModel;
139 
140  double _epoch; // epoch to correct proper motion/parallax to (Julian Epoch year, e.g. J2000.0)
141  double _posError; // constant term on error on position (in pixel unit)
142 
143  void leastSquareDerivativesMeasurement(CcdImage const &ccdImage, TripletList &tripletList,
144  Eigen::VectorXd &grad,
145  MeasuredStarList const *msList = nullptr) const override;
146 
147  void leastSquareDerivativesReference(FittedStarList const &fittedStarList, TripletList &tripletList,
148  Eigen::VectorXd &grad) const override;
149 
150  void accumulateStatImageList(CcdImageList const &ccdImageList, Chi2Accumulator &accum) const override;
151 
152  void accumulateStatRefStars(Chi2Accumulator &accum) const override;
153 
154  void getIndicesOfMeasuredStar(MeasuredStar const &measuredStar, IndexVector &indices) const override;
155 
164  Point transformFittedStar(FittedStar const &fittedStar, AstrometryTransform const &sky2TP,
165  double deltaYears) const;
166 
168  void accumulateStatImage(CcdImage const &ccdImage, Chi2Accumulator &accum) const;
169 };
170 } // namespace jointcal
171 } // namespace lsst
172 #endif // LSST_JOINTCAL_ASTROMETRY_FIT_H
Class that handles the astrometric least squares problem.
Definition: AstrometryFit.h:78
void checkStuff()
DEBUGGING routine.
AstrometryFit(AstrometryFit const &)=delete
No copy or move: there is only ever one fitter of a given type.
void freezeErrorTransform()
The transformations used to propagate errors are freezed to the current state.
std::shared_ptr< AstrometryModel > getModel() const
Return the model being fit.
void saveChi2RefContributions(std::string const &filename) const override
Save a CSV file containing residuals of reference terms.
void offsetParams(Eigen::VectorXd const &delta) override
Offset the parameters by the requested quantities.
AstrometryFit & operator=(AstrometryFit const &)=delete
AstrometryFit(std::shared_ptr< Associations > associations, std::shared_ptr< AstrometryModel > astrometryModel, double posError)
this is the only constructor
void assignIndices(std::string const &whatToFit) override
Set parameters to fit and assign indices in the big matrix.
AstrometryFit & operator=(AstrometryFit &&)=delete
AstrometryFit(AstrometryFit &&)=delete
void saveChi2MeasContributions(std::string const &filename) const override
Save a CSV file containing residuals of measurement terms.
a virtual (interface) class for geometric transformations.
Handler of an actual image from a single CCD.
Definition: CcdImage.h:64
Base class for Chi2Statistic and Chi2List, to allow addEntry inside Fitter for either class.
Definition: Chi2.h:44
FittedStars are objects whose position or flux is going to be fitted, and which come from the associa...
Definition: FittedStar.h:50
A list of FittedStar s. Such a list is typically constructed by Associations.
Definition: FittedStar.h:116
Base class for fitters.
Definition: FitterBase.h:53
Sources measured on images.
Definition: MeasuredStar.h:51
A list of MeasuredStar. They are usually filled in Associations::createCcdImage.
Definition: MeasuredStar.h:151
A point in a plane.
Definition: Point.h:37
A base class for image defects.