Given input LFP data, output iv with nice theta cycles

Possible workflow:

• Filter LFP in theta band (FilterLFP())
• Obtain power of signal
• Threshold LFP to include only data with larger than some threshold theta power
• Get theta cycle times

May be helpful to restrict data by running speed also -- only include data when animal is running so that we exclude SWRs and other things that have theta power in the LFP

Like this:

To be able to average across trials of different lengths, may need to:

1. align to trial start, go length of shortest trial
2. align to trial end, go back length of shortest trial
3. bin time into some fixed number of bins

The metadata has trial start and end times.

Minimal example to reproduce error, found by @lemuelle :

cd('D:\data\R064\R064-2015-04-23');
cfg = []; cfg.fc = ExpKeys.goodTheta(1);
csc = LoadCSC(cfg);
csc_f = FilterLFP([], csc);
plot(csc); hold on; plot(csc_f, 'r'); % note csc_f is all NaNs