LSST Applications g013ef56533+7c9321ec0f,g042eb84c57+c6cfa41bc3,g199a45376c+0ba108daf9,g1fd858c14a+fcad0d0313,g210f2d0738+c0f94c6586,g262e1987ae+a7e710680e,g29ae962dfc+fb55f2edb0,g2ac17093b6+61d6563b1e,g2b1d02342f+df6f932764,g2cef7863aa+aef1011c0b,g2f7ad74990+c0f94c6586,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+53cf87ae69,g47891489e3+4316d04fff,g511e8cfd20+baa56acf6c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+fd7ad03fde,g64539dfbff+c0f94c6586,g67b6fd64d1+4316d04fff,g67fd3c3899+c0f94c6586,g6985122a63+4316d04fff,g74acd417e5+ca833bee28,g786e29fd12+668abc6043,g81db2e9a8d+b2ec8e584f,g87389fa792+8856018cbb,g89139ef638+4316d04fff,g8d7436a09f+0a24083b20,g8ea07a8fe4+760ca7c3fc,g90f42f885a+033b1d468d,g97be763408+11eb8fd5b8,gbf99507273+8c5ae1fdc5,gcdda8b9158+e4c84c9d5c,gce8aa8abaa+8c5ae1fdc5,gd7ef33dd92+4316d04fff,gdab6d2f7ff+ca833bee28,ge410e46f29+4316d04fff,geaed405ab2+c4bbc419c6,gf9a733ac38+8c5ae1fdc5,w.2025.40
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)