Although https://math.stackexchange.com/ would be a better place to ask this question, let me answer it.

Since the characteristic function is basically the Fourier transform, let me explain it in terms of the Fourier transform $\varphi=\hat{f}$ and you can translate it to the language of the characteristic functions. I mean the Fourier transform defined by
$$
\hat{f}(\xi)=\int_{\mathbb{R}^n}f(x)e^{-2\pi ix\cdot\xi}\, dx.
$$
If $g(x)=\bar{f}(-x)$ is the complex conjugte of the "reflection" of $f$, then moving the conjugate under the sign of the integral and applying change of variables $x\mapsto -x$ yields that $\hat{g}(\xi)=\bar{\hat{f}}(\xi)$. Therefore, $\operatorname{re} \varphi=\operatorname{re} \hat{f}=\widehat{\frac{f+g}{2}}$. Now the fact that $\widehat{f*g}=\hat{f}\hat{g}$ implies that
$$
|\varphi|^2=|\hat{f}|^2=\hat{f}\bar{\hat{f}}=\hat{f}\hat{g}=\widehat{f*g}.
$$
In probabilistic terms, if independent random variables have distributions $f$ and $g$, then the sum of random variables has the distribution $f*g$.