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

Timestamp:

Apr 28, 2013, 2:57:07 PM (13 years ago)

Author:

eugene

Message:

psfFlux, psfFluxErr, apFlux from FLUX fields in smf, not MAG fields; set flux <= 0.0 to 1e20

File:

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

trunk/ippToPsps/jython/stackbatch.py

@@ -482 +482 @@
                ,X_PSF_SIG \
                ,Y_PSF_SIG \
-               ,POW(10.0, (-0.4*PSF_INST_MAG)) / " + str(self.expTime) + " \
-               ,ABS((PSF_INST_MAG_SIG*(POW(10.0, (-0.4*PSF_INST_MAG)) / " + str(self.expTime) + ")) / 1.085736) \
+               ,PSF_INST_FLUX / " + str(self.expTime) + " \
+               ,PSF_INST_FLUX_SIG / " + str(self.expTime) + " \
                ,POW(10.0, (-0.4*PEAK_FLUX_AS_MAG)) / " + str(self.expTime) + " \
                ,SKY / " + str(self.expTime) + " \
@@ -504 +504 @@
                ,MOMENTS_R1 \
                ,MOMENTS_RH \
-               ,POW(10.0, -0.4*AP_MAG) / " + str(self.expTime) + " \
+               ,AP_FLUX / " + str(self.expTime) + " \
                ,NULL \
                ,KRON_FLUX / " + str(self.expTime) + " \
@@ -515 +515 @@
                FROM SkyChip_psf"
 
+# these were used above (left over from before PSF_INST_FLUX, PSF_INST_FLUX_SIG were available?
+#               ,POW(10.0, (-0.4*PSF_INST_MAG)) / " + str(self.expTime) + " 
+#               ,ABS((PSF_INST_MAG_SIG*(POW(10.0, (-0.4*PSF_INST_MAG)) / " + str(self.expTime) + ")) / 1.085736) 
+#               ,POW(10.0, -0.4*AP_MAG) / " + str(self.expTime) + " 
+
         # print "sql: ", sql
         # response = raw_input("ready to insert stack det ")
@@ -526 +531 @@
         self.logger.infoPair("Deleting", "entries with StackDetection.objID = 0")
         
+        # XXX EAM : this seems quite inefficient : these commands use updates to set 
+        # fields which are constant across all rows.  Why not just include these in 
+        # the original insert above?
+
         self.scratchDb.updateAllRows("StackDetection", "surveyID", str(self.surveyID))
         
@@ -560 +569 @@
         # response = raw_input("add primary key? ")
   
-        if self.stackType == "DEEP_STACK": 
+        ### XXX this should not be a DEEP_STACK only case...
+        if false and self.stackType == "DEEP_STACK": 
 
             #if deep stack and instFlux = null and err not null
@@ -572 +582 @@
             # response = raw_input("add psf flux ")
             
+        # NOTE : Flux limits : in the current PSPS schema, negative fluxes 
+        # cause problems for sql queries which work in mags as SQL cannot do
+        # something like (f < 0.0) ? -999 : -2.5*log10(f)
+        # as a result, the negative fluxes here result in floating point errors
+        sql = "UPDATE StackDetection SET psfFlux = 1e20 WHERE psfFlux <= 0.0"
+        self.scratchDb.execute(sql)
+        
+        sql = "UPDATE StackDetection SET apFlux = 1e20 WHERE apFlux <= 0.0"
+        self.scratchDb.execute(sql)
+
+        sql = "UPDATE StackDetection SET kronFlux = 1e20 WHERE kronFlux <= 0.0"
+        self.scratchDb.execute(sql)
+
         #leave null instflux in
         #self.scratchDb.reportAndDeleteRowsWithNULLS("StackDetection", "instFlux")