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

Timestamp:

Sep 19, 2009, 8:49:41 AM (17 years ago)

Author:

eugene

Message:

further work on the petrosian analysis & extended sources

File:

• branches/eam_branches/20090715/psphot/src/psphotPetrosianStats.c (modified) (2 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 (pmSource *source, pmPetrosian *petrosian) {
+bool psphotPetrosianStats (pmSource *source) {
+
+    pmSourceRadialProfile *profile = source->extpars->profile;
 
     float petRadius, petFlux;
 
-    psVector *binSB      = petrosian->binSB;
-    psVector *binSBstdev = petrosian->binSBstdev;
-    psVector *binRad     = petrosian->radialBins;
-    psVector *area       = petrosian->area;
+    psVector *binSB      = profile->binSB;
+    psVector *binSBstdev = profile->binSBstdev;
+    psVector *binRad     = profile->radialBins;
+    psVector *area       = profile->area;
 
     psVector *fluxSum     = psVectorAllocEmpty(binSB->n, PS_TYPE_F32);
@@ -139 +141 @@
     }
 
-    // save petRadius, petFlux
+    if (!source->extpars->petrosian_80) {
+        source->extpars->petrosian_80 = pmSourceExtendedFluxAlloc ();
+    }
+    pmSourceExtendedFlux *petrosian = source->extpars->petrosian_80;
+
     // XXX save flags (anyPetro, manyPetro)
-    petrosian->petrosianRadius = petRadius;
-    petrosian->petrosianFlux   = petFlux;
+    petrosian->radius = petRadius;
+    petrosian->flux   = petFlux;
 
-    psphotPetrosianVisualStats (binRad, binSB, refRadius, meanSB, petRatio, petRatioErr, fluxSum, petRadius, PETROSIAN_RATIO, petFlux, apRadius);
+    // psphotPetrosianVisualStats (binRad, binSB, refRadius, meanSB, petRatio, petRatioErr, fluxSum, petRadius, PETROSIAN_RATIO, petFlux, apRadius);
 
     psFree(fluxSum);