Coverage for tests/test_deferredCharge.py: 24%

2# LSST Data Management System

6# This product includes software developed by the

7# LSST Project (http://www.lsst.org/).

9# This program is free software: you can redistribute it and/or modify

10# it under the terms of the GNU General Public License as published by

11# the Free Software Foundation, either version 3 of the License, or

12# (at your option) any later version.

13#

14# This program is distributed in the hope that it will be useful,

15# but WITHOUT ANY WARRANTY; without even the implied warranty of

16# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the

17# GNU General Public License for more details.

18#

19# You should have received a copy of the LSST License Statement and

20# the GNU General Public License along with this program. If not,

21# see <https://www.lsstcorp.org/LegalNotices/>.

22#

23"""Test cases for lsst.cp.pipe.deferredCharge."""

25import unittest

26import numpy

28import lsst.utils

29import lsst.utils.tests

31import lsst.cp.pipe as cpPipe

32from lsst.ip.isr import IsrMock, PhotonTransferCurveDataset

35class CpCtiSolveTaskTestCase(lsst.utils.tests.TestCase):

36 """A test case for the deferred charge/CTI solve task."""

38 def setUp(self):

39 self.camera = IsrMock().getCamera()

41 serialOscanMeansA = [

42 7.18039751e-01, 4.56550479e-01, 4.14261669e-01, 2.88099229e-01, -2.34310962e-02,

43 -4.59854975e-02, -1.14491098e-02, 5.19846082e-02, 2.05635265e-01, 1.25147207e-02,

44 9.00449380e-02, -2.39106059e-01, -1.52413145e-01, 4.63459678e-02, 1.85195580e-01,

45 -1.58051759e-01, -8.76842241e-04, -5.09192124e-02, 2.58496821e-01, 2.54267782e-01,

46 -1.37611866e-01, 3.35322201e-01, 1.04846083e-01, -2.16551736e-01, -8.82746354e-02,

47 -1.00256450e-01, 2.73297966e-01, -4.52805981e-02, 3.40960979e-01, 7.80628920e-02,

48 -2.90697180e-02, -6.99492991e-02, -1.06599867e-01, 6.89002723e-02, 1.46290688e-02,

49 1.19647197e-01, -1.54527843e-01, 9.35689881e-02, -1.06599934e-01, -2.13166289e-02,

50 9.35688764e-02, -1.19286761e-01, 1.18098985e-02, -3.69616691e-03, -6.14914447e-02,

51 -5.81059000e-03, 9.42736641e-02, 3.92978266e-02, -1.55937240e-01, 3.76202404e-01,

52 -1.13648064e-01, 1.71803936e-01, 6.17138995e-03, 8.22918862e-02, -2.84214199e-01,

53 -2.99097435e-03, -1.31973490e-01, -2.84214795e-01, -2.99140741e-03, -3.76546055e-01,

54 5.97376414e-02, -1.91883057e-01, -1.34087920e-01, -3.23684871e-01,

55 ]

56 serialOscanMeansB = [

57 1.50365152e+01, 4.43178511e+00, 2.66550946e+00, 1.67382801e+00, 1.10997069e+00,

58 8.89361799e-01, 4.66469795e-01, 6.10956728e-01, 6.67343795e-01, 5.22854805e-01,

59 -1.15006611e-01, 2.67710119e-01, 2.05686077e-01, 1.84541523e-01, 8.65717679e-02,

60 5.51738311e-03, 2.35288814e-01, 3.45944524e-01, 7.81139359e-02, 1.52119964e-01,

61 2.02162191e-01, 3.44150960e-02, -2.86277920e-01, 1.43662184e-01, 3.21276844e-01,

62 -6.21452965e-02, 8.58670697e-02, -1.63320359e-02, -1.07958235e-01, -1.60820082e-01,

63 -2.19705645e-02, -1.55181482e-01, -2.39055425e-01, -2.75705636e-01, 6.33126274e-02,

64 -5.50971478e-02, -2.42579415e-01, -9.87957790e-02, 1.08421087e-01, -1.12892322e-01,

65 1.89090632e-02, -1.53086300e-03, -2.18615308e-01, -2.19320312e-01, 9.22102109e-02,

66 -4.87535410e-02, -1.81964979e-01, -4.17055413e-02, -4.24422681e-01, -1.96061105e-01,

67 -1.35127297e-02, -1.77031055e-01, -2.30597332e-01, -4.01868790e-01, -4.18784261e-01,

68 -3.75085384e-01, -3.49007100e-01, -1.77735761e-01, -7.41272718e-02, -1.92537069e-01,

69 2.46565759e-01, -3.44777972e-01, -2.85573214e-01, -2.34121397e-01,

70 ]

