Changeset 36934 for branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/mkcmf.c

Timestamp:

Jun 22, 2014, 5:21:17 PM (12 years ago)

Author:

eugene

Message:

handle PS1_SV3 format in addstar

File:

• branches/eam_branches/ipp-20140610/Ohana/src/addstar/src/mkcmf.c (modified) (6 diffs)

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

@@ -21 +21 @@
 void writeStars_PS1_V2 (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars);
 void writeStars_PS1_V1 (FTable *ftable, double *X, double *Y, double *M, int Nstars);
+void writeStars_PS1_SV3 (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars);
 void writeStars_PS1_DEV_1 (FTable *ftable, double *X, double *Y, double *M, int Nstars);
 void writeStars_PS1_DEV_0 (FTable *ftable, double *X, double *Y, double *M, int Nstars);
@@ -339 +340 @@
   if (!strcmp(type, "PS1_V5_Lensing")) {
     writeStars_PS1_V5_Lensing (&ftable, X, Y, M, Flag, Nstars); 
+    found = TRUE;
+  }
+  if (!strcmp(type, "PS1_SV3")) {
+    writeStars_PS1_SV3 (&ftable, X, Y, M, Flag, Nstars); 
     found = TRUE;
   }
@@ -885 +890 @@
 }
 
