Details
-
Type:
Story
-
Status: To Do
-
Resolution: Unresolved
-
Fix Version/s: None
-
Component/s: ap_verify
-
Labels:None
-
Story Points:2
-
Epic Link:
-
Team:Alert Production
-
Urgent?:No
Description
On DM-34068, I noticed that the prep_butler code wasn't loading the refcats shards I thought it should be, based on overplotting the refcat sources on the detector bbox. This might be because the ap_verify_ci_hits2015 dataset doesn't actually include those extra shards, because it didn't use the correct position for the center of the observation. This ticket is to investigate this discrepancy by making some plots of the refcat shards vs. the full focal plane, centered at the exposure boresight, and comparing with the center used in the refcat loading script.
Some plotting code that I was using when investigating
DM-34068, so I don't lose it:xx=raw.wcs.pixelToSky(2100, 4100)plt.scatter(xx.getRa(), xx.getDec(), color='black')xx=raw.wcs.pixelToSky(0, 4100)plt.scatter(xx.getRa(), xx.getDec(), color='black')xx=raw.wcs.pixelToSky(2100, 0)plt.scatter(xx.getRa(), xx.getDec(), color='black')xx=raw.wcs.pixelToSky(0, 0)plt.scatter(xx.getRa(), xx.getDec(), color='black')def plot_shard(id):data=lsst.afw.table.SimpleCatalog.readFits(f'refcats/gen2/gaia/gaia_{id}_refcats_gen2.fits')plt.scatter(data['coord_ra'], data['coord_dec'], s=1.5, label=str(id))for i in [1, 2, 3, 4, 5, 7]:plot_shard(157400+i)plt.legend()