lsst::shapelet::RadialProfile Class Reference

Registry and utility class for multi-Gaussian approximations to radial profiles. More...

#include <RadialProfile.h>

Public Member Functions
std::string const	getName () const
	Return the name of the profile.
virtual double	evaluate (double r) const =0
	Evaluate the profile at the given radius.
ndarray::Array< double, 1, 1 >	evaluate (ndarray::Array< double const, 1, 1 > const &r) const
	Evaluate the profile at the given radius (vectorized).
double	getMomentsRadiusFactor () const
	Return the 2nd-moment radius in units of the half-light radius.
std::shared_ptr< MultiShapeletBasis >	getBasis (int nComponents, int maxRadius=0) const
	Return a multi-Gaussian approximation to the radial profile.
void	registerBasis (std::shared_ptr< MultiShapeletBasis > basis, int nComponents, int maxRadius)
	Register a basis with the profile, making it available to callers of getBasis().

Static Public Member Functions
static RadialProfile &	get (std::string const &name)
	Return a predefined radial profile given its name.

Protected Member Functions
	RadialProfile (std::string const &name, int defaultMaxRadius)

Protected Attributes
double	_momentsRadiusFactor

Detailed Description

Registry and utility class for multi-Gaussian approximations to radial profiles.

RadialProfile provides C++ access to the saved multi-Gaussian approximations stored in the pickle files in the data subdirectory of the shapelet package, by maintaining a singleton registry of these that is populated at import time by tractor.py. It also provides access to the true (exact) profiles, (mostly useful for educational/diagnostic purposes) and information about the relationship between different radii for this profile (see getMomentsRadiusFactor()).

Each RadialProfile object represents a particular "true" profile, and can hold multiple multi-Gaussian (or multi-Shapelet) approximations with different degrees of fidelity, as controlled by the number of components in the approximation and the maximum radius at which the approximation was fit. For more information about these approximations, see tractor.py and references therein. All RadialProfiles are defined in units of the half-light radius of the true profile, even for profiles for which this is not the radius typically used.

Several predefined subclasses of RadialProfile are defined privately, and can be accessed by name through the static get() method:

• "gaussian": a single Gaussian profile
• "exp": Exponential profile (Sersic n=1)
• "ser2": Sersic profile with n=2
• "ser3": Sersic profile with n=3
• "dev": de Vaucouleur profile (Sersic n=4)
• "ser5": Sersic profile with n=5
• "lux": SDSS truncated Exponential profile
• "luv": SDSS truncated and softened de Vaucouleur profile

Definition at line 60 of file RadialProfile.h.

Constructor & Destructor Documentation

RadialProfile()

lsst::shapelet::RadialProfile::RadialProfile ( std::string const & name, int defaultMaxRadius )

explicitprotected

Member Function Documentation

evaluate() [1/2]

virtual double lsst::shapelet::RadialProfile::evaluate ( double r ) const

pure virtual

Evaluate the profile at the given radius.

Radius is given in units of half-light radius, and the profile is normalized to 1 at r=1.

evaluate() [2/2]

ndarray::Array< double, 1, 1 > lsst::shapelet::RadialProfile::evaluate

(

ndarray::Array< double const, 1, 1 > const &

r

)

const

Evaluate the profile at the given radius (vectorized).

Radius is given in units of half-light radius, and the profile is normalized to 1 at r=1.

get()

static RadialProfile & lsst::shapelet::RadialProfile::get ( std::string const & name )

static

Return a predefined radial profile given its name.

New RadialProfile classes are registered when their constructor is invoked, so all subclasses should only be instantiated once (at module scope).

getBasis()

std::shared_ptr< MultiShapeletBasis > lsst::shapelet::RadialProfile::getBasis	(	int	nComponents,
		int	maxRadius = 0 ) const

Return a multi-Gaussian approximation to the radial profile.

Unlike the evaluate() method, the returned basis is normalized that the integrated flux is equal to the (single) amplitude.

Parameters

[in]	nComponents	Number of Gaussians used to approximate the profile.
[in]	maxRadius	Maximum radius used in fitting the approximation. Passing 0 selects the default for that profile.

Gaussian approximations are computed in advance and persisted, and are usually loaded and registered with the RadialProfile object in Python at module-import time.

Throws NotFoundError if the a basis with the given combination of nComponents and maxRadius has not been added to the RadialProfile.

getMomentsRadiusFactor()

double lsst::shapelet::RadialProfile::getMomentsRadiusFactor ( ) const

inline

Return the 2nd-moment radius in units of the half-light radius.

Definition at line 89 of file RadialProfile.h.

{ return _momentsRadiusFactor; }

getName()

std::string const lsst::shapelet::RadialProfile::getName ( ) const

inline

Return the name of the profile.

Definition at line 72 of file RadialProfile.h.

{ return _name; }

registerBasis()

void lsst::shapelet::RadialProfile::registerBasis	(	std::shared_ptr< MultiShapeletBasis >	basis,
		int	nComponents,
		int	maxRadius )

Register a basis with the profile, making it available to callers of getBasis().

Member Data Documentation

_momentsRadiusFactor

double lsst::shapelet::RadialProfile::_momentsRadiusFactor

protected

Definition at line 118 of file RadialProfile.h.

---

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

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-9.0.0/Linux64/shapelet/ge79ae78c31+9634bc57db/include/lsst/shapelet/RadialProfile.h