JointcalStatistics

class lsst.jointcal.utils.JointcalStatistics(match_radius=4.84814e-07 rad, flux_limit=100.0, do_photometry=True, do_astrometry=True, verbose=False)

Bases: object

Compute statistics on jointcal-processed data, and optionally generate plots.

Notes

Instantiate JointcalStatistics and call compute_rms() to get the relevant statistics for e.g. unittests, and call make_plots() to generate a suite of diagnostic plots.

Methods Summary

compute_rms(data_refs, reference)	Match all data_refs to compute the RMS, for all detections above self.flux_limit.
make_plots(data_refs, old_wcs_list[, name, …])	Make plots of various quantites to help with debugging.

Methods Documentation

compute_rms(data_refs, reference)

Match all data_refs to compute the RMS, for all detections above self.flux_limit.

Parameters:

data_refs : list of lsst.daf.persistence.butlerSubset.ButlerDataRef

A list of data refs to do the calculations between.

reference : lsst reference catalog

reference catalog to do absolute matching against.

Returns:

namedtuple:

astropy.Quantity

Post-jointcal relative RMS of the matched sources.

astropy.Quantity

Post-jointcal absolute RMS of matched sources.

float

Post-jointcal photometric repeatability (PA1 from the SRD).

make_plots(data_refs, old_wcs_list, name='', interactive=False, per_ccd_plot=False, outdir='.plots')

Make plots of various quantites to help with debugging. Requires that compute_rms() was run first.

Parameters:

data_refs : list of lsst.daf.persistence.butlerSubset.ButlerDataRef

A list of data refs to do the calculations between.

old_wcs_list : list of lsst.afw.image.wcs.Wcs

A list of the old (pre-jointcal) WCSs, one-to-one corresponding to data_refs.

name : str

Name to include in plot titles and save files.

interactive : bool

Turn on matplotlib interactive mode and drop into a debugger when plotting is finished. Otherwise, use a non-interactive backend.

per_ccd_plot : bool

Plot the WCS per CCD (takes longer and generates many plots for a large camera)

outdir : str

directory to save plots to.