lsst::gauss2d::fit::SersicEllipseData Class Reference

Inheritance diagram for lsst::gauss2d::fit::SersicEllipseData:

Public Member Functions
	SersicEllipseData (const std::shared_ptr< const ReffXParameterD > size_x, const std::shared_ptr< const ReffYParameterD > size_y, const std::shared_ptr< const RhoParameterD > rho, const std::shared_ptr< const SersicMixComponentIndexParameterD > sersicindex, unsigned short index)
double	get_sizeratio () const
double	get_sigma_x () const override
	Get sigma_x.
double	get_sigma_y () const override
	Get sigma_y.
double	get_rho () const override
	Get rho.
double	get_size_x () const override
	Get the size_x value.
double	get_size_y () const override
	Get the size_y value.
std::array< double, 3 >	get_xyr () const override
	Get sigma_x, sigma_y, rho.
void	set (double sigma_x, double sigma_y, double rho) override
	Set sigma_x, sigma_y, rho.
void	set_h (double hwhm_x, double hwhm_y, double rho) override
	Set hwhm_x, hwhm_y, rho (half-width at half-max)
void	set_sigma_x (double sigma_x) override
	Set the x-axis dispersion (sigma)
void	set_sigma_y (double sigma_y) override
	Set the y-axis dispersion (sigma)
void	set_rho (double rho) override
	Set the correlation parameter (rho)
void	set_hxyr (const std::array< double, 3 > &hxyr) override
	Set hwhm_x, hwhm_y, rho from an array.
void	set_xyr (const std::array< double, 3 > &xyr) override
	Set sigma_x, sigma_y, rho from an array.
std::string	repr (bool name_keywords, std::string_view namespace_separator) const override
	Return a full, callable string representation of this.
std::string	str () const override
	Return a brief, human-readable string representation of this.
virtual void	convolve (const Ellipse &ell)
	Convolve this ellipse with another.
virtual double	get_area () const
	Return the area of this ellipse, equal to pi*sigma_major*sigma_minor.
virtual double	get_cov_xy () const
	Return the covariance, equal to sigma_x*sigma_y*rho.
virtual double	get_hwhm_x () const
	Get the x-axis half-width at half-maximum.
virtual double	get_hwhm_y () const
	Get the y-axis half-width at half-maximum.
virtual double	get_radius_trace () const
	Return the trace radius, equal to sqrt(sigma_x^2 + sigma_y^2)
virtual double	get_sigma_x_sq () const
	Get the square of sigma_x.
virtual double	get_sigma_y_sq () const
	Get the square of sigma_y.
virtual double	get_sigma_xy () const
	Return sigma_x*sigma_y.
virtual std::array< double, 3 >	get_hxyr () const
	Get hwhm_x, hwhm_y, rho.
virtual void	set (const Covariance &covar)
	Set values from a Covariance object.
virtual void	set (const EllipseMajor &ellipse)
	Set values from an EllipseMajor object.
virtual void	set_hwhm_x (double hwhm_x)
	Set the x-axis half-width at half-max (FWHM/2)
virtual void	set_hwhm_y (double hwhm_y)
	Set the y-axis half-width at half-max (FWHM/2)
bool	operator== (const EllipseData &other) const
bool	operator!= (const EllipseData &other) const

Static Public Member Functions
static void	check (double size_x, double size_y, double rho, std::string_view error_suffix="")
	Check whether Ellipse parameter values are valid, throwing if not.
static void	check_size (double size, std::string_view error_suffix="")
	Check whether an x- or x-yaxis size value is valid.
static void	check_rho (double rho, std::string_view error_suffix="")
	Check whether a rho value is valid.
static std::string_view	null_str (const std::string_view &namespace_separator)

Static Public Attributes
static constexpr std::string_view	CC_NAMESPACE_SEPARATOR = "::"
	The C++ namespace separator.
static constexpr std::string_view	NULL_STR_GENERAL = "None"
static constexpr std::string_view	PY_NAMESPACE_SEPARATOR = "."

Detailed Description

Definition at line 23 of file sersicmixcomponent.cc.

Constructor & Destructor Documentation

SersicEllipseData()

