function [counts, cbins] = counthist(spikes, tmin, tmax, T, cmax)
tbins = tmin+T/2:T:tmax;
cbins = 0.5:cmax;
counts = zeros(1, length(cbins));
for k = 1:length(spikes)
    times = spikes{k};
    n = hist(times((times>=tmin)&(times<=tmax)), tbins);
    counts = counts + hist(n, cbins);
end
counts = counts / sum(counts);
end