[p, q] = get_products();

A = fft(sum(q'));
A = A .* (1 - (abs(A) >5000) .* (abs(A) < 20000));
plot(sum(q'), '-ob',real(ifft(A)),'-or');