Changeset 38802
- Timestamp:
- Oct 6, 2015, 4:13:12 PM (11 years ago)
- Location:
- trunk/ippToPsps/jython
- Files:
-
- 4 edited
-
detectionbatch.py (modified) (3 diffs)
-
diffbatch.py (modified) (4 diffs)
-
forcedwarpbatch.py (modified) (3 diffs)
-
stackbatch.py (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippToPsps/jython/detectionbatch.py
r38799 r38802 410 410 results['SATDET'] = 0 411 411 412 ex tTimeString = str(self.header['EXPTIME'])412 expTimeString = str(self.header['EXPTIME']) 413 413 414 414 # insert all detections into table : do NOT include the IGNORE unless we have to … … 444 444 sqlLine.group("decErr", "Y_PSF_SIG * 0.257") 445 445 sqlLine.group("extNSigma", "EXT_NSIGMA") 446 sqlLine.group("psfFlux", "PSF_INST_FLUX / " + ex tTimeString)447 sqlLine.group("psfFluxErr", "PSF_INST_FLUX_SIG / " + ex tTimeString)446 sqlLine.group("psfFlux", "PSF_INST_FLUX / " + expTimeString) 447 sqlLine.group("psfFluxErr", "PSF_INST_FLUX_SIG / " + expTimeString) 448 448 if (self.id >= 982483): 449 449 sqlLine.group("psfMajorFWHM", "PSF_FWHM_MAJ * abs(PLTSCALE)") … … 465 465 sqlLine.group("momentM4S", "MOMENTS_M4S * PLTSCALE * PLTSCALE") 466 466 if (self.id >= 982483): 467 sqlLine.group("apFlux", "AP_FLUX / " + ex tTimeString)468 sqlLine.group("apFluxErr", "AP_FLUX_SIG / " + ex tTimeString)467 sqlLine.group("apFlux", "AP_FLUX / " + expTimeString) 468 sqlLine.group("apFluxErr", "AP_FLUX_SIG / " + expTimeString) 469 469 sqlLine.group("apFillF", "AP_NPIX / (3.14159265359 * POW(AP_MAG_RADIUS - 0.5, 2))") 470 470 sqlLine.group("apRadius", "AP_MAG_RADIUS * abs(PLTSCALE)") 471 sqlLine.group("kronFlux", "KRON_FLUX / " + ex tTimeString)472 sqlLine.group("kronFluxErr", "KRON_FLUX_ERR / " + ex tTimeString)471 sqlLine.group("kronFlux", "KRON_FLUX / " + expTimeString) 472 sqlLine.group("kronFluxErr", "KRON_FLUX_ERR / " + expTimeString) 473 473 sqlLine.group("kronRad", "MOMENTS_R1 * 2.5 * abs(PLTSCALE)") 474 sqlLine.group("sky", "SKY / PLTSCALE / PLTSCALE / " + ex tTimeString)475 sqlLine.group("skyErr", "SKY_SIGMA / PLTSCALE / PLTSCALE / " + ex tTimeString)474 sqlLine.group("sky", "SKY / PLTSCALE / PLTSCALE / " + expTimeString) 475 sqlLine.group("skyErr", "SKY_SIGMA / PLTSCALE / PLTSCALE / " + expTimeString) 476 476 sqlLine.group("infoFlag", "FLAGS") 477 477 sqlLine.group("infoFlag2", "FLAGS2") -
trunk/ippToPsps/jython/diffbatch.py
r38783 r38802 264 264 #results['SATDET'] = 0 265 265 266 ex tTimeString = str(self.expTime[num])266 expTimeString = str(self.expTime[num]) 267 267 268 268 # insert all detections into table : do NOT include the IGNORE unless we have to … … 296 296 sqlLine.group("raErr", "X_PSF_SIG * 0.257") 297 297 sqlLine.group("decErr", "Y_PSF_SIG * 0.257") 298 sqlLine.group("DpsfFlux", "PSF_INST_FLUX / " + ex tTimeString)299 sqlLine.group("DpsfFluxErr", "PSF_INST_FLUX_SIG / " + ex tTimeString)298 sqlLine.group("DpsfFlux", "PSF_INST_FLUX / " + expTimeString) 299 sqlLine.group("DpsfFluxErr", "PSF_INST_FLUX_SIG / " + expTimeString) 300 300 sqlLine.group("DpsfMajorFWHM", "PSF_FWHM_MAJ * abs(PLTSCALE)") 301 301 sqlLine.group("DpsfMinorFWHM", "PSF_FWHM_MIN * abs(PLTSCALE)") … … 316 316 sqlLine.group("DmomentM4S", "MOMENTS_M4S * PLTSCALE * PLTSCALE") 317 317 if (self.id >= 982483): 318 sqlLine.group("DapFlux", "AP_FLUX / " + ex tTimeString)319 sqlLine.group("DapFluxErr", "AP_FLUX_SIG / " + ex tTimeString)318 sqlLine.group("DapFlux", "AP_FLUX / " + expTimeString) 319 sqlLine.group("DapFluxErr", "AP_FLUX_SIG / " + expTimeString) 320 320 # sqlLine.group("DapFillF", "AP_NPIX / (3.14159265359 * POW(AP_MAG_RADIUS - 0.5, 2))") 321 321 # sqlLine.group("DapRadius", "AP_MAG_RADIUS") 322 sqlLine.group("DkronFlux", "KRON_FLUX / " + ex tTimeString)323 sqlLine.group("DkronFluxErr", "KRON_FLUX_ERR / " + ex tTimeString)322 sqlLine.group("DkronFlux", "KRON_FLUX / " + expTimeString) 323 sqlLine.group("DkronFluxErr", "KRON_FLUX_ERR / " + expTimeString) 324 324 sqlLine.group("DkronRad", "MOMENTS_R1 * 2.5 * abs(PLTSCALE)") 325 325 sqlLine.group("diffNPos", "DIFF_NPOS") … … 332 332 sqlLine.group("diffPSN", "DIFF_SN_P") 333 333 sqlLine.group("diffNSN", "DIFF_SN_M") 334 sqlLine.group("Dsky", "SKY / PLTSCALE / PLTSCALE / " + ex tTimeString)335 sqlLine.group("DskyErr", "SKY_SIGMA / PLTSCALE / PLTSCALE / " + ex tTimeString)334 sqlLine.group("Dsky", "SKY / PLTSCALE / PLTSCALE / " + expTimeString) 335 sqlLine.group("DskyErr", "SKY_SIGMA / PLTSCALE / PLTSCALE / " + expTimeString) 336 336 sqlLine.group("DinfoFlag", "FLAGS") 337 337 sqlLine.group("DinfoFlag2", "FLAGS2") -
trunk/ippToPsps/jython/forcedwarpbatch.py
r38800 r38802 332 332 333 333 BEFORE = self.scratchDb.getRowCount(ippTableName) 334 ex tTimeString = str(self.expTime[num])334 expTimeString = str(self.expTime[num]) 335 335 336 336 #missing are zp, telluricExt, airmass (not easy to get ?) … … 353 353 print "here4" 354 354 sqlLine.group("obsTime", str(self.obsTime[num])) 355 sqlLine.group("FpsfFlux", "PSF_INST_FLUX / " + ex tTimeString)356 sqlLine.group("FpsfFluxErr", "PSF_INST_FLUX_SIG / " + ex tTimeString)355 sqlLine.group("FpsfFlux", "PSF_INST_FLUX / " + expTimeString) 356 sqlLine.group("FpsfFluxErr", "PSF_INST_FLUX_SIG / " + expTimeString) 357 357 sqlLine.group("FpsfMajorFWHM", "PSF_FWHM_MAJ * abs(PLTSCALE)") 358 358 sqlLine.group("FpsfMinorFWHM", "PSF_FWHM_MIN * abs(PLTSCALE)") … … 372 372 sqlLine.group("FmomentM4C", "MOMENTS_M4C * PLTSCALE * PLTSCALE") 373 373 sqlLine.group("FmomentM4S", "MOMENTS_M4S * PLTSCALE * PLTSCALE") 374 sqlLine.group("FapFlux", "AP_FLUX / " + ex tTimeString)375 sqlLine.group("FapFluxErr", "AP_FLUX_SIG / " + ex tTimeString)374 sqlLine.group("FapFlux", "AP_FLUX / " + expTimeString) 375 sqlLine.group("FapFluxErr", "AP_FLUX_SIG / " + expTimeString) 376 376 sqlLine.group("FapFillF", "AP_NPIX / (3.14159265359 * POW(AP_MAG_RADIUS - 0.5, 2))") 377 377 sqlLine.group("FapRadius", "AP_MAG_RADIUS") 378 sqlLine.group("FkronFlux", "KRON_FLUX / " + ex tTimeString)379 sqlLine.group("FkronFluxErr", "KRON_FLUX_ERR / " + ex tTimeString )378 sqlLine.group("FkronFlux", "KRON_FLUX / " + expTimeString) 379 sqlLine.group("FkronFluxErr", "KRON_FLUX_ERR / " + expTimeString ) 380 380 sqlLine.group("FkronRad", "MOMENTS_R1 * 2.5 * abs(PLTSCALE)") 381 sqlLine.group("Fsky", "SKY / PLTSCALE / PLTSCALE / " + ex tTimeString)382 sqlLine.group("FskyErr", "SKY_SIGMA / PLTSCALE / PLTSCALE / " + ex tTimeString)381 sqlLine.group("Fsky", "SKY / PLTSCALE / PLTSCALE / " + expTimeString) 382 sqlLine.group("FskyErr", "SKY_SIGMA / PLTSCALE / PLTSCALE / " + expTimeString) 383 383 sqlLine.group("FinfoFlag", "FLAGS") 384 384 sqlLine.group("FinfoFlag2", "FLAGS2") -
trunk/ippToPsps/jython/stackbatch.py
r38801 r38802 361 361 raise 362 362 363 # insert values coming from the CMF, not carried by DVO364 header = self.headerSet[filter]365 exptime = self.getKeyFloat(header, "%.5f", "EXPTIME")366 367 363 sqlLine = sqlUtility("UPDATE StackObjectThin AS a , " + filter + "SkyChip_psf AS b SET") 368 364 … … 424 420 header = self.headerSet[filter] 425 421 exptime = self.getKeyFloat(header, "%.5f", "EXPTIME") 422 expTimeString = str(exptime) 426 423 427 424 # insert all the detections … … 447 444 sqlLine.group("a."+filter+"momentR1", "b.MOMENTS_R1 * b.PLTSCALE") 448 445 sqlLine.group("a."+filter+"momentRH", "b.MOMENTS_RH * b.PLTSCALE") 449 sqlLine.group("a."+filter+"PSFFlux", "b.PSF_INST_FLUX / " + str(exptime))450 sqlLine.group("a."+filter+"PSFFluxErr", "b.PSF_INST_FLUX_SIG / " + str(exptime))451 sqlLine.group("a."+filter+"ApFlux", "b.AP_FLUX / " + str(exptime))452 sqlLine.group("a."+filter+"ApFluxErr", "b.AP_FLUX_SIG / " + str(exptime))453 sqlLine.group("a."+filter+"KronFlux", "b.KRON_FLUX / " + str(exptime))454 sqlLine.group("a."+filter+"KronFluxErr", "b.KRON_FLUX_ERR / " + str(exptime))446 sqlLine.group("a."+filter+"PSFFlux", "b.PSF_INST_FLUX / " + expTimeString) 447 sqlLine.group("a."+filter+"PSFFluxErr", "b.PSF_INST_FLUX_SIG / " + expTimeString) 448 sqlLine.group("a."+filter+"ApFlux", "b.AP_FLUX / " + expTimeString) 449 sqlLine.group("a."+filter+"ApFluxErr", "b.AP_FLUX_SIG / " + expTimeString) 450 sqlLine.group("a."+filter+"KronFlux", "b.KRON_FLUX / " + expTimeString) 451 sqlLine.group("a."+filter+"KronFluxErr", "b.KRON_FLUX_ERR / " + expTimeString) 455 452 sqlLine.group("a."+filter+"ApFillFac", "b.AP_NPIX / (3.14159265359 * POW(b.AP_MAG_RADIUS - 0.5, 2))") 456 453 sqlLine.group("a."+filter+"ApRadius", "b.AP_MAG_RADIUS") … … 458 455 sqlLine.group("a."+filter+"PlateScale", "b.PLTSCALE") 459 456 sqlLine.group("a."+filter+"ExtNSigma", "b.EXT_NSIGMA") 460 sqlLine.group("a."+filter+"sky", "b.SKY / " + str(exptime))461 sqlLine.group("a."+filter+"skyErr", "b.SKY_SIGMA / " + str(exptime))457 sqlLine.group("a."+filter+"sky", "b.SKY / " + expTimeString) 458 sqlLine.group("a."+filter+"skyErr", "b.SKY_SIGMA / " + expTimeString) 462 459 463 460 sql = sqlLine.makeEquals("WHERE a." + filter + "ippDetectID = b.IPP_IDET") … … 531 528 continue 532 529 533 header = self.headerSet[filter]534 exptime = self.getKeyFloat(header, "%.5f", "EXPTIME")535 536 530 # insert all the detections 537 531 sqlLine = sqlUtility("UPDATE " + tablename + " AS a , " + filter + "SkyChip_xsrc AS b SET") 538 532 539 sqlLine.group("a."+filter+"haveData", "'1'")540 sqlLine.group("a."+filter+"S2", "b.G_S2")533 sqlLine.group("a."+filter+"haveData", "'1'") 534 sqlLine.group("a."+filter+"S2", "b.G_S2") 541 535 sqlLine.group("a."+filter+"logRT", "b.G_RT") 542 536 sqlLine.group("a."+filter+"logRA", "b.G_RA") 543 537 sqlLine.group("a."+filter+"bumpy", "b.G_BUMPY") 544 sqlLine.group("a."+filter+"halfLightRad", "b.HALF_LIGHT_RADIUS") 538 sqlLine.group("a."+filter+"halfLightRad", "b.HALF_LIGHT_RADIUS") 539 545 540 sql = sqlLine.makeEquals("WHERE a." + filter + "ippDetectID = b.IPP_IDET") 546 541 … … 781 776 # PETRO_MAG is instrumental + 2.5log(exptime) + FPA.ZP 782 777 # we want to apply the DVO zero point, so remove FPA.ZP: 778 ## XXX is this statement true: is FPA.ZP + 2.5log(exptime) applied to PETRO_MAG? 783 779 784 780 hdrZP = "%.5f" % float(header['FPA.ZP']) 785 #magtime = "%.5f" % (2.5*math.log10(float(header['EXPTIME'])))781 magtime = "%.5f" % (2.5*math.log10(float(header['EXPTIME']))) 786 782 787 783 self.logger.infoPair("Adding", "petrosians for extended sources")
Note:
See TracChangeset
for help on using the changeset viewer.
