Changeset 32977 for trunk/psModules/src/objects/pmSourceIO_CMF_PS1_SV1.c

Timestamp:

Dec 19, 2011, 3:12:24 PM (14 years ago)

Author:

bills

Message:

pass readout down to the Read XRAD functions so that we can save the
fwhmValues vector on the analysis metadata

File:

• trunk/psModules/src/objects/pmSourceIO_CMF_PS1_SV1.c (modified) (4 diffs)

trunk/psModules/src/objects/pmSourceIO_CMF_PS1_SV1.c

@@ -1026 +1026 @@
 }
 
-bool pmSourcesRead_CMF_PS1_SV1_XRAD(psFits *fits, psMetadata *hduHeader, psArray *sources, long *sourceIndex)
+bool pmSourcesRead_CMF_PS1_SV1_XRAD(psFits *fits, pmReadout *readout, psMetadata *hduHeader, psArray *sources, long *sourceIndex)
 {
     PS_ASSERT_PTR_NON_NULL(fits, false);
@@ -1037 +1037 @@
         return false;
     }
+
+    long       seq_first = -1;
+    long       seq_last = -1;
+    psVector   *fwhmValues = psVectorAllocEmpty(10, PS_TYPE_F32);
+    long       max_entries = -1;
+    long       num_entries = -1;
 
     for (long i = 0; i < numSources; i++) {
@@ -1061 +1067 @@
             return false;
         }
+        if (seq_first == -1) {
+            seq_first = seq;
+        }
+        if (seq == seq_first) {
+            psF32 value = psMetadataLookupF32(&status, row, "PSF_FWHM");
+            psVectorAppend(fwhmValues, value);
+        }
+        if (seq == seq_last) {
+            num_entries++;
+        } else {
+            num_entries = 1;
+            seq_last = seq;
+        }
+        if (num_entries > max_entries) {
+            max_entries = num_entries;
+        }
+
         if (!source->radialAper) {
             // XXX: where to find the number of models to expect?
@@ -1078 +1101 @@
     }
 
+    // check for consistency between the length of fwhmValues and the maximum number of entries for each row
+    if (fwhmValues->n != max_entries) {
+        psError(PS_ERR_PROGRAMMING, true, "number of PSF_FWHM values found %ld does not match expected number: %ld\n",
+            fwhmValues->n, max_entries);
+        psAssert(0, "fixme");
+    }
+
+    if (!readout->analysis) {
+        readout->analysis = psMetadataAlloc();
+    }
+
+    psMetadataAddVector(readout->analysis, PS_LIST_TAIL, "STACK.PSF.FWHM.VALUES", PS_META_REPLACE, "PSF sizes", fwhmValues);
+    psFree(fwhmValues);
+
     return true;
 }