Changeset 19888

Timestamp:

Oct 3, 2008, 12:09:25 PM (18 years ago)

Author:

eugene

Message:

add plot of CR and EXT nSigma limits

Location:

trunk/psphot/src

Files:

• psphot.h (modified) (1 diff)
• psphotReadout.c (modified) (2 diffs)
• psphotVisual.c (modified) (1 diff)

trunk/psphot/src/psphot.h

@@ -150 +150 @@
 bool psphotVisualShowResidualImage (pmConfig *config, pmReadout *readout);
 bool psphotVisualPlotApResid (pmConfig *config, psArray *sources);
+bool psphotVisualPlotSourceSize (pmConfig *config, psArray *sources);
 
 // structures & functions to support psf-convolved model fitting

trunk/psphot/src/psphotReadout.c

@@ -176 +176 @@
     }
 
+    psphotVisualPlotSourceSize (config, sources);
     psphotVisualShowSourceSize (config, view, sources);
 
@@ -249 +250 @@
     // measure source size for the remaining sources
     psphotSourceSize (config, readout, sources, recipe, 0);
+    psphotVisualPlotSourceSize (config, sources);
     psphotVisualShowSourceSize (config, view, sources);
 

trunk/psphot/src/psphotVisual.c

@@ -1236 +1236 @@
 }
 
+bool psphotVisualPlotSourceSize (pmConfig *config, psArray *sources) {
+
+    bool status;
+    Graphdata graphdata;
+    KapaSection section;
+
+    if (!isVisual) return true;
+
+    if (kapa3 == -1) {
+        kapa3 = KapaOpenNamedSocket ("kapa", "psphot:plots");
+        if (kapa3 == -1) {
+            fprintf (stderr, "failure to open kapa; visual mode disabled\n");
+            isVisual = false;
+            return false;
+        }
+    }  
+
+    // select the current recipe
+    psMetadata *recipe = psMetadataLookupPtr (&status, config->recipes, PSPHOT_RECIPE);
+    if (!recipe) {
+        fprintf (stderr, "missing recipe, skipping moments plot\n");
+        return false;
+    }
+
+    KapaClearPlots (kapa3);
+    KapaInitGraph (&graphdata);
+
+    // first section : mag vs CR nSigma
+    section.dx = 1.0;
+    section.dy = 0.5;
+    section.x = 0.0;
+    section.y = 0.0;
+    section.name = NULL;
+    psStringAppend (&section.name, "a1");
+    KapaSetSection (kapa3, &section);
+    psFree (section.name);
+
+    psVector *x = psVectorAllocEmpty (sources->n, PS_TYPE_F32);
+    psVector *y = psVectorAllocEmpty (sources->n, PS_TYPE_F32);
+
+    graphdata.xmin = +32.0;
+    graphdata.xmax = -32.0;
+    graphdata.ymin = +32.0;
+    graphdata.ymax = -32.0;
+
+    // construct the plot vectors
+    int n = 0;
+    for (int i = 0; i < sources->n; i++) {
+        pmSource *source = sources->data[i];
+        if (!source) continue;
+        if (source->type != PM_SOURCE_TYPE_STAR) continue;
+        if (!isfinite (source->crNsigma)) continue;
+
+        x->data.F32[n] = -2.5*log10(source->peak->flux);
+        y->data.F32[n] = source->crNsigma;
+        graphdata.xmin = PS_MIN(graphdata.xmin, x->data.F32[n]);
+        graphdata.xmax = PS_MAX(graphdata.xmax, x->data.F32[n]);
+        graphdata.ymin = -0.5;
+        graphdata.ymax = 10.0;
+
+        n++;
+    }
+    x->n = y->n = n;
+
+    float range;
+    range = graphdata.xmax - graphdata.xmin;
+    graphdata.xmax += 0.05*range;
+    graphdata.xmin -= 0.05*range;
+
+    // XXX set the plot range to match the image
+    KapaSetLimits (kapa3, &graphdata);
+
+    KapaSetFont (kapa3, "helvetica", 14);
+    KapaBox (kapa3, &graphdata);
+    KapaSendLabel (kapa3, "Peak as Mag", KAPA_LABEL_XM);
+    KapaSendLabel (kapa3, "CR N Sigma", KAPA_LABEL_YM);
+
+    graphdata.color = KapaColorByName ("black");
+    graphdata.ptype = 2;
+    graphdata.size = 0.5;
+    graphdata.style = 2;
+    KapaPrepPlot (kapa3, n, &graphdata);
+    KapaPlotVector (kapa3, n, x->data.F32, "x");
+    KapaPlotVector (kapa3, n, y->data.F32, "y");
+
+    // second section : mag vs EXT nSigma
+    section.dx = 1.0;
+    section.dy = 0.5;
+    section.x = 0.0;
+    section.y = 0.5;
+    section.name = NULL;
+    psStringAppend (&section.name, "a2");
+    KapaSetSection (kapa3, &section);
+    psFree (section.name);
+
+    graphdata.xmin = +32.0;
+    graphdata.xmax = -32.0;
+    graphdata.ymin = +32.0;
+    graphdata.ymax = -32.0;
+
+    // construct the plot vectors
+    n = 0;
+    for (int i = 0; i < sources->n; i++) {
+        pmSource *source = sources->data[i];
+        if (!source) continue;
+        if (source->type != PM_SOURCE_TYPE_STAR) continue;
+        if (!isfinite (source->extNsigma)) continue;
+
+        x->data.F32[n] = -2.5*log10(source->peak->flux);
+        y->data.F32[n] = source->extNsigma;
+        graphdata.xmin = PS_MIN(graphdata.xmin, x->data.F32[n]);
+        graphdata.xmax = PS_MAX(graphdata.xmax, x->data.F32[n]);
+        graphdata.ymin = -0.5;
+        graphdata.ymax = 10.0;
+
+        n++;
+    }
+    x->n = y->n = n;
+
+    range = graphdata.xmax - graphdata.xmin;
+    graphdata.xmax += 0.05*range;
+    graphdata.xmin -= 0.05*range;
+
+    // XXX set the plot range to match the image
+    KapaSetLimits (kapa3, &graphdata);
+
+    KapaSetFont (kapa3, "helvetica", 14);
+    KapaBox (kapa3, &graphdata);
+    KapaSendLabel (kapa3, "EXT N Sigma", KAPA_LABEL_YM);
+
+    graphdata.color = KapaColorByName ("black");
+    graphdata.ptype = 2;
+    graphdata.size = 0.5;
+    graphdata.style = 2;
+    KapaPrepPlot (kapa3, n, &graphdata);
+    KapaPlotVector (kapa3, n, x->data.F32, "x");
+    KapaPlotVector (kapa3, n, y->data.F32, "y");
+
+    psFree (x);
+    psFree (y);
+
+    // pause and wait for user input:
+    // continue, save (provide name), ??
+    char key[10];
+    fprintf (stdout, "[c]ontinue? ");
+    if (!fgets(key, 8, stdin)) {
+        psWarning("Unable to read option");
+    }
+    return true;
+}
+
 bool psphotVisualShowResidualImage (pmConfig *config, pmReadout *readout) {