← Back to branch summary

~plane1/+junk/dissertation

« back to all changes in this revision

Viewing changes to src/Analysis/HPRF/plot_breakdown_decays.py

  • Committer: plane1
  • Date: 2016-04-11 17:09:19 UTC
  • Revision ID: plane1@hawk.iit.edu-20160411170919-6k08ir33ru4lphsv
Added additional Python analysis scripts and NumPy data files

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/Analysis/HPRF/plot_breakdown_decays.py
 
1
# -*- coding: utf-8 -*-
 
2
"""
 
3
Created on Fri Feb 19 16:28:19 2016
 
4
 
 
5
@author: plane
 
6
"""
 
7
 
 
8
import math
 
9
import subprocess as proc
 
10
import sys
 
11
import platform
 
12
import numpy as np  # NumPy (multidimensional arrays, linear algebra, ...)
 
13
import matplotlib as mpl
 
14
import matplotlib.cm as cm
 
15
import matplotlib.collections as mcollections
 
16
import matplotlib.patches as mpatches
 
17
import matplotlib.pyplot as plt  # Matplotlib's pyplot: MATLAB-like syntax
 
18
import scipy.signal as sig
 
19
import scipy.special as special
 
20
import bisect as bi
 
21
import time
 
22
import pgl.comsol as comsol
 
23
import pgl.curve as curve
 
24
import pgl.labview as labview
 
25
import pgl.plot as plot
 
26
import pgl.signal as psig
 
27
import pgl.tektronix as tektronix
 
28
import pgl.mc as mc
 
29
import pgl.progress as prog
 
30
import pgl.cluster as clust
 
31
 
 
32
def reset_plot_params():
 
33
    mpl.rcParams['ytick.labelsize'] = 22
 
34
    mpl.rcParams['xtick.labelsize'] = 22
 
35
    mpl.rcParams['axes.labelsize'] = 26
 
36
    mpl.rcParams['font.size'] = 26
 
37
    mpl.rcParams['mathtext.default'] = 'regular'
 
38
    mpl.rcParams['figure.subplot.left'] = 0.02
 
39
    mpl.rcParams['figure.subplot.right'] = 0.98
 
40
    mpl.rcParams['figure.subplot.top'] = 0.9
 
41
    mpl.rcParams['figure.subplot.bottom'] = 0.1
 
42
    mpl.rcParams['figure.subplot.wspace'] = 0.2
 
43
    mpl.rcParams['figure.subplot.hspace'] = 0.2
 
44
reset_plot_params()
 
45
 
 
46
plt.figure(figsize=(12,9))
 
47
 
 
48
lv_directory = "C:\\Users\\plane\\Desktop\\Data\\HPRF\\HPRF 20120605\sparks_29p5MV-M_15Hz\\"
 
49
lv_file = "spark_66_68.npz"
 
50
(times, raw_signals) = labview.load_data(lv_directory, lv_file, t0=221.324e-3, t1=222.524e-3,
 
51
                                         channels=[0,1,2,6])
 
52
dt = 2.0e-6
 
53
flip_time = 310.e-6
 
54
flip_offset = int(round(flip_time / dt))
 
55
raw_signals[2,flip_offset:] = -raw_signals[2,flip_offset:]
 
56
plt.subplot(211)
 
57
plt.ylim((-1.2,1.2))
 
58
plot.plot_signals(times*1e6, raw_signals, norm=True)  # 22us wavefront delay
 
59
 
 
60
data_dir = "C:\\Users\\plane\\Dropbox\\Research\\MTA\\Analysis\\HPRF\\"
 
61
if platform.system() == 'Linux':
 
62
    data_dir = "/home/lane/Dropbox/Research/MTA/Analysis/HPRF/"
 
63
"""
 
64
comsol_file = "hpc_wall_shock.npy"
 
65
(times, breakdown_signals) = comsol.load_data(data_dir, comsol_file, dt=2.0e-6)
 
66
dt = times[1]
 
67
plt.subplot(212)
 
68
plt.ylim((-1.2,1.2))
 
69
plt.ylabel('Normalized Amplitude')
 
70
plot.plot_signals(times*1e6, breakdown_signals[0], norm=True)  # 22us wavefront delay
 
71
"""
 
72
 
 
73
plt.tight_layout()
 
 
b'\\ No newline at end of file'