Changeset 35456 for trunk/ippToPsps
- Timestamp:
- May 1, 2013, 5:52:45 AM (13 years ago)
- Location:
- trunk/ippToPsps
- Files:
-
- 4 edited
-
doc/notes.txt (modified) (1 diff)
-
jython/detectionbatch.py (modified) (2 diffs)
-
test/fulltest.sh (modified) (3 diffs)
-
test/mkgpc1data.dvo (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippToPsps/doc/notes.txt
r35233 r35456 1 2 2013.04.27 3 4 need to move the conversion between mags and flux into dvopsps. 5 6 * detection table 7 * psfFlux : uses PSF_INST_MAG instead of PSF_INST_FLUX 8 * psfFluxErr : uses PSF_INST_MAG_SIG and PSF_INST_MAG instead of PSF_INST_FLUX_SIG 9 * apFlux : uses AP_MAG instead of AP_FLUX 10 * kronFlux : correctly uses KRON_FLUX 11 * kronFluxErr : correctly uses KRON_FLUX_ERR 12 13 ** psfFlux is forced to 0.0000001 instead of 0.0 -- why? 14 15 * stack detection table (same issues as detection table) 16 * psfFlux : uses PSF_INST_MAG instead of PSF_INST_FLUX 17 * psfFluxErr : uses PSF_INST_MAG_SIG and PSF_INST_MAG instead of PSF_INST_FLUX_SIG 18 * apFlux : uses AP_MAG instead of AP_FLUX 19 * kronFlux : correctly uses KRON_FLUX 20 * kronFluxErr : correctly uses KRON_FLUX_ERR 21 22 * psfFlux is set to 2 sigma if psfFlux is NULL 23 24 * object 25 * Mean PSF Mag : direct from DVO MAG 26 * Mean PSF Mag Err : direct from DVO MAG_ERR 27 * Mean Kron Mag : direct from DVO KRON MAG 28 * Mean Kron Mag Err : direct from DVO KRON MAG_ERR 29 30 * Stack PSF Mag : from DVO stack psf flux 31 * Stack PSF Mag Err : from DVO stack psf flux err 32 * Stack Kron Mag : from DVO stack kron flux 33 * Stack Kron Mag Err : from DVO stack kron flux err 34 35 36 37 38 39 old notes: 1 40 2 41 example of joining ipptopsps and gpc1 tables: -
trunk/ippToPsps/jython/detectionbatch.py
r35452 r35456 483 483 # something like (f < 0.0) ? -999 : -2.5*log10(f) 484 484 # as a result, the negative fluxes here result in floating point errors 485 sql = "UPDATE " + pspsTableName " SET psfFlux = 1e20 WHERE psfFlux <= 0.0"486 self.scratchDb.execute(sql) 487 488 sql = "UPDATE " + pspsTableName " SET apFlux = 1e20 WHERE apFlux <= 0.0"489 self.scratchDb.execute(sql) 490 491 sql = "UPDATE " + pspsTableName " SET kronFlux = 1e20 WHERE kronFlux <= 0.0"485 sql = "UPDATE " + pspsTableName + " SET psfFlux = 1e20 WHERE psfFlux <= 0.0" 486 self.scratchDb.execute(sql) 487 488 sql = "UPDATE " + pspsTableName + " SET apFlux = 1e20 WHERE apFlux <= 0.0" 489 self.scratchDb.execute(sql) 490 491 sql = "UPDATE " + pspsTableName + " SET kronFlux = 1e20 WHERE kronFlux <= 0.0" 492 492 self.scratchDb.execute(sql) 493 493 … … 1035 1035 1036 1036 # XXX EAM NOTE : this is fragile : requires PS1_V4 1037 columns = "IPP_IDET X_PSF Y_PSF X_PSF_SIG Y_PSF_SIG PSF_INST_MAG PSF_INST_MAG_SIG P EAK_FLUX_AS_MAG PSF_MAJOR PSF_MINOR PSF_THETA EXT_NSIGMA PSF_QF MOMENTS_XX MOMENTS_XY MOMENTS_YY AP_MAGKRON_FLUX KRON_FLUX_ERR FLAGS FLAGS2 SKY SKY_SIGMA EXT_NSIGMA PSF_QF_PERFECT PSF_CHISQ"1037 columns = "IPP_IDET X_PSF Y_PSF X_PSF_SIG Y_PSF_SIG PSF_INST_MAG PSF_INST_MAG_SIG PSF_INST_FLUX PSF_INST_FLUX_SIG PEAK_FLUX_AS_MAG PSF_MAJOR PSF_MINOR PSF_THETA EXT_NSIGMA PSF_QF MOMENTS_XX MOMENTS_XY MOMENTS_YY AP_MAG AP_FLUX KRON_FLUX KRON_FLUX_ERR FLAGS FLAGS2 SKY SKY_SIGMA EXT_NSIGMA PSF_QF_PERFECT PSF_CHISQ" 1038 1038 1039 1039 return super(DetectionBatch, self).importIppTables(columns, regex) -
trunk/ippToPsps/test/fulltest.sh
r35417 r35456 1 1 #!/bin/csh -f 2 3 chmod +x *.dvo *.sh 2 4 3 5 # a bit crude : get the datadir path by interpollation of the PATH element for psconfig … … 23 25 set objectbatch = 1 24 26 25 set args = "" 26 27 while ("$1" != "") 28 switch ("$1") 27 set targ = ($argv) 28 set args = "" 29 30 while ($#targ) 31 switch ("$targ[1]") 32 case -h 33 case -help 34 echo "USAGE: fulltest.sh [-XX] [+XX]" 35 echo " -XX : turn off features" 36 echo " +XX : turn on features" 37 echo " features: all, mkgpc1, initdb, initscratch, initbatch, camera, stack, object" 38 exit 2 39 case -all 40 set mkgpc1 = 0 41 set initdb = 0 42 set initscratch = 0 43 set initbatch = 0 44 set camqueue = 0 45 set cambatch = 0 46 set stackqueue = 0 47 set stackbatch = 0 48 set objectqueue = 0 49 set objectbatch = 0 50 breaksw; 29 51 case -mkgpc1 30 52 set mkgpc1 = 0 … … 52 74 breaksw; 53 75 default: 54 set args=($args $ 1);76 set args=($args $targ[1]); 55 77 breaksw; 56 78 endsw 57 shift 79 shift targ 58 80 end 81 82 set targ = ($args) 83 set args = "" 84 85 while ($#targ) 86 switch ("$targ[1]") 87 case +all 88 set mkgpc1 = 1 89 set initdb = 1 90 set initscratch = 1 91 set initbatch = 1 92 set camqueue = 1 93 set cambatch = 1 94 set stackqueue = 1 95 set stackbatch = 1 96 set objectqueue = 1 97 set objectbatch = 1 98 breaksw; 99 case +mkgpc1 100 set mkgpc1 = 1 101 breaksw; 102 case +initdb 103 set initdb = 1 104 breaksw; 105 case +initscratch 106 set initscratch = 1 107 breaksw; 108 case +initbatch 109 set initbatch = 1 110 breaksw; 111 case "+camera" 112 set camqueue = 1 113 set cambatch = 1 114 breaksw; 115 case +stack 116 set stackqueue = 1 117 set stackbatch = 1 118 breaksw; 119 case +object 120 set objectqueue = 1 121 set objectbatch = 1 122 breaksw; 123 default: 124 set args=($args $targ[1]); 125 breaksw; 126 endsw 127 shift targ 128 end 129 130 if ($args[1] != "") then 131 echo "USAGE: fulltest.sh [-XX] [+XX]" 132 echo " -XX : turn off features" 133 echo " +XX : turn on features" 134 echo " features: all, mkgpc1, initdb, initscratch, initbatch, camera, stack, object" 135 exit 2 136 endif 59 137 60 138 set stackqueuename = 0 -
trunk/ippToPsps/test/mkgpc1data.dvo
r35417 r35456 37 37 macro mkfull 38 38 init.db 39 mkcatdir.cam PS1_V3 PS1_V4 40 mkcatdir.stk PS1_V3 PS1_V4 39 # mkcatdir.cam PS1_V3 PS1_V4 40 # mkcatdir.stk PS1_V3 PS1_V4 41 mkcatdir.cam PS1_V4 PS1_V4 42 mkcatdir.stk PS1_V4 PS1_V4 41 43 insert.stack.set 42 44 end
Note:
See TracChangeset
for help on using the changeset viewer.
