In section 9.4 of Peskin & Schroeder's textbook on quantum field theory, when applying the Faddeev Popov procedure to quantize an Abelian gauge theory, they obtain the following functional integral: $$\int \mathcal{D}A e^{iS[A]} = \det\Big(\frac{1}{e}\partial^2\Big)\Big( \int \mathcal{D}\alpha\Big) \int\mathcal{D}Ae^{iS[A]}\delta\big(\partial^\mu A_\mu - \omega(x)\big).\tag{p.296}$$ Here $A$ is the gauge field. They then say:
This equality holds for any $\omega(x)$, so it will also hold if we replace the right-hand side with any properly normalized linear combination involving different functions $\omega(x)$. For our final trick, we will integrate over all $\omega(x)$, with a Gaussian weighting function centered on $\omega = 0$. The above expression is thus equal to $$N(\xi)\int D\omega \exp\Big[-i\int d^4x \frac{\omega^2}{2\xi} \Big] \det\Big(\frac{1}{e}\partial^2\Big) \Big(\int \mathcal{D}\alpha\Big) \int \mathcal{D} A e^{iS[A]} \delta(\partial^\mu A_\mu - \omega(x)) \\ = N(\xi)\det \Big(\frac{1}{e}\partial^2\Big)\Big(\int \mathcal{D}\alpha \Big) \int \mathcal{D}A e^{iS[A]} \exp\Big[-i\int d^4x \frac{1}{2\xi} (\partial^\mu A_\mu)^2\Big], \tag{9.56}$$ where $N(\xi)$ is an unimportant normalization constant and we have used the delta function to perform the integral over $\omega$. We can choose $\xi$ to be any finite constant. Effectively, we have added a new term $-(\partial^\mu A_\mu)^2/2\xi$ to the Lagrangian.
I have understood everything about this derivation except for the Gaussian weighting function and the role of the parameter $\xi$.
Regarding the Gaussian weighting function: How is this a Gaussian weight? A Gaussian integral is usually defined with the spatial integral being outside the exponential, but here the integral is inside the exponential. Also, if it is a Gaussian, why is the argument of the exponential negative?
Regarding $\xi$: What is the purpose of it and why do we introduce it as part of the Gaussian weight? Is it supposed to be interpreted as a variance? Second, one often chooses a value for $\xi$ and this is known as fixing a gauge. For example, $\xi = 0$ is the Landau gauge and $\xi = 1$ is the Feynman gauge. Does this gauge fixing have any relation to the gauge fixing of $A$ (e.g. the Lorenz gauge $\partial_\mu A^\mu = 0$)?