Solution (source code)

= Solution

Condition on $X$. The normal <characteristic function> gives $\mathbb E(e^{iuXY}\mid X)=e^{-u^2X^2/2}$, whose <Gaussian integral> is $(1+u^2)^{-1/2}$. A second <Gaussian integral> gives $\mathbb E(e^{-iuZ^2})=(1+2iu)^{-1/2}$, taking the square-root branch continuous from one at $u=0$. Independence of $(X,Y)$ and $Z$ therefore gives
$$
\boxed{\phi_\eta(u)=\frac1{\sqrt{1+u^2}\sqrt{1+2iu}}.}
$$