Coverage for tests/test_fgcmcal_hsc.py : 15%
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
|
# See COPYRIGHT file at the top of the source tree. # # This file is part of fgcmcal. # # Developed for the LSST Data Management System. # This product includes software developed by the LSST Project # (https://www.lsst.org). # See the COPYRIGHT file at the top-level directory of this distribution # for details of code ownership. # # 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 GNU General Public License # along with this program. If not, see <https://www.gnu.org/licenses/>.
Run test suite on fgcmcal using HSC data from testdata_jointcal. """
def setUpClass(cls): try: cls.dataDir = lsst.utils.getPackageDir('testdata_jointcal') os.environ['ASTROMETRY_NET_DATA_DIR'] = os.path.join(cls.dataDir, 'hsc_and_index') except lsst.pex.exceptions.NotFoundError: raise unittest.SkipTest("testdata_jointcal not setup")
inputDir = os.path.join(self.dataDir, 'hsc')
testDir = tempfile.mkdtemp(dir=ROOT, prefix="TestFgcm-")
self.setUp_base(inputDir=inputDir, testDir=testDir)
lsst.log.setLevel("HscMapper", lsst.log.FATAL)
# Set numpy seed for stability np.random.seed(seed=1000)
visitDataRefName = 'visit' ccdDataRefName = 'ccd'
# First test making the LUT self.config = fgcmcal.FgcmMakeLutConfig() self.config.filterNames = ['r', 'i'] self.config.stdFilterNames = ['r', 'i'] self.config.atmosphereTableName = 'fgcm_atm_subaru2_test' self.otherArgs = []
nBand = 2 i0Std = np.array([0.07877351, 0.06464688]) i10Std = np.array([-0.00061516, -0.00063434]) i0Recon = np.array([0.0689530429, 0.05600673]) i10Recon = np.array([-7.01847144, 3.62675740])
self._testFgcmMakeLut(nBand, i0Std, i0Recon, i10Std, i10Recon)
# Now the star building self.config = fgcmcal.FgcmBuildStarsConfig() self.config.filterToBand = {'r': 'r', 'i': 'i'} self.config.requiredBands = ['r', 'i'] self.config.referenceBands = ['i'] self.config.checkAllCcds = True self.config.visitDataRefName = visitDataRefName self.config.ccdDataRefName = ccdDataRefName self.otherArgs = []
nVisit = 11 nStar = 472 nObs = 5431
self._testFgcmBuildStars(nVisit, nStar, nObs)
# And the fit cycle self.config = fgcmcal.FgcmFitCycleConfig() self.config.outfileBase = 'TestFgcm' self.config.bands = ['r', 'i'] self.config.fitFlag = (1, 1) self.config.requiredFlag = (1, 1) self.config.filterToBand = {'r': 'r', 'i': 'i'} self.config.maxIterBeforeFinalCycle = 1 self.config.nCore = 1 self.config.cycleNumber = 0 self.config.utBoundary = 0.0 self.config.washMjds = (0.0, ) self.config.epochMjds = (0.0, 100000.0) self.config.latitude = 19.8256 self.config.pixelScale = 0.17 self.config.expGrayPhotometricCut = (-0.05, -0.05) self.config.expGrayHighCut = (0.2, 0.2) self.config.aperCorrFitNBins = 0 self.config.sedFudgeFactors = (1.0, 1.0) self.config.starColorCuts = ('r,i,-0.50,2.25',) self.config.freezeStdAtmosphere = True self.config.precomputeSuperStarInitialCycle = False self.config.minExpPerNight = 1 self.config.minStarPerExp = 50 self.config.nStarPerRun = 50 self.config.nExpPerRun = 2 self.config.colorSplitIndices = (0, 1) self.config.superStarSubCcd = True self.config.superStarSubCcdChebyshevOrder = 1 self.config.ccdGraySubCcd = False self.config.ccdGraySubCcdChebyshevOrder = 1 self.config.ccdGraySubCcdTriangular = True self.config.modelMagErrors = False self.config.outputZeropointsBeforeFinalCycle = False self.config.outputStandardsBeforeFinalCycle = False self.config.sigmaCalRange = (0.003, 0.003) self.otherArgs = []
nZp = 1232 nGoodZp = 26 nOkZp = 32 nBadZp = 1200 nStdStars = 391 nPlots = 30
self._testFgcmFitCycle(nZp, nGoodZp, nOkZp, nBadZp, nStdStars, nPlots, skipChecks=True)
# Test the second fit cycle -- need to copy to unfreeze config newConfig = fgcmcal.FgcmFitCycleConfig() newConfig.update(**self.config.toDict()) newConfig.cycleNumber = 1 self.config = newConfig
self._testFgcmFitCycle(nZp, nGoodZp, nOkZp, nBadZp, nStdStars, nPlots, skipChecks=True)
# Test the "final" fit cycle newConfig = fgcmcal.FgcmFitCycleConfig() newConfig.update(**self.config.toDict()) newConfig.cycleNumber = 2 newConfig.ccdGraySubCcd = True newConfig.isFinalCycle = True self.config = newConfig
self._testFgcmFitCycle(nZp, nGoodZp, nOkZp, nBadZp, nStdStars, nPlots)
# And output the products self.config = fgcmcal.FgcmOutputProductsConfig() self.config.cycleNumber = 2 self.config.photoCal.photoCatName = 'sdss-dr9-fink-v5b' self.config.photoCal.colorterms.data = {} self.config.photoCal.colorterms.data['sdss*'] = lsst.pipe.tasks.colorterms.ColortermDict() self.config.photoCal.colorterms.data['sdss*'].data = {} self.config.photoCal.colorterms.data['sdss*'].data['g'] = lsst.pipe.tasks.colorterms.Colorterm() self.config.photoCal.colorterms.data['sdss*'].data['g'].primary = 'g' self.config.photoCal.colorterms.data['sdss*'].data['g'].secondary = 'r' self.config.photoCal.colorterms.data['sdss*'].data['g'].c0 = -0.00816446 self.config.photoCal.colorterms.data['sdss*'].data['g'].c1 = -0.08366937 self.config.photoCal.colorterms.data['sdss*'].data['g'].c2 = -0.00726883 self.config.photoCal.colorterms.data['sdss*'].data['r'] = lsst.pipe.tasks.colorterms.Colorterm() self.config.photoCal.colorterms.data['sdss*'].data['r'].primary = 'r' self.config.photoCal.colorterms.data['sdss*'].data['r'].secondary = 'i' self.config.photoCal.colorterms.data['sdss*'].data['r'].c0 = 0.0013181 self.config.photoCal.colorterms.data['sdss*'].data['r'].c1 = 0.01284177 self.config.photoCal.colorterms.data['sdss*'].data['r'].c2 = -0.03068248 self.config.photoCal.colorterms.data['sdss*'].data['i'] = lsst.pipe.tasks.colorterms.Colorterm() self.config.photoCal.colorterms.data['sdss*'].data['i'].primary = 'i' self.config.photoCal.colorterms.data['sdss*'].data['i'].secondary = 'z' self.config.photoCal.colorterms.data['sdss*'].data['i'].c0 = 0.00130204 self.config.photoCal.colorterms.data['sdss*'].data['i'].c1 = -0.16922042 self.config.photoCal.colorterms.data['sdss*'].data['i'].c2 = -0.01374245 self.config.refObjLoader.ref_dataset_name = "sdss-dr9-fink-v5b"
filterMapping = {'r': 'HSC-R', 'i': 'HSC-I'} # These zeropoint offsets are empirical, and are there # to check if changes in the code are altering the final # output in a measurable way. zpOffsets = np.array([-0.0223396513611, 0.266586363316])
self._testFgcmOutputProducts(visitDataRefName, ccdDataRefName, filterMapping, zpOffsets, 904014, 12, 'i', 1)
lsst.utils.tests.init()
lsst.utils.tests.init() unittest.main() |