LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
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.
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.
void leastSquareDerivativesMeasurement(CcdImage const &ccdImage, TripletList &tripletList, Eigen::VectorXd &grad, MeasuredStarList const *msList=nullptr) const override
Compute the derivatives of the measured stars and model for one CcdImage.
void accumulateStatImageList(CcdImageList const &ccdImageList, Chi2Accumulator &accum) const override
Compute the chi2 (per star or total, depending on which Chi2Accumulator is used) for measurements.
void accumulateStatRefStars(Chi2Accumulator &accum) const override
Compute the chi2 (per star or total, depending on which Chi2Accumulator is used) for RefStars.
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.
void leastSquareDerivativesReference(FittedStarList const &fittedStarList, TripletList &tripletList, Eigen::VectorXd &grad) const override
Compute the derivatives of the reference terms.
void getIndicesOfMeasuredStar(MeasuredStar const &measuredStar, IndexVector &indices) const override
this routine is to be used only in the framework of outlier removal
AstrometryFit(AstrometryFit &&)=delete
void saveChi2MeasContributions(std::string const &filename) const override
Save a CSV file containing residuals of measurement terms.
Base class for fitters.
Definition FitterBase.h:55
std::list< std::shared_ptr< CcdImage > > CcdImageList
Definition CcdImage.h:46