-void writeStars_PS1_V5_Lensing (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars) {
+void writeStars_PS1_SV3 (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars) {
 
   int i;
-  CMF_PS1_V5_Lensing *stars;
+  CMF_PS1_SV3 *stars;
   float flux, fSN;
 
   // XXX add gaussian-distributed noise based on counts
   // this needs to make different output 'stars' entries depending on the desired type
-  ALLOCATE (stars, CMF_PS1_V5_Lensing, Nstars);
+  ALLOCATE (stars, CMF_PS1_SV3, Nstars);
   for (i = 0; i < Nstars; i++) {
 
@@ -964 +969 @@
     stars[i].Mrh       = FX;
 
-    stars[i].X11_sm_obj = FX;
-    stars[i].X12_sm_obj = FX;
-    stars[i].X22_sm_obj = FY;
-    stars[i].E1_sm_obj  = FX;
-    stars[i].E2_sm_obj  = FX;
-
-    stars[i].X11_sh_obj = FX;
-    stars[i].X12_sh_obj = 0.2;
-    stars[i].X22_sh_obj = FY;
-    stars[i].E1_sh_obj  = FX*0.9;
-    stars[i].E2_sh_obj  = FX*0.8;
-
-    stars[i].X11_sm_psf = FX;
-    stars[i].X12_sm_psf = FX;
-    stars[i].X22_sm_psf = FX;
-    stars[i].E1_sm_psf  = FX;
-    stars[i].E2_sm_psf  = FX;
-
-    stars[i].X11_sh_psf = FX;
-    stars[i].X12_sh_psf = FX;
-    stars[i].X22_sh_psf = FX;
-    stars[i].E1_sh_psf  = FX;
-    stars[i].E2_sh_psf  = FX;
-
     stars[i].kronFlux  = flux * 1.25;
     stars[i].kronFluxErr = fSN * flux * 1.25;
@@ -994 +975 @@
     stars[i].kronOuter = fSN * flux * 1.5;
 
-    stars[i].skyLimitRad = 1;
-    stars[i].skyLimitFlux = 2;
-    stars[i].skyLimitSlope = 0.2;
+    // stars[i].skyLimitRad = 1;
+    // stars[i].skyLimitFlux = 2;
+    // stars[i].skyLimitSlope = 0.2;
 
     stars[i].flags     = Flag[i];
@@ -1014 +995 @@
   }
 
+  gfits_table_set_CMF_PS1_SV3 (ftable, stars, Nstars);
+  gfits_modify (ftable->header, "EXTTYPE",   "%s", 1, "PS1_SV3");
+}
+
+void writeStars_PS1_V5_Lensing (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars) {
+
+  int i;
+  CMF_PS1_V5_Lensing *stars;
+  float flux, fSN;
+
+  // XXX add gaussian-distributed noise based on counts
+  // this needs to make different output 'stars' entries depending on the desired type
+  ALLOCATE (stars, CMF_PS1_V5_Lensing, Nstars);
+  for (i = 0; i < Nstars; i++) {
+
+    flux = pow (10.0, -0.4*M[i]);
+    fSN = 1.0 / sqrt(flux);
+
+    stars[i].detID = i;
+    stars[i].X = X[i];
+    stars[i].Y = Y[i];
+    stars[i].dX = FX * fSN;
+    stars[i].dY = FY * fSN;
+
+    stars[i].posangle = 10.0;
+    stars[i].pltscale = 0.25;
+
+    stars[i].M = M[i];
+    stars[i].dM = fSN;
+
+    stars[i].Flux = flux;
+    stars[i].dFlux = flux * fSN;
+
+    stars[i].Map = M[i] - 0.05;
+    stars[i].MapRaw = M[i] - 0.10;
+
+    stars[i].apRadius = 8.0;
+
+    stars[i].apFlux = pow(10.0, -0.4*stars[i].Map);
+    stars[i].apFluxErr = stars[i].apFlux * fSN;
+
+    stars[i].Mcalib = M[i] + ZERO_POINT + 2.5*log10(exptime);
+    stars[i].dMcal = 0.05;
+
+    XY_to_RD (&stars[i].RA, &stars[i].DEC, X[i], Y[i], &coords);
+    stars[i].apNpix = 3.14*8.0*8.0;
+
+    stars[i].Mpeak     = M[i] + 1.0;
+    stars[i].sky       = SKY;
+    stars[i].dSky      = DSKY;
+    stars[i].psfChisq  = PSFCHI;
+    stars[i].crNsigma  = CRN;
+    stars[i].extNsigma = EXTN;
+    stars[i].fx        = FX;
+    stars[i].fy        = FY;
+    stars[i].df        = DF;
+
+    stars[i].k         = 1.0;
+    stars[i].fwhmMaj   = FX*2.8;
+    stars[i].fwhmMin   = FY*2.8;
+
+    // randomly give poor PSFQF values
+    if ((BAD_PSFQF_FRAC > 0.0) && (drand48() < BAD_PSFQF_FRAC)) {
+      stars[i].psfQF     = 0.25;
+      stars[i].psfQFperf = 0.24;
+    } else {
+      stars[i].psfQF     = PSFQUAL;
+      stars[i].psfQFperf = MAX(PSFQUAL - 0.01, 0.0);
+    }
+
+    stars[i].psfNdof   = 1;
+    stars[i].psfNpix   = 2;
+
+    stars[i].Mxx       = FX;
+    stars[i].Mxy       = 0.01;
+    stars[i].Myy       = FX;
+    stars[i].M3c       = FX;
+    stars[i].M3s       = FX;
+    stars[i].M4c       = FX;
+    stars[i].M4s       = FX;
+    stars[i].Mr1       = FX;
+    stars[i].Mrh       = FX;
+
+    stars[i].X11_sm_obj = FX;
+    stars[i].X12_sm_obj = FX;
+    stars[i].X22_sm_obj = FY;
+    stars[i].E1_sm_obj  = FX;
+    stars[i].E2_sm_obj  = FX;
+
+    stars[i].X11_sh_obj = FX;
+    stars[i].X12_sh_obj = 0.2;
+    stars[i].X22_sh_obj = FY;
+    stars[i].E1_sh_obj  = FX*0.9;
+    stars[i].E2_sh_obj  = FX*0.8;
+
+    stars[i].X11_sm_psf = FX;
+    stars[i].X12_sm_psf = FX;
+    stars[i].X22_sm_psf = FX;
+    stars[i].E1_sm_psf  = FX;
+    stars[i].E2_sm_psf  = FX;
+
+    stars[i].X11_sh_psf = FX;
+    stars[i].X12_sh_psf = FX;
+    stars[i].X22_sh_psf = FX;
+    stars[i].E1_sh_psf  = FX;
+    stars[i].E2_sh_psf  = FX;
+
+    stars[i].kronFlux  = flux * 1.25;
+    stars[i].kronFluxErr = fSN * flux * 1.25;
+
+    stars[i].kronInner = fSN * flux * 0.9;
+    stars[i].kronOuter = fSN * flux * 1.5;
+
+    stars[i].skyLimitRad = 1;
+    stars[i].skyLimitFlux = 2;
+    stars[i].skyLimitSlope = 0.2;
+
+    stars[i].flags     = Flag[i];
+    stars[i].flags2    = 0x80;
+
+    stars[i].nFrames   = 1;
+
+    if (ADDNOISE) {
+      stars[i].X += FX * fSN * rnd_gauss(0.0, 1.0);
+      stars[i].Y += FY * fSN * rnd_gauss(0.0, 1.0);
+      float Moff = fSN*rnd_gauss(0.0, 1.0);
+      stars[i].M += Moff;
+      stars[i].Map += Moff;
+      stars[i].MapRaw += Moff;
+      stars[i].Mcalib += Moff;
+    }
+  }
+
   gfits_table_set_CMF_PS1_V5_Lensing (ftable, stars, Nstars);
   gfits_modify (ftable->header, "EXTTYPE",   "%s", 1, "PS1_V5");