;; The following are "pixview" data that I took on Feb 23, 1998 ;; to verify the linearity of the shutter, which broke prior. ;; It is a way to test what the ratio is between the 5B and 5C opals. ;; The number used in the original data reduction for the 5C opal was ;; 8.0e-8 which is a linear interpolation for the POISE wavelength ;; of 620nm... ;; From orw, 4/30/98 caption ;; The spectral band of the image is set by an Andover filter with ;; central wavelength of 620 nm and a bandwidth of 10 nm. ;; From elmore 17 Mar 1998 ;; Alice, ;; doing a linear interpolation between the values at 530 and 656nm, I ;; get ;; 5A 8.39 e-6 ;; 5B 7.76 e-7 ;; 5C 8.00 e-8 ;; david ;; ;; ;; But, working on the spie paper in May of 2000, David believes ;; the 5C opal number has a typo in it, and the interpolated value ;; should be 8.28e-8. ;; I checked the 5B and 5C opal data with code: cmpBtoC_new.pro ;; and this gives a 5C opal of 5B * 0.100351 = 7.78726e-08 ;; ;; BUT! ;; this program ckBC.pro ;; gives a 5C opal of 5B * 0.110048 = 8.53976e-08 ;; SO! ;; I do not know if I can trust the opal exposures used to make the ;; gain images. CTimes = [ 1., 2., 3., 4., 5., 7.5, 10., 15., 17.5, 1.0, 5.0] CData = [ 5358., 7542., 9831., 11902., 14530., 20215., 25523., 36323., 41901., 5404., 14594.] BTimes = [ 0.5, 1.0, 1.5, 2.0, 1.75, 1.825, 1.89 ] BData = [ 14305., 25359., 34991., 44386., 39907., 41365., 43024.] Cbias = 3262.99 Bbias = 4678.49 CData = CData - Cbias BData = BData - Bbias leastsq, CTimes, CData, Caa, Cbb, CRR CFit = CTimes*Cbb + Caa print,Caa,Cbb,CRR leastsq, BTimes, BData, Baa, Bbb, BRR BFit = BTimes*Bbb + Baa print,Baa,Bbb,BRR !p.multi=[0,1,1] plot,[0.,20.],[0.,45000.],xstyle=1,ystyle=1,/nodata oplot,CTimes,Cdata,psym=1 oplot,CTimes,CFit,thick=2 oplot,BTimes,Bdata,psym=1 oplot,BTimes,BFit,thick=2 opal5B=7.76e-7 opal5C=8.00e-8 print,opal5B * Cbb/Bbb ;;; = 8.53976e-08 NEITHER 8e-8 or 8.28e-8 ! print,' ' Bpoint = Bbb * 1.0 +Baa Cpoint = Cbb * 1.0 +Caa print,Cbb/Bbb ;;; = 0.110048 print,Cpoint/Bpoint ;;; = 0.110049 print,' ' plot,[0.,20.],[0.,45000.],xstyle=1,ystyle=1,/nodata oplot,CTimes,Cdata,psym=2 oplot,BTimes*Bbb/Cbb,Bdata,psym=5 oplot,CTimes,CFit,linestyle=2,thick=2 oplot,BTimes*Bbb/Cbb,BFit end