Changeset 30229

Timestamp:

Jan 9, 2011, 12:47:04 PM (15 years ago)

Author:

eugene

Message:

check completeness of a stack (different EXTNAME and zero point)

File:

• branches/eam_branches/ipp-20101205/psphot/test/tap_psphot_deteff.pro (modified) (1 diff)

branches/eam_branches/ipp-20101205/psphot/test/tap_psphot_deteff.pro

@@ -216 +216 @@
  set ldmag = log(PSF_INST_MAG_SIG_raw)
  section overlay 0 0 1 1; lim dMag -5 1.2; box -ypad 5 +ypad 5 -ticks 1011 -labels 1001; plot PSF_INST_MAG_raw ldmag
+ label +y log(S/N)
+
+ line -c red -20 {log(1/25.0)} to 20 {log(1/25.0)}
+ line -c blue -20 {log(1/5.0)} to 20 {log(1/5.0)}
+
+ png -name $3
+end
+
+macro completestack
+ if ($0 != 6)
+   echo "USAGE: completeness (raw) (out) (output) (radius) (zeropt)"
+   break
+ end
+
+ load.cmf $1 Chip.psf raw
+ load.cmf $2 SkyChip.psf out
+ set X_raw = int(X_PSF_raw) + 0.5
+ set Y_raw = int(Y_PSF_raw) + 0.5
+ set Mraw = PSF_INST_MAG_raw + $5
+ match2d X_PSF_raw Y_PSF_raw X_PSF_out Y_PSF_out $4 -index1 index1 -index2 index2 -closest
+
+ histogram Mraw nMag -16.0 -3.0 0.25 -range dMag
+ set fMag = zero(dMag) 
+ for i 0 {dMag[]-1}
+  set inrange = (Mraw > dMag[$i]) && (Mraw <= dMag[$i+1])
+  subset all = index1 if (inrange)
+  subset got = index1 if (inrange) && (index1 >= 0)
+  if (all[] == 0)
+    fMag[$i] = 0
+  else 
+   fMag[$i] = got[] / all[]
+  end
+ end
+
+ device -n complete
+ resize 1000 600
+
+ clear -s
+
+ section default 0 0 1 1
+ lim dMag fMag; clear; box -ypad 5 +ypad 5 -ticks 1110; plot -x 1 dMag fMag
+ label -x mag_inst -y det_frac 
+
+ set found = (index1 >= 0)
+ plot Mraw found -c red
+
+ set ldmag = log(PSF_INST_MAG_SIG_raw)
+ section overlay 0 0 1 1; lim dMag -5 1.2; box -ypad 5 +ypad 5 -ticks 1011 -labels 1001; plot Mraw ldmag
  label +y log(S/N)