Is there a way to characterize all non-constant entire functions such that $|f(z)| = f(|z|)$ ?
The monomials work, but I can't think of any other functions. So far, I have managed to show that $f$ is real on the real line, $f(z) \not= 0 $ if $z \not=0$, $|f'(z)| \geq f'(|z|)$, and $\forall \ \theta \in \mathbb R \ \exists \alpha \in \mathbb R $ s.t $ f(e^{i\theta}z)=e^{i\alpha}f(z) $. It would then suffice to show $\lim_{x \to \infty} f(x)= +\infty$, because $f$ would be a polynomial and have all roots at 0.
Edit: What about functions $f: U \to \mathbb C$ ?
An alternative proof, under the weaker assumption that $f$ is only defined and holomorphic in some disk $B_R(0)$:
If $f$ is not constant and $|f(z)| = f(|z|)$ then $f$ has no zeros except at the origin (because of the identity principle). It follows that $f(z) = z^n e^{g(z)}$ where $n$ is a positive integer and $g$ is holomorphic in $B_R(0)$. Then $$ |f(z)| = r^n e^{\operatorname{Re}(g(z))} $$ is constant on each circle $|z|=r$, so that $\operatorname{Re}(g(z))$ is constant on $|z|=r$. But the real part of $g$ is a harmonic function: If it is constant on a circle then it is constant inside the circle.
This implies that $g$ is constant, so that $f(z) = cz^n$ with some constant $c$ which is necessarily a positive real number.