LSST Applications 24.1.6,g063fba187b+56b85ce14a,g0f08755f38+df8a265115,g12f32b3c4e+891a09f10d,g1524ad2192+7a5d7b3fbd,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g28da252d5a+07cb1400be,g2bbee38e9b+ae03bbfc84,g2bc492864f+ae03bbfc84,g3156d2b45e+6e55a43351,g347aa1857d+ae03bbfc84,g35bb328faa+a8ce1bb630,g3a166c0a6a+ae03bbfc84,g3e281a1b8c+c5dd892a6c,g414038480c+6b9177ef31,g41af890bb2+8f257c4c0b,g781aacb6e4+a8ce1bb630,g7af13505b9+7137b3b17d,g80478fca09+6df6903293,g82479be7b0+091ce1d07f,g858d7b2824+df8a265115,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,g9726552aa6+414189b318,ga5288a1d22+4a2bca08d7,gacef1a1666+c9a8ff65f4,gb58c049af0+d64f4d3760,gbcfae0f0a0+de1d42d831,gc28159a63d+ae03bbfc84,gcf0d15dbbd+72117bf34e,gda6a2b7d83+72117bf34e,gdaeeff99f8+1711a396fd,ge500cccec5+c8c9c9af63,ge79ae78c31+ae03bbfc84,gf0baf85859+c1f95f4921,gfa517265be+df8a265115,gfa999e8aa5+17cd334064,gfb92a5be7c+df8a265115
LSST Data Management Base Package
Loading...
Searching...
No Matches
Public Member Functions | List of all members
lsst::meas::base::SafeCentroidExtractor Class Reference

Utility class for measurement algorithms that extracts a position from the Centroid slot and handles errors in a safe and consistent way. More...

#include <InputUtilities.h>

Public Member Functions

 SafeCentroidExtractor (afw::table::Schema &schema, std::string const &name, bool isCentroider=false)
 Construct the extractor, creating a flag alias that indicates failure in the input centroid by linking to the slot centroid flag.
 
geom::Point2D operator() (afw::table::SourceRecord &record, FlagHandler const &flags) const
 Extract a position from the given record.
 

Detailed Description

Utility class for measurement algorithms that extracts a position from the Centroid slot and handles errors in a safe and consistent way.

Definition at line 39 of file InputUtilities.h.

Constructor & Destructor Documentation

◆ SafeCentroidExtractor()

lsst::meas::base::SafeCentroidExtractor::SafeCentroidExtractor ( afw::table::Schema & schema,
std::string const & name,
bool isCentroider = false )

Construct the extractor, creating a flag alias that indicates failure in the input centroid by linking to the slot centroid flag.

Parameters
[out]schemaSchema to which the alias should be added. The "slot_Centroid" alias must already be present in the Schema's AliasMap.
[in]nameThe name of the algorithm; the flag alias added will be "<name>_flag_badCentroid", or "<name>_flag_badInitialCentroid" if isCentroider=true.
[in]isCentroiderIndicates whether the calling algorithm is itself a centroid measurement algorithm. If true,, falling back to the Peak because there was no previous centroider or a previous centroider failed will not cause the general failure flag of the current algorithm to be set.

Definition at line 35 of file InputUtilities.cc.

