for allCells = 1:size(spikes.times,2)
    
    spikes0 = spikes.times{1,allCells}*sessionInfo.rates.wideband;    %all of the spikes of current cell in samples
 recordingLength = size(thermistorC,2)/10;
    FRmatrix = [];
    for i = 1:recordingLength
        edges = (i-1)*2e4*10+1:(i-1)*2e4*10+2e5;
        spikeok = histc(spikes0,edges);
        FRmatrix(:,i) = spikeok(1:length(edges));
    end
    FR = sum(FRmatrix)/10;
FRkell = FR(FR > 0);
tempFR = thermistorC(3,1:10:86950);
tempFR2 = tempFR(FR > 0);
figure
plot(tempFR2,FRkell,'o')
end


