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
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"
32 #include "lsst/afw/image/Calib.h"
33 
34 namespace lsst {
35 namespace afw {
36 
37 namespace image {
38 
40 template <typename T>
41 ndarray::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 
50 template <typename T>
51 ndarray::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 
66 template <typename T>
67 ndarray::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 
76 template <typename T>
77 ndarray::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 
100 INSTANTIATE(float);
101 INSTANTIATE(double);
102 
103 } // namespace image
104 } // namespace afw
105 } // namespace lsst
#define INSTANTIATE(TYPE)
Definition: Calib.cc:92
#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
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
double abMagErrFromFluxErr(double fluxErr, double flux)
Compute AB magnitude error from flux and flux error in Janskys.
Definition: Calib.h:52
double fluxErrFromABMagErr(double magErr, double mag) noexcept
Compute flux error in Janskys from AB magnitude error and AB magnitude.
Definition: Calib.h:60
double abMagFromFlux(double flux)
Compute AB magnitude from flux in Janskys.
Definition: Calib.h:49
double fluxFromABMag(double mag) noexcept
Compute flux in Janskys from AB magnitude.
Definition: Calib.h:57
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:174
A base class for image defects.