37 : _name(name), _isCentroider(isCentroider) {
38 // Instead of aliasing e.g. MyAlgorithm_flag_badCentroid->slot_Centroid_flag, we actually
39 // look up the target of slot_Centroid_flag, and alias that to MyAlgorithm_flag_badCentroid.
40 // That way, if someone changes the slots later, after we've already done the measurement,
41 // this alias still points to the right thing.
42 std::string aliasedFlagName = schema.join("slot", "Centroid", "flag");
43 std::string slotFlagName = schema.getAliasMap()->apply(aliasedFlagName);
44 if (_isCentroider) {
45 if (slotFlagName != schema.join(name, "flag")) {
46 // only setup the alias if this isn't the slot algorithm itself (otherwise it'd be circular)
47 schema.getAliasMap()->set(schema.join(name, "flag", "badInitialCentroid"), slotFlagName);
48 }
49 } else {
50 if (aliasedFlagName == slotFlagName) {
51 throw LSST_EXCEPT(
52 pex::exceptions::LogicError,
53 (boost::format("Alias for '%s' must be defined before initializing '%s' plugin.") %
54 aliasedFlagName % name)
55 .str());
56 }
57 schema.getAliasMap()->set(schema.join(name, "flag", "badCentroid"), slotFlagName);
58 }
59}
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
table::Schema schema
Definition python.h:134

Member Function Documentation

◆ operator()()

geom::Point2D lsst::meas::base::SafeCentroidExtractor::operator() ( afw::table::SourceRecord & record,
FlagHandler const & flags ) const

Extract a position from the given record.

We use the Centroid slot if it is not NaN, and fall back to the Peak on the Footprint otherwise.

If the Centroid slot is not defined, we throw FatalAlgorithmError, as this indicates a configuration problem.

If the Centroid slot value is NaN and is not flagged (or there is no Centroid slot flag), we throw RuntimeError, which should cause the measurement framework to log a warning and set the current algorithm's general failure flag if it is allowed to propagate out of the algorithm implementation.

If the Centroid slot is NaN and there is no Footprint or Peak (even if the centroid is flagged), we also throw RuntimeError, as this indicates something wrong in another stage of the pipeline that should be addressed before measurement is run.

If the Centroid slot is flagged and we nevertheless obtain a usable position (either from a the Centroid slot itself or from a successful fall-back to the Peak), we set the current algorithm's general failure flag, but return the position as well, allowing it to continue while indicating that the result may not be reliable.

Definition at line 84 of file InputUtilities.cc.

85 {
86 if (!record.getTable()->getCentroidSlot().getMeasKey().isValid()) {
87 if (_isCentroider) {
88 return extractPeak(record, _name);
89 } else {
90 throw LSST_EXCEPT(
91 FatalAlgorithmError,
92 (boost::format("%s requires a centroid, but the centroid slot is not defined") % _name)
93 .str());
94 }
95 }
96 geom::Point2D result = record.getCentroid();
97 if (std::isnan(result.getX()) || std::isnan(result.getY())) {
98 if (!record.getTable()->getCentroidSlot().getFlagKey().isValid()) {
99 if (_isCentroider) {
100 return extractPeak(record, _name);
101 } else {
102 throw LSST_EXCEPT(
103 pex::exceptions::RuntimeError,
104 (boost::format(
105 "%s: Centroid slot value is NaN, but there is no Centroid slot flag "
106 "(is the executionOrder for %s lower than that of the slot Centroid?)") %
107 _name % _name)
108 .str());
109 }
110 }
111 if (!record.getCentroidFlag() && !_isCentroider) {
112 throw LSST_EXCEPT(
113 pex::exceptions::RuntimeError,
114 (boost::format("%s: Centroid slot value is NaN, but the Centroid slot flag is not set "
115 "(is the executionOrder for %s lower than that of the slot Centroid?)") %
116 _name % _name)
117 .str());
118 }
119 result = extractPeak(record, _name);
120 if (!_isCentroider) {
121 // set the general flag, because using the Peak might affect the current measurement
122 flags.setValue(record, flags.getFailureFlagNumber(), true);
123 }
124 } else if (!_isCentroider && record.getTable()->getCentroidSlot().getFlagKey().isValid() &&
125 record.getCentroidFlag()) {
126 // we got a usable value, but the centroid flag is still be set, and that might affect
127 // the current measurement
128 flags.setValue(record, flags.getFailureFlagNumber(), true);
129 }
130 return result;
131}
py::object result
Definition _schema.cc:429
T isnan(T... args)

The documentation for this class was generated from the following files: