Changeset 39097 for trunk/ippToPsps/jython/stackbatch.py

Timestamp:

Nov 10, 2015, 3:09:38 PM (11 years ago)

Author:

watersc1

Message:

DetEff/ detectionThreshold/ nInjected fixes.

File:

• trunk/ippToPsps/jython/stackbatch.py (modified) (3 diffs)

trunk/ippToPsps/jython/stackbatch.py

@@ -265 +265 @@
         
         # Convert detectionThreshold to appropriate magnitudes
-        detectionThreshold = self.getKeyFloat(header, "%.8f", 'DETEFF.MAGREF')
+        # Grab this from the appropriate extension.
+        deteffHeader = self.fits[filter].findAndReadHeader("SkyChip.deteff", self.config.test)
+        if not deteffHeader:
+            if self.config.test: self.logger.errorPair("No header found for Skychip.deteff for", filter)
+            return False
+
+        detectionThreshold = self.getKeyFloat(deteffHeader, "%.8f", 'DETEFF.MAGREF')
         expTime = self.getKeyFloat(header, "%.8f", 'EXPTIME')
         zp      = self.getKeyFloat(header, "%.8f", 'ZPT_OBS')
@@ -1168 +1174 @@
 
             header = self.headerSet[filter]
+
+            # Convert detectionThreshold to appropriate magnitudes
+            # Grab this from the appropriate extension.
+            deteffHeader = self.fits[filter].findAndReadHeader("SkyChip.deteff", self.config.test)
+            if not deteffHeader:
+                if self.config.test: self.logger.errorPair("No header found for Skychip.deteff for", filter)
+                return False
+
+            detectionThreshold = self.getKeyFloat(deteffHeader, "%.8f", 'DETEFF.MAGREF')
+            expTime = self.getKeyFloat(header, "%.8f", 'EXPTIME')
+            zp      = self.getKeyFloat(header, "%.8f", 'ZPT_OBS')
+            zpErr   = self.getKeyFloat(header, "%.8f", 'ZPT_ERR')
+
+        # zp correction should comes from DVO (unless image is not in DVO)
+            zpImage = self.scratchDb.getImageZeroPoint(stackID)
+            if zpImage < 0.0:
+                zpImage = zp
+
+            detectionThreshold = detectionThreshold + zpImage + 2.5 * math.log10(expTime)
+            nInjected = int(self.getKeyInt(deteffHeader, 0, 'DETEFF.NUM'))
+
             zpImage = self.scratchDb.getImageZeroPoint(stackImageID)
             if zpImage < 0.0:
@@ -1176 +1203 @@
             sqlLine.group("stackImageID", str(stackImageID))
             sqlLine.group("magref",       magref)
+            sqlLine.group("nInjected",    nInjected)
             sql = sqlLine.make(") VALUES ( ", ")")