lsst::gauss2d::fit::SersicEllipseData::SersicEllipseData ( const std::shared_ptr< const ReffXParameterD > size_x, const std::shared_ptr< const ReffYParameterD > size_y, const std::shared_ptr< const RhoParameterD > rho, const std::shared_ptr< const SersicMixComponentIndexParameterD > sersicindex, unsigned short index )

inline

Definition at line 25 of file sersicmixcomponent.cc.

            : _size_x(std::move(size_x)),
              _size_y(std::move(size_y)),
              _rho(std::move(rho)),
              _sersicindex(std::move(sersicindex)),
              _index(index) {
        if ((_size_x == nullptr) || (_size_y == nullptr) || (_rho == nullptr) || (_sersicindex == nullptr)) {
            throw std::invalid_argument("SersicEllipseData args must not be nullptr");
        }
        if (!(_index < _sersicindex->get_order())) {
            throw std::invalid_argument("index=" + std::to_string(_index) + "!< sersicindex->get_order()="
                                        + std::to_string(_sersicindex->get_order()));
        }
    }

Member Function Documentation

check()

static void lsst::gauss2d::EllipseData::check ( double size_x, double size_y, double rho, std::string_view error_suffix = "" )

inlinestaticinherited

Check whether Ellipse parameter values are valid, throwing if not.

Definition at line 137 of file ellipse.h.

                                                                                                  {
        if (!(size_x >= 0) || !(size_y >= 0) || !(rho >= -1 && rho <= 1)) {
            throw std::invalid_argument("Invalid size_x, size_y, rho=" + to_string_float(size_x) + ","
                                        + to_string_float(size_y) + "," + to_string_float(rho)
                                        + "; sigma_x,y >= 0 and 1 >= rho >= -1 required."
                                        + std::string(error_suffix));
        }
    }

check_rho()

static void lsst::gauss2d::EllipseData::check_rho ( double rho, std::string_view error_suffix = "" )

inlinestaticinherited

Check whether a rho value is valid.

Definition at line 153 of file ellipse.h.

                                                                        {
        if (!(rho >= -1 && rho <= 1)) {
            throw std::invalid_argument("Invalid rho=" + to_string_float(rho) + "; 1 >= rho >= -1 required."
                                        + std::string(error_suffix));
        }
    }

check_size()

static void lsst::gauss2d::EllipseData::check_size ( double size, std::string_view error_suffix = "" )

inlinestaticinherited

Check whether an x- or x-yaxis size value is valid.

Definition at line 146 of file ellipse.h.

                                                                          {
        if (!(size >= 0)) {
            throw std::invalid_argument("Invalid size=" + to_string_float(size) + "; size >= 0 required."
                                        + std::string(error_suffix));
        }
    }

convolve()

void lsst::gauss2d::EllipseData::convolve ( const Ellipse & ell )

virtualinherited

Convolve this ellipse with another.

Definition at line 149 of file ellipse.cc.

                                             {
    double sigma_x_ell = ell.get_sigma_x();
    double sigma_y_ell = ell.get_sigma_y();
    double sigma_x = this->get_sigma_x();
    sigma_x = sqrt(sigma_x * sigma_x + sigma_x_ell * sigma_x_ell);
    double sigma_y = this->get_sigma_y();
    sigma_y = sqrt(sigma_y * sigma_y + sigma_y_ell * sigma_y_ell);
    double rho = (this->get_cov_xy() + ell.get_cov_xy()) / (sigma_x * sigma_y);
    this->set(sigma_x, sigma_y, rho);
}

get_area()

double lsst::gauss2d::EllipseData::get_area ( ) const

virtualinherited

Return the area of this ellipse, equal to pi*sigma_major*sigma_minor.

Definition at line 142 of file ellipse.cc.

                                   {
    double rho = this->get_rho();
    return M_PI * this->get_sigma_xy() * sqrt(1 - rho * rho);
}

get_cov_xy()

double lsst::gauss2d::EllipseData::get_cov_xy ( ) const

virtualinherited

Return the covariance, equal to sigma_x*sigma_y*rho.

Definition at line 160 of file ellipse.cc.

{ return this->get_sigma_x() * this->get_sigma_y() * this->get_rho(); }

