Changeset 9587

Timestamp:

Oct 16, 2006, 2:51:46 PM (20 years ago)

Author:

eugene

Message:

adding refstar subset operation using magnitude limits

Location:

trunk/psastro/src

Files:

• Makefile.am (modified) (2 diffs)
• psastro.c (modified) (2 diffs)
• psastro.h (modified) (1 diff)
• psastroChooseRefstars.c (added)
• psastroLoadReferences.c (deleted)
• psastroLoadRefstars.c (added)
• psastroLuminosityFunction.c (added)
• psastroMosaicGetRefstars.c (deleted)
• psastroRefstarSubset.c (added)

trunk/psastro/src/Makefile.am

@@ -17 +17 @@
 psastroDataSave.c           \
 psastroAstromGuess.c        \
-psastroLoadReferences.c     \
+psastroLoadRefstars.c     \
+psastroChooseRefstars.c  \
 psastroConvert.c            \
 psastroChipAstrom.c         \
@@ -23 +24 @@
 psastroUtils.c              \
 psastroTestFuncs.c          \
+psastroLuminosityFunction.c \
+psastroRefstarSubset.c      \
 psastroMosaicAstrom.c       \
 psastroMosaicGetGrads.c     \
-psastroMosaicGetRefstars.c  \
 psastroMosaicChipAstrom.c   \
 psastroMosaicSetAstrom.c    \

trunk/psastro/src/psastro.c

@@ -39 +39 @@
 
     // load the reference stars overlapping the data stars
-    psArray *refs = psastroLoadReferences (config);
+    psArray *refs = psastroLoadRefstars (config);
     if (!refs) {
         psErrorStackPrint(stderr, "failed to load reference data");
@@ -45 +45 @@
     }
 
-    // XXX ?? what does this do ??
-    psastroMosaicGetRefstars (config, refs);
+    // choose reference stars corresponding to the selected chips
+    if (psastroChooseRefstars (config, refs)) {
+        psErrorStackPrint(stderr, "failed to select reference data for chips");
+        exit (1);
+    }
 
     char *mosastro = psMetadataLookupStr (NULL, config->arguments, "MOSASTRO");

trunk/psastro/src/psastro.h

@@ -20 +20 @@
 # define toTPA toTangentPlane
 # define toSky projection
+
+// this structure represents a fit to the logN / logS curve for a set of stars
+// logN = offset + slope * logS 
+typedef struct {
+  double mMin;
+  double mMax;
+  double offset;
+  double slope;
+} pmLumFunc;
 
 pmConfig         *psastroArguments (int argc, char **argv);