LSST Applications g013ef56533+b8d55c8942,g083dd6704c+a047e97985,g199a45376c+0ba108daf9,g1fd858c14a+7a3b874d60,g210f2d0738+7416ca6900,g262e1987ae+1d557ba9a3,g29ae962dfc+519d34895e,g2cef7863aa+aef1011c0b,g30d7c61c20+36d16ea71a,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+cb326ad149,g47891489e3+f459a6810c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+890c8e1e5d,g5a60e81ecd+6240c63dbc,g64539dfbff+7416ca6900,g67b6fd64d1+f459a6810c,g6ebf1fc0d4+8c5ae1fdc5,g74acd417e5+0bae3c876a,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+f459a6810c,g8d7436a09f+dee7680868,g8ea07a8fe4+81eaaadc04,g90f42f885a+34c0557caf,g97be763408+14b8164b5b,g98a1a72a9c+8389601a76,g98df359435+fff771c62d,gb8cb2b794d+6728931916,gbf99507273+8c5ae1fdc5,gc2a301910b+7416ca6900,gca7fc764a6+f459a6810c,gd7ef33dd92+f459a6810c,gdab6d2f7ff+0bae3c876a,ge410e46f29+f459a6810c,ge41e95a9f2+7416ca6900,geaed405ab2+e3b4b2a692,gf9a733ac38+8c5ae1fdc5,w.2025.43
LSST Data Management Base Package
Loading...
Searching...
No Matches
photometryMappings.cc
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#include "pybind11/pybind11.h"
27#include "pybind11/eigen.h"
28#include "pybind11/stl.h"
29#include "Eigen/Core"
30
33
34namespace py = pybind11;
35using namespace pybind11::literals;
36
37namespace lsst {
38namespace jointcal {
39namespace {
40
41void declarePhotometryMappingBase(lsst::cpputils::python::WrapperCollection &wrappers) {
42 using PyPhotometryMappingBase = py::classh<PhotometryMappingBase>;
43
44 wrappers.wrapType(PyPhotometryMappingBase(wrappers.module, "PhotometryMappingBase"), [](auto &mod, auto &cls) {
45 cls.def("getNpar", &PhotometryMappingBase::getNpar);
46 cls.def("transform", &PhotometryMappingBase::transform);
47 cls.def("transformError", &PhotometryMappingBase::transformError);
48 cls.def("computeParameterDerivatives",
49 [](PhotometryMappingBase const &self, MeasuredStar const &star, double instFlux) {
50 Eigen::VectorXd derivatives(self.getNpar());
51 self.computeParameterDerivatives(star, instFlux, derivatives);
52 return derivatives;
53 });
54
55 cls.def("getParameters", &PhotometryMappingBase::getParameters);
56
57 cls.def("getMappingIndices", [](PhotometryMappingBase const &self) {
58 IndexVector indices(0);
59 self.getMappingIndices(indices);
60 return indices;
61 });
62
63 cls.def("getIndex", &PhotometryMappingBase::getIndex);
64 cls.def("setIndex", &PhotometryMappingBase::setIndex);
65 });
66}
67
68void declarePhotometryMapping(lsst::cpputils::python::WrapperCollection &wrappers) {
69 using PyPhotometryMapping = py::classh<PhotometryMapping, PhotometryMappingBase>;
70
71 wrappers.wrapType(PyPhotometryMapping(wrappers.module, "PhotometryMapping"), [](auto &mod, auto &cls) {
72 cls.def(py::init<std::shared_ptr<PhotometryTransform>>(), "transform"_a);
73
74 cls.def("offsetParams", &PhotometryMapping::offsetParams);
75 cls.def("getTransform", &PhotometryMapping::getTransform);
76 cls.def("getTransformErrors", &PhotometryMapping::getTransformErrors);
77 });
78}
79
80void declareChipVisitPhotometryMapping(lsst::cpputils::python::WrapperCollection &wrappers) {
81 using PyChipVisitPhotometryMapping =
82 py::classh<ChipVisitPhotometryMapping, PhotometryMappingBase>;
83
84 wrappers.wrapType(
85 PyChipVisitPhotometryMapping(wrappers.module, "ChipVisitPhotometryMapping"), [](auto &mod, auto &cls) {
86 cls.def("setWhatToFit", &ChipVisitPhotometryMapping::setWhatToFit);
87 cls.def("getChipMapping", &ChipVisitPhotometryMapping::getChipMapping);
88 cls.def("getVisitMapping", &ChipVisitPhotometryMapping::getVisitMapping);
89 cls.def("getNParChip", &ChipVisitPhotometryMapping::getNParChip);
90 cls.def("getNParVisit", &ChipVisitPhotometryMapping::getNParVisit);
91 });
92}
93
94void declareChipVisitFluxMapping(lsst::cpputils::python::WrapperCollection &wrappers) {
95 using PyChipVisitFluxMapping =
96 py::classh<ChipVisitFluxMapping, ChipVisitPhotometryMapping>;
97
98 wrappers.wrapType(PyChipVisitFluxMapping(wrappers.module, "ChipVisitFluxMapping"), [](auto &mod, auto &cls) {
99 cls.def(py::init<std::shared_ptr<PhotometryMapping>, std::shared_ptr<PhotometryMapping>>(),
100 "chipMapping"_a, "visitMapping"_a);
101 });
102}
103
104void declareChipVisitMagnitudeMapping(lsst::cpputils::python::WrapperCollection &wrappers) {
105 using PyChipVisitMagnitudeMapping = py::classh<ChipVisitMagnitudeMapping,
107
108 wrappers.wrapType(PyChipVisitMagnitudeMapping(wrappers.module, "ChipVisitMagnitudeMapping"), [](auto &mod, auto &cls) {
109 cls.def(py::init<std::shared_ptr<PhotometryMapping>, std::shared_ptr<PhotometryMapping>>(),
110 "chipMapping"_a, "visitMapping"_a);
111 });
112}
113} // namespace
114
116 declarePhotometryMappingBase(wrappers);
117 declarePhotometryMapping(wrappers);
118 declareChipVisitPhotometryMapping(wrappers);
119 declareChipVisitFluxMapping(wrappers);
120 declareChipVisitMagnitudeMapping(wrappers);
121}
122
123} // namespace jointcal
124} // namespace lsst
std::vector< std::ptrdiff_t > IndexVector
Definition Eigenstuff.h:38
A helper class for subdividing pybind11 module across multiple translation units (i....
Definition python.h:242
PyType wrapType(PyType cls, ClassWrapperCallback function, bool setModuleName=true)
Add a type (class or enum) wrapper, deferring method and other attribute definitions until finish() i...
Definition python.h:391
pybind11::module module
The module object passed to the PYBIND11_MODULE block that contains this WrapperCollection.
Definition python.h:448
A two-level photometric transform: one for the ccd and one for the visit.
Relates transform(s) to their position in the fitting matrix and allows interaction with the transfor...
virtual Eigen::VectorXd getParameters()=0
void setIndex(Eigen::Index i)
Set the index of this mapping in the grand fit.
Eigen::Index getIndex() const
Get the index of this mapping in the grand fit.
void wrapPhotometryMappings(WrapperCollection &wrappers)