get_hwhm_x()

double lsst::gauss2d::EllipseData::get_hwhm_x ( ) const

virtualinherited

Get the x-axis half-width at half-maximum.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 162 of file ellipse.cc.

{ return M_SIGMA_HWHM * this->get_sigma_x(); };

get_hwhm_y()

double lsst::gauss2d::EllipseData::get_hwhm_y ( ) const

virtualinherited

Get the y-axis half-width at half-maximum.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 164 of file ellipse.cc.

{ return M_SIGMA_HWHM * this->get_sigma_y(); };

get_hxyr()

std::array< double, 3 > lsst::gauss2d::EllipseData::get_hxyr ( ) const

virtualinherited

Get hwhm_x, hwhm_y, rho.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 166 of file ellipse.cc.

                                                {
    return {this->get_hwhm_x(), this->get_hwhm_y(), this->get_rho()};
}

get_radius_trace()

double lsst::gauss2d::EllipseData::get_radius_trace ( ) const

virtualinherited

Return the trace radius, equal to sqrt(sigma_x^2 + sigma_y^2)

Definition at line 173 of file ellipse.cc.

{ return sqrt(this->get_sigma_x_sq() + this->get_sigma_y_sq()); }

get_rho()

double lsst::gauss2d::fit::SersicEllipseData::get_rho ( ) const

inlineoverridevirtual

Get rho.

Implements lsst::gauss2d::EllipseData.

Definition at line 48 of file sersicmixcomponent.cc.

{ return _rho->get_value(); }

get_sigma_x()

double lsst::gauss2d::fit::SersicEllipseData::get_sigma_x ( ) const

inlineoverridevirtual

Get sigma_x.

Implements lsst::gauss2d::EllipseData.

Definition at line 46 of file sersicmixcomponent.cc.

{ return get_sizeratio() * _size_x->get_value(); }

get_sigma_x_sq()

double lsst::gauss2d::EllipseData::get_sigma_x_sq ( ) const

virtualinherited

Get the square of sigma_x.

Definition at line 175 of file ellipse.cc.

                                         {
    double sigma = this->get_sigma_x();
    return sigma * sigma;
}

get_sigma_xy()

double lsst::gauss2d::EllipseData::get_sigma_xy ( ) const

virtualinherited

Return sigma_x*sigma_y.

Definition at line 147 of file ellipse.cc.

{ return this->get_sigma_x() * this->get_sigma_y(); }

get_sigma_y()

double lsst::gauss2d::fit::SersicEllipseData::get_sigma_y ( ) const

inlineoverridevirtual

Get sigma_y.

Implements lsst::gauss2d::EllipseData.

Definition at line 47 of file sersicmixcomponent.cc.

{ return get_sizeratio() * _size_y->get_value(); }

get_sigma_y_sq()

double lsst::gauss2d::EllipseData::get_sigma_y_sq ( ) const

virtualinherited

Get the square of sigma_y.

Definition at line 180 of file ellipse.cc.

                                         {
    double sigma = this->get_sigma_y();
    return sigma * sigma;
}

get_size_x()

double lsst::gauss2d::fit::SersicEllipseData::get_size_x ( ) const

inlineoverridevirtual

Get the size_x value.

Implements lsst::gauss2d::fit::QuasiEllipse.

Definition at line 50 of file sersicmixcomponent.cc.

{ return get_sigma_x(); };

get_size_y()

double lsst::gauss2d::fit::SersicEllipseData::get_size_y ( ) const

inlineoverridevirtual

Get the size_y value.

Implements lsst::gauss2d::fit::QuasiEllipse.

Definition at line 51 of file sersicmixcomponent.cc.

{ return get_sigma_y(); };

get_sizeratio()

double lsst::gauss2d::fit::SersicEllipseData::get_sizeratio ( ) const

inline

Definition at line 44 of file sersicmixcomponent.cc.

{ return _sersicindex->get_sizeratio(_index); }

get_xyr()

std::array< double, 3 > lsst::gauss2d::fit::SersicEllipseData::get_xyr ( ) const

inlineoverridevirtual

Get sigma_x, sigma_y, rho.

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 52 of file sersicmixcomponent.cc.

