from burst import Burst import numpy as np import matplotlib.pyplot as plt import sys import numpy.ma as ma # 00143: Norris Fit out_file = '/home/thomas/git/MSDataSci/Thesis/out_files/00143_b536_3_0.out' burst_file = './sample_data/cat64ms.00143' burst = Burst(burst_file) burst.parse_batse_file() data = burst.burst_data burst_meta = burst.header_data four_channel_data = np.sum(data, axis=1) trig = burst_meta['trig#'] nlasc = burst_meta['nlasc'] npts = burst_meta['npts'] bg_height = 1047.2220 bg_slope = 0.025935641 start = 44.965676 ampl = 2823.6814 tau1 = 8.1735936 tau2 = 2.2154609 taupk = 48.006471 time = (np.linspace((-nlasc), (npts - nlasc), num=npts) * 0.064) - 2