CalcRhoStatistics¶

class lsst.analysis.tools.actions.vector.CalcRhoStatistics(*args, **kw)¶

Bases: KeyedDataAction

Calculate rho statistics.

Rho statistics refer to a collection of correlation functions involving PSF ellipticity and size residuals. They quantify the contribution from PSF leakage due to errors in PSF modeling to the weak lensing shear correlation functions.

The exact definitions of rho statistics as defined in [1] are given below.

\[ \begin{align}\begin{aligned}\rho_1(\theta) &= \left\langle \delta e^*_{PSF}(x) \delta e_{PSF}(x+\theta) \right\rangle\\\rho_2(\theta) &= \left\langle e^*_{PSF}(x) \delta e_{PSF}(x+\theta \right\rangle\\\rho_3(\theta) &= \left\langle (e^*_{PSF}\frac{\delta T_{PSF}}{T_{PSF}}(x)) (e_{PSF}\frac{\delta T_{PSF}}{T_{PSF}})(x+\theta) \right\rangle\\\rho_4(\theta) &= \left\langle \delta e^*_{PSF}(x) (e_{PSF}\frac{\delta T_{PSF}}{T_{PSF}})(x+\theta) \right\rangle\\\rho_5(\theta) &= \left\langle e^*_{PSF}(x) (e_{PSF}\frac{\delta T_{PSF}}{T_{PSF}})(x+\theta) \right\rangle\end{aligned}\end{align} \]

In addition to these five, we also compute the auto-correlation function of the fractional size residuals and call it as the \(\rho'_3( \theta )\), as referred to in Melchior et al. (2015) [2].

\[\rho'_3(\theta) = \left\langle\frac{\delta T_{PSF}}{T_{PSF}}(x) \frac{\delta T_{PSF}}{T_{PSF}}(x+\theta) \right\rangle\]

The definition of ellipticity used in [1] correspond to shear-type, which is typically smaller by a factor of 4 than using distortion-type.

References

[1] (1,2)

Jarvis, M., Sheldon, E., Zuntz, J., Kacprzak, T., Bridle, S. L., et. al (2016). The DES Science Verification weak lensing shear catalogues MNRAS, 460, 2245–2281. https://doi.org/10.1093/mnras/stw990; https://arxiv.org/abs/1507.05603

[2]

Melchior, P., et. al (2015) Mass and galaxy distributions of four massive galaxy clusters from Dark Energy Survey Science Verification data MNRAS, 449, no. 3, pp. 2219–2238. https://doi:10.1093/mnras/stv398 https://arxiv.org/abs/1405.4285

Attributes Summary

`applyContext`	Apply a `Context` to an `AnalysisAction` recursively.
`colDec`	Dec column (`str`, default `'coord_dec'`)
`colPsfXx`	The column name to get the PSF xx shape component from.
`colPsfXy`	The column name to get the PSF xy shape component from.
`colPsfYy`	The column name to get the PSF yy shape component from.
`colRa`	RA column (`str`, default `'coord_ra'`)
`colXx`	The column name to get the xx shape component from.
`colXy`	The column name to get the xy shape component from.
`colYy`	The column name to get the yy shape component from.
`ellipticityType`	The type of ellipticity to calculate (`str`, default `'distortion'`)
`history`
`identity`	If a configurable action is assigned to a `ConfigurableActionField`, or a `ConfigurableActionStructField` the name of the field will be bound to this variable when it is retrieved.
`sizeType`	The type of size to calculate (`str`, default `'trace'`)
`treecorr`	TreeCorr configuration (`TreecorrConfig`, default `<class 'lsst.meas.algorithms.treecorrUtils.TreecorrConfig'>`)

Methods Summary

`__call__`(data, **kwargs)	Call self as a function.
`addInputSchema`(inputSchema)	Add the supplied inputSchema argument to the class such that it will be returned along side any other arguments in a call to `getInputSchema`.
`compare`(other[, shortcut, rtol, atol, output])	Compare this configuration to another `Config` for equality.
`formatHistory`(name, **kwargs)	Format a configuration field's history to a human-readable format.
`freeze`()	Make this config, and all subconfigs, read-only.
`getFormattedInputSchema`(**kwargs)	Return input schema, with keys formatted with any arguments supplied by kwargs passed to this method.
`getInputSchema`()	Return the schema an `AnalysisAction` expects to be present in the arguments supplied to the __call__ method.
`getOutputSchema`()	Return the schema an `AnalysisAction` will produce, if the `__call__` method returns `KeyedData`, otherwise this may return None.
`items`()	Get configurations as `(field name, field value)` pairs.
`keys`()	Get field names.
`load`(filename[, root])	Modify this config in place by executing the Python code in a configuration file.
`loadFromStream`(stream[, root, filename, ...])	Modify this Config in place by executing the Python code in the provided stream.
`loadFromString`(code[, root, filename, ...])	Modify this Config in place by executing the Python code in the provided string.
`names`()	Get all the field names in the config, recursively.
`save`(filename[, root])	Save a Python script to the named file, which, when loaded, reproduces this config.
`saveToStream`(outfile[, root, skipImports])	Save a configuration file to a stream, which, when loaded, reproduces this config.
`saveToString`([skipImports])	Return the Python script form of this configuration as an executable string.
`setDefaults`()	Subclass hook for computing defaults.
`toDict`()	Make a dictionary of field names and their values.
`update`(**kw)	Update values of fields specified by the keyword arguments.
`validate`()	Validate the Config, raising an exception if invalid.
`values`()	Get field values.

Attributes Documentation

applyContext¶

Apply a Context to an AnalysisAction recursively.

Generally this method is called from within an AnalysisTool to configure all AnalysisActions at one time to make sure that they all are consistently configured. However, it is permitted to call this method if you are aware of the effects, or from within a specific execution environment like a python shell or notebook.

Parameters:

contextContext: The specific execution context, this may be a single context or a joint context, see Context for more info.

colDec¶: Dec column (str, default 'coord_dec')

colPsfXx¶: The column name to get the PSF xx shape component from. (str, default '{band}_ixxPSF')

colPsfXy¶: The column name to get the PSF xy shape component from. (str, default '{band}_ixyPSF')

colPsfYy¶: The column name to get the PSF yy shape component from. (str, default '{band}_iyyPSF')

colRa¶: RA column (str, default 'coord_ra')

colXx¶: The column name to get the xx shape component from. (str, default '{band}_ixx')

colXy¶: The column name to get the xy shape component from. (str, default '{band}_ixy')

colYy¶: The column name to get the yy shape component from. (str, default '{band}_iyy')

ellipticityType¶

The type of ellipticity to calculate (str, default 'distortion')

Allowed values:

'distortion': Distortion, measured as \((I_{xx}-I_{yy})/(I_{xx}+I_{yy})\)
'shear': Shear, measured as \((I_{xx}-I_{yy})/(I_{xx}+I_{yy}+2\sqrt{I_{xx}I_{yy}-I_{xy}^2})\)

history¶: Read-only history.

identity: str | None = None¶: If a configurable action is assigned to a ConfigurableActionField, or a ConfigurableActionStructField the name of the field will be bound to this variable when it is retrieved.

sizeType¶

The type of size to calculate (str, default 'trace')

Allowed values:

'trace': trace radius
'determinant': determinant radius
'None': Field is optional

treecorr¶: TreeCorr configuration (TreecorrConfig, default <class 'lsst.meas.algorithms.treecorrUtils.TreecorrConfig'>)

Methods Documentation

__call__(data: MutableMapping[str, ndarray[Any, dtype[_ScalarType_co]] | Scalar | HealSparseMap | Tensor | Mapping], **kwargs) → MutableMapping[str, ndarray[Any, dtype[_ScalarType_co]] | Scalar | HealSparseMap | Tensor | Mapping]¶: Call self as a function.

addInputSchema(inputSchema: Mapping]]]) → None¶

