Changeset 38058 for branches/eam_branches/ipp-20150112/psModules/src/objects

Timestamp:

Mar 29, 2015, 6:08:31 AM (11 years ago)

Author:

eugene

Message:

merge changes from trunk

Location:

branches/eam_branches/ipp-20150112/psModules/src/objects

Files:

• . (modified) (2 props)
• Makefile.am (modified) (3 diffs)
• mksource.pl (modified) (1 diff)
• pmFootprintAssignPeaks.c (modified) (1 diff)
• pmModelFuncs.h (modified) (1 diff)
• pmSourceIO.c (modified) (6 diffs)
• pmSourceIO.h (modified) (1 diff)
• pmSourceIO_CMF.c.in (modified) (8 diffs, 1 prop)
• pmSourceIO_PS1_CAL_0.c (modified) (1 diff)
• pmSourceIO_PS1_DEV_1.c (modified) (1 diff)
• pmSourceLensing.c (modified) (2 diffs)
• pmSourceLensing.h (modified) (1 diff)
• pmSourceOutputs.c (modified) (1 diff)

branches/eam_branches/ipp-20150112/psModules/src/objects

@@ -21 +21 @@
 pmSourceIO_CMF_PS1_DV3.c
 pmSourceIO_CMF_PS1_DV4.c
-
+pmSourceIO_CMF_PS1_DV5.c

@@ -21 +21 @@
 pmSourceIO_CMF_PS1_DV3.c
 pmSourceIO_CMF_PS1_DV4.c
-
+pmSourceIO_CMF_PS1_DV5.c

branches/eam_branches/ipp-20150112/psModules/src/objects/Makefile.am

@@ -58 +58 @@
         pmSourceIO_CMF_PS1_DV3.c \
         pmSourceIO_CMF_PS1_DV4.c \
+        pmSourceIO_CMF_PS1_DV5.c \
         pmSourceIO_MatchedRefs.c \
         pmSourcePlots.c \
@@ -156 +157 @@
 pmSourceIO_CMF_PS1_DV3.c \
 pmSourceIO_CMF_PS1_DV4.c \
+pmSourceIO_CMF_PS1_DV5.c \
 pmSourceIO_CMF_PS1_SV1.c \
 pmSourceIO_CMF_PS1_SV2.c \
@@ -188 +190 @@
         mksource.pl pmSourceIO_CMF.c.in PS1_DV4 pmSourceIO_CMF_PS1_DV4.c
 
+pmSourceIO_CMF_PS1_DV5.c : pmSourceIO_CMF.c.in mksource.pl
+        mksource.pl pmSourceIO_CMF.c.in PS1_DV5 pmSourceIO_CMF_PS1_DV5.c
+
 pmSourceIO_CMF_PS1_SV1.c : pmSourceIO_CMF.c.in mksource.pl
         mksource.pl pmSourceIO_CMF.c.in PS1_SV1 pmSourceIO_CMF_PS1_SV1.c

branches/eam_branches/ipp-20150112/psModules/src/objects/mksource.pl

@@ -25 +25 @@
                 "PS1_DV3", 3,
                 "PS1_DV4", 4,
