From 940a4619782e87f95521338bb233b3f5d79c7fe8 Mon Sep 17 00:00:00 2001
From: Diana <diana.stoll@stundent.uni-tuebingen.de>
Date: Wed, 23 Oct 2024 13:31:58 +0200
Subject: [PATCH] power_spectrum binary statt rate

---
 code/useful_functions.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/code/useful_functions.py b/code/useful_functions.py
index 796ccb5..7153e09 100644
--- a/code/useful_functions.py
+++ b/code/useful_functions.py
@@ -273,7 +273,7 @@ def power_spectrum(stimulus):
     # computes firing rates
     rate = firing_rate(binary, dt = dt)
     # creates power spectrum
-    freq, power = welch(rate, fs = 1/dt, nperseg = 2**16, noverlap = 2**15)
+    freq, power = welch(binary, fs = 1/dt, nperseg = 2**16, noverlap = 2**15)
     return freq, power
 
 def prepare_harmonic(frequency, category, num_harmonics, color):