Add the supplied inputSchema argument to the class such that it will be returned along side any other arguments in a call to getInputSchema.

Parameters:

inputSchemaKeyedDataSchema: A schema that is to be merged in with any existing schema when a call to getInputSchema is made.

compare(other, shortcut=True, rtol=1e-08, atol=1e-08, output=None)¶

Compare this configuration to another Config for equality.

Parameters:

otherlsst.pex.config.Config: Other Config object to compare against this config.
shortcutbool, optional: If True, return as soon as an inequality is found. Default is True.
rtolfloat, optional: Relative tolerance for floating point comparisons.
atolfloat, optional: Absolute tolerance for floating point comparisons.
outputcallable, optional: A callable that takes a string, used (possibly repeatedly) to report inequalities.

Returns:

isEqualbool: True when the two lsst.pex.config.Config instances are equal. False if there is an inequality.

See also

lsst.pex.config.compareConfigs

Notes

Unselected targets of RegistryField fields and unselected choices of ConfigChoiceField fields are not considered by this method.

Floating point comparisons are performed by numpy.allclose.

formatHistory(name, **kwargs)¶

Format a configuration field’s history to a human-readable format.

Parameters:

namestr: Name of a Field in this config.
**kwargs: Keyword arguments passed to lsst.pex.config.history.format.

