Changeset 33376

Timestamp:

Feb 27, 2012, 11:38:59 AM (14 years ago)

Author:

eugene

Message:

relphot tests

Location:

branches/eam_branches/ipp-20111122/Ohana/src/addstar

Files:

• src/mkcmf.c (modified) (10 diffs)
• test/relphot.reject.dvo (added)

branches/eam_branches/ipp-20111122/Ohana/src/addstar/src/mkcmf.c

@@ -14 +14 @@
 void gauss_init (int Nbin);
 double rnd_gauss (double mean, double sigma);
-void writeStars_PS1_V2 (FTable *ftable, double *X, double *Y, double *M, int Nstars);
+void writeStars_PS1_V2 (FTable *ftable, double *X, double *Y, double *M, int *Flag, int Nstars);
 void writeStars_PS1_V1 (FTable *ftable, double *X, double *Y, double *M, int Nstars);
 void writeStars_PS1_DEV_1 (FTable *ftable, double *X, double *Y, double *M, int Nstars);
@@ -20 +20 @@
 
 int ADDNOISE = TRUE;
+float BAD_PSFQF_FRAC = 0.0;
 
 int main (int argc, char **argv) {
@@ -26 +27 @@
   // load a text table with X,Y,Mag (instrumental?)
 
-  int N, Nstars, NSTARS, found;
+  int N, Nstars, NSTARS, found, *Flag;
   double *X, *Y, *M, Xmax, Ymax;
 
@@ -125 +126 @@
   }
 
+  // expect a field with the photom flags for each detection
+  int READ_FLAGS = FALSE;
+  if ((N = get_argument (argc, argv, "-flags"))) {
+    remove_argument (N, &argc, argv);
+    READ_FLAGS = TRUE;
+  }
+
   // add support for all cmf types
   ADDNOISE = TRUE;
@@ -130 +138 @@
     remove_argument (N, &argc, argv);
     ADDNOISE = FALSE;
+  }
+
+  // random bad PSF_QF values
+  if ((N = get_argument (argc, argv, "-bad-psfqf-frac"))) {
+    remove_argument (N, &argc, argv);
+    BAD_PSFQF_FRAC = atof(argv[N]);
+    remove_argument (N, &argc, argv);
   }
 
@@ -166 +181 @@
   ALLOCATE (Y, double, NSTARS);
   ALLOCATE (M, double, NSTARS);
+  ALLOCATE (Flag, int, NSTARS);
 
   Xmax = Ymax = 0;
   while (1) {
     int status;
+    double ra, dec, mag, xobs, yobs, xraw, yraw, mraw;
+    int flags;
     if (FROM_COORDS) {
-      double ra, dec, mag, xraw, yraw, mraw;
       status = fscanf (f, "%lf %lf %lf %lf %lf %lf", &ra, &dec, &mag, &xraw, &yraw, &mraw);
-      if (status != EOF) {
-        RD_to_XY (&X[Nstars], &Y[Nstars], ra, dec, &coords);
+      RD_to_XY (&xobs, &yobs, ra, dec, &coords);
+    } else {
+      status = fscanf (f, "%lf %lf %lf", &xobs, &yobs, &mraw);
+    }
+    if (status == EOF) break;
+    if (READ_FLAGS) {
+      status = fscanf (f, "%d", &flags);
+      if (status == EOF) {
+          fprintf (stderr, "error: missing flag for last star?\n");
+          exit (1);
       }
-      M[Nstars] = mraw;
+      Flag[Nstars] = mraw;
     } else {
-      status = fscanf (f, "%lf %lf %lf", &X[Nstars], &Y[Nstars], &M[Nstars]);
-    }
-    if (status == EOF) break;
-    
+      Flag[Nstars] = 0;
+    }
+
+    X[Nstars] = xobs;
+    Y[Nstars] = yobs;
+    M[Nstars] = mraw;
+
     Xmax  = MAX(Xmax, X[Nstars]);
     Ymax  = MAX(Ymax, Y[Nstars]);
@@ -190 +218 @@
       REALLOCATE (Y, double, NSTARS);
       REALLOCATE (M, double, NSTARS);
+      REALLOCATE (Flag, int, NSTARS);
     }
     Nstars ++;
@@ -488 +517 @@
 }
 
-void writeStars_PS1_V2 (FTable *ftable, double *X, double *Y, double *M, int Nstars) {
+void writeStars_PS1_V2 (FTable *ftable, double *X, double *Y, double *M, int *Flag, int Nstars) {
 
   int i;
@@ -514 +543 @@
     }
 
+    // randomly give poor PSFQF values
+    if ((BAD_PSFQF_FRAC > 0.0) && (drand48() < BAD_PSFQF_FRAC*Nstars)) {
+      stars[i].psfQual   = 0.25;
+    } else {
+      stars[i].psfQual   = PSFQUAL;
+    }
+    
     stars[i].dX = FX * fSN;
     stars[i].dY = FY * fSN;
@@ -527 +563 @@
     stars[i].fy        = FY;
     stars[i].df        = DF;
-    stars[i].psfQual   = PSFQUAL;
     stars[i].nFrames   = 1;
-    stars[i].flags     = FLAGS;
+    stars[i].flags     = Flag[i];
   }