Changeset 28939

Timestamp:

Aug 17, 2010, 10:42:55 AM (16 years ago)

Author:

bills

Message:

minor fixes to addstar for sdss files

File:

• trunk/Ohana/src/addstar/src/ReadStarsSDSS.c (modified) (3 diffs)

trunk/Ohana/src/addstar/src/ReadStarsSDSS.c

@@ -11 +11 @@
   NAME = (TYPE *) gfits_get_bintable_column_data (table.header, &table, #NAME, type, &Nrow, &Ncol); \
   assert (NAME); assert (!strcmp (type, #TYPE)); assert (Nrow == Nstars); assert (Ncol == NFILTER);
+
+# define GET_COLUMN_5_NOASSERT(NAME,TYPE) \
+  TYPE *NAME; \
+  NAME = (TYPE *) gfits_get_bintable_column_data (table.header, &table, #NAME, type, &Nrow, &Ncol);
 
 # define GET_COLUMN_1(NAME,TYPE) \
@@ -123 +127 @@
   GET_COLUMN_5 (flags2, int);
 
+#ifdef notyet
   GET_COLUMN_5 (prob_psf, float);
+#else
+  GET_COLUMN_5_NOASSERT (prob_psf, float);
+#endif
 
   GET_COLUMN_1 (ra, double);
@@ -160 +168 @@
       stars[N].measure.FWx       = ShortPixels(seeing[j]); // reported in arcsec?
       stars[N].measure.FWy       = ShortPixels(seeing[j]);
-      stars[N].measure.psfChisq  = prob_psf[N]; // XXX not really the correct value...
+      if (prob_psf) {
+          stars[N].measure.psfChisq  = prob_psf[N]; // XXX not really the correct value...
+      } else {
+          stars[N].measure.psfChisq  = NAN;
+      }
       stars[N].measure.detID     = N;
       stars[N].measure.t         = tzero[j] + clockRate*rowc[N]; // time since row 0
-      stars[N].measure.dt        = 53.907456; // is this 2048*clockRate ?
+      stars[N].measure.dt        = 4.32912209; // 2.5 * log(53.907456) the sdss exposure time // old comment is 53907456 is this 2048*clockRate ?
 
       SetSDSSFlags (&stars[N], flags[N], flags2[N]);