Changeset 10880 for trunk/psastro/src/psastroMosaicSetAstrom.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/psastroMosaicSetAstrom.c (modified) (1 diff)

trunk/psastro/src/psastroMosaicSetAstrom.c

@@ -47 +47 @@
         }
     }
+    psFree (view);
     return true;
 }
-
-bool psastroMosaicSetAstrom_tmp (pmFPA *fpa) {
-
-    pmChip *chip = NULL;
-    pmCell *cell = NULL;
-    pmReadout *readout = NULL;
-    pmFPAview *view = pmFPAviewAlloc (0);
-
-    FILE *f1 = fopen ("raw.tmp.dat", "w");
-    FILE *f2 = fopen ("ref.tmp.dat", "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; }
-
-                for (int i = 0; i < rawstars->n; i++) {
-                    pmAstromObj *raw = rawstars->data[i];
-        
-                    fprintf (f1, "%f %f  %f %f  %f %f  ", raw->sky->r, raw->sky->d, raw->TP->x, raw->TP->y, raw->FP->x, raw->FP->y);
-
-                    psPlaneTransformApply (raw->FP, chip->toFPA, raw->chip);
-                    psPlaneTransformApply (raw->TP, fpa->toTPA, raw->FP);
-                    psDeproject (raw->sky, raw->TP, fpa->toSky);
-
-                    fprintf (f1, "|  %f %f  %f %f  %f %f   %f\n", raw->sky->r, raw->sky->d, raw->TP->x, raw->TP->y, raw->FP->x, raw->FP->y, raw->Mag);
-                }
-
-                psArray *refstars = psMetadataLookupPtr (NULL, readout->analysis, "PSASTRO.REFSTARS");
-                if (refstars == NULL) { continue; }
-
-                for (int i = 0; i < refstars->n; i++) {
-                    pmAstromObj *ref = refstars->data[i];
-        
-                    fprintf (f2, "%f %f  %f %f  %f %f  ", ref->sky->r, ref->sky->d, ref->TP->x, ref->TP->y, ref->FP->x, ref->FP->y);
-
-                    psProject (ref->TP, ref->sky, fpa->toSky);
-                    psPlaneTransformApply (ref->FP, fpa->fromTPA, ref->TP);
-                    psPlaneTransformApply (ref->chip, chip->fromFPA, ref->FP);
-
-                    fprintf (f2, "|  %f %f  %f %f  %f %f   %f\n", ref->sky->r, ref->sky->d, ref->TP->x, ref->TP->y, ref->FP->x, ref->FP->y, ref->Mag);
-                }
-            }
-        }
-    }
-    fclose (f1);
-    fclose (f2);
-    return true;
-}