lsst.jointcal 22.0.1-26-g6000b2f+22c1dfde1c
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
34#include "lsst/jointcal/Chi2.h"
39
40namespace lsst {
41namespace jointcal {
42
78class AstrometryFit : public FitterBase {
79public:
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
129protected:
131 void saveChi2MeasContributions(std::string const &filename) const override;
132
134 void saveChi2RefContributions(std::string const &filename) const override;
135
136private:
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.
std::shared_ptr< AstrometryModel > getModel() const
Return the model being fit.
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.
AstrometryFit & operator=(AstrometryFit &&)=delete
AstrometryFit & operator=(AstrometryFit const &)=delete
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(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(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
Class for a simple mapping implementing a generic AstrometryTransform.