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()
