Loading [MathJax]/extensions/tex2jax.js
LSST Applications 28.0.0,g1653933729+a8ce1bb630,g1a997c3884+a8ce1bb630,g28da252d5a+5bd70b7e6d,g2bbee38e9b+638fca75ac,g2bc492864f+638fca75ac,g3156d2b45e+07302053f8,g347aa1857d+638fca75ac,g35bb328faa+a8ce1bb630,g3a166c0a6a+638fca75ac,g3e281a1b8c+7bbb0b2507,g4005a62e65+17cd334064,g414038480c+5b5cd4fff3,g41af890bb2+4ffae9de63,g4e1a3235cc+0f1912dca3,g6249c6f860+3c3976f90c,g80478fca09+46aba80bd6,g82479be7b0+77990446f6,g858d7b2824+78ba4d1ce1,g89c8672015+f667a5183b,g9125e01d80+a8ce1bb630,ga5288a1d22+2a6264e9ca,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc22bb204ba+78ba4d1ce1,gc28159a63d+638fca75ac,gcf0d15dbbd+32ddb6096f,gd6b7c0dfd1+3e339405e9,gda3e153d99+78ba4d1ce1,gda6a2b7d83+32ddb6096f,gdaeeff99f8+1711a396fd,gdd5a9049c5+b18c39e5e3,ge2409df99d+a5e4577cdc,ge33fd446bb+78ba4d1ce1,ge79ae78c31+638fca75ac,gf0baf85859+64e8883e75,gf5289d68f6+e1b046a8d7,gfa443fc69c+91d9ed1ecf,gfda6b12a05+8419469a56
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
Calib.cc
Go to the documentation of this file.
1// -*- lsst-c++ -*-
2
3/*
4 * LSST Data Management System
5 * Copyright 2008-2016 LSST Corporation.
6 *
7 * This product includes software developed by the
8 * LSST Project (http://www.lsst.org/).
9 *
10 * This program is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 3 of the License, or
13 * (at your option) any later version.
14 *
15 * This program is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the LSST License Statement and
21 * the GNU General Public License along with this program. If not,
22 * see <http://www.lsstcorp.org/LegalNotices/>.
23 */
24
25/*
26 * Classes to support calibration (e.g. photometric zero points, exposure times)
27 */
28#include <cstdint>
29
30#include "ndarray.h"
31#include "lsst/pex/exceptions.h"
33
34namespace lsst {
35namespace afw {
36
37namespace image {
38
40template <typename T>
41ndarray::Array<T, 1> abMagFromFlux(ndarray::Array<T const, 1> const& flux) {
42 ndarray::Array<T, 1> out = ndarray::allocate(flux.getShape());
43 for (std::size_t ii = 0; ii < flux.getNumElements(); ++ii) {
44 out[ii] = abMagFromFlux(flux[ii]);
45 }
46 return out;
47}
48
50template <typename T>
51ndarray::Array<T, 1> abMagErrFromFluxErr(ndarray::Array<T const, 1> const& fluxErr,
52 ndarray::Array<T const, 1> const& flux) {
53 if (flux.getNumElements() != fluxErr.getNumElements()) {
54 throw LSST_EXCEPT(pex::exceptions::LengthError, (boost::format("Length mismatch: %d vs %d") %
55 flux.getNumElements() % fluxErr.getNumElements())
56 .str());
57 }
58 ndarray::Array<T, 1> out = ndarray::allocate(flux.getShape());
59 for (std::size_t ii = 0; ii < flux.getNumElements(); ++ii) {
60 out[ii] = abMagErrFromFluxErr(fluxErr[ii], flux[ii]);
61 }
62 return out;
63}
64
66template <typename T>
67ndarray::Array<T, 1> fluxFromABMag(ndarray::Array<T const, 1> const& mag) {
68 ndarray::Array<T, 1> out = ndarray::allocate(mag.getShape());
69 for (std::size_t ii = 0; ii < mag.getNumElements(); ++ii) {
70 out[ii] = fluxFromABMag(mag[ii]);
71 }
72 return out;
73}
74
76template <typename T>
77ndarray::Array<T, 1> fluxErrFromABMagErr(ndarray::Array<T const, 1> const& magErr,
78 ndarray::Array<T const, 1> const& mag) {
79 if (mag.getNumElements() != magErr.getNumElements()) {
80 throw LSST_EXCEPT(pex::exceptions::LengthError, (boost::format("Length mismatch: %d vs %d") %
81 mag.getNumElements() % magErr.getNumElements())
82 .str());
83 }
84 ndarray::Array<T, 1> out = ndarray::allocate(mag.getShape());
85 for (std::size_t ii = 0; ii < mag.getNumElements(); ++ii) {
86 out[ii] = fluxErrFromABMagErr(magErr[ii], mag[ii]);
87 }
88 return out;
89}
90
91// Explicit instantiation
92#define INSTANTIATE(TYPE) \
93 template ndarray::Array<TYPE, 1> abMagFromFlux(ndarray::Array<TYPE const, 1> const& flux); \
94 template ndarray::Array<TYPE, 1> abMagErrFromFluxErr(ndarray::Array<TYPE const, 1> const& fluxErr, \
95 ndarray::Array<TYPE const, 1> const& flux); \
96 template ndarray::Array<TYPE, 1> fluxFromABMag(ndarray::Array<TYPE const, 1> const& mag); \
97 template ndarray::Array<TYPE, 1> fluxErrFromABMagErr(ndarray::Array<TYPE const, 1> const& magErr, \
98 ndarray::Array<TYPE const, 1> const& mag);
99
102
103} // namespace image
104} // namespace afw
105} // namespace lsst
#define INSTANTIATE(FROMSYS, TOSYS)
Definition Detector.cc:509
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
Reports attempts to exceed implementation-defined length limits for some classes.
Definition Runtime.h:76
double abMagErrFromFluxErr(double fluxErr, double flux)
Compute AB magnitude error from flux and flux error in Janskys.
Definition Calib.h:60
double fluxErrFromABMagErr(double magErr, double mag) noexcept
Compute flux error in Janskys from AB magnitude error and AB magnitude.
Definition Calib.h:68
double abMagFromFlux(double flux)
Compute AB magnitude from flux in Janskys.
Definition Calib.h:57
double fluxFromABMag(double mag) noexcept
Compute flux in Janskys from AB magnitude.
Definition Calib.h:65