lsst.jointcal 22.0.1-26-g6000b2f+918e6122b5
|
Public Member Functions | |
def | __init__ (self, match_radius=0.1 *arcseconds, flux_limit=100.0, do_photometry=True, do_astrometry=True, verbose=False) |
def | compute_rms (self, data_refs, reference) |
def | make_plots (self, data_refs, old_wcs_list, name='', interactive=False, per_ccd_plot=False, outdir='.plots') |
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.
def lsst.jointcal.utils.JointcalStatistics.__init__ | ( | self, | |
match_radius = 0.1*arcseconds , |
|||
flux_limit = 100.0 , |
|||
do_photometry = True , |
|||
do_astrometry = True , |
|||
verbose = False |
|||
) |
Parameters ---------- match_radius : lsst.geom.Angle match sources within this radius for RMS statistics flux_limit : float Signal/Noise (flux/fluxErr) for sources to be included in the RMS cross-match. 100 is a balance between good centroids and enough sources. do_photometry : bool, optional Perform calculations/make plots for photometric metrics. do_astrometry : bool, optional Perform calculations/make plots for astrometric metrics. verbose : bool, optional Print extra things
def lsst.jointcal.utils.JointcalStatistics.compute_rms | ( | self, | |
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).
def lsst.jointcal.utils.JointcalStatistics.make_plots | ( | self, | |
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.