displayingStackDiffs

"""
Short script for displaying difference Image, from Simon's instructions

The basic philosophy is to *NOT* just look at the files using ds9, but rather
data products created by the DM stack in a repository which knows how to fetch
these data products.
"""
import lsst.afw.display
import lsst.daf.persistence as dp

# Display the differences produced by the difference imaging algorithm
# The argument is a data repository (which has a `_parent`)
butler = dp.Butler('Run3')
diff_exp = butler.get('deepDiff_differenceExp', dataId=dict(visit=250), immediate=True)
display.mtv(diff_exp)

# Find the difference image sources that are above a threshold and show their positions.
srcs = butler.get('deepDiff_diaSrc', dataId=dict(visit=250), immediate=True)
for src in srcs:
    display.dot('X', src.getX(), src.getY())