LSSTApplications  18.1.0
LSSTDataManagementBasePackage
slots.cc
Go to the documentation of this file.
1 #include "lsst/pex/exceptions.h"
2 #include "lsst/afw/table/slots.h"
3 
4 namespace lsst {
5 namespace afw {
6 namespace table {
7 
8 namespace {
9 
10 // Return true if 'a' starts with 'b'
11 bool startsWith(std::string const &a, std::string const b) { return a.compare(0, b.size(), b) == 0; }
12 
13 // helper class that resolves aliases for a slot's main measurement field, while
14 // distiguishing between two cases where the field can't be found:
15 // - if the alias points to an invalid field, we should throw an exception (which is actually
16 // done by the calling setKeys() method, by allowing a lower-level exception to propagate up).
17 // - if the alias simply isn't defined, we should just reset the slot with invalid keys
18 class MeasFieldNameGetter {
19 public:
20  MeasFieldNameGetter(SubSchema const &s, Schema const &schema)
21  : replaced(schema[schema.getAliasMap()->apply(s.getPrefix())]),
22  defined(replaced.getPrefix() !=
23  s.getPrefix()) // slot is defined if applying alias wasn't a no-op
24  {}
25 
26  SubSchema replaced; // a SubSchema that includes alias replacement
27  bool defined; // whether the slot is defined at all
28 };
29 
30 } // namespace
31 
33  SubSchema s = schema["slot"][_name];
34  if (!alias.empty() && !startsWith(alias, s.getPrefix())) return;
35  _measKey = MeasKey();
36  _errKey = ErrKey();
37  _flagKey = Key<Flag>();
38  MeasFieldNameGetter helper(s["instFlux"], schema);
39  if (!helper.defined) {
40  return;
41  }
42  _measKey = helper.replaced;
43  try {
44  _errKey = s["instFluxErr"];
45  } catch (pex::exceptions::NotFoundError &) {
46  }
47  try {
48  _flagKey = s["flag"];
49  } catch (pex::exceptions::NotFoundError &) {
50  }
51 }
52 
53 namespace {
54 
55 CentroidSlotDefinition::ErrKey::NameArray makeCentroidNameArray() {
57  v.push_back("x");
58  v.push_back("y");
59  return v;
60 }
61 
62 } // namespace
63 
65  SubSchema s = schema["slot"][_name];
66  if (!alias.empty() && !startsWith(alias, s.getPrefix())) return;
67  static ErrKey::NameArray names = makeCentroidNameArray();
68  _measKey = MeasKey();
69  _errKey = ErrKey();
70  _flagKey = Key<Flag>();
71  MeasFieldNameGetter helper(s, schema);
72  if (!helper.defined) return;
73  _measKey = helper.replaced;
74  try {
75  _errKey = ErrKey(s, names);
76  } catch (pex::exceptions::NotFoundError &) {
77  }
78  try {
79  _flagKey = s["flag"];
80  } catch (pex::exceptions::NotFoundError &) {
81  }
82 }
83 
84 namespace {
85 
86 ShapeSlotDefinition::ErrKey::NameArray makeShapeNameArray() {
88  v.push_back("xx");
89  v.push_back("yy");
90  v.push_back("xy");
91  return v;
92 }
93 
94 } // namespace
95 
97  SubSchema s = schema["slot"][_name];
98  if (!alias.empty() && !startsWith(alias, s.getPrefix())) return;
99  static ErrKey::NameArray names = makeShapeNameArray();
100  _measKey = MeasKey();
101  _errKey = ErrKey();
102  _flagKey = Key<Flag>();
103  MeasFieldNameGetter helper(s, schema);
104  if (!helper.defined) return;
105  _measKey = helper.replaced;
106  try {
107  _errKey = ErrKey(s, names);
108  } catch (pex::exceptions::NotFoundError &) {
109  }
110  try {
111  _flagKey = s["flag"];
112  } catch (pex::exceptions::NotFoundError &) {
113  }
114 }
115 
117  defPsfFlux.setKeys(alias, schema);
118  defApFlux.setKeys(alias, schema);
119  defGaussianFlux.setKeys(alias, schema);
120  defModelFlux.setKeys(alias, schema);
121  defCalibFlux.setKeys(alias, schema);
122  defCentroid.setKeys(alias, schema);
123  defShape.setKeys(alias, schema);
124 }
125 
127  : defPsfFlux("PsfFlux"),
128  defApFlux("ApFlux"),
129  defGaussianFlux("GaussianFlux"),
130  defModelFlux("ModelFlux"),
131  defCalibFlux("CalibFlux"),
132  defCentroid("Centroid"),
133  defShape("Shape") {
134  defPsfFlux.setKeys("", schema);
135  defApFlux.setKeys("", schema);
136  defGaussianFlux.setKeys("", schema);
137  defModelFlux.setKeys("", schema);
138  defCalibFlux.setKeys("", schema);
139  defCentroid.setKeys("", schema);
140  defShape.setKeys("", schema);
141 }
142 } // namespace table
143 } // namespace afw
144 } // namespace lsst
Defines the fields and offsets for a table.
Definition: Schema.h:50
T empty(T... args)
A proxy type for name lookups in a Schema.
Definition: Schema.h:360
FluxSlotDefinition defCalibFlux
Definition: slots.h:224
table::Key< int > b
FluxSlotDefinition defPsfFlux
Definition: slots.h:220
table::Key< int > a
CentroidSlotDefinition defCentroid
Definition: slots.h:225
void setKeys(std::string const &alias, Schema const &schema)
Update the cached Keys following an change of aliases in the given Schema.
Definition: slots.cc:96
std::string const & getPrefix() const
Return the prefix that defines this SubSchema relative to its parent Schema.
Definition: Schema.h:413
STL class.
Reports attempts to access elements using an invalid key.
Definition: Runtime.h:151
T push_back(T... args)
A base class for image defects.
ShapeSlotDefinition defShape
Definition: slots.h:226
FluxSlotDefinition defModelFlux
Definition: slots.h:223
table::Schema schema
Definition: Camera.cc:161
solver_t * s
SlotSuite(Schema const &schema)
Initialize the slots.
Definition: slots.cc:126
T size(T... args)
SubSchema replaced
Definition: slots.cc:26
bool defined
Definition: slots.cc:27
Key specialization for Flag.
Definition: Flag.h:94
void setKeys(std::string const &alias, Schema const &schema)
Update the cached Keys following an change of aliases in the given Schema.
Definition: slots.cc:32
void setKeys(std::string const &alias, Schema const &schema)
Update the cached Keys following an change of aliases in the given Schema.
Definition: slots.cc:64
A FunctorKey used to get or set a geom::ellipses::Quadrupole from a tuple of constituent Keys...
Definition: aggregates.h:282
FluxSlotDefinition defGaussianFlux
Definition: slots.h:222
T compare(T... args)
FluxSlotDefinition defApFlux
Definition: slots.h:221
void handleAliasChange(std::string const &alias, Schema const &schema)
Handle a callback from an AliasMap informing the table that an alias has changed. ...
Definition: slots.cc:116