PRO plotbw1 ;*** Make tables of exit-status counts. The variables which should be ;*** in memory are ;*** infile: file to be read, such as '../bbscan_feb2305.fit' ;*** fname: name of file, w/o path ;*** npt: number of points in vectors ;*** exits: exit status ;*** index: index in the input array ;*** f0 central frequency ;*** alpha negative of slope ;*** phi0 phase at center of call ;*** res2 residuals per degree of freedom ;*** nchi2 normalized chi-squared ;*** rinfile: Roger's file of triggers from scanning the "c" file, such as ;*** /usr/data/bland03/pioneer/2091/284/bbscan.bwc ;*** rfname: file name without path, such as bbscan.bwc ;*** npeaks: number of triggers from Roger ;*** cfile[npeaks]: name of "c" file that event came from (NB: this ;*** variable is a string 50 characters long, possibly with ;*** trailing blanks ;*** iyear=intarr(npeaks) year number (0 to 9, for 2000 through 2009) ;*** cday=dblarr(npeaks) time of "c" peak, in ms ;*** cdb=fltarr(npeaks) height of peak ;*** ccen=fltarr(npeaks) center of noise peak ;*** csig=fltarr(npeaks) standard deviation of noise peak ;*** cthnsig=fltarr(npeaks) threshold in terms of number of sigmas ;*** cthfx=fltarr(npeaks) fixed threshold str1=''; input buffer chrszsav=!p.charsize ;*** Read the data into memory. rdfit restore ihtml=0 read,'Make HTML file? (y or return) ',str1 IF str1 eq 'y' THEN BEGIN ihtml=1 openw,lunhtml,'plotbw1.html',/get_lun printf,lunhtml,'' printf,lunhtml,'Exit Status Counts printf,lunhtml,'' printf,lunhtml,'
Exit Status Counts. Bits 0-4 indicate', $ ' convergence in five LSQ fits, with 200, 600, 1000, 3000, and ', $ ' 6000 points; bit 5 indicates that the "cheating fit" converged, ', $ 'and had a distinctly better chi-squared than the fifth fit. Central ', $ 'means f0 in [15.7,16.3] and alpha in [-.1,+.2].' ENDIF IF ihtml THEN printf,lunhtml,'
' ;*** Exit status, all fits. nval=64 ibin=histogram(exits,min=0,nbins=nval) print,'Exit status, all fits ('+string(n_elements(exits),format= $ '(i6)')+' events).' IF ihtml THEN BEGIN printf,lunhtml,'' printf,lunhtml,'', $ '' ENDIF nline=nval/8 FOR iline=0,nline-1 DO BEGIN print,ibin(iline*8:iline*8+7),total(ibin(iline*8:iline*8+7)), $ format='(8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,' FOR icol=0,7 DO printf,lunhtml,' printf,lunhtml,'' printf,lunhtml,' ENDIF ENDFOR itots=intarr(8) FOR icol=0,7 DO itots[icol]=fix(total(ibin[indgen(nline)*8+icol])) print,itots,fix(total(itots)),format='(/8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,'', $ '' printf,lunhtml,'' FOR icol=0,7 DO printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'
Exit status, all fits ('+ $ string(n_elements(exits),format='(i6)')+' events).' printf,lunhtml,'
0001020304050607total
',ibin(iline*8+icol),'*',fix(total(ibin(iline*8:iline*8+7))),'
**********
',itots(icol),'*',fix(total(itots)),'
' ENDIF good=nchi2 lt .7 central=f0 gt 15.7 and f0 lt 16.3 and alpha gt -.1 and alpha lt 0.2 IF ihtml THEN printf,lunhtml,'		' ;*** Exit status, nchi2 < .7. nval=64 ix=where(good) ibin=histogram(exits[ix],min=0,nbins=nval) print,'Exit status, nchi2 < 0.7('+string(n_elements(ix),format= $ '(i6)')+' events).' IF ihtml THEN BEGIN printf,lunhtml,'' printf,lunhtml,'', $ '' ENDIF nline=nval/8 FOR iline=0,nline-1 DO BEGIN print,ibin(iline*8:iline*8+7),total(ibin(iline*8:iline*8+7)), $ format='(8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,' FOR icol=0,7 DO printf,lunhtml,' printf,lunhtml,'' printf,lunhtml,' ENDIF ENDFOR itots=intarr(8) FOR icol=0,7 DO itots[icol]=fix(total(ibin[indgen(nline)*8+icol])) print,itots,fix(total(itots)),format='(/8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,'', $ '' printf,lunhtml,'' FOR icol=0,7 DO printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'
Exit status, nchi2 < .7 ('+ $ string(n_elements(ix),format='(i6)')+' events).' printf,lunhtml,'
0001020304050607total
',ibin(iline*8+icol),'*',fix(total(ibin(iline*8:iline*8+7))),'
**********
',itots(icol),'*',fix(total(itots)),'
' ENDIF IF ihtml THEN printf,lunhtml,'
' ;*** Exit status, central events. nval=64 ix=where(central) ibin=histogram(exits[ix],min=0,nbins=nval) print,'Exit status, central events ('+string(n_elements(ix),format= $ '(i6)')+' events).' IF ihtml THEN BEGIN printf,lunhtml,'' printf,lunhtml,'', $ '' ENDIF nline=nval/8 FOR iline=0,nline-1 DO BEGIN print,ibin(iline*8:iline*8+7),total(ibin(iline*8:iline*8+7)), $ format='(8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,' FOR icol=0,7 DO printf,lunhtml,' printf,lunhtml,'' printf,lunhtml,' ENDIF ENDFOR itots=intarr(8) FOR icol=0,7 DO itots[icol]=fix(total(ibin[indgen(nline)*8+icol])) print,itots,fix(total(itots)),format='(/8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,'', $ '' printf,lunhtml,'' FOR icol=0,7 DO printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'
Exit status, central events ('+ $ string(n_elements(ix),format='(i6)')+' events).' printf,lunhtml,'
0001020304050607total
',ibin(iline*8+icol),'*',fix(total(ibin(iline*8:iline*8+7))),'
**********
',itots(icol),'*',fix(total(itots)),'
' ENDIF IF ihtml THEN printf,lunhtml,'		' ;*** Exit status, good, central events. nval=64 ix=where(good and central) ibin=histogram(exits[ix],min=0,nbins=nval) print,'Exit status, good central events ('+string(n_elements(ix),format= $ '(i6)')+' events).' IF ihtml THEN BEGIN printf,lunhtml,'' printf,lunhtml,'', $ '' ENDIF nline=nval/8 FOR iline=0,nline-1 DO BEGIN print,ibin(iline*8:iline*8+7),total(ibin(iline*8:iline*8+7)), $ format='(8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,' FOR icol=0,7 DO printf,lunhtml,' printf,lunhtml,'' printf,lunhtml,' ENDIF ENDFOR itots=intarr(8) FOR icol=0,7 DO itots[icol]=fix(total(ibin[indgen(nline)*8+icol])) print,itots,fix(total(itots)),format='(/8i6,i9)' IF ihtml THEN BEGIN printf,lunhtml,'', $ '' printf,lunhtml,'' FOR icol=0,7 DO printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'' printf,lunhtml,'
Exit status, good, central events ('+ $ string(n_elements(ix),format='(i6)')+' events).' printf,lunhtml,'
0001020304050607total
',ibin(iline*8+icol),'*',fix(total(ibin(iline*8:iline*8+7))),'
**********
',itots(icol),'*',fix(total(itots)),'
' ENDIF IF ihtml THEN BEGIN printf,lunhtml,'
' printf,lunhtml,'

Table made with ', $ 'plotbw1.pro' printf,lunhtml,'' close,lunhtml free_lun,lunhtml ENDIF stop END