+                "PS1_DV5", 5,
     );
 %cmfmodes_sv = ("PS1_SV1", 1,

branches/eam_branches/ipp-20150112/psModules/src/objects/pmFootprintAssignPeaks.c

@@ -56 +56 @@
         const int y = peak->y - row0;
         
-        if (ids) { assert (x >= 0 && x < numCols && y >= 0 && y < numRows);}
+        if (ids) { psAssert (x >= 0 && x < numCols && y >= 0 && y < numRows, "out of range");}
         int id = ids ? ids->data.S32[y][x] : 0;
 

branches/eam_branches/ipp-20150112/psModules/src/objects/pmModelFuncs.h

@@ -36 +36 @@
 
 typedef enum {
-    PM_MODEL_STATUS_NONE           = 0x000, ///< model fit not yet attempted, no other info
-    PM_MODEL_STATUS_FITTED         = 0x001, ///< model fit completed
-    PM_MODEL_STATUS_NONCONVERGE    = 0x002, ///< model fit did not converge
-    PM_MODEL_STATUS_OFFIMAGE       = 0x004, ///< model fit drove out of range
-    PM_MODEL_STATUS_BADARGS        = 0x008, ///< model fit called with invalid args
-    PM_MODEL_STATUS_LIMITS         = 0x010, ///< model parameters hit limits
-    PM_MODEL_STATUS_WEAK_FIT       = 0x020, ///< model fit met loose tolerance, but not tight tolerance
-    PM_MODEL_STATUS_NAN_CHISQ      = 0x040, ///< model fit failed with a NAN chisq 
-    PM_MODEL_SERSIC_PCM_FAIL_GUESS = 0x080, ///< sersic model fit failed on the initial moments-based guess
-    PM_MODEL_SERSIC_PCM_FAIL_GRID  = 0x100, ///< sersic model fit failed on the grid search
-    PM_MODEL_PCM_FAIL_GUESS        = 0x200, ///< non-sersic model fit failed on the initial moments-based guess
-    PM_MODEL_BEST_FIT              = 0x400, ///< this model was the best fit and was subtracted
+    PM_MODEL_STATUS_NONE           = 0x0000, ///< model fit not yet attempted, no other info
+    PM_MODEL_STATUS_FITTED         = 0x0001, ///< model fit completed
+    PM_MODEL_STATUS_NONCONVERGE    = 0x0002, ///< model fit did not converge
+    PM_MODEL_STATUS_OFFIMAGE       = 0x0004, ///< model fit drove out of range
+    PM_MODEL_STATUS_BADARGS        = 0x0008, ///< model fit called with invalid args
+    PM_MODEL_STATUS_LIMITS         = 0x0010, ///< model parameters hit limits
+    PM_MODEL_STATUS_WEAK_FIT       = 0x0020, ///< model fit met loose tolerance, but not tight tolerance
+    PM_MODEL_STATUS_NAN_CHISQ      = 0x0040, ///< model fit failed with a NAN chisq 
+    PM_MODEL_SERSIC_PCM_FAIL_GUESS = 0x0080, ///< sersic model fit failed on the initial moments-based guess
+    PM_MODEL_SERSIC_PCM_FAIL_GRID  = 0x0100, ///< sersic model fit failed on the grid search
+    PM_MODEL_PCM_FAIL_GUESS        = 0x0200, ///< non-sersic model fit failed on the initial moments-based guess
+    PM_MODEL_BEST_FIT              = 0x0400, ///< this model was the best fit and was subtracted
+    PM_MODEL_STATUS_NAN_SHAPE      = 0x0800, ///< model ellipse parameters do not transform to valid ellipse
+    PM_MODEL_STATUS_NAN_SHAPE_ERR  = 0x1000, ///< model ellipse parameters errors do not transform to valid ellipse
+    PM_MODEL_STATUS_FAIL_SHAPE_ERR = 0x2000, ///< could not find an MC solution for ellipse parameter errors
 } pmModelStatus;
 

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceIO.c

@@ -607 +607 @@
             PM_SOURCES_WRITE("PS1_DV3",   CMF_PS1_DV3);
             PM_SOURCES_WRITE("PS1_DV4",   CMF_PS1_DV4);
+            PM_SOURCES_WRITE("PS1_DV5",   CMF_PS1_DV5);
 
             psFree (outhead);
@@ -1130 +1131 @@
             PM_SOURCES_READ_PSF("PS1_DV3",   CMF_PS1_DV3);
             PM_SOURCES_READ_PSF("PS1_DV4",   CMF_PS1_DV4);
+            PM_SOURCES_READ_PSF("PS1_DV5",   CMF_PS1_DV5);
 
             if (!sources) {
@@ -1419 +1421 @@
         PM_SOURCES_READ_XSRC("PS1_DV3",   CMF_PS1_DV3);
         PM_SOURCES_READ_XSRC("PS1_DV4",   CMF_PS1_DV4);
+        PM_SOURCES_READ_XSRC("PS1_DV5",   CMF_PS1_DV5);
     }
     psFree(tableHeader);
@@ -1462 +1465 @@
         PM_SOURCES_READ_XFIT("PS1_DV3",   CMF_PS1_DV3);
         PM_SOURCES_READ_XFIT("PS1_DV4",   CMF_PS1_DV4);
+        PM_SOURCES_READ_XFIT("PS1_DV5",   CMF_PS1_DV5);
     }
     psFree(tableHeader);
@@ -1504 +1508 @@
         PM_SOURCES_READ_XRAD("PS1_DV3",   CMF_PS1_DV3);
         PM_SOURCES_READ_XRAD("PS1_DV4",   CMF_PS1_DV4);
+        PM_SOURCES_READ_XRAD("PS1_DV5",   CMF_PS1_DV5);
     }
     psFree(tableHeader);
@@ -1546 +1551 @@
         PM_SOURCES_READ_XGAL("PS1_DV3",   CMF_PS1_DV3);
         PM_SOURCES_READ_XGAL("PS1_DV4",   CMF_PS1_DV4);
+        PM_SOURCES_READ_XGAL("PS1_DV5",   CMF_PS1_DV5);
     }
     psFree(tableHeader);

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceIO.h

