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
ChipVisitAstrometryMapping.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_CHIP_VISIT_ASTROMETRY_MAPPING_H
26 #define LSST_JOINTCAL_CHIP_VISIT_ASTROMETRY_MAPPING_H
27 
28 #include "memory"
29 
33 
34 namespace lsst {
35 namespace jointcal {
36 
39 public:
43 
49 
50  std::size_t getNpar() const override;
51 
52  void getMappingIndices(IndexVector &indices) const override;
53 
55  void computeTransformAndDerivatives(FatPoint const &where, FatPoint &outPoint,
56  Eigen::MatrixX2d &H) const override;
58  void transformPosAndErrors(FatPoint const &where, FatPoint &outPoint) const override;
59 
66  void offsetParams(Eigen::VectorXd const &delta) override {
67  _m1->offsetParams(delta.segment(_m1->getIndex(), _m1->getNpar()));
68  _m2->offsetParams(delta.segment(_m2->getIndex() + _m1->getNpar(), _m2->getNpar()));
69  }
70 
72  AstrometryTransform const &getTransform1() const { return _m1->getTransform(); }
73 
75  AstrometryTransform const &getTransform2() const { return _m2->getTransform(); }
76 
77  void positionDerivative(Point const &where, Eigen::Matrix2d &derivative, double epsilon) const override;
78 
80  void freezeErrorTransform();
81 
82  void print(std::ostream &out) const override;
83 
84 private:
87  void setWhatToFit(const bool fittingT1, const bool fittingT2);
88 
90  Eigen::Index _nPar1, _nPar2;
91  struct tmpVars // just there to get around constness issues
92  {
93  Eigen::MatrixX2d h1, h2;
94  Eigen::Matrix2d dt2dx;
95  };
96 
98 };
99 } // namespace jointcal
100 } // namespace lsst
101 
102 #endif // LSST_JOINTCAL_CHIP_VISIT_ASTROMETRY_MAPPING_H
Eigen::Matrix< double, Eigen::Dynamic, 2 > MatrixX2d
Definition: Eigenstuff.h:33
virtual class needed in the abstraction of the distortion model
a virtual (interface) class for geometric transformations.
The mapping with two transforms in a row.
std::size_t getNpar() const override
Number of parameters in total.
void getMappingIndices(IndexVector &indices) const override
Sets how this set of parameters (of length Npar()) map into the "grand" fit Expects that indices has ...
ChipVisitAstrometryMapping & operator=(ChipVisitAstrometryMapping const &)=delete
void freezeErrorTransform()
Currently not implemented.
void offsetParams(Eigen::VectorXd const &delta) override
Remember the error scale and freeze it.
void computeTransformAndDerivatives(FatPoint const &where, FatPoint &outPoint, Eigen::MatrixX2d &H) const override
Actually applies the AstrometryMapping and evaluates the derivatives w.r.t the fitted parameters.
AstrometryTransform const & getTransform1() const
access to transforms
ChipVisitAstrometryMapping(ChipVisitAstrometryMapping &&)=delete
ChipVisitAstrometryMapping & operator=(ChipVisitAstrometryMapping &&)=delete
void positionDerivative(Point const &where, Eigen::Matrix2d &derivative, double epsilon) const override
The derivative w.r.t. position.
void print(std::ostream &out) const override
Print a string representation of the contents of this mapping, for debugging.
AstrometryTransform const & getTransform2() const
access to transforms
ChipVisitAstrometryMapping(ChipVisitAstrometryMapping const &)=delete
No copy or move: there is only ever one instance of a given model (i.e.. per ccd+visit)
ChipVisitAstrometryMapping(std::shared_ptr< SimpleAstrometryMapping > chipMapping, std::shared_ptr< SimpleAstrometryMapping > visitMapping)
void transformPosAndErrors(FatPoint const &where, FatPoint &outPoint) const override
The same as above but without the parameter derivatives (used to evaluate chi^2)
A multi-component model, fitting mappings for sensors and visits simultaneously.
A Point with uncertainties.
Definition: FatPoint.h:34
A point in a plane.
Definition: Point.h:37
A base class for image defects.