Changeset 34308

Timestamp:

Aug 14, 2012, 1:26:39 PM (14 years ago)

Author:

bills

Message:

get sigma values for moments window measurement from recipe. Omit the 12 and 18 pixel
sigmas for GPC1

Location:

tags/ipp-20120802

Files:

• ippconfig/gpc1/psphot.config (modified) (1 diff)
• ippconfig/recipes/psphot.config (modified) (1 diff)
• psphot/src/psphotSourceStats.c (modified) (2 diffs)

tags/ipp-20120802/ippconfig/gpc1/psphot.config

@@ -19 +19 @@
 # PSF model parameters : choose the PSF model desired
 PSF_MODEL                           STR   PS_MODEL_PS1_V1
+# vector of simga values to use in when setting the moments window (this is shorter for gpc1 than
+# the psphot default)
+@PSF.SIGMA.VALUES                   F32   1.0  2.0  3.0  4.5  6.0  9.0 # comment due to ticket #1476
+
 PSPHOT.CR.NSIGMA.SOFTEN             F32   0.0025        # Softening parameter for weights in CR stat
-
 PSPHOT.CR.MAX.SIZE                  F32   0.8             # guess for limit on M_minor in pix^2
 PSPHOT.CR.MAX.MAG                   F32  -8.0             # guess for limit on Kron Mag

tags/ipp-20120802/ippconfig/recipes/psphot.config

@@ -122 +122 @@
 PSF_MODEL                           STR   PS_MODEL_GAUSS
 # PSF_MODEL                         STR   PS_MODEL_PGAUSS
+# sigma values used by psphotSetMomentsWindow. Maximum length of this array is 8
+@PSF.SIGMA.VALUES                   F32   1.0  2.0  3.0  4.5  6.0  9.0  12.0  18.0 # comment due to ticket #1476
 
 # PSF.TREND.MASK must be a 2D polynomial

tags/ipp-20120802/psphot/src/psphotSourceStats.c

@@ -521 +521 @@
     float PSF_SN_LIM = 2.0*psMetadataLookupF32(&status, recipe, "PSF_SN_LIM"); psAssert (status, "missing PSF_SN_LIM");
 
-    // XXX move this to a config file?
+
     # define NSIGMA 8
-    float sigma[NSIGMA] = {1.0, 2.0, 3.0, 4.5, 6.0, 9.0, 12.0, 18.0};
+    // moved to config file
+    // float sigma[NSIGMA]  = {1.0, 2.0, 3.0, 4.5, 6.0, 9.0, 12.0, 18.0};
+    psVector  *sigmavec = psMetadataLookupPtr (&status, recipe, "PSF.SIGMA.VALUES"); psAssert(status, "missing PSF.SIGMA.VALUES");
+    psAssert(sigmavec->n <= NSIGMA, "too many sigma values in recipe %ld maximum is %d", sigmavec->n, NSIGMA);
+
+    float sigma[NSIGMA];
     float Sout[NSIGMA];
     float Rmin[NSIGMA];
     int   Nout[NSIGMA]; // number of stars found in clump : use this to control the number of regions measured by psphotRoughClass
 
+    // copy the data from vector to local array to keep the code below easier to read
+    for (int i = 0 ; i < sigmavec->n; i++) {
+        sigma[i] = sigmavec->data.F32[i];
+    }
+
     // this sorts by peak->rawFlux
     sources = psArraySort (sources, pmSourceSortByFlux);
 
     // loop over radii:
-    for (int i = 0; i < NSIGMA; i++) {
+    for (int i = 0; i < sigmavec->n; i++) {
 
         // XXX move max source number to config
@@ -589 +599 @@
     float minS = Sout[0];
     float maxS = Sout[0];
-    for (int i = 0; i < NSIGMA; i++) {
+    for (int i = 0; i < sigmavec->n; i++) {
         minS = PS_MIN(Sout[i], minS);
         maxS = PS_MAX(Sout[i], maxS);
     }
-    if (minS > 0.65) Sigma = sigma[NSIGMA-1];
+    if (minS > 0.65) Sigma = sigma[sigmavec->n-1];
     if (maxS < 0.65) Sigma = sigma[0];
 
-    for (int i = 0; i < NSIGMA - 1 && isnan(Sigma); i++) {
+    for (int i = 0; i < sigmavec->n - 1 && isnan(Sigma); i++) {
         if (!isfinite(Sout[i]) || !isfinite(Sout[i+1])) continue;
         if ((Sout[i] > 0.65) && (Sout[i+1] > 0.65)) continue;