;made by H.Watanabe 2007.3.13;pro dark_countspawn,'ls capture200?-??-??-??-??-??-???.tif',fl,count=nnprint,nn/16+' set'if((nn mod 16) ne 0) then begin  print,'Data set is not correct!!'  stopendif;-;-; calculate average and standard deviation ;-;-;heikin=fltarr(nn)hensa=fltarr(nn)for i=0,nn-1 do begin  a=read_tiff(fl[i]) & a=a/16  aa=moment(a)  heikin[i]=aa(0)  hensa[i]=sqrt(aa(1))endfor;-;-; draw graph ;-;-;window,1,xs=600,ys=500plot,heikin,psym=4,thick=2,chars=2,$  background=255,color=0,pos=[0.,0.1,1.,1.],$  xticks=8,xtickname=[0,8,0,8,0,8,0,8,0,8,0,8,0,8,0,8,0,8,0,8,0,8,0,8,0,8],$  xtitle="Frame No.(time)",ytitle="Dark Count"xxx=fltarr(2)error=fltarr(2);gr=intarr(2);gr(0)=0 & gr(1)=1000for i=0,nn-1 do beginxxx(0)=i & xxx(1)=ierror(0)=heikin[i]+hensa[i]error(1)=heikin[i]-hensa[i]oplot,xxx,error,color=0;if(i mod 16 eq 0) then begin;  oplot,xxx,gr,linestyle=2,color=126;endifendfor;-;-; draw times ;-;-;;xyouts,'Dark Count', , ,/dev,color=255,chars=2;xyouts,'Frame No.(time)', , ,/dev,color=255,chars=2plots,[,30],[,470],/dev,linestyle=1,color=125xyouts,'15:00', , ,/dev,color=255,chars=2xyouts,'15:10', , ,/dev,color=255,chars=2xyouts,'15:20', , ,/dev,color=255,chars=2xyouts,'15:30', , ,/dev,color=255,chars=2xyouts,'15:40', , ,/dev,color=255,chars=2xyouts,'15:50', , ,/dev,color=255,chars=2xyouts,'16:00', , ,/dev,color=255,chars=2xyouts,'16:10', , ,/dev,color=255,chars=2xyouts,'16:20', , ,/dev,color=255,chars=2xyouts,'16:30', , ,/dev,color=255,chars=2xyouts,'16:40', , ,/dev,color=255,chars=2xyouts,'16:50', , ,/dev,color=255,chars=2xyouts,'17:00', , ,/dev,color=255,chars=2end</pre></body></html>