  Quasi-Monte Carlo methods are used for numerically integrating multivariate
functions. However, the error bounds for these methods typically rely on a
priori knowledge of some semi-norm of the integrand, not on the sampled
function values. In this article, we propose an error bound based on the
discrete Fourier coefficients of the integrand. If these Fourier coefficients
decay more quickly, the integrand has less fine scale structure, and the
accuracy is higher. We focus on rank-1 lattices because they are a commonly
used quasi-Monte Carlo design and because their algebraic structure facilitates
an error analysis based on a Fourier decomposition of the integrand. This leads
to a guaranteed adaptive cubature algorithm with computational cost $O(mb^m)$,
where $b$ is some fixed prime number and $b^m$ is the number of data points.
