lsst.jointcal 22.0.1-26-g6000b2f+cc219f2c90
PhotometryFit.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_PHOTOMETRY_FIT_H
26#define LSST_JOINTCAL_PHOTOMETRY_FIT_H
27
28#include <string>
29#include <iostream>
30#include <sstream>
31
32#include "lsst/log/Log.h"
35#include "lsst/jointcal/Chi2.h"
40
41namespace lsst {
42namespace jointcal {
43
45class PhotometryFit : public FitterBase {
46public:
55 : FitterBase(associations),
56 _fittingModel(false),
57 _fittingFluxes(false),
58 _photometryModel(photometryModel) {
59 _log = LOG_GET("lsst.jointcal.PhotometryFit");
60 }
61
63 PhotometryFit(PhotometryFit const &) = delete;
67
81 void assignIndices(std::string const &whatToFit) override;
82
83 void offsetParams(Eigen::VectorXd const &delta) override;
84
86 std::shared_ptr<PhotometryModel> getModel() const { return _photometryModel; }
87
88protected:
90 void saveChi2MeasContributions(std::string const &filename) const override;
91
93 void saveChi2RefContributions(std::string const &filename) const override;
94
95private:
96 bool _fittingModel, _fittingFluxes;
97 std::shared_ptr<PhotometryModel> _photometryModel;
98
99 void accumulateStatImageList(CcdImageList const &ccdImageList, Chi2Accumulator &accum) const override;
100
101 void accumulateStatRefStars(Chi2Accumulator &accum) const override;
102
103 void getIndicesOfMeasuredStar(MeasuredStar const &measuredStar, IndexVector &indices) const override;
104
105 void leastSquareDerivativesMeasurement(CcdImage const &ccdImage, TripletList &tripletList,
106 Eigen::VectorXd &grad,
107 MeasuredStarList const *measuredStarList = nullptr) const override;
108
110 void leastSquareDerivativesReference(FittedStarList const &fittedStarList, TripletList &tripletList,
111 Eigen::VectorXd &grad) const override;
112};
113} // namespace jointcal
114} // namespace lsst
115#endif // LSST_JOINTCAL_PHOTOMETRY_FIT_H
#define LOG_GET(logger)
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
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
Class that handles the photometric least squares problem.
Definition: PhotometryFit.h:45
void assignIndices(std::string const &whatToFit) override
Set parameters to fit and assign indices in the big matrix.
PhotometryFit & operator=(PhotometryFit &&)=delete
PhotometryFit(PhotometryFit const &)=delete
No copy or move: there is only ever one fitter of a given type.
std::shared_ptr< PhotometryModel > getModel() const
Return the model being fit.
Definition: PhotometryFit.h:86
PhotometryFit & operator=(PhotometryFit const &)=delete
PhotometryFit(PhotometryFit &&)=delete
PhotometryFit(std::shared_ptr< Associations > associations, std::shared_ptr< PhotometryModel > photometryModel)
Construct a photometry fitter.
Definition: PhotometryFit.h:53
void saveChi2RefContributions(std::string const &filename) const override
Save a CSV file containing residuals of reference terms.
void saveChi2MeasContributions(std::string const &filename) const override
Save a CSV file containing residuals of measurement terms.
void offsetParams(Eigen::VectorXd const &delta) override
Offset the parameters by the requested quantities.
Class for a simple mapping implementing a generic AstrometryTransform.