Changeset 18004

Timestamp:

Jun 8, 2008, 2:47:14 PM (18 years ago)

Author:

eugene

Message:

some temporary test outputs

File:

• trunk/psphot/src/psphotReadout.c (modified) (3 diffs)

trunk/psphot/src/psphotReadout.c

@@ -17 +17 @@
 
     // set the photcode for this image
-    if (!psphotAddPhotcode (recipe, config, view)) {
+    if (!psphotAddPhotcode (recipe, config, view, "PSPHOT.INPUT")) {
         psError (PSPHOT_ERR_CONFIG, false, "trouble defining the photcode");
         return false;
@@ -118 +118 @@
 
     if (0) { 
-        FILE *out = fopen ("out.dat", "w");
+        FILE *out = fopen ("out.pass1.dat", "w");
         
         for (int i = 0; i < sources->n; i++) {
@@ -194 +194 @@
     psphotFitSourcesLinear (readout, sources, recipe, psf, TRUE);
 
+    if (0) { 
+        FILE *out = fopen ("out.pass2.dat", "w");
+        
+        for (int i = 0; i < sources->n; i++) {
+            pmSource *source = sources->data[i];
+            if (!source) continue;
+            pmPeak *peak = source->peak;
+            if (!peak) continue;
+            pmModel *model = source->modelPSF;
+            if (!model) continue;
+            if (!model->params) continue;
+            
+            fprintf (out, "%d %f %f  %f %f  %f %f\n", 
+                     i, peak->xf, peak->yf, peak->flux, peak->SN,
+                     model->params->data.F32[PM_PAR_I0],
+                     model->dparams->data.F32[PM_PAR_I0]);
+        }
+        fclose (out);
+    }
+
 pass1finish: