From 9952a13bf937b4de4813d601c6fbf13d10146be3 Mon Sep 17 00:00:00 2001
From: Ramona <efish@verliernix.com>
Date: Fri, 16 Nov 2018 16:20:01 +0100
Subject: [PATCH] gauss function

---
 code/spikes_analysis.py | 22 +++++++++++++++++++---
 code/utility.py         |  6 ++++++
 2 files changed, 25 insertions(+), 3 deletions(-)

diff --git a/code/spikes_analysis.py b/code/spikes_analysis.py
index ffe8086..d0b1867 100644
--- a/code/spikes_analysis.py
+++ b/code/spikes_analysis.py
@@ -18,8 +18,24 @@ for k in spikes.keys():
     else:
         df_map[df] = [k]
 
-embed()
-exit()
+spikes_mat = {}
+for deltaf in df_map.keys():
+    for rep in df_map[deltaf]:
+        for phase in spikes[rep]:
+            #print(phase)
+            spikes_one_chirp = spikes[rep][phase]
+            if deltaf == '-50Hz' and phase == (9, 0.54):
+                spikes_mat[deltaf, rep, phase] = spikes_one_chirp
+
+#embed()
+#exit()
+
 plot_spikes = spikes[(0, '-50Hz', '20%', '100Hz')][(0, 0.789)]
-plt.scatter(plot_spikes, np.ones(len(plot_spikes)))
+#hist_data = plt.hist(plot_spikes, bins=np.arange(-200, 400, 20))
+fig, ax = plt.subplots()
+ax.scatter(plot_spikes, np.ones(len(plot_spikes))*10, marker='|', color='k')
+#ax.plot(hist_data[1][:-1], hist_data[0])
 plt.show()
+
+#embed()
+#exit()
diff --git a/code/utility.py b/code/utility.py
index 1ad9601..68102d3 100644
--- a/code/utility.py
+++ b/code/utility.py
@@ -17,3 +17,9 @@ def vector_strength(spike_times, eod_durations):
 		phase_times[i] = (spike_times[i] / eod_durations[i]) * 2 * np.pi
 	vs = np.sqrt((1/n*sum(np.cos(phase_times)))**2 + (1/n*sum(np.sin(phase_times)))**2)
 	return vs
+
+
+def gaussian(x, mu, sig):
+	y = np.exp(-np.power(x - mu, 2.) / (2 * np.power(sig, 2.)))
+	return y
+