Changeset 10880 for trunk/psastro/src/psastroDemoDump.c

Timestamp:

Jan 2, 2007, 3:26:52 PM (20 years ago)

Author:

eugene

Message:

fixed mosaic analysis; cleaned up test code (some in traces, some removed), general output cleanups

File:

• trunk/psastro/src/psastroDemoDump.c (modified) (1 diff)

trunk/psastro/src/psastroDemoDump.c

@@ -58 +58 @@
     return true;
 }
+
+bool psastroDumpMatches (pmFPA *fpa, char *filename) {
+
+    pmChip *chip = NULL;
+    pmCell *cell = NULL;
+    pmReadout *readout = NULL;
+    pmFPAview *view = pmFPAviewAlloc (0);
+
+    FILE *f = fopen (filename, "w");
+
+    // this loop selects the matched stars for all chips
+    while ((chip = pmFPAviewNextChip (view, fpa, 1)) != NULL) {
+        psTrace ("psastro", 4, "Chip %d: %x %x\n", view->chip, chip->file_exists, chip->process);
+        if (!chip->process || !chip->file_exists) continue;
+        
+        while ((cell = pmFPAviewNextCell (view, fpa, 1)) != NULL) {
+            psTrace ("psastro", 4, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
+            if (!cell->process || !cell->file_exists) continue;
+
+            // process each of the readouts
+            // XXX there can only be one readout per chip, right?
+            while ((readout = pmFPAviewNextReadout (view, fpa, 1)) != NULL) {
+                if (! readout->data_exists) continue;
+
+                // select the raw objects for this readout
+                psArray *rawstars = psMetadataLookupPtr (NULL, readout->analysis, "PSASTRO.RAWSTARS");
+                if (rawstars == NULL) continue;
+
+                // select the raw objects for this readout
+                psArray *refstars = psMetadataLookupPtr (NULL, readout->analysis, "PSASTRO.REFSTARS");
+                if (refstars == NULL) continue;
+                psTrace ("psastro", 4, "Trying %ld refstars\n", refstars->n);
+
+                psArray *matches = psMetadataLookupPtr (NULL, readout->analysis, "PSASTRO.MATCH");
+                if (matches == NULL) continue;
+
+                for (int i = 0; i < matches->n; i++) {
+                    pmAstromMatch *match = matches->data[i];
+
+                    pmAstromObj *raw = rawstars->data[match->raw];
+                    fprintf (f, "%f %f  %f %f  %f %f  %f %f  %f   |   ",  
+                             DEG_RAD*raw->sky->r, DEG_RAD*raw->sky->d, 
+                             raw->TP->x, raw->TP->y, 
+                             raw->FP->x, raw->FP->y, 
+                             raw->chip->x, raw->chip->y, raw->Mag);
+
+                    pmAstromObj *ref = refstars->data[match->ref];
+                    fprintf (f, "%f %f  %f %f  %f %f  %f %f  %f\n", 
+                             DEG_RAD*ref->sky->r, DEG_RAD*ref->sky->d, 
+                             ref->TP->x, ref->TP->y, 
+                             ref->FP->x, ref->FP->y, 
+                             ref->chip->x, ref->chip->y, ref->Mag);
+                }
+            }
+        }
+    }
+    fclose (f);
+    psFree (view);
+    return true;
+}