From 5aa0369e69a3546b3099dc14ceb367fd7b903cbd Mon Sep 17 00:00:00 2001
From: nmtimme <38138179+nmtimme@users.noreply.github.com>
Date: Thu, 7 Dec 2023 10:28:28 -0500
Subject: [PATCH] Remove spikes near edge of recording to prevent index error

---
 analysis/getWaveForms.m | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/analysis/getWaveForms.m b/analysis/getWaveForms.m
index 2ecbdeb..2b4d3e0 100644
--- a/analysis/getWaveForms.m
+++ b/analysis/getWaveForms.m
@@ -47,6 +47,8 @@
 for curUnitInd=1:numUnits
     curUnitID = unitIDs(curUnitInd);
     curSpikeTimes = gwfparams.spikeTimes(gwfparams.spikeClusters==curUnitID);
+    curSpikeTimes((curSpikeTimes + gwfparams.wfWin(1)) < 1) = [];
+    curSpikeTimes((curSpikeTimes + gwfparams.wfWin(end)) > size(mmf.Data.x,2)) = [];
     curUnitnSpikes = size(curSpikeTimes,1);
     spikeTimesRP = curSpikeTimes(randperm(curUnitnSpikes));
     spikeTimeKeeps(curUnitInd,1:min([gwfparams.nWf curUnitnSpikes])) = sort(spikeTimesRP(1:min([gwfparams.nWf curUnitnSpikes])));