parent = '/home/jaime/Documents/Python Code/Data Sets/DataFiles/' files = None for root, dir, file in os.walk(parent): if root == parent: files = file for file in files: print(file) data = pickle.load(open(parent+file, 'rb'))['II_HI'] t = data.dataTime * 1E-6 i = data.data __, di, __ = diff.differentiate(t, i) i = i / np.max(i) di = di / np.max(di) di_filtered = filters.moving_average(di) fig = plt.figure() ax = fig.add_subplot(111) ax.plot(t, i, t, di_filtered) ax.set_title(file[0:6] + '_'+ file[-8:-5] + ' - Normalized and filtered dI/dt (5 pt MA)') p = df.pickerPlot(fig, ax) p.plot() fig1 = plt.figure() ax1 = fig1.add_subplot(111)
"Raw Data/081713/Scope18/C4AC00001.trc" f = lc.lecroy_data(dataFile) wantOffset = False trace = 1 traceUnit = 'V/m' traceUnitsFactor = 1.0 calFactor = 1 df32_tStart = 9.99993682e-001 df32_tStop = 1.00003088e+000 f.get_seg_time() f.get_segment(trace, calFactor) t, d_dt, wave = diff.differentiate(f.dataTime-df32_tStart, f.data) d_dt *= 1E9 fig, (ax1, ax2) = plt.subplots(2) ax1.plot(t*1E6, wave) ax2.plot(t*1E6, d_dt) ax1.set_ylabel('E (V/m)') ax2.set_ylabel('dE/dt (kV/m/$\mu$s)') ax2.set_xlabel('Time ($\mu$s)') plt.show()