Changeset 16496

Timestamp:

Feb 15, 2008, 9:32:14 AM (18 years ago)

Author:

eugene

Message:

allow skyRate to come from mags/square asec

Location:

trunk/ppSim/src

Files:

• ppSimInsertGalaxies.c (modified) (2 diffs)
• ppSimInsertStars.c (modified) (2 diffs)

trunk/ppSim/src/ppSimInsertGalaxies.c

@@ -23 +23 @@
 
     float expTime = psMetadataLookupF32(NULL, config->arguments, "EXPTIME"); // Exposure time
-    float skyRate = psMetadataLookupF32(NULL, config->arguments, "SKY.RATE"); // Sky rate
     float darkRate = psMetadataLookupF32(NULL, config->arguments, "DARK.RATE"); // Dark rate
     float readnoise = psMetadataLookupF32(NULL, cell->concepts, "CELL.READNOISE");// CCD read noise, e
@@ -35 +34 @@
     }
 
+    float skyRate = psMetadataLookupF32(NULL, config->arguments, "SKY.RATE"); // Sky rate
+    if (isnan(skyRate)) {
+        float zp      = psMetadataLookupF32(&mdok, config->arguments, "ZEROPOINT"); assert (mdok);
+        float scale   = psMetadataLookupF32(&mdok, config->arguments, "SCALE");     assert (mdok);
+        float skyMags = psMetadataLookupF32(&mdok, config->arguments, "SKY.MAGS");  assert (mdok);
+        skyRate = scale * scale * pow (10.0, -0.4*(skyMags - zp));
+    }
+    
     // Rough noise estimate, appropriate for entire cell (use for source radius?)
     float roughNoise = sqrtf(PS_SQR(readnoise) + (darkRate + skyRate) * expTime);

trunk/ppSim/src/ppSimInsertStars.c

@@ -28 +28 @@
 
     float expTime = psMetadataLookupF32(NULL, config->arguments, "EXPTIME"); // Exposure time
-    float skyRate = psMetadataLookupF32(NULL, config->arguments, "SKY.RATE"); // Sky rate
     float darkRate = psMetadataLookupF32(NULL, config->arguments, "DARK.RATE"); // Dark rate
     float readnoise = psMetadataLookupF32(NULL, cell->concepts, "CELL.READNOISE");// CCD read noise, e
@@ -38 +37 @@
             return false;
         }
+    }
+
+    float skyRate = psMetadataLookupF32(NULL, config->arguments, "SKY.RATE"); // Sky rate
+    if (isnan(skyRate)) {
+        float zp      = psMetadataLookupF32(&mdok, config->arguments, "ZEROPOINT"); assert (mdok);
+        float scale   = psMetadataLookupF32(&mdok, config->arguments, "SCALE");     assert (mdok);
+        float skyMags = psMetadataLookupF32(&mdok, config->arguments, "SKY.MAGS");  assert (mdok);
+        skyRate = scale * scale * pow (10.0, -0.4*(skyMags - zp));
     }