@@ -46 +46 @@
 MK_PROTO(CMF_PS1_DV3);
 MK_PROTO(CMF_PS1_DV4);
+MK_PROTO(CMF_PS1_DV5);
 
 int pmSourceGetDophotType (pmSource *source);

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceIO_CMF.c.in

@@ -197 +197 @@
           @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SM_OBJ",   PS_DATA_F32, "smear polarizability element (objects)",     source->lensingOBJ->smear->e2);
         }
-
         if (source->lensingOBJ && source->lensingOBJ->shear) {
           @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_X11_SH_OBJ",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->shear->X11);
@@ -205 +204 @@
           @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SH_OBJ",   PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->shear->e2); 
         }
+        if (false && source->lensingOBJ) {
+          // do not bother to save these as they are equivalent to Mxx,Mxy,Myy above
+          @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E1_PSF",      PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->e1); 
+          @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_PSF",      PS_DATA_F32, "shear polarizability element (objects)",     source->lensingOBJ->e2); 
+        }
 
         if (source->lensingPSF && source->lensingPSF->smear) {
@@ -213 +217 @@
           @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SM_PSF",   PS_DATA_F32, "smear polarizability element (objects)",     source->lensingPSF->smear->e2); 
         }
-
         if (source->lensingPSF && source->lensingPSF->shear) {
           @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_X11_SH_PSF",  PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->shear->X11);
@@ -221 +224 @@
           @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_SH_PSF",   PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->shear->e2); 
         }
+        if (source->lensingPSF) {
+          @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E1_PSF",      PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->e1); 
+          @>PS1_V4@               psMetadataAdd (row, PS_LIST_TAIL, "LENS_E2_PSF",      PS_DATA_F32, "shear polarizability element (objects)",     source->lensingPSF->e2); 
+        }
 
         // if lensing params exist also include the backmapped chipID and chip coordinates
@@ -240 +247 @@
         @>PS1_V3@                 psMetadataAdd (row, PS_LIST_TAIL, "SKY_LIMIT_FLUX",   PS_DATA_F32, "Flux / pix where object hits sky",           source->skyFlux);
         @>PS1_V3@                 psMetadataAdd (row, PS_LIST_TAIL, "SKY_LIMIT_SLOPE",  PS_DATA_F32, "d(Flux/pix)/dRadius where object hits sky",  source->skySlope);
+
+        @>PS1_DV4@                psMetadataAdd (row, PS_LIST_TAIL, "SRC_CHIP_NUM",     PS_DATA_S16, "id of warp input chip",                      source->chipNum); 
+        @>PS1_DV4@                psMetadataAdd (row, PS_LIST_TAIL, "SRC_CHIP_X",       PS_DATA_S16, "x coord in warp input chip",                 source->chipX); 
+        @>PS1_DV4@                psMetadataAdd (row, PS_LIST_TAIL, "SRC_CHIP_Y",       PS_DATA_S16, "y coord in warp input chip",                 source->chipY); 
+        @>PS1_DV4@                psMetadataAdd (row, PS_LIST_TAIL, "PADDING3",         PS_DATA_S16, "more padding", 0);
 
         @PS1_DV?@                 psMetadataAdd (row, PS_LIST_TAIL, "DIFF_NPOS",        PS_DATA_S32, "nPos (n pix > 3 sigma)",                     diffStats.nGood);
@@ -481 +493 @@
           @>PS1_V4@ source->lensingPSF->shear->e1  = psMetadataLookupF32 (&status, row, "LENS_E1_SH_PSF");
           @>PS1_V4@ source->lensingPSF->shear->e2  = psMetadataLookupF32 (&status, row, "LENS_E2_SH_PSF");
+          @>PS1_V4@ source->lensingPSF->e1         = psMetadataLookupF32 (&status, row, "LENS_E1_PSF");
+          @>PS1_V4@ source->lensingPSF->e2         = psMetadataLookupF32 (&status, row, "LENS_E2_PSF");
         }
 
@@ -1238 +1252 @@
         }
 