Returns:

historystr: A string containing the formatted history.

See also

lsst.pex.config.history.format

freeze()¶: Make this config, and all subconfigs, read-only.

getFormattedInputSchema(**kwargs) → Mapping]]]¶

Return input schema, with keys formatted with any arguments supplied by kwargs passed to this method.

Returns:

resultKeyedDataSchema: The schema this action requires to be present when calling this action, formatted with any input arguments (e.g. band=’i’)

getInputSchema() → KeyedDataSchema¶

Return the schema an AnalysisAction expects to be present in the arguments supplied to the __call__ method.

Returns:

resultKeyedDataSchema: The schema this action requires to be present when calling this action, keys are unformatted.

getOutputSchema() → Mapping]]] | None¶

Return the schema an AnalysisAction will produce, if the __call__ method returns KeyedData, otherwise this may return None.

Returns:

resultKeyedDataSchema or None: The schema this action will produce when returning from call. This will be unformatted if any templates are present. Should return None if action does not return KeyedData.

items()¶

Get configurations as (field name, field value) pairs.

Returns:

itemsItemsView

Iterator of tuples for each configuration. Tuple items are:

Field name.
Field value.

keys()¶

Get field names.

Returns:

namesKeysView: List of lsst.pex.config.Field names.

load(filename, root='config')¶

Modify this config in place by executing the Python code in a configuration file.

Parameters:

filenamestr

Name of the configuration file. A configuration file is Python module.

rootstr, optional

Name of the variable in file that refers to the config being overridden.

For example, the value of root is "config" and the file contains:

config.myField = 5

Then this config’s field myField is set to 5.

See also

lsst.pex.config.Config.loadFromStream
lsst.pex.config.Config.loadFromString
lsst.pex.config.Config.save
lsst.pex.config.Config.saveToStream
lsst.pex.config.Config.saveToString

loadFromStream(stream, root='config', filename=None, extraLocals=None)¶

Modify this Config in place by executing the Python code in the provided stream.

Parameters:

streamfile-like object, str, bytes, or CodeType

Stream containing configuration override code. If this is a code object, it should be compiled with mode="exec".

rootstr, optional

Name of the variable in file that refers to the config being overridden.

For example, the value of root is "config" and the file contains:

config.myField = 5

Then this config’s field myField is set to 5.

filenamestr, optional

Name of the configuration file, or None if unknown or contained in the stream. Used for error reporting.

extraLocalsdict of str to object, optional

Any extra variables to include in local scope when loading.

See also

lsst.pex.config.Config.load
lsst.pex.config.Config.loadFromString
lsst.pex.config.Config.save
lsst.pex.config.Config.saveToStream
lsst.pex.config.Config.saveToString

Notes

For backwards compatibility reasons, this method accepts strings, bytes and code objects as well as file-like objects. New code should use loadFromString instead for most of these types.

loadFromString(code, root='config', filename=None, extraLocals=None)¶

Modify this Config in place by executing the Python code in the provided string.

Parameters:

codestr, bytes, or CodeType

Stream containing configuration override code.

rootstr, optional

Name of the variable in file that refers to the config being overridden.

For example, the value of root is "config" and the file contains:

config.myField = 5

Then this config’s field myField is set to 5.

filenamestr, optional

Name of the configuration file, or None if unknown or contained in the stream. Used for error reporting.

extraLocalsdict of str to object, optional

Any extra variables to include in local scope when loading.

Raises:

ValueError: Raised if a key in extraLocals is the same value as the value of the root argument.

See also

