Changeset 25433 for branches/eam_branches/20090715/psphot/src/psphotPetrosianStats.c

Timestamp:

Sep 17, 2009, 2:26:32 PM (17 years ago)

Author:

eugene

Message:

change radius for extended sources to use footprint; clean up some of the visualizations; plug some leaks

File:

• branches/eam_branches/20090715/psphot/src/psphotPetrosianStats.c (modified) (3 diffs)

branches/eam_branches/20090715/psphot/src/psphotPetrosianStats.c

@@ -8 +8 @@
 float InterpolateValues (float X0, float Y0, float X1, float Y1, float X);
 
-bool psphotPetrosianStats (pmPetrosian *petrosian) {
+bool psphotPetrosianStats (pmSource *source, pmPetrosian *petrosian) {
 
     float petRadius, petFlux;
@@ -77 +77 @@
         psVectorAppend(refRadius, binRad->data.F32[i]);
 
-        if (1) {
-            fprintf (stderr, "%3d : %5.2f : %5.3f %5.3f : %5.3f %5.3f : %5.3f %5.3f : %5.3f %5.3f : %5.1f %5.1f\n", 
-                     i, refRadius->data.F32[nOut], 
-                     binSB->data.F32[i], binSBstdev->data.F32[i], 
-                     meanSB->data.F32[nOut], meanSBerr, 
-                     petRatio->data.F32[nOut], petRatioErr->data.F32[nOut], 
-                     fluxSum->data.F32[nOut], sqrt(fluxSumErr2->data.F32[nOut]), areaSum->data.F32[nOut], areaInner);
-        }
+        psTrace ("psphot", 4, "%3d : %5.2f : %5.3f %5.3f : %5.3f %5.3f : %5.3f %5.3f : %5.3f %5.3f : %5.1f %5.1f\n", 
+                 i, refRadius->data.F32[nOut], 
+                 binSB->data.F32[i], binSBstdev->data.F32[i], 
+                 meanSB->data.F32[nOut], meanSBerr, 
+                 petRatio->data.F32[nOut], petRatioErr->data.F32[nOut], 
+                 fluxSum->data.F32[nOut], sqrt(fluxSumErr2->data.F32[nOut]), areaSum->data.F32[nOut], areaInner);
     
         // anytime we transition below the PETROSIAN_RATIO, calculate the radius and flux
@@ -149 +147 @@
 
     psFree(fluxSum);
+    psFree(fluxSumErr2);
+    psFree(refRadius);
     psFree(petRatio);
+    psFree(petRatioErr);
     psFree(meanSB);
     psFree(areaSum);