-        // XXX : make this depend on what is in the cmf
-        if (0) {
+        // NOTE: we no longer write out the covariance matrix
+        if (false) {
             // read the covariance matrix
             int nparams = model->params->n;
@@ -1252 +1266 @@
             model->covar = covar;
         }
+
+        // we are only saving the values stored in dPAR[SXX,etc]
+        dPAR[PM_PAR_SXX] = psMetadataLookupF32(&status, row, "EXT_WIDTH_MAJ_ERR");
+        dPAR[PM_PAR_SYY] = psMetadataLookupF32(&status, row, "EXT_WIDTH_MIN_ERR");
+        dPAR[PM_PAR_SXY] = psMetadataLookupF32(&status, row, "EXT_THETA_ERR");
+
+        // other parameters that we need to read
+        PAR[PM_PAR_SKY] = psMetadataLookupF32(&status, row, "SKY_EXT");
 
         if (modelType == extModelType) {

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceIO_PS1_CAL_0.c

@@ -376 +376 @@
     psVector *radialBinsLower = psMetadataLookupPtr (&status, recipe, "RADIAL.ANNULAR.BINS.LOWER");
     psVector *radialBinsUpper = psMetadataLookupPtr (&status, recipe, "RADIAL.ANNULAR.BINS.UPPER");
-    assert (radialBinsLower->n == radialBinsUpper->n);
+    psAssert (radialBinsLower->n == radialBinsUpper->n, "upper and lower bins must match");
 
     // we write out all sources, regardless of quality.  the source flags tell us the state

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceIO_PS1_DEV_1.c

@@ -322 +322 @@
     psVector *radialBinsLower = psMetadataLookupPtr (&status, recipe, "RADIAL.ANNULAR.BINS.LOWER");
     psVector *radialBinsUpper = psMetadataLookupPtr (&status, recipe, "RADIAL.ANNULAR.BINS.UPPER");
-    assert (radialBinsLower->n == radialBinsUpper->n);
+    psAssert (radialBinsLower->n == radialBinsUpper->n, "upper and lower bins must match");
 
     // we write out all sources, regardless of quality.  the source flags tell us the state

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceLensing.c

@@ -55 +55 @@
   tmp->smear = NULL;
   tmp->shear = NULL;
+  tmp->e1 = NAN;
+  tmp->e2 = NAN;
 
   return tmp;
@@ -105 +107 @@
   smear->X12 = R*0.5*s4*(moments->Mxxxy - moments->Mxyyy);
 
-  smear->e1  = R*(s2*(moments->Myy - moments->Mxx) + 0.25*s4 * (moments->Myyyy - moments->Mxxxx));
+  smear->e1  = R*(s2*(moments->Myy - moments->Mxx) + 0.25*s4 * (moments->Mxxxx - moments->Myyyy));
 
   smear->e2  = R*(0.5*s4*(moments->Mxxxy + moments->Mxyyy) - 2.0*s2*moments->Mxy);

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceLensing.h

@@ -25 +25 @@
   pmLensingPars *smear;
   pmLensingPars *shear;
+  float e1;
+  float e2;
 } pmSourceLensing; 
 

branches/eam_branches/ipp-20150112/psModules/src/objects/pmSourceOutputs.c

@@ -60 +60 @@
     // we need a measure of the image quality (FWHM) for this image, in order to get the positional errors
     *fwhmMajor = psMetadataLookupF32(&status1, readout->analysis, "FWHM_MAJ");
-    if (!status1) {
-        *fwhmMajor = psMetadataLookupF32(&status1, readout->analysis, "IQ_FW1");
-        if (!status1) {
-            *fwhmMajor = 5.0; // XXX just a guess!
-        }
-    }
+    if (status1) goto got_major;
+
+    *fwhmMajor = psMetadataLookupF32(&status1, header, "FWHM_MAJ");
+    if (status1) goto got_major;
+
+    *fwhmMajor = psMetadataLookupF32(&status1, readout->analysis, "IQ_FW1");
+    if (status1) goto got_major;
+
+    *fwhmMajor = psMetadataLookupF32(&status1, header, "IQ_FW1");
+    if (status1) goto got_major;
+
+    *fwhmMajor = 5.0; // XXX just a guess!
+
+    got_major:
+
     *fwhmMinor = psMetadataLookupF32(&status1, readout->analysis, "FWHM_MIN");
-    if (!status1) {
-        *fwhmMinor = psMetadataLookupF32(&status1, readout->analysis, "IQ_FW2");
-        if (!status1) {
-            *fwhmMinor = 5.0; // XXX just a guess!
-        }
-    }
-
+    if (status1) goto got_minor;
+    *fwhmMinor = psMetadataLookupF32(&status1, header, "FWHM_MIN");
+    if (status1) goto got_minor;
+
+    *fwhmMinor = psMetadataLookupF32(&status1, readout->analysis, "IQ_FW2");
+    if (status1) goto got_minor;
+    *fwhmMinor = psMetadataLookupF32(&status1, header, "IQ_FW2");
+    if (status1) goto got_minor;
+
+    *fwhmMinor = 5.0; // XXX just a guess!
+
+    got_minor:
     return true;
 

@@ -21 +21 @@
 pmSourceIO_CMF_PS1_DV3.c
 pmSourceIO_CMF_PS1_DV4.c
-
+pmSourceIO_CMF_PS1_DV5.c