71 serialOscanMeansC = [

72 0.212578110, 0.107817403, -0.122200218, -0.0089812368, -0.067990060, 0.040077099,

73 -0.021402006, 0.090923088, -0.099587158, 0.274797124, -0.016930788, 0.045007070,

74 -0.00379911056, -0.16088248, 0.055911896, 0.0601755001, -0.046872945, 0.210018355,

75 0.081641635, -0.046147249, -0.0059020276, 0.108368757, -0.033966731, -0.0058644798,

76 -0.075746922, -0.203826510, 0.12620401, -0.0156685544, -0.09631182, 0.089754454,

77 0.03789926, 0.0304515115, -0.082173715, -0.061332140, -0.24894494, -0.155137551,

78 -0.073825312, 0.24538413, -0.069597074, 0.192338801, -0.0539746876, -0.184556000,

79 -0.173069382, -0.209975778, 0.086679191, 0.016299034, -0.0094125706, -0.100099911,

80 0.061981365, 0.086250364, 0.209128404, -0.0067993622, 0.171072270, -0.29266333,

81 0.075172274, -0.29375612, -0.13377650, 0.0125964781, -0.124991264, 0.226516831,

82 0.128244484, -0.05019844, -0.149249925, -0.1557398,

83 ]

84 serialOscanMeansD = [

85 4.0867248, 1.43194193, 0.95319573, 0.43219185, 0.53112239, 0.28648, 0.323903486,

86 0.27622156, 0.26031138, 0.144442975, 0.0149878587, 0.062969929, 0.018541051,

87 -0.237687056, 0.22804558, 0.0600504708, 0.140250022, -0.137477808, 0.119911710,

88 0.03770870, -0.20021377, 0.187175400, 0.0168790129, -0.110724371, 0.099311580,

89 0.0079969534, -0.157593577, -0.178876067, -0.214948580, -0.11354382, 0.148154530,

90 -0.056012520, 0.11851939, 0.067902033, 0.18970736, -0.181487703, -0.0101017127,

91 0.100998570, -0.0309096733, -0.034450136, -0.066357072, -0.058662959, 0.146185921,

92 -0.218474021, -0.173691633, 0.055349625, -0.178158524, -0.012917378, -0.166576555,

93 -0.063862754, 0.113169933, -0.33518338, -0.074239500, 0.22262230, -0.066653975,

94 -0.200271016, -0.013275277, 0.100596499, -0.092528954, 0.0339541714, 0.113119135,

95 -0.150720824, 0.038237873, 0.17603852613429813,

96 ]

97 parallelOscanMeans = [

98 7.18039751e-01, 4.56550479e-01, 4.14261669e-01, 2.88099229e-01, -2.34310962e-02,

99 -4.59854975e-02, -1.14491098e-02, 5.19846082e-02, 2.05635265e-01, 1.25147207e-02,

100 ]

101

102 self.inputMeasurements = [

103 {'CTI': {'C:0,0': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

104 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

105 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

106 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

107 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

108 'C:0,1': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

109 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

110 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

111 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

112 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

113 'C:0,2': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

114 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

115 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

116 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

117 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

118 'C:0,3': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

119 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

120 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

121 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

122 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

123 'C:1,0': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

124 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

125 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

126 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

127 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

128 'C:1,1': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

129 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

130 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

131 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

132 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

133 'C:1,2': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

134 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

135 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

136 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

137 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

138 'C:1,3': {'FIRST_COLUMN_MEAN': 117.810165, 'LAST_COLUMN_MEAN': 1.09791130e+02,

139 'IMAGE_MEAN': 117.810165, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansA,

140 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

141 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

142 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]}}},

143

144 {'CTI': {'C:0,0': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

145 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

146 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

147 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

148 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

149 'C:0,1': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

150 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

151 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

152 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

153 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

154 'C:0,2': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

155 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

156 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

157 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

158 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

159 'C:0,3': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

160 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

161 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

162 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

163 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

164 'C:1,0': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

165 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

166 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

167 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

168 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

169 'C:1,1': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

170 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

171 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

172 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

173 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

174 'C:1,2': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

175 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

176 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

177 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

178 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

179 'C:1,3': {'FIRST_COLUMN_MEAN': 36562.082, 'LAST_COLUMN_MEAN': 3.45901172e+04,

180 'IMAGE_MEAN': 36562.082, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansB,

181 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

182 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

183 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]}}},

184

