DipoleFitPlugin¶
- class lsst.ip.diffim.DipoleFitPlugin(config, name, schema, metadata, logName=None)¶
Bases:
SingleFramePlugin
A single frame measurement plugin that fits dipoles to all merged (two-peak)
diaSources
.This measurement plugin accepts up to three input images in its
measure
method. If these are provided, it includes data from the pre-subtraction posImage (science image) and optionally negImage (template image) to constrain the fit. The meat of the fitting routines are in the classDipoleFitAlgorithm
.Notes
The motivation behind this plugin and the necessity for including more than one exposure are documented in DMTN-007 (http://dmtn-007.lsst.io).
This class is named
ip_diffim_DipoleFit
so that it may be used alongside the existingip_diffim_DipoleMeasurement
classes until such a time as those are deemed to be replaceable by this.Attributes Summary
Order for algorithms which require shape, centroid and flux (
float
).Order for algorithms which require only Footprint and Peaks (
float
).Order for catalog calculation plugins.
Order for algorithms which require a shape and a centroid (
float
).Order for algorithms which require a centroid (
float
).Registry of subclasses of
SingleFramePlugin
(PluginRegistry
).Methods Summary
doClassify
(measRecord, chi2val)Classify a source as a dipole.
fail
(measRecord[, error])Catch failures and set the correct flags.
This algorithm simultaneously fits the centroid and flux, and does not require any previous centroid fit.
Get the measurement transformation appropriate to this plugin.
measure
(measRecord, exposure)Measure the properties of a source on a single image.
measureDipoles
(measRecord, exposure[, ...])Perform the non-linear least squares minimization on the putative dipole source.
measureN
(measCat, exposure)Measure the properties of blended sources on a single image.
Attributes Documentation
- APCORR_ORDER = 3.0¶
Order for algorithms which require shape, centroid and flux (
float
).Notes
These algorithms may assume that
getCentroid
andgetShape
will return good values, that flux has been measured, and that and that a Footprint and its Peaks are available.
- CENTROID_ORDER = 0.0¶
Order for algorithms which require only Footprint and Peaks (
float
).Notes
Algorithms with this execution order include centroids.
- DEFAULT_CATALOGCALCULATION = 4.0¶
Order for catalog calculation plugins.
Notes
These plugins only operate on catalogs; they may not access pixel values.
- FAILURE_EDGE = 1¶
- FAILURE_FIT = 2¶
- FAILURE_NOT_DIPOLE = 4¶
- FAILURE_TOO_LARGE = 8¶
- FLUX_ORDER = 2.0¶
Order for algorithms which require a shape and a centroid (
float
).Notes
These algorithms may assume that both
getCentroid
andgetShape
will return good values, and that a Footprint and its Peaks are available.
- SHAPE_ORDER = 1.0¶
Order for algorithms which require a centroid (
float
).Notes
These algorithms may refer assume that
getCentroid
will return a good centroid, and that a Footprint and its Peaks are available.
- registry = <lsst.meas.base.pluginRegistry.PluginRegistry object>¶
Registry of subclasses of
SingleFramePlugin
(PluginRegistry
).
Methods Documentation
- doClassify(measRecord, chi2val)¶
Classify a source as a dipole.
- Parameters:
- measRecordTODO: DM-17458
TODO: DM-17458
- chi2valTODO: DM-17458
TODO: DM-17458
Notes
Sources are classified as dipoles, or not, according to three criteria:
Does the total signal-to-noise surpass the
minSn
?Are the pos/neg fluxes greater than 1.0 and no more than 0.65 (
maxFluxRatio
) of the total flux? By default this will never happen sinceposFlux == negFlux
.Is it a good fit (
chi2dof
< 1)? (Currently not used.)
- fail(measRecord, error=None)¶
Catch failures and set the correct flags.
Fallback on the current slot centroid positions, but set the dipole failure flag, since we attempted to fit the source.
- classmethod getExecutionOrder()¶
This algorithm simultaneously fits the centroid and flux, and does not require any previous centroid fit.
- getLogName()¶
- static getTransformClass()¶
Get the measurement transformation appropriate to this plugin.
This returns a subclass of
transforms.MeasurementTransform
, which may be instantiated with details of the algorithm configuration and then called with information about calibration and WCS to convert from raw measurement quantities to calibrated units. Calibrated data is then provided in a separate output table.Notes
By default, we copy everything from the input to the output without transformation.
- measure(measRecord, exposure)¶
Measure the properties of a source on a single image.
The image may be from a single epoch, or it may be a coadd.
- Parameters:
- measRecord
lsst.afw.table.SourceRecord
Record describing the object being measured. Previously-measured quantities may be retrieved from here, and it will be updated in-place tih the outputs of this plugin.
- exposure
lsst.afw.image.ExposureF
The pixel data to be measured, together with the associated PSF, WCS, etc. All other sources in the image should have been replaced by noise according to deblender outputs.
- measRecord
- measureDipoles(measRecord, exposure, posExp=None, negExp=None)¶
Perform the non-linear least squares minimization on the putative dipole source.
- Parameters:
- measRecord
lsst.afw.table.SourceRecord
diaSources that will be measured using dipole measurement
- exposure
lsst.afw.image.Exposure
Difference exposure on which the diaSources were detected;
exposure = posExp-negExp
If bothposExp
andnegExp
areNone
, will attempt to fit the dipole to just theexposure
with no constraint.- posExp
lsst.afw.image.Exposure
, optional “Positive” exposure, typically a science exposure, or None if unavailable When
posExp
isNone
, will computeposImage = exposure + negExp
.- negExp
lsst.afw.image.Exposure
, optional “Negative” exposure, typically a template exposure, or None if unavailable When
negExp
isNone
, will computenegImage = posExp - exposure
.
- measRecord
Notes
The main functionality of this routine was placed outside of this plugin (into
DipoleFitAlgorithm.fitDipole()
) so thatDipoleFitAlgorithm.fitDipole()
can be called separately for testing (@seetests/testDipoleFitter.py
)
- measureN(measCat, exposure)¶
Measure the properties of blended sources on a single image.
This operates on all members of a blend family at once. The image may be from a single epoch, or it may be a coadd.
- Parameters:
- measCat
lsst.afw.table.SourceCatalog
Catalog describing the objects (and only those objects) being measured. Previously-measured quantities will be retrieved from here, and it will be updated in-place with the outputs of this plugin.
- exposure
lsst.afw.image.ExposureF
The pixel data to be measured, together with the associated PSF, WCS, etc. All other sources in the image should have been replaced by noise according to deblender outputs.
- measCat
Notes
Derived classes that do not implement
measureN
should just inherit this disabled version. Derived classes that do implementmeasureN
should additionally add a bool doMeasureN config field to their config class to signal that measureN-mode is available.