| | 7 | |
| | 8 | It looks like I've fixed the ppSub / magic masking problems. The last |
| | 9 | two issues I found this week where (a) I needed to turn on SVD dynamic range |
| | 10 | clipping in the linear algebra solution to avoid introducing noise in |
| | 11 | cases where the solution was particularly poorly conditioned and (b) I |
| | 12 | discovered that the covariance averaging process was weighting by the |
| | 13 | image variance (not the inverse variance). As far as I can tell from |
| | 14 | the math and the results, a straight average (unweighted) is actually the correct |
| | 15 | answer. This latter issue meant that the covariance of the difference |
| | 16 | image was dominated by the covariance of the less convolved image, |
| | 17 | resulting in a somewhat too small flux for a given detection limit. This |
| | 18 | explains the magic issue, but it also explains the excess sources in |
| | 19 | the difference images which correspond to empty sky. This weighting |
| | 20 | choice has actually been in place since long before the 8/23 software |
| | 21 | change. I believe the reason it did not impact us earlier was that, before |
| | 22 | the 8/23 release, the dual convolution process penalties tended to |
| | 23 | result in a similar amount of smoothing on both images (with excessive overall smoothing), so the |
| | 24 | different weighting did not have too large of an impact. The code |
| | 25 | since 8/23 has been trying harder to get the minimal convolution, with |
| | 26 | the result that the the mean variance of the 2 images are more likely |
| | 27 | to be rather different from one another. |
| | 28 | |
| | 29 | I ran the 3pi test set that Heather has been running. Here are the |
| | 30 | numbers of magic streaks from the various test runs: |
| | 31 | |
| | 32 | || label || max(magic_id) || avg(magicMask.streaks) || |
| | 33 | || magictest.3Pi.trunk.default.20101230 ||108422.5833 || 122.5833 || |
| | 34 | || magictest.3Pi.trunk.modified.20100105 ||109485.9667 || 170.3000 || |
| | 35 | || magictest.3Pi.ipp20100701.20110106 ||109801.5000 || 42.9000 || |
| | 36 | || magictest.r28950 ||111028.7000 || 36.7250 || |
| | 37 | || magictest.ippeam.r30312 ||113452.5000 || 147.2500 || |
| | 38 | || magictest.ippeam.r30312.v2 ||113930.5000 || 134.9500 || |
| | 39 | || magictest.ippeam.r30312.v3 ||114296.5000 || 156.5000 || |
| | 40 | || magictest.ippeam.r30312.v4 ||114692.5000 || 222.9737 || |
| | 41 | || magictest.ippeam.r30312.v5 ||114730.5000 || 98.3611 || |
| | 42 | || magictest.ippeam.r30312.v8 ||118059.5000 || 45.4412 || |
| | 43 | |
| | 44 | At this point, we can get back to the task of running a good-sized mask / magic test suite to check on the quality of the new static masking. |
| | 45 | |
| | 46 | Meanwhile, I've started to work on the reference database while that is processing. |