Changeset 30707 for branches/eam_branches/ipp-20110213/psphot/src/psphotPetrosianRadialBins.c

Timestamp:

Feb 19, 2011, 10:40:50 AM (15 years ago)

Author:

eugene

Message:

remove invalid references to isophot and kron from extended source analysis; use the measured sky stdev instead of the user guess; optionally allow stars as well as galaxies to get petrosians; fix target vs input seeing in psphotStack; allow radial apertures for single-image analysis as well as stack; calculate the radial aperture flux errors; calculate the petrosian parameter errors; turn on threads for ext source fits; save covar matrix for ext source fits; calculate petrosian fill factor

Location:

branches/eam_branches/ipp-20110213/psphot

Files:

• . (modified) (1 prop)
• src/psphotPetrosianRadialBins.c (modified) (4 diffs)

branches/eam_branches/ipp-20110213/psphot/src/psphotPetrosianRadialBins.c

@@ -51 +51 @@
     psVector *binRad     = psVectorAllocEmpty(nMax, PS_TYPE_F32); // mean radius of radial bin
     psVector *binArea    = psVectorAllocEmpty(nMax, PS_TYPE_F32); // area of radial bin (contiguous, non-overlapping)
+    psVector *binFill    = psVectorAllocEmpty(nMax, PS_TYPE_F32); // fraction of radial bin with valid pixels
 
     psVectorInit (binSB, 0.0);
@@ -144 +145 @@
             binSB->data.F32[nOut] = value;
             binSBstdev->data.F32[nOut] = sqrt(PS_SQR(dvalue) / values->n + skyModelErrorSQ);
+            binFill->data.F32[nOut] = values->n / binArea->data.F32[nOut];
 
             // error in the SB is the stdev per bin / sqrt (number of pixels) 
@@ -176 +178 @@
         psFree(binRad);
         psFree(binArea);
+        psFree(binFill);
         psFree(radMin);
         psFree(radMax);
@@ -202 +205 @@
     psFree(profile->radialBins);
     psFree(profile->area);
+    psFree(profile->binFill);
 
     // save the vectors
     profile->radialBins = binRad;
     profile->area       = binArea;
+    profile->binFill    = binFill;
     profile->binSB      = binSB;
     profile->binSBstdev = binSBstdev;