{ return {get_size_x(), get_size_y(), get_rho()}; };

null_str()

static std::string_view lsst::gauss2d::Object::null_str ( const std::string_view & namespace_separator )

inlinestaticinherited

Definition at line 49 of file object.h.

                                                                              {
        return namespace_separator == CC_NAMESPACE_SEPARATOR ? "nullptr" : NULL_STR_GENERAL;
    }

operator!=()

bool lsst::gauss2d::EllipseData::operator!= ( const EllipseData & other ) const

inlineinherited

Definition at line 217 of file ellipse.h.

{ return !(*this == other); };

operator==()

bool lsst::gauss2d::EllipseData::operator== ( const EllipseData & other ) const

inlineinherited

Definition at line 216 of file ellipse.h.

{ return get_xyr() == other.get_xyr(); };

repr()

std::string lsst::gauss2d::fit::SersicEllipseData::repr ( bool name_keywords, std::string_view namespace_separator ) const

inlineoverridevirtual

Return a full, callable string representation of this.

Parameters

name_keywords	Whether to prefix arguments with "{name}=", where name is the arg name in the header (as with keyword arguments in Python).
namespace_separator	The string to use to delimit namespaces, i.e. :: in C++ and . in Python.

Returns

A callable string representation of this, which should return an an identical object to this.

Note

The representation with name_keywords=false must be callable in C++. The representation with name_keywords=true should be callable in Python, if there are any bindings.

Implements lsst::gauss2d::EllipseData.

Definition at line 70 of file sersicmixcomponent.cc.

                                                                                          {
        return type_name_str<SersicEllipseData>(false, namespace_separator) + "("
               + (name_keywords ? "size_x=" : "") + _size_x->repr(name_keywords, namespace_separator) + ", "
               + (name_keywords ? "size_y=" : "") + _size_y->repr(name_keywords, namespace_separator) + ", "
               + (name_keywords ? "rho=" : "") + _rho->repr(name_keywords, namespace_separator) + ")";
    }

set() [1/3]

void lsst::gauss2d::EllipseData::set ( const Covariance & covar )

virtualinherited

Set values from a Covariance object.

Definition at line 263 of file ellipse.cc.

                                             {
    double sigma_x = covar.get_sigma_x_sq();
    double sigma_y = covar.get_sigma_y_sq();
    sigma_x = sqrt(sigma_x);
    sigma_y = sqrt(sigma_y);
    double rho = (sigma_x == 0 || sigma_y == 0) ? 0 : covar.get_cov_xy() / (sigma_x * sigma_y);
    this->set(sigma_x, sigma_y, rho);
}

set() [2/3]

void lsst::gauss2d::EllipseData::set ( const EllipseMajor & ellipse )

virtualinherited

Set values from an EllipseMajor object.

Definition at line 272 of file ellipse.cc.

                                                 {
    const double r_major = ellipse.get_r_major();
    if (r_major == 0) {
        this->set(0, 0, 0);
        return;
    }
    const double axrat = ellipse.get_axrat();
    if (axrat == 1) {
        this->set(r_major, r_major, 0);
        return;
    }
    const auto [sin_th, cos_th] = sincos(ellipse.get_angle_radians());
    const double sin_th_sq = sin_th * sin_th;
    const double cos_th_sq = cos_th * cos_th;
 
    const double r_major_sq = r_major * r_major;
    const double r_minor_sq = r_major_sq * axrat * axrat;
    const double sigma_x = sqrt(cos_th_sq * r_major_sq + sin_th_sq * r_minor_sq);
    const double sigma_y = sqrt(sin_th_sq * r_major_sq + cos_th_sq * r_minor_sq);
    double rho = (sigma_x == 0 || sigma_y == 0)
                         ? 0
                         : sin_th * cos_th * (r_major_sq - r_minor_sq) / (sigma_x * sigma_y);
    this->set(sigma_x, sigma_y, rho);
}

set() [3/3]

void lsst::gauss2d::fit::SersicEllipseData::set ( double sigma_x, double sigma_y, double rho )

inlineoverridevirtual

Set sigma_x, sigma_y, rho.

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 54 of file sersicmixcomponent.cc.

                                                                  {
        throw std::runtime_error("Can't set on SersicEllipseData");
    }

