This is a first attempt at a list of the LSST SQuaRE metrics. These metrics are for quality checks of both the data and the software.
Feel free to comment or suggest other metrics!
Alert Production (AP) Metrics
We need QA metrics for the single-frame real-time processing (AP) of the images to make sure that the data isn't useless. Here are some ideas:
- Background/sky level: Don't want it too high, maybe we are looking at the moon, or it's too early to be observing.
- Structure in background: If there's a lot of structure in the background level across a chip this can be caused by clouds or electronics problems (at least on DECam chip has this problem).
- Seeing: Measure the seeing, i.e. the PSF FWHM, if it's terrible then maybe stop the rest of the processing. If it's 6 arc sec seeing then we probably won't get useful difference images or transient detection.
- Check for transparency/clouds: use the magnitude zeropoint from known stars, if there are clouds then maybe we should stop the processing. Also check how much the zero point is varying across the field.
- Other basic QA checks to make sure that something bad hasn't happened to the telescope/camera system, i.e. garbage/noise in the images due to an electronics problem. "Health" checks.
- Tracking problems: look for streaks of stars in the images.
- Automated check for "linear" features in the image that could be satellites, airplanes, meteors, etc.
- Check for residuals in the difference images around bright, isolated stars that might indicate a problem with the difference images
- Check the PSF FWHM of the image and the template to make sure they are actually similar
We need QA metrics for the calibration products (e.g., flats, biases, etc.) to make sure that the individual ones are good enough and also track how they are changing with time.
Individual calibration frame checks:
- Domeflats: Need decent S/N, maybe 5000 minimum and shouldn't be too high to get near saturation or nonlinear regime. Also, make sure that the median signal for each chip isn't dramatically different from the last few times the dome flats were taken.
- Bias: Make sure it's consistent with biases from the last few days.
- Darks: Make sure they are long enough to detect the dark current. Also make sure that they frames are consistent with the darks from the last few days.
Track with time:
- Number of bad pixels per chip
- Bias level in each chip
- Median domeflat counts per chip/filter
- Median dark count rate per chip
- Readnoise for each amplifier
- Gain per amplifier?
Data Release Products (DRP) Metrics
Instrument Signature Removal (ISR)
- Cross-talk: Double-check that cross-talk has actually been removed. See if there is still a signature of bright stars in one chip in other chips.
- Bias correction:
- Mask bad pixels:
- Check saturation and bleed problems:
- Check the noise in background of the image. Is it similar to expectations
- Issues from flat fielding / illumination corrections
- Issues from fringe corrections: Check if there are any fringe patterns left in the image, maybe use fourier transform.
- Check the background level across the field. If gain corrections are applied then the background level should be pretty smooth once ISR is finished.
- If there's a sky/background subtraction stage (i.e. setting the background close to zero), then check that this hasn't created artifacts
- Check for leftover cosmic rays that got through the "repair" step.
Point Spread Function (PSF)
- Check seeing, is it consistent across the chip/focal plane. It should be. If not, then there might be a problem.
- Check shape of PSF, is it elliptical, weird, double/bimodal? If the telescope is out of focus then you can get bimodal PSFs or even donuts.
- Check that enough the programs are allowing enough spatial variation across chip/focal plane. If the problem is saying the PSF needs to be fixed across a chip and it's actually varying then this is a problem.
- Check for brighter-fatter effects, i.e. the PSF FWHM as a function of magnitude.
World Coordinate System (WCS)
- Check the astrometric scatter around the WCS fit for reference stars. See how this varies across the focal plane. If it varies a lot then there's a problem. For example, large scatter at the edges might indicate that there's a problem with the distortion terms.
- Check for lots of sources detected around bright, saturated stars.
- Check that the aperture corrections have the right sign, and see how they vary across the field. Should be pretty uniform.
- Check the S/N in the background and the stars, is it what you expect
- Check that regions around bright/saturated stars are okay
Multi-fit / Forced Photometry
- Check depth of the photometry, is it as deep as expected
- Check resulting chi/sharp or equivalent PSF/shape metrics. Are the median/sigma values as expected? How do they vary with magnitude?
- Check aperture corrections
- Check calibrated photometry against our master/reference LSST catalog
- Is it as expected?
- Any color, airmass or magnitude trends?
- Is the scatter wrt the reference catalog as expected and similar to the photometry noise?
- Are we getting the S/N that we should be getting according to the model?
- At the catalog level we can:
- Compare stack-reduced image to a reference catalog (e.g. SDSS Stripe82)
- Compare reductions of images with the stack and other software packages (e.g. DAOPHOT, DoPHOT, SExtractor, etc.)
- Compare stack reductions of simulated images to truth tables.
- Check the calibration by measuring the scatter in the stellar locus across a field. Compare this to the stellar locus scatter in the instrumental photometry. This isn't an absolute calibration check, but it's a check that the calibration terms are actually improving the relative calibration over a larger field.
- Issues due to differential refraction not being dealt with properly
- Issues in crowded regions?
- Resampling/warping problems?
- Calibrating/rescaling of image problems