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
shapeprior.h
Go to the documentation of this file.
1#ifndef LSST_GAUSS2D_FIT_SHAPEPRIOR_H
2#define LSST_GAUSS2D_FIT_SHAPEPRIOR_H
3
4#include <memory>
5
7
8#include "parametricellipse.h"
10#include "prior.h"
11#include "transforms.h"
12
13namespace lsst::gauss2d::fit {
14
18class ShapePriorOptions : public Object {
19public:
20 static inline const double delta_jacobian_default = 1e-5;
21 static inline const double size_maj_floor_default = 1e-3;
22 static inline const double axrat_floor_default = 1e-3;
23
31 explicit ShapePriorOptions(double delta_jacobian = delta_jacobian_default,
32 double size_maj_floor = size_maj_floor_default,
33 double axrat_floor = axrat_floor_default);
34
35 bool check_delta_jacobian(double delta_jacobian, bool do_throw = false);
36 bool check_size_maj_floor(double size_maj_floor, bool do_throw = false);
37 bool check_axrat_floor(double axrat_floor, bool do_throw = false);
38
39 double get_delta_jacobian() const;
40 double get_size_maj_floor() const;
41 double get_axrat_floor() const;
42
43 void set_delta_jacobian(double delta_jacobian);
44 void set_size_maj_floor(double size_maj_floor);
45 void set_axrat_floor(double axrat_floor);
46
47 std::string repr(bool name_keywords = false,
48 std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
49 std::string str() const override;
50
51private:
52 double _delta_jacobian;
53 double _size_maj_floor;
54 double _axrat_floor;
55};
56
63class ShapePrior : public Prior {
64public:
73 std::shared_ptr<ParametricGaussian1D> prior_size = nullptr,
74 std::shared_ptr<ParametricGaussian1D> prior_axrat = nullptr,
77
78 PriorEvaluation evaluate(bool calc_jacobians = false, bool normalize_loglike = false) const override;
79
82
84
87
88 size_t size() const override;
89
90 std::string repr(bool name_keywords = false,
91 std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
92 std::string str() const override;
93
94private:
99};
100} // namespace lsst::gauss2d::fit
101
102#endif
A generic object from the gauss2d library.
Definition object.h:40
static constexpr std::string_view CC_NAMESPACE_SEPARATOR
The C++ namespace separator.
Definition object.h:45
Results from the evaluation of a prior probability function.
Definition prior.h:16
Interface for a prior probability function.
Definition prior.h:38
A two-part prior on the shape of a parametric ellipse.
Definition shapeprior.h:63
std::string repr(bool name_keywords=false, std::string_view namespace_separator=Object::CC_NAMESPACE_SEPARATOR) const override
Return a full, callable string representation of this.
size_t size() const override
PriorEvaluation evaluate(bool calc_jacobians=false, bool normalize_loglike=false) const override
Evaluate the log likelihood and residual-dependent terms.
std::string str() const override
Return a brief, human-readable string representation of this.
std::shared_ptr< ParametricGaussian1D > get_prior_axrat() const
ShapePrior(std::shared_ptr< const ParametricEllipse > ellipse, std::shared_ptr< ParametricGaussian1D > prior_size=nullptr, std::shared_ptr< ParametricGaussian1D > prior_axrat=nullptr, std::shared_ptr< ShapePriorOptions > options=nullptr)
Construct a ShapePrior from a Parameter and mean_size/std.
void set_prior_size(std::shared_ptr< ParametricGaussian1D > prior_size)
std::shared_ptr< ParametricGaussian1D > get_prior_size() const
std::vector< double > get_loglike_const_terms() const override
Return the constant terms of the log likelihood (dependent on stddevs only)
void set_prior_axrat(std::shared_ptr< ParametricGaussian1D > prior_axrat)
Options for a ShapePrior.
Definition shapeprior.h:18
static const double delta_jacobian_default
Definition shapeprior.h:20
void set_size_maj_floor(double size_maj_floor)
Definition shapeprior.cc:86
bool check_axrat_floor(double axrat_floor, bool do_throw=false)
Definition shapeprior.cc:64
bool check_size_maj_floor(double size_maj_floor, bool do_throw=false)
Definition shapeprior.cc:53
static const double size_maj_floor_default
Definition shapeprior.h:21
std::string str() const override
Return a brief, human-readable string representation of this.
void set_delta_jacobian(double delta_jacobian)
Definition shapeprior.cc:81
void set_axrat_floor(double axrat_floor)
Definition shapeprior.cc:91
ShapePriorOptions(double delta_jacobian=delta_jacobian_default, double size_maj_floor=size_maj_floor_default, double axrat_floor=axrat_floor_default)
Construct a ShapePriorOptions from values.
Definition shapeprior.cc:35
bool check_delta_jacobian(double delta_jacobian, bool do_throw=false)
Definition shapeprior.cc:42
std::string repr(bool name_keywords=false, std::string_view namespace_separator=Object::CC_NAMESPACE_SEPARATOR) const override
Return a full, callable string representation of this.
Definition shapeprior.cc:96
static const double axrat_floor_default
Definition shapeprior.h:22