set_h()

void lsst::gauss2d::fit::SersicEllipseData::set_h ( double hwhm_x, double hwhm_y, double rho )

inlineoverridevirtual

Set hwhm_x, hwhm_y, rho (half-width at half-max)

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 57 of file sersicmixcomponent.cc.

                                                                  {
        throw std::runtime_error("Can't set on SersicEllipseData");
    }

set_hwhm_x()

void lsst::gauss2d::EllipseData::set_hwhm_x ( double hwhm_x )

virtualinherited

Set the x-axis half-width at half-max (FWHM/2)

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 303 of file ellipse.cc.

{ this->set_sigma_x(M_HWHM_SIGMA * hwhm_x); }

set_hwhm_y()

void lsst::gauss2d::EllipseData::set_hwhm_y ( double hwhm_y )

virtualinherited

Set the y-axis half-width at half-max (FWHM/2)

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 304 of file ellipse.cc.

{ this->set_sigma_y(M_HWHM_SIGMA * hwhm_y); }

set_hxyr()

void lsst::gauss2d::fit::SersicEllipseData::set_hxyr ( const std::array< double, 3 > & hxyr )

inlineoverridevirtual

Set hwhm_x, hwhm_y, rho from an array.

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 63 of file sersicmixcomponent.cc.

                                                            {
        throw std::runtime_error("Can't set on SersicEllipseData");
    }

set_rho()

void lsst::gauss2d::fit::SersicEllipseData::set_rho ( double rho )

inlineoverridevirtual

Set the correlation parameter (rho)

Implements lsst::gauss2d::EllipseData.

Definition at line 62 of file sersicmixcomponent.cc.

{ throw std::runtime_error("Can't set on SersicEllipseData"); }

set_sigma_x()

void lsst::gauss2d::fit::SersicEllipseData::set_sigma_x ( double sigma_x )

inlineoverridevirtual

Set the x-axis dispersion (sigma)

Implements lsst::gauss2d::EllipseData.

Definition at line 60 of file sersicmixcomponent.cc.

{ throw std::runtime_error("Can't set on SersicEllipseData"); }

set_sigma_y()

void lsst::gauss2d::fit::SersicEllipseData::set_sigma_y ( double sigma_y )

inlineoverridevirtual

Set the y-axis dispersion (sigma)

Implements lsst::gauss2d::EllipseData.

Definition at line 61 of file sersicmixcomponent.cc.

{ throw std::runtime_error("Can't set on SersicEllipseData"); }

set_xyr()

void lsst::gauss2d::fit::SersicEllipseData::set_xyr ( const std::array< double, 3 > & xyr )

inlineoverridevirtual

Set sigma_x, sigma_y, rho from an array.

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 66 of file sersicmixcomponent.cc.

                                                          {
        throw std::runtime_error("Can't set on SersicEllipseData");
    }

str()

std::string lsst::gauss2d::fit::SersicEllipseData::str ( ) const

inlineoverridevirtual

Return a brief, human-readable string representation of this.

Implements lsst::gauss2d::EllipseData.

Definition at line 77 of file sersicmixcomponent.cc.

                                   {
        return type_name_str<SersicEllipseData>(true) + "(size_x=" + _size_x->str()
               + ", size_y=" + _size_y->str() + ", rho=" + _rho->str() + +")";
    }

Member Data Documentation

CC_NAMESPACE_SEPARATOR

std::string_view lsst::gauss2d::Object::CC_NAMESPACE_SEPARATOR = "::"

staticconstexprinherited

The C++ namespace separator.

Definition at line 45 of file object.h.

NULL_STR_GENERAL

std::string_view lsst::gauss2d::Object::NULL_STR_GENERAL = "None"

staticconstexprinherited

Definition at line 46 of file object.h.

PY_NAMESPACE_SEPARATOR

std::string_view lsst::gauss2d::Object::PY_NAMESPACE_SEPARATOR = "."

staticconstexprinherited

Definition at line 47 of file object.h.

---

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

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-12.1.0/Linux/gauss2d_fit/g199a45376c+5137f08352/src/sersicmixcomponent.cc