def test_load_reference_catalog(self):
refObjLoader, center, radius, filterName, fakeRefCat = self._make_fake_refcat()
config = lsst.jointcal.jointcal.JointcalConfig()
jointcal = lsst.jointcal.JointcalTask(config=config, butler=self.butler)
refCat, fluxField = jointcal._load_reference_catalog(refObjLoader,
jointcal.astrometryReferenceSelector,
center,
radius,
filterName)
# operator== isn't implemented for Catalogs, so we have to check like
# this, in case the records are copied during load.
self.assertEqual(len(refCat), len(fakeRefCat))
for r1, r2 in zip(refCat, fakeRefCat):
self.assertEqual(r1, r2)
def test_load_reference_catalog(self):
refObjLoader, center, radius, filterName, fakeRefCat = self._make_fake_refcat(
)
config = lsst.jointcal.jointcal.JointcalConfig()
config.astrometryReferenceErr = 0.1 # our test refcats don't have coord errors
jointcal = lsst.jointcal.JointcalTask(config=config,
butler=self.butler)
refCat, fluxField = jointcal._load_reference_catalog(
refObjLoader, jointcal.astrometryReferenceSelector, center, radius,
filterName)
# operator== isn't implemented for Catalogs, so we have to check like
# this, in case the records are copied during load.
self.assertEqual(len(refCat), len(fakeRefCat))
for r1, r2 in zip(refCat, fakeRefCat):
self.assertEqual(r1, r2)
def test_load_reference_catalog_subselect(self):
"""Test that we can select out the one source in the fake refcat
with a ridiculous S/N cut.
"""
refObjLoader, center, radius, filterName, fakeRefCat = self._make_fake_refcat()
config = lsst.jointcal.jointcal.JointcalConfig()
config.astrometryReferenceSelector.doSignalToNoise = True
config.astrometryReferenceSelector.signalToNoise.minimum = 1e10
config.astrometryReferenceSelector.signalToNoise.fluxField = "fake_flux"
config.astrometryReferenceSelector.signalToNoise.errField = "fake_fluxErr"
jointcal = lsst.jointcal.JointcalTask(config=config, butler=self.butler)
refCat, fluxField = jointcal._load_reference_catalog(refObjLoader,
jointcal.astrometryReferenceSelector,
center,
radius,
filterName)
self.assertEqual(len(refCat), 0)
def test_load_reference_catalog_subselect(self):
"""Test that we can select out the one source in the fake refcat
with a ridiculous S/N cut.
"""
refObjLoader, center, radius, filterName, fakeRefCat = self._make_fake_refcat(
)
config = lsst.jointcal.jointcal.JointcalConfig()
config.astrometryReferenceErr = 0.1 # our test refcats don't have coord errors
config.astrometryReferenceSelector.doSignalToNoise = True
config.astrometryReferenceSelector.signalToNoise.minimum = 1e10
config.astrometryReferenceSelector.signalToNoise.fluxField = "fake_flux"
config.astrometryReferenceSelector.signalToNoise.errField = "fake_fluxErr"
jointcal = lsst.jointcal.JointcalTask(config=config,
butler=self.butler)
refCat, fluxField = jointcal._load_reference_catalog(
refObjLoader, jointcal.astrometryReferenceSelector, center, radius,
filterName)
self.assertEqual(len(refCat), 0)