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
FluxUtilities.cc
Go to the documentation of this file.
1// -*- lsst-c++ -*-
2/*
3 * LSST Data Management System
4 * Copyright 2008-2015 AURA/LSST.
5 *
6 * This product includes software developed by the
7 * LSST Project (http://www.lsst.org/).
8 *
9 * This program is free software: you can redistribute it and/or modify
10 * it under the terms of the GNU General Public License as published by
11 * the Free Software Foundation, either version 3 of the License, or
12 * (at your option) any later version.
13 *
14 * This program is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 * GNU General Public License for more details.
18 *
19 * You should have received a copy of the LSST License Statement and
20 * the GNU General Public License along with this program. If not,
21 * see <http://www.lsstcorp.org/LegalNotices/>.
22 */
23
27
28namespace lsst {
29namespace meas {
30namespace base {
31
33 : instFlux(std::numeric_limits<meas::base::Flux>::quiet_NaN()),
34 instFluxErr(std::numeric_limits<meas::base::FluxErrElement>::quiet_NaN()) {}
35
37 std::string const& doc) {
39 result._instFlux = schema.addField<meas::base::Flux>(schema.join(name, "instFlux"), doc, "count");
40 result._instFluxErr = schema.addField<meas::base::FluxErrElement>(
41 schema.join(name, "instFluxErr"), "1-sigma instFlux uncertainty", "count");
42 return result;
43}
44
46 FluxResult r;
47 r.instFlux = record.get(_instFlux);
48 r.instFluxErr = record.get(_instFluxErr);
49 return r;
50}
51
52void FluxResultKey::set(afw::table::BaseRecord& record, FluxResult const& value) const {
53 record.set(_instFlux, value.instFlux);
54 record.set(_instFluxErr, value.instFluxErr);
55}
56
59 result._magKey = schema.addField<Mag>(schema.join(name, "mag"), "Magnitude");
60 result._magErrKey = schema.addField<MagErrElement>(schema.join(name, "magErr"), "Error on magnitude");
61 return result;
62}
63
65 MagResult result = {record.get(_magKey), record.get(_magErrKey)};
66 return result;
67}
68
69void MagResultKey::set(afw::table::BaseRecord& record, MagResult const& magResult) const {
70 record.set(_magKey, magResult.mag);
71 record.set(_magErrKey, magResult.magErr);
72}
73
75 record.set(_magKey, magResult.value);
76 record.set(_magErrKey, magResult.error);
77}
78
80 : BaseTransform{name} {
81 // Map the flag through to the output
82 mapper.addMapping(mapper.getInputSchema().find<afw::table::Flag>(name + "_flag").key);
83
84 // Add keys for the magnitude and associated error
85 _magKey = MagResultKey::addFields(mapper.editOutputSchema(), name);
86}
87
89 afw::table::BaseCatalog& outputCatalog, afw::geom::SkyWcs const& wcs,
90 afw::image::PhotoCalib const& photoCalib) const {
91 checkCatalogSize(inputCatalog, outputCatalog);
92 FluxResultKey instFluxKey(inputCatalog.getSchema()[_name]);
93 afw::table::SourceCatalog::const_iterator inSrc = inputCatalog.begin();
94 afw::table::BaseCatalog::iterator outSrc = outputCatalog.begin();
95 {
96 for (; inSrc != inputCatalog.end() && outSrc != outputCatalog.end(); ++inSrc, ++outSrc) {
97 FluxResult instFluxResult = instFluxKey.get(*inSrc);
98 _magKey.set(*outSrc,
99 photoCalib.instFluxToMagnitude(instFluxResult.instFlux, instFluxResult.instFluxErr));
100 }
101 }
102}
103
104} // namespace base
105} // namespace meas
106} // namespace lsst
SchemaMapper * mapper
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition SkyWcs.h:117
The photometric calibration of an exposure.
Definition PhotoCalib.h:114
Base class for all records.
Definition BaseRecord.h:31
Field< T >::Value get(Key< T > const &key) const
Return the value of a field for the given key.
Definition BaseRecord.h:151
void set(Key< T > const &key, U const &value)
Set value of a field for the given key.
Definition BaseRecord.h:164
Iterator class for CatalogT.
Definition Catalog.h:40
iterator begin()
Iterator access.
Definition Catalog.h:400
Defines the fields and offsets for a table.
Definition Schema.h:51
A mapping between the keys of two Schemas, used to copy data between them.
typename Base::const_iterator const_iterator
Abstract base class for all C++ measurement transformations.
Definition Transform.h:86
void checkCatalogSize(afw::table::BaseCatalog const &cat1, afw::table::BaseCatalog const &cat2) const
Ensure that catalogs have the same size.
Definition Transform.h:102
A FunctorKey for FluxResult.
static FluxResultKey addFields(afw::table::Schema &schema, std::string const &name, std::string const &doc)
Add a pair of _instFlux, _instFluxErr fields to a Schema, and return a FluxResultKey that points to t...
virtual FluxResult get(afw::table::BaseRecord const &record) const
Get a FluxResult from the given record.
virtual void set(afw::table::BaseRecord &record, FluxResult const &other) const
Set a FluxResult in the given record.
virtual void operator()(afw::table::SourceCatalog const &inputCatalog, afw::table::BaseCatalog &outputCatalog, afw::geom::SkyWcs const &wcs, afw::image::PhotoCalib const &photoCalib) const
FluxTransform(std::string const &name, afw::table::SchemaMapper &mapper)
A FunctorKey for MagResult.
virtual void set(afw::table::BaseRecord &record, MagResult const &magResult) const
Set a MagResult in the given record.
static MagResultKey addFields(afw::table::Schema &schema, std::string const &name)
Add a pair of _mag, _magErr fields to a Schema, and return a MagResultKey that points to them.
virtual MagResult get(afw::table::BaseRecord const &record) const
Get a MagResult from the given record.
This defines the base of measurement transformations.
double MagErrElement
Definition constants.h:54
double FluxErrElement
Definition constants.h:52
STL namespace.
A value and its error.
Definition PhotoCalib.h:51
A reusable result struct for instFlux measurements.
meas::base::Flux instFlux
Measured instFlux in DN.
meas::base::FluxErrElement instFluxErr
Standard deviation of instFlux in DN.
FluxResult()
Default constructor; initializes everything to NaN.
A reusable result struct for magnitudes.