In case anyone may find it useful, we’ve written a driver for forced photometry. It’s nothing special - it basically wraps ForcedPhotCcdTask
in a BatchParallelTask, which enables forced photometry to be pooled out to multiple cores. Like ForcedPhotCcdTask
, it performs forced photometry on calexp images using the catalogues produced by multiBandDriver on coadds as references. It can be found in our obs package here: https://github.com/GOTO-OBS/obs_goto/tree/master/python/lsst/obs/goto
It’s in two parts: forcedPhotDriver.py
which calls forcedPhotCalexp.py
. The latter is a small modification of forcedPhotCcd.py
which identifies the tract covered by the calexp.
Comments welcome.
J