Coverage for tests/test_defaultZeroPoint.py : 41%

Hot-keys on this page
r m x p toggle line displays
j k next/prev highlighted chunk
0 (zero) top of page
1 (one) first highlighted chunk
# # LSST Data Management System # # Copyright 2008-2016 AURA/LSST. # # This product includes software developed by the # LSST Project (http://www.lsst.org/). # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by # the Free Software Foundation, either version 3 of the License, or # (at your option) any later version. # # This program is distributed in the hope that it will be useful, # but WITHOUT ANY WARRANTY; without even the implied warranty of # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the # GNU General Public License for more details. # # You should have received a copy of the LSST License Statement and # the GNU General Public License along with this program. If not, # see <https://www.lsstcorp.org/LegalNotices/>. #
'''Test that photoCalib objects containing invalid calibrations are replaced with a default having a zero-point of magnitude 27. '''
'''Create two calibs, one with a valid zero-point and one without. Use these to create two UnitSystem objects. ''' self.mag2Flux = lambda m: 10.0**(m/2.5) self.flux2Mag = lambda f: 2.5*np.log10(f)
photoCalibNoZero = afwImage.PhotoCalib() photoCalibWithZero = afwImage.makePhotoCalibFromCalibZeroPoint(self.mag2Flux(25))
scale = 0.2 * geom.arcseconds wcs = afwGeom.makeSkyWcs(crpix=geom.Point2D(), crval=geom.SpherePoint(45.0, 45.0, geom.degrees), cdMatrix=afwGeom.makeCdMatrix(scale=scale))
self.unitNoZero = measModel.UnitSystem(wcs, photoCalibNoZero) self.unitWithZero = measModel.UnitSystem(wcs, photoCalibWithZero)
'''Verify that the UnitSystem with an invalid photoCalib gets populated with a valid photoCalib, and that nothing happens to the UnitSystem with a valid photoCalib. ''' magNoZero = self.flux2Mag(self.unitNoZero.photoCalib.getInstFluxAtZeroMagnitude()) magWithZero = self.flux2Mag(self.unitWithZero.photoCalib.getInstFluxAtZeroMagnitude())
print((magNoZero, magWithZero))
self.assertAlmostEqual(magNoZero, 27, 6) self.assertAlmostEqual(magWithZero, 25, 6)
del self.unitNoZero del self.unitWithZero
lsst.utils.tests.init()
lsst.utils.tests.init() unittest.main() |