Changeset 16498

Timestamp:

Feb 15, 2008, 9:33:37 AM (18 years ago)

Author:

eugene

Message:

allow sky rate from sky mags

Location:

trunk/ppSim/src

Files:

• ppSimMakeGalaxies.c (modified) (2 diffs)
• ppSimMakeStars.c (modified) (2 diffs)

trunk/ppSim/src/ppSimMakeGalaxies.c

@@ -32 +32 @@
 
     // XXX push these into the recipe...
-    float skyRate = psMetadataLookupF32(&mdok, config->arguments, "SKY.RATE"); // Sky rate
     float darkRate = psMetadataLookupF32(&mdok, config->arguments, "DARK.RATE"); // Dark rate
     float expTime = psMetadataLookupF32(&mdok, config->arguments, "EXPTIME"); // Exposure time
@@ -39 +38 @@
     float seeing = psMetadataLookupF32(&mdok, config->arguments, "SEEING"); // Seeing sigma (pix)
     float scale = psMetadataLookupF32(&mdok, config->arguments, "SCALE") * M_PI / 3600.0 / 180.0; // Plate scale (radians/pixel)
+
+    float skyRate = psMetadataLookupF32(&mdok, config->arguments, "SKY.RATE"); // Sky rate
+    if (isnan(skyRate)) {
+        float skyMags = psMetadataLookupF32(&mdok, config->arguments, "SKY.MAGS");  assert (mdok);
+        skyRate = scale * scale * pow (10.0, -0.4*(skyMags - zp));
+    }
 
     if (galaxyDensity <= 0) return NULL;

trunk/ppSim/src/ppSimMakeStars.c

@@ -18 +18 @@
     float starsDensity = psMetadataLookupF32(NULL, config->arguments, "STARS.DENSITY"); // Density of fakes
 
-    float skyRate = psMetadataLookupF32(NULL, config->arguments, "SKY.RATE"); // Sky rate
     float darkRate = psMetadataLookupF32(NULL, config->arguments, "DARK.RATE"); // Dark rate
     float expTime = psMetadataLookupF32(NULL, config->arguments, "EXPTIME"); // Exposure time
@@ -25 +24 @@
     float seeing = psMetadataLookupF32(NULL, config->arguments, "SEEING"); // Seeing sigma (pix)
     float scale = psMetadataLookupF32(NULL, config->arguments, "SCALE") * M_PI / 3600.0 / 180.0; // Plate scale (radians/pixel)
+
+    float skyRate = psMetadataLookupF32(&mdok, config->arguments, "SKY.RATE"); // Sky rate
+    if (isnan(skyRate)) {
+        float skyMags = psMetadataLookupF32(&mdok, config->arguments, "SKY.MAGS");  assert (mdok);
+        skyRate = scale * scale * pow (10.0, -0.4*(skyMags - zp));
+    }
 
     if (starsDensity <= 0) return true;