;+ ; NAME: bfits_burst ; ; PURPOSE: read in BATSE bfits time history for bursts ; ; CALLING SEQUENCE: ; bfits_burst, filename, t, rate, erate=erate, edges=edges, head=head, burst=burst ; ; INPUTS: ; filename - bfits filename ; ; OPTIONAL INPUTS: ; ; OUTPUTS: ; t - 2xnsamples, times relative to burst trigger time in seconds ; rate - count rate (/s) in each channel, (nchans x nsamples) ; ; OPTIONAL OUTPUTS: ; erate - uncertainty on rate, NaN for bad samples ; edges - nominal energy channels (2 x nchans) in keV ; head - composite fits headers ; burst - burst trigger numbers ; PROCEDURE: ; ; CALLS: ; extensive calls to astronomy library binary table reader ; fits_spectra ; COMMON BLOCKS: ; ; RESTRICTIONS: ; ; MODIFICATION HISTORY: ; ras, 13-feb-95 ;- ; pro bfits_burst, filename, t, rate, edges=edges, erate=erate, head=head, burst=burst file = (findfile(filename, count=count))(0) if count eq 0 then message, 'Could not find '+filename fits_spectra, file, t, rate, erate=erate, live=live, edges=edges, head=head ;edges = edges(*,0:3) t= t - utime(tjd2ymd(fxpar(head,'basetime'))) burst=fxpar(head,'batse_tr') print,'This is the data for burst number ',burst help, t, rate, edges end