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
ScaledApertureFlux.cc
Go to the documentation of this file.
1// -*- lsst-c++ -*-
2/*
3 * LSST Data Management System
4 * Copyright 2008-2016 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
29
30namespace lsst {
31namespace meas {
32namespace base {
33
35 afw::table::Schema& schema)
36 : _ctrl(ctrl),
37 _instFluxResultKey(
38 FluxResultKey::addFields(schema, name, "instFlux derived from PSF-scaled aperture")),
39 _centroidExtractor(schema, name) {
41}
42
44 afw::image::Exposure<float> const& exposure) const {
45 geom::Point2D const center = _centroidExtractor(measRecord, _flagHandler);
46 double const radius = exposure.getPsf()->computeShape(center).getDeterminantRadius();
47 double const fwhm = 2.0 * std::sqrt(2.0 * std::log(2)) * radius;
48 double const size = _ctrl.scale * fwhm;
49 afw::geom::ellipses::Axes const axes(size, size);
50
51 // ApertureFluxAlgorithm::computeSincFlux requires an ApertureFluxControl as an
52 // argument. All that it uses it for is to read the type of warping kernel.
54 apCtrl.shiftKernel = _ctrl.shiftKernel;
55
57 exposure.getMaskedImage(), afw::geom::ellipses::Ellipse(axes, center), apCtrl);
58 measRecord.set(_instFluxResultKey, result);
59
62 if (result.getFlag(iter.number)) {
63 _flagHandler.setValue(measRecord, iter.number, true);
64 }
65 }
66}
67
69 _flagHandler.handleFailure(measRecord, error);
70}
71
74 : FluxTransform{name, mapper} {
78 mapper.getInputSchema().find<afw::table::Flag>(name + "_" + flagName).key;
79 if (key.isValid()) {
80 mapper.addMapping(key);
81 }
82 }
83}
84
85} // namespace base
86} // namespace meas
87} // namespace lsst
An ellipse core for the semimajor/semiminor axis and position angle parametrization (a,...
Definition Axes.h:47
An ellipse defined by an arbitrary BaseCore and a center point.
Definition Ellipse.h:51
A class to contain the data, WCS, and other information needed to describe an image of the sky.
Definition Exposure.h:72
void set(Key< T > const &key, U const &value)
Set value of a field for the given key.
Definition BaseRecord.h:164
A class used as a handle to a particular field in a table.
Definition Key.h:53
bool isValid() const noexcept
Return true if the key was initialized to valid offset.
Definition Key.h:97
Defines the fields and offsets for a table.
Definition Schema.h:51
SchemaItem< T > find(std::string const &name) const
Find a SchemaItem in the Schema by name.
Definition Schema.cc:467
A mapping between the keys of two Schemas, used to copy data between them.
Key< T > addMapping(Key< T > const &inputKey, bool doReplace=false)
Add a new field to the output Schema that is a copy of a field in the input Schema.
Schema const getInputSchema() const
Return the input schema (copy-on-write).
Record class that contains measurements made on a single exposure.
Definition Source.h:78
static FlagDefinitionList const & getFlagDefinitions()
static Result computeSincFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using Sinc photometry.
Configuration object for multiple-aperture flux algorithms.
std::string shiftKernel
"Warping kernel used to shift Sinc photometry coefficients to different center positions" ;
std::size_t size() const
return the current size (number of defined elements) of the collection
static FlagHandler addFields(afw::table::Schema &schema, std::string const &prefix, FlagDefinitionList const &flagDefs, FlagDefinitionList const &exclDefs=FlagDefinitionList::getEmptyList())
Add Flag fields to a schema, creating a FlagHandler object to manage them.
A FunctorKey for FluxResult.
FluxTransform(std::string const &name, afw::table::SchemaMapper &mapper)
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Definition exceptions.h:48
ScaledApertureFluxAlgorithm(Control const &control, std::string const &name, afw::table::Schema &schema)
virtual void measure(afw::table::SourceRecord &measRecord, afw::image::Exposure< float > const &exposure) const override
Measure the scaled aperture instFlux on the given image.
virtual void fail(afw::table::SourceRecord &measRecord, MeasurementError *error=nullptr) const override
Handle an exception thrown by the current algorithm by setting flags in the given record.
ScaledApertureFluxTransform(Control const &ctrl, std::string const &name, afw::table::SchemaMapper &mapper)
T log(T... args)
Point< double, 2 > Point2D
Definition Point.h:324
T sqrt(T... args)
Simple class used to define and document flags The name and doc constitute the identity of the FlagDe...
Definition FlagHandler.h:40