- Timestamp:
- Jun 23, 2012, 6:56:38 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20120601/psphot/test/tap_psphot_varmodel.pro
r34055 r34056 23 23 # variance mode: CONSTANT, IMAGE_VAR, MODEL_VAR 24 24 25 macro ppsimtest 26 data $1 27 # sf1 : measured flux inserted in image 28 # sm2 : theoretical magnitude for star 29 read x 1 y 2 sf1 3 sm2 5 Io 12 30 set sm1 = -2.5*log(sf1) 31 set dsm = sm1 - sm2 32 lim -n 2 sm1 -0.2 0.2; clear; box; plot sm1 dsm 33 end 34 25 35 macro go.one 26 36 mkdir test … … 43 53 sprintf name "test/test.%02d.%02d.%s" $sky {10*$fwhm} $psf_in 44 54 mkexp $name $sky $fwhm $psf_in 55 data $name.dat 56 read fluxIn 3 MagPredIn 5 57 set MagRealIn = -2.5*log(fluxIn) 58 set dMagIn = MagRealIn - MagPredIn 45 59 46 60 foreach psf_out GAUSS PS1_V1 … … 256 270 section a1 0.0 0.50 1.0 0.25 257 271 lim rv -0.025 0.075; label -fn courier 14; box; 272 plot -c black -pt 0 -sz 0.3 MagPredIn dMagIn 258 273 plot -c red -pt 7 -sz 2.0 imag_V dmag_V 259 274 plot -c darkgreen -pt 3 -sz 2.0 imag_V Dmag_V … … 264 279 sprintf line "OFFSET: %6.3f" $offset 265 280 textline -frac 0.1 0.8 -fn courier 24 "$line" 281 266 282 267 283 set lChiNorm = log(PSF_CHISQ_out / PSF_NDOF_out)
Note:
See TracChangeset
for help on using the changeset viewer.