185 {'CTI': {'C:0,0': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

186 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

187 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

188 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

189 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

190 'C:0,1': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

191 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

192 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

193 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

194 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

195 'C:0,2': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

196 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

197 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

198 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

199 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

200 'C:0,3': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

201 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

202 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

203 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

204 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

205 'C:1,0': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

206 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

207 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

208 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

209 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

210 'C:1,1': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

211 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

212 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

213 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

214 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

215 'C:1,2': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

216 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

217 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

218 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

219 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

220 'C:1,3': {'FIRST_COLUMN_MEAN': 994.811, 'LAST_COLUMN_MEAN': 936.415,

221 'IMAGE_MEAN': 994.811, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansC,

222 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

223 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

224 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]}}},

225

226 {'CTI': {'C:0,0': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

227 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

228 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

229 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

230 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

231 'C:0,1': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

232 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

233 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

234 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

235 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

236 'C:0,2': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

237 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

238 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

239 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

240 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

241 'C:0,3': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

242 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

243 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

244 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

245 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

246 'C:1,0': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

247 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

248 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

249 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

250 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

251 'C:1,1': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

252 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

253 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

254 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

255 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

256 'C:1,2': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

257 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

258 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

259 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

260 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]},

261 'C:1,3': {'FIRST_COLUMN_MEAN': 12215.778, 'LAST_COLUMN_MEAN': 11536.875,

262 'IMAGE_MEAN': 12215.778, 'SERIAL_OVERSCAN_VALUES': serialOscanMeansD,

263 'SERIAL_OVERSCAN_COLUMNS': [x for x in range(0, 64)],

264 'PARALLEL_OVERSCAN_VALUES': parallelOscanMeans,

265 'PARALLEL_OVERSCAN_ROWS': [x for x in range(0, 10)]}}}]

266 self.inputDims = [{'detector': 20, 'instrument': 'IsrMock', 'exposure': 2019101200433},

267 {'detector': 20, 'instrument': 'IsrMock', 'exposure': 2019101300154},

268 {'detector': 20, 'instrument': 'IsrMock', 'exposure': 2019101300004},

269 {'detector': 20, 'instrument': 'IsrMock', 'exposure': 2019101200333}]

270

271 self.task = cpPipe.CpCtiSolveTask()

272

273 self.ptc = PhotonTransferCurveDataset()

274 self.ptc.updateMetadata(camera=self.camera, detector=self.camera[20])

275 self.sequencerMetadata = {

276 "SEQNAME": "a_sequencer",

277 "SEQFILE": "a_sequencer_file",

278 "SEQCKSUM": "deadbeef",

279 }

280 self.ptc.updateMetadata(**self.sequencerMetadata, setCalibInfo=True)

281

282 def test_task(self):

283 """A test for the main CpCtiSolveTask.

284

285 This should excercise most of the new code.

286 """

287 results = self.task.run(self.inputMeasurements, self.camera, self.inputDims, self.ptc)

288

289 calib = results.outputCalib

290

291 for key, value in self.sequencerMetadata.items():

292 self.assertEqual(calib.metadata[key], value)

293

294 for key in ["INSTRUME", "DETECTOR", "DET_NAME", "DET_SER"]:

295 self.assertEqual(calib.metadata[key], self.ptc.metadata[key])

296

297 # Check that the signals array is sorted

298 for ampName in calib.signals.keys():

299 assert numpy.all(numpy.diff(calib.signals[ampName]) > 0)

300

301 # Check that the result matches expectation.

302 self.assertAlmostEqual(calib.globalCti['C:0,0'], 1.0e-7, 4)

303 self.assertAlmostEqual(calib.driftScale['C:0,0'], 1.8105e-4, 4)

304 self.assertAlmostEqual(calib.decayTime['C:0,0'], 3.08095, 4)

305

306 # Check that all amps are equal.

307 for ampName in calib.globalCti.keys():

308 self.assertEqual(calib.globalCti['C:0,0'], calib.globalCti[ampName])

309 self.assertEqual(calib.driftScale['C:0,0'], calib.driftScale[ampName])

310 self.assertEqual(calib.decayTime['C:0,0'], calib.decayTime[ampName])

311 self.assertEqual(calib.signals[ampName][-2], calib.serialCtiTurnoff[ampName])

312 self.assertEqual(calib.signals[ampName][-1], calib.parallelCtiTurnoff[ampName])

313 self.assertEqual(

314 (calib.signals[ampName][-1] - calib.signals[ampName][-3]) / 2.0,

315 calib.serialCtiTurnoffSamplingErr[ampName],

316 )

317 self.assertEqual(

318 calib.signals[ampName][-1] - calib.signals[ampName][-2],

319 calib.parallelCtiTurnoffSamplingErr[ampName],

320 )

321

322

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

324 pass

325

326

327def setup_module(module):

328 lsst.utils.tests.init()

329

330

331if __name__ == "__main__": 331 ↛ 332line 331 didn't jump to line 332 because the condition on line 331 was never true

332 lsst.utils.tests.init()

333 unittest.main()