;
;
;

othertemp=strarr(100)
itime2=[0.299,0.599,0.999,1.8,3.0,6.0,10.,18.,30.,60.]

;goto,plot2
;goto,plot1
;flist=file_search('20050306/lt*q1.fit',count=numfiles)
;flist=file_search('20050313/lt_dt/ltb*q1.fit',count=numfiles)
;flist=file_search('20050313/lt_dt/lt.0*q1.fit',count=numfiles)
flist=file_search('20050325/lt6_dt/lt6_full*q1.fit',count=numfiles)
window,0,xs=300,ys=300,retain=2
;
filetimes=fltarr(numfiles)
counts=fltarr(numfiles)
itime=fltarr(numfiles)
data=fltarr(256,256,numfiles)

for i=0,numfiles-1 do begin
 print,flist(i)
 r=fitshline(file=flist(i),error=e,KEY='ITIME')
 itime(i)=float(strmid(r,11,8))
 r=fitshline(file=flist(i),error=e,KEY='OTHERT')
 othertemp(i)=r
 temp=readfits(flist(i))
 data(*,*,i)=temp
; tvhist,rebin(temp,256,256)
; wait,1.0
 counts(i)=mean(temp)
endfor

counts2=fltarr(256,256,10)
noise=fltarr(256,256,10)

for m=0,9 do begin
 indx1=m*5
 indx2=(m+1)*5-1
 for i=0,255 do begin
  for j=0,255 do begin
   noise(i,j,m)=stdev([data(i,j,indx1:indx2)])
   counts2(i,j,m)=mean([data(i,j,indx1:indx2)])
  endfor
 endfor
endfor

n_noise=fltarr(10)
s_signal=fltarr(10)

plot2:

for i=0,9 do begin
 n_noise(i)=median(noise(128:255,*,i))
 s_signal(i)=median(counts2(128:255,*,i))
endfor

lognoise=alog(n_noise)/alog(10.)
logsignal=alog(s_signal)/alog(10.)

cgain=poly_fit(logsignal(7:8),lognoise(7:8),1,ygain)
cread=poly_fit(logsignal(0:1),lognoise(0:1),1,yread)

fitgain=poly(logsignal,cgain)
fitgain2=poly(findgen(400)/100.,[-0.6021,0.5])
fitread=poly(logsignal,cread)
;fitread2=poly(logsignal,[lognoise(0),0.0])
fitread2=poly(findgen(400)/100.,[0.514692,0.0])

set_plot,'ps'
device,/encapsulated,filename='transfer_test_mar25.eps'
plot,logsignal,lognoise,$
     xran=[0,4.5],xstyl=1,$
     yran=[0.4,2.0],ystyl=1,$
     psym=1,$
     xtit='Log Signal',ytit='Log Noise',tit='NAC Light Transfer',$
     chars=1.6
;oplot,logsignal,fitgain,linestyl=2
oplot,findgen(400)/100.,fitgain2,linestyl=2
;oplot,logsignal,fitread,linestyl=1
oplot,findgen(400)/100.,fitread2,linestyl=1

n_noise1=n_noise
s_signal1=s_signal
;openr,1,'light_trans_20050313_dark.txt'
;readf,1,n_noise1,s_signal1
;close,1
;oplot,alog(s_signal1)/alog(10.),alog(n_noise1)/alog(10.),psym=4
device,/close
set_plot,'x'
;openw,1,'light_trans_20050313_dark.txt'
;printf,1,n_noise,s_signal
;close,1

goto,the_end

;plot1:
;set_plot,'ps'
;device,/encapsulated,filename='transfer_test_mar05.eps'
;plot,alog([counts2])/alog(10.),alog([noise])/alog(10.),psym=3,$
;     yran=[-1,2.],ystyl=1,xran=[1.,4.],xstyl=1,$
;     xtit='Log Signal',ytit='Log Noise',tit='NAC Light Transfer',$
;     chars=1.6
;device,/close
;set_plot,'x'

the_end:
end