Hide keyboard shortcuts

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

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

# 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/>. 

"""Test the fgcmcal fgcmLoadReferenceCatalog code with testdata_jointcal/hsc. 

 

""" 

 

import unittest 

import os 

import numpy as np 

import healpy as hp 

import esutil 

 

import lsst.utils 

import lsst.pipe.tasks 

import lsst.daf.persistence as dafPersist 

 

import lsst.fgcmcal as fgcmcal 

 

 

class FgcmLoadReferenceTestHSC(lsst.utils.tests.TestCase): 

@classmethod 

def setUpClass(cls): 

try: 

cls.dataDir = lsst.utils.getPackageDir('testdata_jointcal') 

except LookupError: 

raise unittest.SkipTest("testdata_jointcal not setup") 

 

def setUp(self): 

self.inputDir = os.path.join(self.dataDir, 'hsc') 

 

lsst.log.setLevel("HscMapper", lsst.log.FATAL) 

 

def test_fgcmLoadReference(self): 

""" 

Test loading of the fgcm reference catalogs. 

""" 

 

filterList = ['r', 'i'] 

 

config = fgcmcal.FgcmLoadReferenceCatalogConfig() 

config.applyColorTerms = True 

config.refObjLoader.ref_dataset_name = 'sdss-dr9-fink-v5b' 

config.colorterms.data = {} 

config.colorterms.data['sdss*'] = lsst.pipe.tasks.colorterms.ColortermDict() 

config.colorterms.data['sdss*'].data = {} 

config.colorterms.data['sdss*'].data['g'] = lsst.pipe.tasks.colorterms.Colorterm() 

config.colorterms.data['sdss*'].data['g'].primary = 'g' 

config.colorterms.data['sdss*'].data['g'].secondary = 'r' 

config.colorterms.data['sdss*'].data['g'].c0 = -0.00816446 

config.colorterms.data['sdss*'].data['g'].c1 = -0.08366937 

config.colorterms.data['sdss*'].data['g'].c2 = -0.00726883 

config.colorterms.data['sdss*'].data['r'] = lsst.pipe.tasks.colorterms.Colorterm() 

config.colorterms.data['sdss*'].data['r'].primary = 'r' 

config.colorterms.data['sdss*'].data['r'].secondary = 'i' 

config.colorterms.data['sdss*'].data['r'].c0 = 0.0013181 

config.colorterms.data['sdss*'].data['r'].c1 = 0.01284177 

config.colorterms.data['sdss*'].data['r'].c2 = -0.03068248 

config.colorterms.data['sdss*'].data['i'] = lsst.pipe.tasks.colorterms.Colorterm() 

config.colorterms.data['sdss*'].data['i'].primary = 'i' 

config.colorterms.data['sdss*'].data['i'].secondary = 'z' 

config.colorterms.data['sdss*'].data['i'].c0 = 0.00130204 

config.colorterms.data['sdss*'].data['i'].c1 = -0.16922042 

config.colorterms.data['sdss*'].data['i'].c2 = -0.01374245 

 

butler = dafPersist.Butler(self.inputDir) 

loadCat = fgcmcal.FgcmLoadReferenceCatalogTask(butler, config=config) 

 

ra = 320.0 

dec = 0.0 

rad = 0.1 

 

refCat = loadCat.getFgcmReferenceStarsSkyCircle(ra, dec, rad, filterList) 

 

# Check the number of mags and ranges 

self.assertEqual(len(filterList), refCat['refMag'].shape[1]) 

self.assertEqual(len(filterList), refCat['refMagErr'].shape[1]) 

self.assertLess(np.max(refCat['refMag'][:, 0]), 99.1) 

self.assertLess(np.max(refCat['refMag'][:, 1]), 99.1) 

self.assertLess(np.max(refCat['refMagErr'][:, 0]), 99.1) 

self.assertLess(np.max(refCat['refMagErr'][:, 1]), 99.1) 

 

# Check the separations from the center 

self.assertLess(np.max(esutil.coords.sphdist(ra, dec, refCat['ra'], refCat['dec'])), rad) 

 

# And load a healpixel 

nside = 256 

pixel = 393614 

 

refCat = loadCat.getFgcmReferenceStarsHealpix(nside, pixel, filterList) 

 

ipring = hp.ang2pix(nside, np.radians(90.0 - refCat['dec']), np.radians(refCat['ra'])) 

self.assertEqual(pixel, np.max(ipring)) 

self.assertEqual(pixel, np.min(ipring)) 

 

 

class TestMemory(lsst.utils.tests.MemoryTestCase): 

pass 

 

 

def setup_module(module): 

lsst.utils.tests.init() 

 

 

124 ↛ 125line 124 didn't jump to line 125, because the condition on line 124 was never trueif __name__ == "__main__": 

lsst.utils.tests.init() 

unittest.main()