lsst.pex.config.Config.load
lsst.pex.config.Config.loadFromStream
lsst.pex.config.Config.save
lsst.pex.config.Config.saveToStream
lsst.pex.config.Config.saveToString

names()¶

Get all the field names in the config, recursively.

Returns:

nameslist of str: Field names.

save(filename, root='config')¶

Save a Python script to the named file, which, when loaded, reproduces this config.

Parameters:

filenamestr: Desination filename of this configuration.
rootstr, optional: Name to use for the root config variable. The same value must be used when loading (see lsst.pex.config.Config.load).

See also

lsst.pex.config.Config.saveToStream
lsst.pex.config.Config.saveToString
lsst.pex.config.Config.load
lsst.pex.config.Config.loadFromStream
lsst.pex.config.Config.loadFromString

saveToStream(outfile, root='config', skipImports=False)¶

Save a configuration file to a stream, which, when loaded, reproduces this config.

Parameters:

outfilefile-like object: Destination file object write the config into. Accepts strings not bytes.
rootstr, optional: Name to use for the root config variable. The same value must be used when loading (see lsst.pex.config.Config.load).
skipImportsbool, optional: If True then do not include import statements in output, this is to support human-oriented output from pipetask where additional clutter is not useful.

See also

lsst.pex.config.Config.save
lsst.pex.config.Config.saveToString
lsst.pex.config.Config.load
lsst.pex.config.Config.loadFromStream
lsst.pex.config.Config.loadFromString

saveToString(skipImports=False)¶

Return the Python script form of this configuration as an executable string.

Parameters:

skipImportsbool, optional: If True then do not include import statements in output, this is to support human-oriented output from pipetask where additional clutter is not useful.

Returns:

codestr: A code string readable by loadFromString.

See also

lsst.pex.config.Config.save
lsst.pex.config.Config.saveToStream
lsst.pex.config.Config.load
lsst.pex.config.Config.loadFromStream
lsst.pex.config.Config.loadFromString

setDefaults()¶

Subclass hook for computing defaults.

Notes

Derived Config classes that must compute defaults rather than using the Field instances’s defaults should do so here. To correctly use inherited defaults, implementations of setDefaults must call their base class’s setDefaults.

toDict()¶

Make a dictionary of field names and their values.

Returns:

dict_dict: Dictionary with keys that are Field names. Values are Field values.

See also

lsst.pex.config.Field.toDict

Notes

This method uses the toDict method of individual fields. Subclasses of Field may need to implement a toDict method for this method to work.

update(**kw)¶

Update values of fields specified by the keyword arguments.

Parameters:

**kw: Keywords are configuration field names. Values are configuration field values.

Notes

The __at and __label keyword arguments are special internal keywords. They are used to strip out any internal steps from the history tracebacks of the config. Do not modify these keywords to subvert a Config instance’s history.

Examples

This is a config with three fields:

>>> from lsst.pex.config import Config, Field
>>> class DemoConfig(Config):
...     fieldA = Field(doc='Field A', dtype=int, default=42)
...     fieldB = Field(doc='Field B', dtype=bool, default=True)
...     fieldC = Field(doc='Field C', dtype=str, default='Hello world')
...
>>> config = DemoConfig()

These are the default values of each field:

>>> for name, value in config.iteritems():
...     print(f"{name}: {value}")
...
fieldA: 42
fieldB: True
fieldC: 'Hello world'

Using this method to update fieldA and fieldC:

>>> config.update(fieldA=13, fieldC='Updated!')

Now the values of each field are:

>>> for name, value in config.iteritems():
...     print(f"{name}: {value}")
...
fieldA: 13
fieldB: True
fieldC: 'Updated!'

validate()¶

Validate the Config, raising an exception if invalid.

Raises:

lsst.pex.config.FieldValidationError: Raised if verification fails.

Notes

The base class implementation performs type checks on all fields by calling their validate methods.

Complex single-field validation can be defined by deriving new Field types. For convenience, some derived lsst.pex.config.Field-types (ConfigField and ConfigChoiceField) are defined in lsst.pex.config that handle recursing into subconfigs.

Inter-field relationships should only be checked in derived Config classes after calling this method, and base validation is complete.

values()¶

Get field values.

Returns:

valuesValuesView: Iterator of field values.

Navigation

CalcRhoStatistics¶