Changeset 34615
- Timestamp:
- Oct 26, 2012, 5:15:52 PM (14 years ago)
- Location:
- trunk/ippToPsps
- Files:
-
- 6 edited
-
config/OB/tables.vot (modified) (8 diffs)
-
config/P2/tables.vot (modified) (5 diffs)
-
config/ST/tables.vot (modified) (10 diffs)
-
jython/detectionbatch.py (modified) (14 diffs)
-
jython/objectbatch.py (modified) (3 diffs)
-
jython/stackbatch.py (modified) (18 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippToPsps/config/OB/tables.vot
r34275 r34615 45 45 <DESCRIPTION> Declination error</DESCRIPTION> 46 46 </FIELD> 47 <FIELD name="cx" arraysize="1" datatype="double" unit="unknown">48 <DESCRIPTION>Cartesian x</DESCRIPTION>49 </FIELD>50 <FIELD name="cy" arraysize="1" datatype="double" unit="unknown">51 <DESCRIPTION>Cartesian y</DESCRIPTION>52 </FIELD>53 <FIELD name="cz" arraysize="1" datatype="double" unit="unknown">54 <DESCRIPTION>Cartesian z</DESCRIPTION>55 </FIELD>56 47 <FIELD name="lambda" arraysize="1" datatype="double" unit="degree"> 57 48 <DESCRIPTION> ecliptic longitude</DESCRIPTION> … … 66 57 <DESCRIPTION> galactic latitude</DESCRIPTION> 67 58 </FIELD> 68 <FIELD name="lsg" arraysize="1" datatype="double" unit="unknown">69 <DESCRIPTION>supergalactic longitude</DESCRIPTION>70 </FIELD>71 <FIELD name="bsg" arraysize="1" datatype="double" unit="unknown">72 <DESCRIPTION>supergalactic latitude</DESCRIPTION>73 </FIELD>74 59 <FIELD name="nDetections" arraysize="1" datatype="short" unit="unknown"> 75 60 <DESCRIPTION>total number of detection measurements in all filters</DESCRIPTION> … … 201 186 <DESCRIPTION> 80 pct w magnitude</DESCRIPTION> 202 187 </FIELD> 203 <FIELD name="gStackPSF Flux" arraysize="1" datatype="float" unit="mag">188 <FIELD name="gStackPSFMag" arraysize="1" datatype="float" unit="mag"> 204 189 <DESCRIPTION> g stack psf flux </DESCRIPTION> 205 190 </FIELD> 206 <FIELD name="gStackPSF FluxErr" arraysize="1" datatype="double" unit="mag">191 <FIELD name="gStackPSFMagErr" arraysize="1" datatype="double" unit="mag"> 207 192 <DESCRIPTION> estimated error in g stack psf flux</DESCRIPTION> 208 193 </FIELD> … … 216 201 <DESCRIPTION>flags for g stack</DESCRIPTION> 217 202 </FIELD> 218 <FIELD name="rStackPSFFlux" arraysize="1" datatype="float" unit="mag"> 203 <FIELD name="gStackDetectID" arraysize="1" datatype="long" unit="unknown"> 204 <DESCRIPTION>ODM detection identifier</DESCRIPTION> 205 </FIELD> 206 <FIELD name="rStackPSFMag" arraysize="1" datatype="float" unit="mag"> 219 207 <DESCRIPTION> r stack psf flux </DESCRIPTION> 220 208 </FIELD> 221 <FIELD name="rStackPSF FluxErr" arraysize="1" datatype="double" unit="mag">209 <FIELD name="rStackPSFMagErr" arraysize="1" datatype="double" unit="mag"> 222 210 <DESCRIPTION> estimated error in r stack psf flux</DESCRIPTION> 223 211 </FIELD> … … 231 219 <DESCRIPTION>flags for r stack</DESCRIPTION> 232 220 </FIELD> 233 <FIELD name="iStackPSFFlux" arraysize="1" datatype="float" unit="mag"> 221 <FIELD name="rStackDetectID" arraysize="1" datatype="long" unit="unknown"> 222 <DESCRIPTION>ODM detection identifier</DESCRIPTION> 223 </FIELD> 224 <FIELD name="iStackPSFMag" arraysize="1" datatype="float" unit="mag"> 234 225 <DESCRIPTION> i stack psf flux </DESCRIPTION> 235 226 </FIELD> 236 <FIELD name="iStackPSF FluxErr" arraysize="1" datatype="float" unit="mag">227 <FIELD name="iStackPSFMagErr" arraysize="1" datatype="float" unit="mag"> 237 228 <DESCRIPTION> estimated error in i stack psf flux</DESCRIPTION> 238 229 </FIELD> … … 246 237 <DESCRIPTION>flags for i stack</DESCRIPTION> 247 238 </FIELD> 248 <FIELD name="zStackPSFFlux" arraysize="1" datatype="float" unit="mag"> 239 <FIELD name="iStackDetectID" arraysize="1" datatype="long" unit="unknown"> 240 <DESCRIPTION>ODM detection identifier</DESCRIPTION> 241 </FIELD> 242 <FIELD name="zStackPSFMag" arraysize="1" datatype="float" unit="mag"> 249 243 <DESCRIPTION> z stack psf flux </DESCRIPTION> 250 244 </FIELD> 251 <FIELD name="zStackPSF FluxErr" arraysize="1" datatype="float" unit="mag">245 <FIELD name="zStackPSFMagErr" arraysize="1" datatype="float" unit="mag"> 252 246 <DESCRIPTION> estimated error in z stack psf flux</DESCRIPTION> 253 247 </FIELD> … … 261 255 <DESCRIPTION>flags for z stack</DESCRIPTION> 262 256 </FIELD> 263 <FIELD name="yStackPSFFlux" arraysize="1" datatype="float" unit="unknown"> 257 <FIELD name="zStackDetectID" arraysize="1" datatype="long" unit="unknown"> 258 <DESCRIPTION>ODM detection identifier</DESCRIPTION> 259 </FIELD> 260 <FIELD name="yStackPSFMag" arraysize="1" datatype="float" unit="unknown"> 264 261 <DESCRIPTION>y band stack PSF flux</DESCRIPTION> 265 262 </FIELD> 266 <FIELD name="yStackPSF FluxErr" arraysize="1" datatype="float" unit="unknown">263 <FIELD name="yStackPSFMagErr" arraysize="1" datatype="float" unit="unknown"> 267 264 <DESCRIPTION>y band stack psf flux err</DESCRIPTION> 268 265 </FIELD> … … 276 273 <DESCRIPTION>flags for y stack</DESCRIPTION> 277 274 </FIELD> 278 <FIELD name="extra" arraysize="1" datatype="float" unit="unknown"> 279 <DESCRIPTION>currently unused</DESCRIPTION> 275 <FIELD name="yStackDetectID" arraysize="1" datatype="long" unit="unknown"> 276 <DESCRIPTION>ODM detection identifier</DESCRIPTION> 277 </FIELD> 278 <FIELD name="qfPerfect" arraysize="1" datatype="float" unit="unknown"> 279 <DESCRIPTION>qf perfect</DESCRIPTION> 280 280 </FIELD> 281 281 <FIELD name="sgSep" arraysize="1" datatype="float" unit="unknown"> -
trunk/ippToPsps/config/P2/tables.vot
r32031 r34615 404 404 <DESCRIPTION> estimated error in y</DESCRIPTION> 405 405 </FIELD> 406 <FIELD name=" instFlux" arraysize="1" datatype="float" unit="adu/s">406 <FIELD name="psfFlux" arraysize="1" datatype="float" unit="adu/s"> 407 407 <DESCRIPTION> PSF instrumental flux</DESCRIPTION> 408 408 </FIELD> 409 <FIELD name=" instFluxErr" arraysize="1" datatype="float" unit="adu/s">409 <FIELD name="psfFluxErr" arraysize="1" datatype="float" unit="adu/s"> 410 410 <DESCRIPTION> estimated error in flux</DESCRIPTION> 411 411 </FIELD> … … 425 425 <DESCRIPTION>PSF likelihood</DESCRIPTION> 426 426 </FIELD> 427 <FIELD name="psf Cf" arraysize="1" datatype="float" unit="unknown">428 <DESCRIPTION>PSF coverage factor</DESCRIPTION>427 <FIELD name="psfQf" arraysize="1" datatype="float" unit="unknown"> 428 <DESCRIPTION>PSF qf</DESCRIPTION> 429 429 </FIELD> 430 430 <FIELD name="momentXX" arraysize="1" datatype="float" unit="arcsec"> … … 437 437 <DESCRIPTION> momeny YY</DESCRIPTION> 438 438 </FIELD> 439 <FIELD name="ap Mag" arraysize="1" datatype="float" unit="magnitude">440 <DESCRIPTION> Aperture magnitude</DESCRIPTION>441 </FIELD> 442 <FIELD name="ap MagErr" arraysize="1" datatype="float" unit="magnitude">443 <DESCRIPTION> Aperture magnitudeerror</DESCRIPTION>444 </FIELD> 445 <FIELD name="kronFlux" arraysize="1" datatype="float" unit=" magnitude">439 <FIELD name="apFlux" arraysize="1" datatype="float" unit="adu/s"> 440 <DESCRIPTION> Aperture Flux</DESCRIPTION> 441 </FIELD> 442 <FIELD name="apFluxErr" arraysize="1" datatype="float" unit="adu/s"> 443 <DESCRIPTION> Aperture Flux error</DESCRIPTION> 444 </FIELD> 445 <FIELD name="kronFlux" arraysize="1" datatype="float" unit="adu/s"> 446 446 <DESCRIPTION> Kron flux</DESCRIPTION> 447 447 </FIELD> 448 <FIELD name="kronFluxErr" arraysize="1" datatype="float" unit=" magnitude">448 <FIELD name="kronFluxErr" arraysize="1" datatype="float" unit="adu/s"> 449 449 <DESCRIPTION> Kron flux error</DESCRIPTION> 450 450 </FIELD> … … 455 455 <DESCRIPTION> Kron radius error</DESCRIPTION> 456 456 </FIELD> 457 <FIELD name=" crLikelihood" arraysize="1" datatype="float" unit="unknown">458 <DESCRIPTION> Likelihood the source is a cosmic ray</DESCRIPTION>459 </FIELD> 460 <FIELD name=" extendedLikelihood" arraysize="1" datatype="float" unit="unknown">461 <DESCRIPTION> Likelihood the source is extended</DESCRIPTION>457 <FIELD name="psfQfPerfect" arraysize="1" datatype="float" unit="unknown"> 458 <DESCRIPTION>psf QF Perfect</DESCRIPTION> 459 </FIELD> 460 <FIELD name="psfChiSq" arraysize="1" datatype="float" unit="unknown"> 461 <DESCRIPTION>PSF Chi sq</DESCRIPTION> 462 462 </FIELD> 463 463 <FIELD name="infoFlag" arraysize="1" datatype="long" unit="unknown"> 464 464 <DESCRIPTION>flag indicating provenance information</DESCRIPTION> 465 465 </FIELD> 466 <FIELD name="sky" arraysize="1" datatype="float" unit="adu ">466 <FIELD name="sky" arraysize="1" datatype="float" unit="adu/s/pixel"> 467 467 <DESCRIPTION> PSF sky level at source</DESCRIPTION> 468 468 </FIELD> 469 <FIELD name="skyErr" arraysize="1" datatype="float" unit="adu ">469 <FIELD name="skyErr" arraysize="1" datatype="float" unit="adu/s/pixel"> 470 470 <DESCRIPTION> estimated error in sky</DESCRIPTION> 471 471 </FIELD> … … 539 539 </DATA> 540 540 </TABLE> 541 <TABLE name="DetectionCalib"> 542 <DESCRIPTION>VOTable description of PSPS table DetectionCalib</DESCRIPTION> 543 <PARAM arraysize="1" datatype="char" ucd="meta.bib.author" name="Author" value="PSPS"></PARAM> 544 <FIELD name="objID" arraysize="1" datatype="long" unit="unknown"> 545 <DESCRIPTION>ODM object identifier index</DESCRIPTION> 546 </FIELD> 547 <FIELD name="detectID" arraysize="1" datatype="long" unit="unknown"> 548 <DESCRIPTION>ODM detection identifier index</DESCRIPTION> 549 </FIELD> 550 <FIELD name="ippObjID" arraysize="1" datatype="long" unit="unknown"> 551 <DESCRIPTION>IPP object number</DESCRIPTION> 552 </FIELD> 553 <FIELD name="ippDetectID" arraysize="1" datatype="long" unit="unknown"> 554 <DESCRIPTION>detection ID generated by IPP</DESCRIPTION> 555 </FIELD> 556 <FIELD name="filterID" arraysize="1" datatype="unsignedByte" unit="unknown"> 557 <DESCRIPTION>filter identifier</DESCRIPTION> 558 </FIELD> 559 <FIELD name="surveyID" arraysize="1" datatype="unsignedByte" unit="unknown"> 560 <DESCRIPTION>survey identifier</DESCRIPTION> 561 </FIELD> 562 <FIELD name="ra" arraysize="1" datatype="double" unit="degrees"> 563 <DESCRIPTION>right ascension</DESCRIPTION> 564 </FIELD> 565 <FIELD name="dec" arraysize="1" datatype="double" unit="degrees"> 566 <DESCRIPTION>declination</DESCRIPTION> 567 </FIELD> 568 <FIELD name="raErr" arraysize="1" datatype="float" unit="arcsec"> 569 <DESCRIPTION>right ascension error</DESCRIPTION> 570 </FIELD> 571 <FIELD name="decErr" arraysize="1" datatype="float" unit="arcsec"> 572 <DESCRIPTION>declination error</DESCRIPTION> 573 </FIELD> 574 <FIELD name="zp" arraysize="1" datatype="float" unit="magnitude"> 575 <DESCRIPTION>zero point offset</DESCRIPTION> 576 </FIELD> 577 <FIELD name="zpErr" arraysize="1" datatype="float" unit="magnitude"> 578 <DESCRIPTION>zerop point offset error</DESCRIPTION> 579 </FIELD> 580 <FIELD name="expTime" arraysize="1" datatype="float" unit="s"> 581 <DESCRIPTION>exposure time</DESCRIPTION> 582 </FIELD> 583 <FIELD name="airMass" arraysize="1" datatype="float" unit="unknown"> 584 <DESCRIPTION>airmass</DESCRIPTION> 585 </FIELD> 586 <FIELD name="dataRelease" arraysize="1" datatype="unsignedByte" unit="unknown"> 587 <DESCRIPTION>Data release to propagate to the object</DESCRIPTION> 588 </FIELD> 589 <DATA> 590 <TABLEDATA></TABLEDATA> 591 </DATA> 592 </TABLE> 541 593 </RESOURCE> 542 594 </VOTABLE> -
trunk/ippToPsps/config/ST/tables.vot
r32031 r34615 148 148 <DESCRIPTION>projection cell identifier</DESCRIPTION> 149 149 </FIELD> 150 <FIELD name="obsTime" arraysize="1" datatype="double">151 <DESCRIPTION> Time of mid observation (unit = day)</DESCRIPTION>152 </FIELD>153 150 <FIELD name="xPos" arraysize="1" datatype="float"> 154 151 <DESCRIPTION>measured x on CCD from PSF fit</DESCRIPTION> … … 163 160 <DESCRIPTION>estimated error in y</DESCRIPTION> 164 161 </FIELD> 165 <FIELD name=" instFlux" arraysize="1" datatype="float">162 <FIELD name="psfFlux" arraysize="1" datatype="float"> 166 163 <DESCRIPTION>PSF instrumental flux</DESCRIPTION> 167 164 </FIELD> 168 <FIELD name=" instFluxErr" arraysize="1" datatype="float">165 <FIELD name="psfFluxErr" arraysize="1" datatype="float"> 169 166 <DESCRIPTION>estimated error in instrumental flux</DESCRIPTION> 170 167 </FIELD> … … 193 190 <DESCRIPTION>PSF likelihood</DESCRIPTION> 194 191 </FIELD> 195 <FIELD name="psfCf" arraysize="1" datatype="float"> 196 <DESCRIPTION>PSF coverage factor</DESCRIPTION> 197 </FIELD> 192 <FIELD name="psfQf" arraysize="1" datatype="float"> 193 <DESCRIPTION>PSF Q factor</DESCRIPTION> 194 </FIELD> 195 <FIELD name="psfQfPerfect" arraysize="1" datatype="float"> 196 <DESCRIPTION>PSF QF Perfect</DESCRIPTION> 197 </FIELD> 198 198 <FIELD name="momentXX" arraysize="1" datatype="float"> 199 199 <DESCRIPTION> moment XX (unit = arcsec)</DESCRIPTION> … … 223 223 <DESCRIPTION> moment RH (unit = arcsec)</DESCRIPTION> 224 224 </FIELD> 225 <FIELD name="ap Mag" arraysize="1" datatype="float">226 <DESCRIPTION> Aperture magnitude (unit = magnitude)</DESCRIPTION>227 </FIELD> 228 <FIELD name="ap MagErr" arraysize="1" datatype="float">229 <DESCRIPTION> Aperture magnitude error (unit = magnitude)</DESCRIPTION>225 <FIELD name="apFlux" arraysize="1" datatype="float"> 226 <DESCRIPTION> Aperture flux (unit = adu/s)</DESCRIPTION> 227 </FIELD> 228 <FIELD name="apFluxErr" arraysize="1" datatype="float"> 229 <DESCRIPTION> Aperture flux error (unit = adu/s)</DESCRIPTION> 230 230 </FIELD> 231 231 <FIELD name="kronFlux" arraysize="1" datatype="float"> 232 <DESCRIPTION> Kron flux (unit = magnitude)</DESCRIPTION>232 <DESCRIPTION> Kron flux (unit = adu/s)</DESCRIPTION> 233 233 </FIELD> 234 234 <FIELD name="kronFluxErr" arraysize="1" datatype="float"> 235 <DESCRIPTION> Kron flux error (unit = magnitude)</DESCRIPTION>235 <DESCRIPTION> Kron flux error (unit = adu/s)</DESCRIPTION> 236 236 </FIELD> 237 237 <FIELD name="kronRad" arraysize="1" datatype="float"> 238 238 <DESCRIPTION> Kron radius (unit = arcsecs)</DESCRIPTION> 239 239 </FIELD> 240 <FIELD name=" kronRadErr" arraysize="1" datatype="float">241 <DESCRIPTION> Kron radius error (unit = arcsecs)</DESCRIPTION>240 <FIELD name="psfChiSq" arraysize="1" datatype="float"> 241 <DESCRIPTION> psf chi sq</DESCRIPTION> 242 242 </FIELD> 243 243 <FIELD name="infoFlag" arraysize="1" datatype="long"> … … 293 293 <DESCRIPTION>stack identifier</DESCRIPTION> 294 294 </FIELD> 295 <FIELD name="isophot Mag" arraysize="1" datatype="float">296 <DESCRIPTION>isophotal magnitude</DESCRIPTION>297 </FIELD> 298 <FIELD name="isophot MagErr" arraysize="1" datatype="float">299 <DESCRIPTION>estimated error in isophotal magnitude</DESCRIPTION>295 <FIELD name="isophotFlux" arraysize="1" datatype="float"> 296 <DESCRIPTION>isophotal flux</DESCRIPTION> 297 </FIELD> 298 <FIELD name="isophotFluxErr" arraysize="1" datatype="float"> 299 <DESCRIPTION>estimated error in isophotal flux</DESCRIPTION> 300 300 </FIELD> 301 301 <FIELD name="isophotMajAxis" arraysize="1" datatype="float"> … … 332 332 <DESCRIPTION>estimated error inPetrosian radius</DESCRIPTION> 333 333 </FIELD> 334 <FIELD name="pet Mag" arraysize="1" datatype="float">335 <DESCRIPTION>Petrosian magntiude</DESCRIPTION>336 </FIELD> 337 <FIELD name="pet MagErr" arraysize="1" datatype="float">338 <DESCRIPTION>estimated error in pet Mag</DESCRIPTION>334 <FIELD name="petFlux" arraysize="1" datatype="float"> 335 <DESCRIPTION>Petrosian Flux</DESCRIPTION> 336 </FIELD> 337 <FIELD name="petFluxErr" arraysize="1" datatype="float"> 338 <DESCRIPTION>estimated error in petFlux</DESCRIPTION> 339 339 </FIELD> 340 340 <FIELD name="petR50" arraysize="1" datatype="float"> … … 765 765 <DESCRIPTION>estimated error in deVaucouleurs radius</DESCRIPTION> 766 766 </FIELD> 767 <FIELD name="deV Mag" arraysize="1" datatype="float">767 <FIELD name="deVFlux" arraysize="1" datatype="float"> 768 768 <DESCRIPTION>deVaucouleurs magntiude</DESCRIPTION> 769 769 </FIELD> 770 <FIELD name="deV MagErr" arraysize="1" datatype="float">770 <FIELD name="deVFluxErr" arraysize="1" datatype="float"> 771 771 <DESCRIPTION>estimated error in deV_mag</DESCRIPTION> 772 772 </FIELD> … … 891 891 <DESCRIPTION>estimated error in Exponential fit radius</DESCRIPTION> 892 892 </FIELD> 893 <FIELD name="exp Mag" arraysize="1" datatype="float">894 <DESCRIPTION>Exponential fit magntiude</DESCRIPTION>895 </FIELD> 896 <FIELD name="exp MagErr" arraysize="1" datatype="float">897 <DESCRIPTION>estimated error in exp Mag</DESCRIPTION>893 <FIELD name="expFlux" arraysize="1" datatype="float"> 894 <DESCRIPTION>Exponential fit flux</DESCRIPTION> 895 </FIELD> 896 <FIELD name="expFluxErr" arraysize="1" datatype="float"> 897 <DESCRIPTION>estimated error in expFlux</DESCRIPTION> 898 898 </FIELD> 899 899 <FIELD name="expAb" arraysize="1" datatype="float"> … … 1017 1017 <DESCRIPTION>estimated error in Sersic radius</DESCRIPTION> 1018 1018 </FIELD> 1019 <FIELD name="ser Mag" arraysize="1" datatype="float">1020 <DESCRIPTION>Sersic magntiude</DESCRIPTION>1021 </FIELD> 1022 <FIELD name="ser MagErr" arraysize="1" datatype="float">1023 <DESCRIPTION>estimated error in ser Mag</DESCRIPTION>1019 <FIELD name="serFlux" arraysize="1" datatype="float"> 1020 <DESCRIPTION>Sersic Flux</DESCRIPTION> 1021 </FIELD> 1022 <FIELD name="serFluxErr" arraysize="1" datatype="float"> 1023 <DESCRIPTION>estimated error in serFlux</DESCRIPTION> 1024 1024 </FIELD> 1025 1025 <FIELD name="serAb" arraysize="1" datatype="float"> … … 1240 1240 </DATA> 1241 1241 </TABLE> 1242 <TABLE name="StackDetectionCalib"> 1243 <DESCRIPTION>VOTable description of PSPS table StackDetectionCalib</DESCRIPTION> 1244 <PARAM arraysize="1" datatype="char" ucd="meta.bib.author" name="Author" value="PSPS"></PARAM> 1245 <FIELD name="objID" arraysize="1" datatype="long" unit="unknown"> 1246 <DESCRIPTION>ODM object identifier index</DESCRIPTION> 1247 </FIELD> 1248 <FIELD name="stackDetectID" arraysize="1" datatype="long" unit="unknown"> 1249 <DESCRIPTION>ODM detection identifier index</DESCRIPTION> 1250 </FIELD> 1251 <FIELD name="ippObjID" arraysize="1" datatype="long" unit="unknown"> 1252 <DESCRIPTION>IPP object number</DESCRIPTION> 1253 </FIELD> 1254 <FIELD name="ippDetectID" arraysize="1" datatype="long" unit="unknown"> 1255 <DESCRIPTION>detection ID generated by IPP</DESCRIPTION> 1256 </FIELD> 1257 <FIELD name="filterID" arraysize="1" datatype="unsignedByte" unit="unknown"> 1258 <DESCRIPTION>filter identifier</DESCRIPTION> 1259 </FIELD> 1260 <FIELD name="surveyID" arraysize="1" datatype="unsignedByte" unit="unknown"> 1261 <DESCRIPTION>survey identifier</DESCRIPTION> 1262 </FIELD> 1263 <FIELD name="ra" arraysize="1" datatype="double" unit="degrees"> 1264 <DESCRIPTION>right ascension</DESCRIPTION> 1265 </FIELD> 1266 <FIELD name="dec" arraysize="1" datatype="double" unit="degrees"> 1267 <DESCRIPTION>declination</DESCRIPTION> 1268 </FIELD> 1269 <FIELD name="raErr" arraysize="1" datatype="float" unit="arcsec"> 1270 <DESCRIPTION>right ascension error</DESCRIPTION> 1271 </FIELD> 1272 <FIELD name="decErr" arraysize="1" datatype="float" unit="arcsec"> 1273 <DESCRIPTION>declination error</DESCRIPTION> 1274 </FIELD> 1275 <FIELD name="zp" arraysize="1" datatype="float" unit="magnitude"> 1276 <DESCRIPTION>zero point offset</DESCRIPTION> 1277 </FIELD> 1278 <FIELD name="zpErr" arraysize="1" datatype="float" unit="magnitude"> 1279 <DESCRIPTION>zerop point offset error</DESCRIPTION> 1280 </FIELD> 1281 <FIELD name="expTime" arraysize="1" datatype="float" unit="s"> 1282 <DESCRIPTION>exposure time</DESCRIPTION> 1283 </FIELD> 1284 <FIELD name="airMass" arraysize="1" datatype="float" unit="unknown"> 1285 <DESCRIPTION>airmass</DESCRIPTION> 1286 </FIELD> 1287 <FIELD name="dataRelease" arraysize="1" datatype="unsignedByte" unit="unknown"> 1288 <DESCRIPTION>Data release to propagate to the object</DESCRIPTION> 1289 </FIELD> 1290 <DATA> 1291 <TABLEDATA></TABLEDATA> 1292 </DATA> 1293 </TABLE> 1242 1294 </RESOURCE> 1243 1295 </VOTABLE> -
trunk/ippToPsps/jython/detectionbatch.py
r33830 r34615 388 388 # delete all detections with PSF_INST_MAG < 17.5, as decreed by Gene 389 389 BEFORE = self.scratchDb.getRowCount(ippTableName) 390 sql = "DELETE FROM " + ippTableName + " WHERE PSF_INST_MAG < -17.5" 391 self.scratchDb.execute(sql) 390 #don't do this anymore 391 #sql = "DELETE FROM " + ippTableName + " WHERE PSF_INST_MAG < -17.5" 392 #self.scratchDb.execute(sql) 392 393 #DELETED = BEFORE - self.scratchDb.getRowCount(ippTableName) 393 results['SATDET'] = BEFORE - self.scratchDb.getRowCount(ippTableName) 394 #results['SATDET'] = BEFORE - self.scratchDb.getRowCount(ippTableName) 395 #we aren't deletng these now 396 results['SATDET'] = 0 394 397 #self.logger.infoPair("Saturated detections", "%d deleted" % DELETED) 395 398 … … 404 407 ,xPosErr \ 405 408 ,yPosErr \ 406 , instFlux \407 , instFluxErr \409 ,psfFlux \ 410 ,psfFluxErr \ 408 411 ,peakADU \ 409 412 ,psfWidMajor \ … … 411 414 ,psfTheta \ 412 415 ,psfLikelihood \ 413 ,psf Cf \416 ,psfQf \ 414 417 ,momentXX \ 415 418 ,momentXY \ 416 419 ,momentYY \ 417 ,ap Mag\420 ,apFlux \ 418 421 ,kronFlux \ 419 422 ,kronFluxErr \ 423 ,psfQfPerfect \ 424 ,psfChiSq \ 420 425 ,infoFlag \ 421 426 ,sky \ … … 450 455 ,KRON_FLUX / " + self.header['EXPTIME'] + " \ 451 456 ,KRON_FLUX_ERR / " + self.header['EXPTIME'] + " \ 452 ,FLAGS\ 457 ,PSF_QF_PERFECT \ 458 ,PSF_CHISQ \ 459 ,FLAGS2 << 32 | FLAGS\ 453 460 ,SKY \ 454 461 ,SKY_SIGMA \ … … 459 466 , " + str(self.config.dataRelease) + "\ 460 467 FROM " + ippTableName 468 # SA10 todo 469 # need to correct formulas for psfFlux/psfFluxErr 470 # need to correct apMag apFlux 471 # need to correct unites kronFlux/Err to correct units 472 # need to change psfqfperfect -> psfqfperfect PSF_QF_PERFECT? 473 # need to change psfchisq PSF_CHI_SQ? 474 # need to change unites of sky/skyErr to adu/s/pixel 475 # need formula for sgsep/psflikelhood 476 # self.logger.info(sql) 461 477 self.scratchDb.execute(sql) 462 478 463 479 # now delete bad flux and bad chip positions 464 sql="DELETE FROM " + pspsTableName + " WHERE instFlux < 0.0000001" # TODO clearly a hack = 0 not allowed for instFlux 465 self.scratchDb.execute(sql) 466 480 #sql="DELETE FROM " + pspsTableName + " WHERE instFlux < 0.0000001" # TODO clearly a hack = 0 not allowed for instFlux 481 #self.scratchDb.execute(sql) 482 483 # add a instFlux = 0.0 -> 0.000001 484 sql="UPDATE " + pspsTableName + " SET psfFlux = 0.0000001 WHERE psfFlux = 0" 485 self.scratchDb.execute(sql) 486 467 487 # update cosmic ray and extended likelihoods 468 sql="UPDATE " + pspsTableName + " SET extendedLikelihood = 0, crLikelihood = 1.0 - psfLikelihood WHERE sgSep <= 0"469 self.scratchDb.execute(sql)470 sql="UPDATE " + pspsTableName + " SET crLikelihood = 0, extendedLikelihood = 1.0 - psfLikelihood WHERE sgSep > 0"471 self.scratchDb.execute(sql)488 #sql="UPDATE " + pspsTableName + " SET extendedLikelihood = 0, crLikelihood = 1.0 - psfLikelihood WHERE sgSep <= 0" 489 #self.scratchDb.execute(sql) 490 #sql="UPDATE " + pspsTableName + " SET crLikelihood = 0, extendedLikelihood = 1.0 - psfLikelihood WHERE sgSep > 0" 491 #self.scratchDb.execute(sql) 472 492 473 493 # remove detections will NULL inst flux or NULL peak ADU 474 results['NULLINSTFLUX'] = self.scratchDb.reportAndDeleteRowsWithNULLS(pspsTableName, "instFlux") 475 results['NULLPEAKADU'] = self.scratchDb.reportAndDeleteRowsWithNULLS(pspsTableName, "peakADU") 494 # leave them in at this stage- don't delete them. 495 #results['NULLINSTFLUX'] = self.scratchDb.reportAndDeleteRowsWithNULLS(pspsTableName, "instFlux") 496 #results['NULLPEAKADU'] = self.scratchDb.reportAndDeleteRowsWithNULLS(pspsTableName, "peakADU") 497 #we don't delete thos anymore 498 results['NULLINSTFLUX'] = 0; 499 results['NULLPEAKADU'] = 0; 476 500 477 501 ''' … … 504 528 505 529 ''' 530 Populates the DetectionCalib table for this OTA 531 ''' 532 def populateDetectionCalibTable(self, ota): 533 534 tableName = "DetectionCalib_" + ota 535 # drop then re-create table 536 self.scratchDb.dropTable(tableName) 537 sql = "CREATE TABLE " + tableName + " LIKE DetectionCalib" 538 try: self.scratchDb.execute(sql) 539 except: pass 540 541 # insert all detections into table 542 sql = "INSERT INTO " + tableName + " ( \ 543 objID \ 544 ,detectID \ 545 ,ippObjID \ 546 ,ippDetectID \ 547 ,filterID \ 548 ,surveyID \ 549 ) \ 550 SELECT \ 551 objID \ 552 ,detectID \ 553 ,ippObjID \ 554 ,ippDetectID \ 555 ,filterID \ 556 ,surveyID \ 557 FROM Detection_" + ota 558 self.scratchDb.execute(sql) 559 #SA10 TODO 560 # insert ra/dec/calstuff into Table 561 # heather discovered dec is called 'dec_' 562 sql = "UPDATE " + tableName + " SET ra = -999, dec_ =-999, raErr = -999, decErr = -999, zp = -999, zpErr = -999, expTime = -999, airMass = -999" 563 self.scratchDb.execute(sql) 564 self.scratchDb.updateAllRows(tableName, "dataRelease", str(self.config.dataRelease)) 565 566 ''' 506 567 Populates the Detection table for this OTA 507 568 ''' … … 608 669 extension = "XY%d%d_psf" % (x, y) 609 670 self.scratchDb.createIndex(extension, "IPP_IDET") 610 671 self.logger.infoPair("created indexes on", "IPP tables") 611 672 ''' 612 673 Updates provided table with DVO IDs from DVO table … … 620 681 a.detectID = b.detectID, \ 621 682 a.objID = b.objID, \ 622 a.infoFlag = b.flags << 32| a.infoFlag \683 a.infoFlag = b.flags << 46 | a.infoFlag \ 623 684 WHERE a.ippDetectID = b.ippDetectID \ 624 685 AND b.imageID = " + str(imageID) … … 658 719 self.logger.info("| %5s | ------------- Bad astrometric solution : rejecting ------------ |" % ota) 659 720 continue 660 721 # self.logger.info("populate stuff "); 661 722 # populate remainder of tables 662 723 self.populateDetectionTable(ota, results) 663 724 # self.logger.info("successful populate "); 664 725 # now add DVO IDs 665 726 self.updateDvoIDs("Detection_" + ota, self.imageIDs[ota]) 727 # self.logger.info("updated dvoids") 666 728 results['NULLOBJID'] = self.scratchDb.reportAndDeleteRowsWithNULLS("Detection_" + ota, "objID") 667 self.updateImageID("Detection_" + ota, x, y) 729 #self.logger.info("deleted nulls") 730 self.updateImageID("Detection_" + ota, x, y) 731 # self.logger.info("updateImageId") 668 732 rowCount = self.scratchDb.getRowCount("Detection_" + ota) 733 # self.logger.info("got row count") 669 734 self.logger.info("| %5s | %13d | %13d | %13d | %13d | %13d | %13d |", 670 735 ota, … … 675 740 results['NULLOBJID'], 676 741 rowCount) 677 678 742 totalOriginal = totalOriginal + results['ORIGINALTOTAL'] 679 743 totalSatDet = totalSatDet + results['SATDET'] … … 682 746 totalNullObjID = totalNullObjID + results['NULLOBJID'] 683 747 totalDetections = totalDetections + rowCount 684 748 #self.logger.info("updated totals") 685 749 # check we have something in this Detection table TODO add this to table above 686 750 if rowCount < 1: … … 692 756 SET nDetect = %d, photoCalID = %d" % (self.scratchDb.getRowCount("Detection_" + ota), self.scratchDb.getPhotoCalID(self.imageIDs[ota])) 693 757 self.scratchDb.execute(sql) 694 758 759 #self.logger.info("updated imagedata") 695 760 self.populateSkinnyObjectTable(ota) 761 #self.logger.info("updated skinnyobject") 696 762 self.populateObjectCalColorTable(ota) 697 763 #self.logger.info("updated objectcalcolor") 764 self.populateDetectionCalibTable(ota) 765 #self.logger.info("updated detectioncalibtable") 766 698 767 # add these to list of tables to export later 699 768 self.tablesToExport.append("ImageMeta_" + ota) … … 701 770 self.tablesToExport.append("SkinnyObject_" + ota) 702 771 self.tablesToExport.append("ObjectCalColor_" + ota) 772 self.tablesToExport.append("DetectionCalib_" + ota) 703 773 tables.append("Detection_" + ota) 704 774 705 775 otaCount = otaCount + 1 706 776 #self.logger.info("end of ota") 707 777 # print totals 708 778 self.logger.info("+-------+---------------+---------------+---------------+---------------+---------------+---------------+") … … 757 827 else : regex = ".*.psf" 758 828 759 columns = "IPP_IDET X_PSF Y_PSF X_PSF_SIG Y_PSF_SIG PSF_INST_MAG PSF_INST_MAG_SIG PEAK_FLUX_AS_MAG PSF_MAJOR PSF_MINOR PSF_THETA EXT_NSIGMA PSF_QF MOMENTS_XX MOMENTS_XY MOMENTS_YY AP_MAG KRON_FLUX KRON_FLUX_ERR FLAGS SKY SKY_SIGMA EXT_NSIGMA"829 columns = "IPP_IDET X_PSF Y_PSF X_PSF_SIG Y_PSF_SIG PSF_INST_MAG PSF_INST_MAG_SIG PEAK_FLUX_AS_MAG PSF_MAJOR PSF_MINOR PSF_THETA EXT_NSIGMA PSF_QF MOMENTS_XX MOMENTS_XY MOMENTS_YY AP_MAG KRON_FLUX KRON_FLUX_ERR FLAGS FLAGS2 SKY SKY_SIGMA EXT_NSIGMA PSF_QF_PERFECT PSF_CHISQ" 760 830 761 831 return super(DetectionBatch, self).importIppTables(columns, regex) -
trunk/ippToPsps/jython/objectbatch.py
r34441 r34615 114 114 for filter in filters: 115 115 116 sql = "UPDATE Object \117 JOIN " + cpsTable + " AS cps ON (cps.row = (Object.row*" + str(filterCount) + ")-(" + str(filterCount) + " - " + str(filter[0]) + ")) \116 sql = "UPDATE Object JOIN \ 117 " + cpsTable + " AS cps ON (cps.row = (Object.row* " + str(filterCount) + ")-(" + str(filterCount) + " - " + str(filter[0]) + ")) \ 118 118 SET \ 119 119 n" + filter[1] + " = NCODE \ … … 122 122 ," + filter[1] + "MeanKronMag = MAG_KRON \ 123 123 ," + filter[1] + "MeanKronMagErr = MAG_KRON_ERR \ 124 ," + filter[1] + "StackPSF Flux = FLUX_PSF\125 ," + filter[1] + "StackPSF FluxErr = FLUX_PSF_ERR\124 ," + filter[1] + "StackPSFMag = -2.5 * log10(FLUX_PSF) + 8.9 \ 125 ," + filter[1] + "StackPSFMagErr = FLUX_PSF_ERR / FLUX_PSF \ 126 126 ," + filter[1] + "StackKronMag = -2.5 * log10(FLUX_KRON) + 8.9 \ 127 127 ," + filter[1] + "StackKronMagErr = FLUX_KRON_ERR / FLUX_KRON \ 128 128 ," + filter[1] + "20pct = MAG_20/1000 \ 129 129 ," + filter[1] + "80pct = MAG_80/1000 \ 130 ," + filter[1] + "Flags = FLAGS " 131 130 ," + filter[1] + "Flags = FLAGS \ 131 ," + filter[1] + "StackDetectID = STACK_ID" 132 133 #stackDetectID needs to be filled still for SA10 - currently at STACK_ID for test purposes 134 # everthing should be correctly calculated as mags for SA10 135 #self.logger.info(sql) 136 self.scratchDb.execute(sql) 137 138 self.logger.infoPair("Calculating nDetections from", "n[filters]") 139 for filter in filters: 140 # now do a sum of n[filters], but do not include the ones with -999 141 sql = "UPDATE Object \ 142 SET nDetections = nDetections + n" + filter[1] + " \ 143 WHERE n" + filter[1] + " != -999" 144 self.scratchDb.execute(sql) 132 145 133 self.scratchDb.execute(sql) 146 147 134 148 135 149 # now set to null all MeanMagErr values > 0.5 (cut set by Gene, 2012-04-12) … … 204 218 ,decErr \ 205 219 ,nDetections \ 206 , extra\220 ,qfPerfect \ 207 221 ,sgSep \ 208 222 ,dataRelease \ -
trunk/ippToPsps/jython/stackbatch.py
r34441 r34615 89 89 self.scratchDb.dropTable("SkinnyObject") 90 90 self.scratchDb.dropTable("ObjectCalColor") 91 self.scratchDb.dropTable("StackDetectionCalib") 91 92 92 93 # delete IPP tables … … 140 141 141 142 sql = "UPDATE StackApFlx AS a, SkyChip_xrad AS b SET \ 142 "+prefix+"flxR1=b.APER_FLUX_1 , \143 "+prefix+"flxR1Err=b.APER_FLUX_ERR_1 , \144 "+prefix+"flxR1Std=b.APER_FLUX_STDEV_1 , \143 "+prefix+"flxR1=b.APER_FLUX_1 / " + str(self.expTime) + ", \ 144 "+prefix+"flxR1Err=b.APER_FLUX_ERR_1 / " + str(self.expTime) + ", \ 145 "+prefix+"flxR1Std=b.APER_FLUX_STDEV_1 / " + str(self.expTime) + ", \ 145 146 "+prefix+"flxR1Fill=b.APER_FILL_1, \ 146 "+prefix+"flxR2=b.APER_FLUX_2 , \147 "+prefix+"flxR2Err=b.APER_FLUX_ERR_2 , \148 "+prefix+"flxR2Std=b.APER_FLUX_STDEV_2 , \147 "+prefix+"flxR2=b.APER_FLUX_2 / " + str(self.expTime) + ", \ 148 "+prefix+"flxR2Err=b.APER_FLUX_ERR_2 / " + str(self.expTime) + ", \ 149 "+prefix+"flxR2Std=b.APER_FLUX_STDEV_2 / " + str(self.expTime) + ", \ 149 150 "+prefix+"flxR2Fill=b.APER_FILL_2, \ 150 "+prefix+"flxR3=b.APER_FLUX_3 , \151 "+prefix+"flxR3Err=b.APER_FLUX_ERR_3 , \152 "+prefix+"flxR3Std=b.APER_FLUX_STDEV_3 , \151 "+prefix+"flxR3=b.APER_FLUX_3 / " + str(self.expTime) + ", \ 152 "+prefix+"flxR3Err=b.APER_FLUX_ERR_3 / " + str(self.expTime) + ", \ 153 "+prefix+"flxR3Std=b.APER_FLUX_STDEV_3 / " + str(self.expTime) + ", \ 153 154 "+prefix+"flxR3Fill=b.APER_FILL_3, \ 154 "+prefix+"flxR4=b.APER_FLUX_4 , \155 "+prefix+"flxR4Err=b.APER_FLUX_ERR_4 , \156 "+prefix+"flxR4Std=b.APER_FLUX_STDEV_4 , \155 "+prefix+"flxR4=b.APER_FLUX_4 / " + str(self.expTime) + ", \ 156 "+prefix+"flxR4Err=b.APER_FLUX_ERR_4 / " + str(self.expTime) + ", \ 157 "+prefix+"flxR4Std=b.APER_FLUX_STDEV_4 / " + str(self.expTime) + ", \ 157 158 "+prefix+"flxR4Fill=b.APER_FILL_4, \ 158 "+prefix+"flxR5=b.APER_FLUX_5 , \159 "+prefix+"flxR5Err=b.APER_FLUX_ERR_5 , \160 "+prefix+"flxR5Std=b.APER_FLUX_STDEV_5 , \159 "+prefix+"flxR5=b.APER_FLUX_5 / " + str(self.expTime) + ", \ 160 "+prefix+"flxR5Err=b.APER_FLUX_ERR_5 / " + str(self.expTime) + ", \ 161 "+prefix+"flxR5Std=b.APER_FLUX_STDEV_5 / " + str(self.expTime) + ", \ 161 162 "+prefix+"flxR5Fill=b.APER_FILL_5, \ 162 "+prefix+"flxR6=b.APER_FLUX_6 , \163 "+prefix+"flxR6Err=b.APER_FLUX_ERR_6 , \164 "+prefix+"flxR6Std=b.APER_FLUX_STDEV_6 , \163 "+prefix+"flxR6=b.APER_FLUX_6 / " + str(self.expTime) + ", \ 164 "+prefix+"flxR6Err=b.APER_FLUX_ERR_6 / " + str(self.expTime) + ", \ 165 "+prefix+"flxR6Std=b.APER_FLUX_STDEV_6 / " + str(self.expTime) + ", \ 165 166 "+prefix+"flxR6Fill=b.APER_FILL_6, \ 166 "+prefix+"flxR7=b.APER_FLUX_7 , \167 "+prefix+"flxR7Err=b.APER_FLUX_ERR_7 , \168 "+prefix+"flxR7Std=b.APER_FLUX_STDEV_7 , \167 "+prefix+"flxR7=b.APER_FLUX_7 / " + str(self.expTime) + ", \ 168 "+prefix+"flxR7Err=b.APER_FLUX_ERR_7 / " + str(self.expTime) + ", \ 169 "+prefix+"flxR7Std=b.APER_FLUX_STDEV_7 / " + str(self.expTime) + ", \ 169 170 "+prefix+"flxR7Fill=b.APER_FILL_7, \ 170 "+prefix+"flxR8=b.APER_FLUX_8 , \171 "+prefix+"flxR8Err=b.APER_FLUX_ERR_8 , \172 "+prefix+"flxR8Std=b.APER_FLUX_STDEV_8 , \171 "+prefix+"flxR8=b.APER_FLUX_8 / " + str(self.expTime) + ", \ 172 "+prefix+"flxR8Err=b.APER_FLUX_ERR_8 / " + str(self.expTime) + ", \ 173 "+prefix+"flxR8Std=b.APER_FLUX_STDEV_8 / " + str(self.expTime) + ", \ 173 174 "+prefix+"flxR8Fill=b.APER_FILL_8, \ 174 "+prefix+"flxR9=b.APER_FLUX_9 , \175 "+prefix+"flxR9Err=b.APER_FLUX_ERR_9 , \176 "+prefix+"flxR9Std=b.APER_FLUX_STDEV_9 , \175 "+prefix+"flxR9=b.APER_FLUX_9 / " + str(self.expTime) + ", \ 176 "+prefix+"flxR9Err=b.APER_FLUX_ERR_9 / " + str(self.expTime) + ", \ 177 "+prefix+"flxR9Std=b.APER_FLUX_STDEV_9 / " + str(self.expTime) + ", \ 177 178 "+prefix+"flxR9Fill=b.APER_FILL_9, \ 178 "+prefix+"flxR10=b.APER_FLUX_10 , \179 "+prefix+"flxR10Err=b.APER_FLUX_ERR_10 , \180 "+prefix+"flxR10Std=b.APER_FLUX_STDEV_10 , \179 "+prefix+"flxR10=b.APER_FLUX_10 / " + str(self.expTime) + ", \ 180 "+prefix+"flxR10Err=b.APER_FLUX_ERR_10 / " + str(self.expTime) + ", \ 181 "+prefix+"flxR10Std=b.APER_FLUX_STDEV_10 / " + str(self.expTime) + ", \ 181 182 "+prefix+"flxR10Fill=b.APER_FILL_10 \ 182 183 WHERE a.ippDetectID=b.IPP_IDET AND b.PSF_FWHM "+psfCondition … … 196 197 sql = "UPDATE StackModelFit AS a, SkyChip_xfit AS b SET \ 197 198 "+model+"Radius=b.EXT_WIDTH_MAJ, \ 198 "+model+" Mag=b.EXT_INST_MAG, \199 "+model+" MagErr=b.EXT_INST_MAG_SIG, \199 "+model+"Flux=b.EXT_INST_MAG, \ 200 "+model+"FluxErr=b.EXT_INST_MAG_SIG, \ 200 201 "+model+"Ab=b.EXT_WIDTH_MAJ/b.EXT_WIDTH_MIN, \ 201 202 "+model+"Phi=b.EXT_THETA, \ … … 229 230 "+modelLong+"Covar77=b.EXT_COVAR_06_06 \ 230 231 WHERE a.ippDetectID=b.IPP_IDET AND b.MODEL_TYPE = '"+ippModelType+"'" 232 233 #SA10 convert mag to flux in calculation above 231 234 232 235 self.scratchDb.execute(sql) … … 319 322 ,skyCellID \ 320 323 ,projectionCellID \ 321 ,obsTime \322 324 ,xPos \ 323 325 ,yPos \ 324 326 ,xPosErr \ 325 327 ,yPosErr \ 326 , instFlux \327 , instFluxErr \328 ,psfFlux \ 329 ,psfFluxErr \ 328 330 ,peakFlux \ 329 331 ,sky \ … … 335 337 ,psfLikelihood \ 336 338 ,infoFlag \ 337 ,psfCf \ 339 ,psfQf \ 340 ,psfQfPerfect \ 338 341 ,momentXX \ 339 342 ,momentXY \ … … 345 348 ,momentR1 \ 346 349 ,momentRH \ 347 ,ap Mag\348 ,ap MagErr \350 ,apFlux \ 351 ,apFluxErr \ 349 352 ,kronFlux \ 350 353 ,kronFluxErr \ 351 354 ,kronRad \ 352 , kronRadErr\355 ,psfChiSq \ 353 356 ,nFrames \ 354 357 ,assocDate \ … … 359 362 ," + self.skycell + " \ 360 363 ," + self.projectioncell + " \ 361 ," + self.header['MJD-OBS'] + " \362 364 ,X_PSF \ 363 365 ,Y_PSF \ … … 374 376 ,PSF_THETA \ 375 377 ,psfLikelihood(EXT_NSIGMA) \ 376 , FLAGS2 << 32 | FLAGS \378 , FLAGS2 << 32 | FLAGS \ 377 379 ,PSF_QF \ 380 ,PSF_QF_PERFECT \ 378 381 ,MOMENTS_XX \ 379 382 ,MOMENTS_XY \ … … 389 392 ,KRON_FLUX / " + str(self.expTime) + " \ 390 393 ,KRON_FLUX_ERR / " + str(self.expTime) + " \ 391 , MOMENTS_R1 * 2.5 \392 , NULL\394 ,MOMENTS_R1 * 2.5 \ 395 ,PSF_CHISQ \ 393 396 ,N_FRAMES \ 394 397 , '" + self.dateStr + "' \ 395 398 ," + self.historyModNum + " \ 396 399 FROM SkyChip_psf" 400 401 #SA10 to do 402 # put in correct formula for apFlux, psfFlux, psfLikelihood, sgsep 403 # need correct units of anything in Flux 397 404 self.scratchDb.execute(sql) 398 405 … … 417 424 self.updateStackTypeID("StackDetection") 418 425 419 self.updateDvoIDs ("StackDetection")426 self.updateDvoIDsAndFlags("StackDetection") 420 427 421 428 sql = "ALTER IGNORE TABLE StackDetection ADD PRIMARY KEY (objID)" … … 427 434 #if deep stack and instFlux = null and err not null 428 435 sql = "UPDATE StackDetection AS a, SkyChip_psf AS b \ 429 SET instFlux = 2*b.PSF_INST_FLUX_SIG / " + str(self.expTime) + " \430 WHERE instFlux IS NULL \436 SET psfFlux = 2*b.PSF_INST_FLUX_SIG / " + str(self.expTime) + " \ 437 WHERE psfFlux IS NULL \ 431 438 AND a.ippDetectID = b.IPP_IDET \ 432 439 AND b.PSF_INST_FLUX_SIG IS NOT NULL" … … 434 441 self.scratchDb.execute(sql) 435 442 436 437 self.scratchDb.reportAndDeleteRowsWithNULLS("StackDetection", "instFlux")443 #leave null instflux in 444 #self.scratchDb.reportAndDeleteRowsWithNULLS("StackDetection", "instFlux") 438 445 439 446 self.scratchDb.reportAndDeleteRowsWithNULLS("StackDetection", "objID") … … 474 481 petRadius=b.PETRO_RADIUS \ 475 482 ,petRadiusErr=b.PETRO_RADIUS_ERR \ 476 ,pet Mag=b.PETRO_MAG \477 ,pet MagErr=b.PETRO_MAG_ERR \483 ,petFlux=b.PETRO_MAG \ 484 ,petFluxErr=b.PETRO_MAG_ERR \ 478 485 ,petR50=b.PETRO_RADIUS_50 \ 479 486 ,petR50Err=b.PETRO_RADIUS_50_ERR \ … … 482 489 ,petCf=b.PETRO_FILL \ 483 490 WHERE a.ippDetectID=b.IPP_IDET" 491 #SA10 Put correct calculation for PetFlux 484 492 self.scratchDb.execute(sql) 485 493 … … 586 594 self.scratchDb.updateAllRows("SkinnyObject", "surveyID", str(self.surveyID)) 587 595 self.scratchDb.updateAllRows("SkinnyObject", "dataRelease", str(self.config.dataRelease)) 596 597 ''' 598 Populates the StackDetectionCalib table 599 ''' 600 def populateStackDetectionCalib(self): 601 self.logger.infoPair("Processing table", "StackDetectionCalib") 602 tableName = "StackDetectionCalib" 603 # insert all detections into table 604 sql = "INSERT INTO " + tableName + " ( \ 605 objID \ 606 ,stackDetectID \ 607 ,ippObjID \ 608 ,ippDetectID \ 609 ,filterID \ 610 ,surveyID \ 611 ) \ 612 SELECT \ 613 objID \ 614 ,stackDetectID \ 615 ,ippObjID \ 616 ,ippDetectID \ 617 ,filterID \ 618 ,surveyID \ 619 FROM StackDetection" 620 self.logger.info(sql) 621 self.scratchDb.execute(sql) 622 #SA10 TODO 623 # insert ra/dec/calstuff into Table 624 # heather discovered dec is called 'dec_' 625 sql = "UPDATE " + tableName + " SET ra = -999, dec_ =-999, raErr = -999, decErr = -999, zp = -999, zpErr = -999, expTime = -999, airMass = -999" 626 self.logger.info(sql) 627 self.scratchDb.execute(sql) 628 self.scratchDb.updateAllRows(tableName, "dataRelease", str(self.config.dataRelease)) 629 630 631 588 632 589 633 ''' … … 649 693 WHERE a.ippDetectID = b.ippDetectID \ 650 694 AND b.imageID = " + str(imageID) 651 652 self.scratchDb.execute(sql) 653 695 self.scratchDb.execute(sql) 696 697 ''' 698 Updates Flags 699 ''' 700 def updateDvoIDsAndFlags(self, table): 701 self.logger.infoPair("getting","imageID") 702 imageID = self.scratchDb.getImageIDFromExternID(self.header['IMAGEID']) 703 self.logger.infoPair("obtained","imageID") 704 self.logger.debug("Updating table '" + table + "' with DVO IDs...") 705 sql = "UPDATE " + table + " AS a, " + self.scratchDb.dvoDetectionTable + " AS b SET \ 706 a.ippObjID = b.ippObjID, \ 707 a.stackDetectID = 10 * b.detectID + " + str(self.filterID) + ", \ 708 a.objID = b.objID, \ 709 a.infoFlag = b.flags << 46 | a.infoFlag \ 710 WHERE a.ippDetectID = b.ippDetectID \ 711 AND b.imageID = " + str(imageID) 712 self.scratchDb.execute(sql) 713 714 654 715 ''' 655 716 Does the processing, i.e. pulling stuff from IPP tables into PSPS tables … … 660 721 if not self.getIDsFromDVO(): 661 722 return False 662 723 #self.logger.infoPair("populating","StackMeta") 663 724 self.populateStackMeta() 725 #self.logger.infoPair("populating","StackDetection") 664 726 self.populateStackDetection() 665 727 666 728 if self.stackType != "NIGHTLY_STACK": 729 #self.logger.infoPair("populating","StackModelFit") 667 730 self.populateStackModelFit() 731 #self.logger.infoPair("populating","StackApFlx") 668 732 self.populateStackApFlx() 669 733 #self.logger.infoPair("populating","StackToImage") 670 734 self.populateStackToImage() 735 #self.logger.infoPair("populating","SkinnyObject") 671 736 self.populateSkinnyObject() 737 #self.logger.infoPair("populating","ObjectCalColor") 672 738 self.populateObjectCalColor() 739 #self.logger.infoPair("populating","StackDetectionCalib") 740 self.populateStackDetectionCalib() 673 741 674 742 self.setMinMaxObjID(["StackDetection"])
Note:
See TracChangeset
for help on using the changeset viewer.
