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

Timestamp:

Nov 12, 2015, 4:16:35 PM (11 years ago)

Author:

watersc1

Message:

Add update sql call to set XpetRadius and XpetRadiusErr to -999 when
XpetRadius is less than -1e10. This will exclude the -1.56496e+18
-1.56496e+17 values that are showing up in the sastest. The bad
radius measurement is found in the skycal CMF, converted to janskies,
and the error set using the 10% fix.

From old batch:
gippDetectID(7) gstackImageID(9) gpetRadius(10) gpetRadiusErr(11) gpetMag(12) gpetMagErr(13) gpetR50(14)
616 3942577 -1.56496e+18 -1.56496e+17 -999 -999 -999

From CMF:
IPP_IDET(1) PETRO_RADIUS(10) PETRO_RADIUS_ERR(11) PETRO_MAG(8) PETRO_MAG_ERR(9) PETRO_RADIUS_50(12)
616 -6.25985e+18 nan nan nan nan

File:

• trunk/ippToPsps/jython/stackbatch.py (modified) (1 diff)

trunk/ippToPsps/jython/stackbatch.py

@@ -890 +890 @@
 
         # modify petRadiusErr if it is not set correctlyu
-
         sqlLine = sqlUtility("UPDATE " + tablename + " AS a SET")
         sqlLine.group("a." + filter + "petRadiusErr", "a." + filter + "petRadius * 0.1")
         sql = sqlLine.makeEquals("WHERE a." + filter + "petRadiusErr IS NULL")
 
+        self.scratchDb.execute(sql)
+
+        # modify petRadius/petRadiusErr if they have the very wrong values.
+        # Be conservative and assume all values smaller than -1e10 are very wrong values.
+        sqlLine = sqlUtility("UPDATE " + tablename + " AS a SET")
+        sqlLine.group("a." + filter + "petRadius", "-999")
+        sqlLine.group("a." + filter + "petRadiusErr", "-999")
+        sql = sqlLine.makeEquals("WHERE a." + filter + "petRadius <= -1e+10")
+        
         self.scratchDb.execute(sql)