One Normalization of the Fourier Transform
There are a several ways to normalize the Fourier Transform. The way I find nicest is to define
$$
\hat{f}(\xi)=\int_{-\infty}^\infty f(x)\,e^{-2\pi ix\cdot\xi}\,\mathrm{d}x\tag{1}
$$
Define
$$
g(x)=\sum_{k=-\infty}^\infty f(x+k)\tag{2}
$$
Then $g$ has period $1$. Plugging $x=0$ into $(2)$ gives
$$
g(0)=\sum_{k=-\infty}^\infty f(k)\tag{3}
$$
For $n\in\mathbb{Z}$, the coefficients of the Fourier Series are
$$
\begin{align}
\hat{g}(n)
&=\int_0^1g(x)\,e^{-2\pi inx}\,\mathrm{d}x\\
&=\sum_{k=-\infty}^\infty\int_0^1f(x+k)\,e^{-2\pi inx}\,\mathrm{d}x\\
&=\int_{-\infty}^\infty f(x)\,e^{-2\pi inx}\,\mathrm{d}x\\
&=\hat{f}(n)\tag{4}
\end{align}
$$
Inverting the Fourier Series, we get
$$
g(x)=\sum_{n=-\infty}^\infty\hat{g}(n)e^{2\pi inx}\tag{5}
$$
Plugging $x=0$ into $(5)$ and applying $(4)$ yields
$$
\begin{align}
g(0)
&=\sum_{n=-\infty}^\infty\hat{g}(n)\\
&=\sum_{n=-\infty}^\infty\hat{f}(n)\tag{6}
\end{align}
$$
$(3)$ and $(6)$ give the Poisson Summation Formula:
$$
\sum_{n=-\infty}^\infty f(n)=\sum_{n=-\infty}^\infty\hat{f}(n)\tag{7}
$$
Let $f(x)=g(ax)$. Then we have
$$
\begin{align}
\hat{f}(n)
&=\int_{-\infty}^\infty f(x) e^{-2\pi inx}\,\mathrm{d}x\\
&=\int_{-\infty}^\infty g(ax) e^{-2\pi inx}\,\mathrm{d}x\\
&=\frac1a\int_{-\infty}^\infty g(x) e^{-2\pi inx/a}\,\mathrm{d}x\\
&=\frac1a\hat{g}(n/a)\tag{8}
\end{align}
$$
Applying $(8)$ gives the scaled version of $(7)$:
$$
\sqrt{a}\sum_{n=-\infty}^\infty g(an)=\frac1{\sqrt{a}}\sum_{n=-\infty}^\infty\hat{g}(n/a)\tag{9}
$$
Equation $(9)$ is the equivalent of the equation in the question for the Fourier Transform defined in $(1)$. To get the equation in the question, we need to use a different normalization of the Fourier Transform.
Another Normalization of the Fourier Transform
If we define the Fourier Transform as
$$
\hat{f}(\xi)=\int_{-\infty}^\infty f(x)\,e^{-ix\cdot\xi}\,\mathrm{d}x\tag{1a}
$$
and define
$$
g(x)=\sum_{k=-\infty}^\infty f(x+2\pi k)\tag{2a}
$$
then $g$ has period $2\pi$ and plugging $x=0$ into $\mathrm{(2a)}$ gives
$$
g(0)=\sum_{k=-\infty}^\infty f(2\pi k)\tag{3a}
$$
For $n\in\mathbb{Z}$, the coefficients of the Fourier Series are
$$
\begin{align}
\hat{g}(n)
&=\int_0^{2\pi}g(x)\,e^{-inx}\,\mathrm{d}x\\
&=\sum_{k=-\infty}^\infty\int_0^{2\pi}f(x+2\pi k)\,e^{-inx}\,\mathrm{d}x\\
&=\int_{-\infty}^\infty f(x)\,e^{-inx}\,\mathrm{d}x\\
&=\hat{f}(n)\tag{4a}
\end{align}
$$
Inverting the Fourier Series, we get
$$
g(x)=\frac1{2\pi}\sum_{n=-\infty}^\infty\hat{g}(n)e^{inx}\tag{5a}
$$
Plugging $x=0$ into $\mathrm{(5a)}$ and applying $\mathrm{(4a)}$ yields
$$
\begin{align}
g(0)
&=\frac1{2\pi}\sum_{n=-\infty}^\infty\hat{g}(n)\\
&=\frac1{2\pi}\sum_{n=-\infty}^\infty\hat{f}(n)\tag{6a}
\end{align}
$$
$\mathrm{(3a)}$ and $\mathrm{(6a)}$ give the Poisson Summation Formula:
$$
\sum_{n=-\infty}^\infty f(2\pi n)=\frac1{2\pi}\sum_{n=-\infty}^\infty\hat{f}(n)\tag{7a}
$$
Let $f(x)=g(ax)$. Then we have
$$
\begin{align}
\hat{f}(n)
&=\int_{-\infty}^\infty f(x) e^{-inx}\,\mathrm{d}x\\
&=\int_{-\infty}^\infty g(ax) e^{-inx}\,\mathrm{d}x\\
&=\frac1a\int_{-\infty}^\infty g(x) e^{-inx/a}\,\mathrm{d}x\\
&=\frac1a\hat{g}(n/a)\tag{8a}
\end{align}
$$
Applying $\mathrm{(8a)}$ gives the scaled version of $\mathrm{(7a)}$:
$$
\sqrt{2\pi a}\sum_{n=-\infty}^\infty g(2\pi an)=\frac1{\sqrt{a}}\sum_{n=-\infty}^\infty\hat{g}(n/a)\tag{9a}
$$
$\mathrm{(9a)}$ is the equation in the question.