From 67353f5c0b4cb2de4be1f2d9f7530996a68d8a3c Mon Sep 17 00:00:00 2001 From: efish Date: Thu, 15 Nov 2018 14:33:11 +0100 Subject: [PATCH] bug fixes --- code/read_baseline_data.py | 16 ++++++++++++++-- code/read_chirp_data.py | 6 +++--- 2 files changed, 17 insertions(+), 5 deletions(-) diff --git a/code/read_baseline_data.py b/code/read_baseline_data.py index 6b3cd29..9e951a9 100644 --- a/code/read_baseline_data.py +++ b/code/read_baseline_data.py @@ -7,7 +7,13 @@ def read_baseline_eod(dataset): base = dataset.split(os.path.sep)[-1] + ".nix" nix_file = nix.File.open(os.path.join(dataset, base), nix.FileMode.ReadOnly) b = nix_file.blocks[0] - t = b.tags["BaselineActivity_1"] + if 'BaselineActivity_1' in b.tags: + t = b.tags["BaselineActivity_1"] + elif "BaselineActivity_2" in b.tags: + t = b.tags["BaselineActivity_2"] + else: + f.close() + return [],[] eod_da = b.data_arrays["LocalEOD-1"] eod = t.retrieve_data("LocalEOD-1")[:] time = np.asarray(eod_da.dimensions[0].axis(len(eod))) @@ -19,7 +25,13 @@ def read_baseline_spikes(dataset): base = dataset.split(os.path.sep)[-1] + ".nix" nix_file = nix.File.open(os.path.join(dataset, base), nix.FileMode.ReadOnly) b = nix_file.blocks[0] - t = b.tags["BaselineActivity_1"] + if 'BaselineActivity_1' in b.tags: + t = b.tags["BaselineActivity_1"] + elif "BaselineActivity_2" in b.tags: + t = b.tags["BaselineActivity_2"] + else: + f.close() + return [],[] spikes_da = b.data_arrays["Spikes-1"] spike_times = spikes_da[:spikes_da.shape[0]-5000] baseline_spikes = spike_times[(spike_times > t.position[0]) & (spike_times < (t.position[0] + t.extent[0]))] diff --git a/code/read_chirp_data.py b/code/read_chirp_data.py index 4283c48..4b84139 100644 --- a/code/read_chirp_data.py +++ b/code/read_chirp_data.py @@ -2,7 +2,7 @@ import numpy as np import os -def load_chirp_spikes(dataset): +def read_chirp_spikes(dataset): spikes_file = os.path.join(dataset, "chirpspikess1.dat") if not os.path.exists(spikes_file): print("found no chirps!") @@ -32,7 +32,7 @@ def load_chirp_spikes(dataset): return spikes -def load_chirp_eod(dataset): +def read_chirp_eod(dataset): eod_file = os.path.join(dataset, "chirpeodampls.dat") if not os.path.exists(eod_file): print("found no chirpeodampls.dat!") @@ -60,7 +60,7 @@ def load_chirp_eod(dataset): return chirp_eod -def load_chirp_times(dataset): +def read_chirp_times(dataset): chirp_times_file = os.path.join(dataset, "chirpss.dat") if not os.path.exists(chirp_times_file): print("found no chirpss.dat!")