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

# 

# LSST Data Management System 

# Copyright 2017 LSST/AURA. 

# 

# 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 <http://www.lsstcorp.org/LegalNotices/>. 

# 

from __future__ import division, print_function, absolute_import 

 

import lsst.meas.base 

 

__all__ = ["FilterFractionPlugin"] 

 

 

class FilterFractionPluginConfig(lsst.meas.base.SingleFramePluginConfig): 

pass 

 

 

@lsst.meas.base.register("subaru_FilterFraction") 

class FilterFractionPlugin(lsst.meas.base.SingleFramePlugin): 

 

ConfigClass = FilterFractionPluginConfig 

 

@classmethod 

def getExecutionOrder(cls): 

return cls.FLUX_ORDER 

 

def __init__(self, config, name, schema, metadata): 

lsst.meas.base.SingleFramePlugin.__init__(self, config, name, schema, metadata) 

self.unweightedKey = schema.addField( 

schema.join(name, "unweighted"), type=float, 

doc=("Fraction of observations taken with the preferred version of this filter. " 

"For filters with a single version, this will always be 1; for HSC i or r, " 

"the preferred filter is the replacement filter (i2 or r2)."), 

doReplace=True 

) 

self.weightedKey = schema.addField( 

schema.join(name, "weighted"), type=float, 

doc=("Contribution-weighted fraction of observations taken with the preferred " 

"version of this filter. " 

"For filters with a single version, this will always be 1; for HSC i or r, " 

"the preferred filter is the replacement filter (i2 or r2)."), 

doReplace=True 

) 

 

def measure(self, measRecord, exposure): 

try: 

ccds = exposure.getInfo().getCoaddInputs().ccds 

except AttributeError: 

raise lsst.meas.base.FatalAlgorithmError("FilterFraction can only be run on coadds.") 

overlapping = ccds.subsetContaining(measRecord.getCentroid(), exposure.getWcs(), 

includeValidPolygon=True) 

67 ↛ 68line 67 didn't jump to line 68, because the condition on line 67 was never true if not overlapping: 

measRecord.set(self.key, float("NaN")) # no inputs in any filter 

return 

counts = {} 

weights = {} 

filterKey = overlapping.schema.find("filter").key 

try: 

weightKey = overlapping.schema.find("weight").key 

except KeyError: 

weightKey = None 

for ccd in overlapping: 

filterName = ccd.get(filterKey) 

counts[filterName] = counts.get(filterName, 0) + 1 

80 ↛ 77line 80 didn't jump to line 77, because the condition on line 80 was never false if weightKey is not None: 

weight = ccd.get(weightKey) 

weights[filterName] = weights.get(filterName, 0.0) + weight 

if "i" in counts: 

weird = set(counts.keys()) - set(["i", "i2"]) 

if weird: 

raise lsst.meas.base.FatalAlgorithmError( 

"Unexpected filter combination found in coadd: %s" % counts.keys() 

) 

measRecord.set(self.unweightedKey, counts.get("i2", 0.0)/len(overlapping)) 

90 ↛ exitline 90 didn't return from function 'measure', because the condition on line 90 was never false if weightKey is not None: 

measRecord.set(self.weightedKey, weights.get("i2", 0.0)/sum(weights.values())) 

elif "r" in counts: 

weird = set(counts.keys()) - set(["r", "r2"]) 

94 ↛ 95line 94 didn't jump to line 95, because the condition on line 94 was never true if weird: 

raise lsst.meas.base.FatalAlgorithmError( 

"Unexpected filter combination found in coadd: %s" % counts.keys() 

) 

measRecord.set(self.unweightedKey, counts.get("r2", 0.0)/len(overlapping)) 

99 ↛ exitline 99 didn't return from function 'measure', because the condition on line 99 was never false if weightKey is not None: 

measRecord.set(self.weightedKey, weights.get("r2", 0.0)/sum(weights.values())) 

101 ↛ 102line 101 didn't jump to line 102, because the condition on line 101 was never true elif len(counts) > 1: 

raise lsst.meas.base.FatalAlgorithmError( 

"Unexpected filter combination found in coadd: %s" % counts.keys() 

) 

else: 

measRecord.set(self.unweightedKey, 1.0) 

107 ↛ exitline 107 didn't return from function 'measure', because the condition on line 107 was never false if weightKey is not None: 

measRecord.set(self.weightedKey, 1.0) 

 

def fail(self, measRecord, error=None): 

# this plugin can only raise fatal exceptions 

pass