off-motion-receptive-fields/receptiveFields/crossCorrShiftRF.m

function [rf, lagDiff] = crossCorrShiftRF(rf, refRF)
    [acor, lag] = xcorr(rf, refRF);
    [~, I] = max(abs(acor));
    lagDiff = lag(I);
    rf = circshift(rf, -lagDiff);
end