Changeset 25178 for branches/eam_branches/20090715/psphot/src/psphotPetrosianVisual.c

Timestamp:

Aug 24, 2009, 8:40:34 AM (17 years ago)

Author:

eugene

Message:

various petrosian / radial profile analysis improvements

File:

• branches/eam_branches/20090715/psphot/src/psphotPetrosianVisual.c (modified) (5 diffs)

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

@@ -161 +161 @@
 bool psphotPetrosianVisualStats (psVector *radBin, psVector *fluxBin, 
                                  psVector *refRadius, psVector *meanSB, 
-                                 psVector *petRatio, psVector *fluxSum, 
-                                 float petRadius, float petFlux)
+                                 psVector *petRatio, psVector *petRatioErr,
+                                 psVector *fluxSum, 
+                                 float petRadius, float ratioForRadius,
+                                 float petFlux, float radiusForFlux)
 {
     Graphdata graphdata;
@@ -229 +231 @@
 
     graphdata.color = KapaColorByName ("black");
-    pmVisualLimitsFromVectors (&graphdata, radBin, petRatio);
+    graphdata.ymax = +1.05;
+    graphdata.ymin = -0.05;
+    pmVisualLimitsFromVectors (&graphdata, radBin, NULL);
     KapaSetLimits (kapa, &graphdata);
 
@@ -240 +244 @@
     graphdata.ptype = 0;
     graphdata.size = 1.0;
+    graphdata.etype = 0x01;
     KapaPrepPlot (kapa, refRadius->n, &graphdata);
     KapaPlotVector (kapa, refRadius->n, refRadius->data.F32, "x");
     KapaPlotVector (kapa, refRadius->n, petRatio->data.F32, "y");
-
-    float ratio_90 = 0.1;
+    KapaPlotVector (kapa, refRadius->n, petRatioErr->data.F32, "dym");
+    KapaPlotVector (kapa, refRadius->n, petRatioErr->data.F32, "dyp");
+    graphdata.etype = 0;
+
     graphdata.color = KapaColorByName ("red");
     graphdata.style = 2;
@@ -251 +258 @@
     KapaPrepPlot   (kapa, 1, &graphdata);
     KapaPlotVector (kapa, 1, &petRadius, "x");
-    KapaPlotVector (kapa, 1, &ratio_90, "y");
+    KapaPlotVector (kapa, 1, &ratioForRadius, "y");
 
     // *** section 3: radius vs integrated flux
@@ -283 +290 @@
     graphdata.size = 2.0;
     KapaPrepPlot   (kapa, 1, &graphdata);
-    KapaPlotVector (kapa, 1, &petRadius, "x");
+    KapaPlotVector (kapa, 1, &radiusForFlux, "x");
     KapaPlotVector (kapa, 1, &petFlux, "y");