| | 1 | = Gene 1 = |
| | 2 | |
| | 3 | {{{ |
| | 4 | |
| | 5 | Oct 5 |
| | 6 | Here are some notes related to the Great Reprocessing which will be |
| | 7 | useful for discussion today. |
| | 8 | |
| | 9 | In terms of IPP issues, I believe the following are needed before the |
| | 10 | reprocessing can start: |
| | 11 | |
| | 12 | 1) new hardware delivered and installed (shipped late last week, |
| | 13 | expected in 2-3 weeks according to Freddie) |
| | 14 | |
| | 15 | 2) star/galaxy/CR separation finalized. While I was in Harvard, I got |
| | 16 | Eddie Schlafly up to speed on this issue. I gave him a collection of |
| | 17 | SMF files in the latest format that includes all of the parameters |
| | 18 | which I am currently trying to use to make the star / galaxy / CR |
| | 19 | separation. These exposures overlap SDSS, so he is able to use that |
| | 20 | as ground truth. He is looking at the parameters I measure (Mxx, Myy, |
| | 21 | PSF mag, Kron mag, other aperture-like mags) and attempting to define |
| | 22 | the cuts in this multi-dimensional space that define stars, galaxies, |
| | 23 | and CRs. He is first working to define the a-priori probability that |
| | 24 | a given object is consistent with a PSF. Next, I've asked him to work |
| | 25 | with Shaun Cole on defining the a-posteriori probabilities (star, |
| | 26 | galaxy) given some other information (specifically galactic latitude |
| | 27 | and longitude). He has already made some headway in better defining |
| | 28 | the CRs from the faint stars. There is some subtlety here that I've |
| | 29 | asked him to explore: is it enough to define the probability of |
| | 30 | something being a star (PSF) using on the locations of stars in the |
| | 31 | space, or is it necessary to include some description of the shape |
| | 32 | distribution of the galaxies themselves. |
| | 33 | |
| | 34 | 3) decisions / implications to difference image analysis from |
| | 35 | photfest. The big advance over last week was that the CfA group now |
| | 36 | agree that ppSub SINGLE direction convolution produces results |
| | 37 | indistinguishable from the photpipe difference image analysis. They |
| | 38 | also convinced me that, for the stack-stack diffs where we have a deep |
| | 39 | reference stack, it is acceptable to always choose the deep reference |
| | 40 | stack as the convolution target -- the occasional instances of small |
| | 41 | amounts of deconvolutions are acceptable given the high S/N of the |
| | 42 | reference stack and the over-sampling. We also agreed that the DUAL |
| | 43 | convolution is at worst equivalent to an optimal SINGLE convolution |
| | 44 | plus an additional smoothing, and it may be the case that the DUAL |
| | 45 | convolution is acceptable as well. In the near term, we are going to |
| | 46 | switch the MD stack-stack diff processing to use SINGLE convolution of |
| | 47 | the refstack, and the CfA group will start to use at least the IPP |
| | 48 | pixels. We (mostly they -- Mark Huber, Micheal Wood-Vasey, and Ken |
| | 49 | Smith) are looking at the detections produced by the IPP, along with |
| | 50 | the feedback I gave them on the flags and exclusions, to see if (with |
| | 51 | appropriate filtering) the false positives and false negatives are |
| | 52 | acceptable and if the photometry is is consistent with the photpipe |
| | 53 | equivalent. The quick look at the last point was positive, but they |
| | 54 | are doing more checks on this. |
| | 55 | |
| | 56 | As far as the Great Reprocessing goes, the main lesson & outstanding |
| | 57 | point from the photfest was the remaining concern that DUAL |
| | 58 | convolution STILL results in more smoothing than one would expect or |
| | 59 | hope. Since DUAL is required for the warp-warp diffs (ie, 3pi and |
| | 60 | SS), it is important to continue to push on this. Two possible |
| | 61 | avenues: 1) the range of kernels which we use is not as large as CfA: |
| | 62 | we have a set of 3 Gaussians which we scale by the FWHM, while Armin |
| | 63 | adds an additional large Gaussian if the FWHM gets too large. 2) the |
| | 64 | DUAL convolution could potentially be made to fall-over to SINGLE |
| | 65 | convolution if only one of the two images is convolved. This is |
| | 66 | tricky to determine, but a bit of work on that may be useful. |
| | 67 | |
| | 68 | 4) address recent large variations in processing throughput (perhaps |
| | 69 | due to nebulous mysql fragmentation?). Serge & Roy have learned that |
| | 70 | running the mysql optimization process may help to speed up the |
| | 71 | nebulous interactions; Serge is investigating this and will also look |
| | 72 | into partitioning the big nebulous table. The throughput variations |
| | 73 | can be a factor of 3.5! |
| | 74 | |
| | 75 | 5) updates to non-linear correction, dark correction, and resulting |
| | 76 | mask modifcations. |
| | 77 | |
| | 78 | 6) a full processing strategy (what area first, exactly what |
| | 79 | processing is done) and any resulting scheduling code changes (ie, |
| | 80 | updates to the nightly science code to manage the processing). |
| | 81 | |
| | 82 | 7) astrometric and photometric reference catalog for the portion of |
| | 83 | the sky that is to be first processed. |
| | 84 | |
| | 85 | It would be nice, but is not required, that the 3pi static sky |
| | 86 | analysis be finalized, and that the IPP->PSPS interface be finished. |
| | 87 | It is also important to have a plan for the storage upgrades, but we |
| | 88 | have enough space to start. |
| | 89 | |
| | 90 | ]]] |