Given $\epsilon>0$, define $$g(x) = \frac{1}{\sqrt{2\pi}}e^{-x^2/2}$$ and $$g_\epsilon(x) = \frac{1}{\epsilon}g\left(\frac{x}{\epsilon}\right) = \frac{1}{\epsilon\sqrt{2\pi}}e^{-x^2/(2\epsilon^2)}.$$
Then for $f\in L^1(\mathbb{R})\cap L^2(\mathbb{R})$, show that the integral
$$\int_{-\infty}^{\infty}\hat{g_\epsilon}(\xi)\vert \hat{f}(\xi)\vert^2 d\xi = 2\pi\int_{-\infty}^{\infty}\int_{-\infty}^{\infty} f(x)\overline{f(y+x)}g_\epsilon(y) dydx$$
using Fubini's Theorem and the Fourier Inversion Theorem.
Here is my attempt:
\begin{align} \int_{-\infty}^{\infty}\hat{g_\epsilon}(\xi)\vert \hat{f}(\xi)\vert^2 d\xi &= \int_{-\infty}^{\infty}\hat{g_\epsilon}(\xi) \hat{f}(\xi)\overline{\hat{f}(\xi)} d\xi\\ &= \int_{-\infty}^{\infty} \left(\int_{-\infty}^{\infty}g_\epsilon(y)e^{-i\xi y} dy\right) \left(\int_{-\infty}^{\infty}f(x)e^{-i\xi x}dx\right) \left(\int_{-\infty}^{\infty}f(t)e^{-i\xi t}dt\right)d\xi \end{align}
Beyond this, I haven't made any progress. Any hints on how to apply both theorems?
Begin with $$ |\hat{f}(\xi)|^2 = \hat{f}(\xi) \overline{\hat{f}(\xi)} = \left( \int f(x) \, e^{-i\xi x} \, dx \right) \left( \int \overline{f(z)} \, e^{i\xi z} \, dz \right) . $$ Then change to the variable $y=z-x$ in the second integral: $$ |\hat{f}(\xi)|^2 = \left( \int f(x) \, e^{-i\xi x} \, dx \right) \left( \int \overline{f(y+x)} \, e^{i\xi (y+x)} \, dy \right) = \iint f(x) \overline{f(y+x)} \, e^{i\xi y} \, dy \, dx . $$
Now multiply with $\hat{g}_\epsilon(\xi),$ integrate over $\xi$ and swap the integrals: $$ \int |\hat{f}(\xi)|^2 \, \hat{g}_\epsilon(\xi) \, d\xi = \int \left( \iint f(x) \overline{f(y+x)} \, e^{i\xi y} \, dy \, dx \right) \, \hat{g}_\epsilon(\xi) \, d\xi \\ = \iint f(x) \overline{f(y+x)} \left( \int \hat{g}_\epsilon(\xi) \, e^{i\xi y} \, d\xi \right) \, dy \, dx \\ = \iint f(x) \overline{f(y+x)} \, 2\pi \, g_\epsilon(y) \, dy \, dx \\ = 2\pi \iint f(x) \overline{f(y+x)} \, g_\epsilon(y) \, dy \, dx . $$