Changeset 17901 for branches/eam_branch_20080511/ppSim/src/ppSimDetections.c

Timestamp:

Jun 4, 2008, 1:52:32 PM (18 years ago)

Author:

eugene

Message:

significant work on the I/O pmFPAfiles, especially to get the format / metaCamera correct for the input psf vs the input image

File:

• branches/eam_branch_20080511/ppSim/src/ppSimDetections.c (modified) (1 diff)

branches/eam_branch_20080511/ppSim/src/ppSimDetections.c

@@ -2 +2 @@
 
 bool ppSimDetections (psImage *significance, psMetadata *recipe, psArray *sources) {
+    psAssert (sources, "programming error: ppSimDetections passed NULL sources");
 
-  bool status;
+    bool status;
 
-  // for each source, measure the significance of the peak at the given coordinate
-  // where does this get stored?
+    // for each source, measure the significance of the peak at the given coordinate
+    // where does this get stored?
 
-  // XXX need to get the effective Area from the PSF sigma
-  float SIGMA_SMTH  = psMetadataLookupF32 (&status, recipe, "SIGMA_SMOOTH"); 
-  psAssert (status, "SIGMA_SMOOTH missing: call psphotSignificanceImage first"); 
-  float effArea = 4.0*M_PI*PS_SQR(SIGMA_SMTH);
+    // XXX need to get the effective Area from the PSF sigma
+    float SIGMA_SMTH  = psMetadataLookupF32 (&status, recipe, "SIGMA_SMOOTH"); 
+    psAssert (status, "SIGMA_SMOOTH missing: call psphotSignificanceImage first"); 
+    float effArea = 4.0*M_PI*PS_SQR(SIGMA_SMTH);
 
-  int row0 = significance->row0;
-  int col0 = significance->col0;
+    int row0 = significance->row0;
+    int col0 = significance->col0;
 
-  for (int i = 0; i < sources->n; i++) {
+    for (int i = 0; i < sources->n; i++) {
 
-    pmSource *source = sources->data[i];
-    pmPeak *peak = source->peak;
-    psAssert (peak, "peak is not defined for the source");
+        pmSource *source = sources->data[i];
+        pmPeak *peak = source->peak;
+        psAssert (peak, "peak is not defined for the source");
 
-    peak->value = significance->data.F32[peak->y-row0][peak->x-col0];
-    peak->SN = sqrt(peak->value*effArea);
+        peak->value = significance->data.F32[peak->y-row0][peak->x-col0];
+        peak->SN = sqrt(peak->value*effArea);
 
-  }
-  return true;
+    }
+    return true;
 }