;cycles 5
;itime=[0.3, 0.6, 1.0, 1.8, 3.0, 6.0, 10.0, 18.0, 30.0, 60.0, 100.0, 180.0]
;cycles 1
;itime 300.0
;
; data only goes up to 100 sec integration!!!!
;
;

;goto,plot1
flist=file_search('20050306/lt*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))
 temp=readfits(flist(i))
 data(*,*,i)=temp
; tvhist,rebin(temp,256,256)
 counts(i)=mean(temp)
endfor

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

for m=0,10 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


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

;device,/close
;set_plot,'x'

the_end:
end