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
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::class_<PhotometryMappingBase, std::shared_ptr<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::class_<PhotometryMapping, std::shared_ptr<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::class_<ChipVisitPhotometryMapping, std::shared_ptr<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::class_<ChipVisitFluxMapping, std::shared_ptr<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::class_<ChipVisitMagnitudeMapping, std::shared_ptr<ChipVisitMagnitudeMapping>,
106 ChipVisitPhotometryMapping>;
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
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
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)