function [time, rates] = getFiringRates(all_spikes, min_t, max_t, dt)


time = min_t:dt:max_t-dt;
rates = zeros(length(time), length(all_spikes));
for i = 1:length(all_spikes)
    spike_times = all_spikes{i}/1000 - min_t;
    spike_times(spike_times == 0.0) = dt;
    [t, rates(:, i)] = convolutionRate(spike_times, 0.01, dt, max_t-min_t);
end