I am trying to understand the application of the operator product expansion to calculate the radially ordered product in the complex plain of $T_{zz}(z)\partial_w X^{\rho}(w)$ which should result in

$$ \langle R(T_{zz}(z)\partial_w X^{\rho}(w))\rangle = -l_s^2\frac{1}{(z-w)^2}\partial_w X^{\rho}(w) - l_s^2\frac{1}{(z-w)}\partial_z^2 X^{\rho}(z) + \cdots $$

but embarassingly I encounter a stumbling block right at the beginning. After inserting $T_{zz}(z) \doteq \, :\eta_{\mu\nu}\partial_z X^{\mu}\partial_zX^{\nu}:$ one has

$$ \langle R(T_{zz}(z)\partial_w X^{\rho}(w))\rangle = R(:\eta_{\mu\nu}\partial_z X^{\mu}(z)\partial_zX^{\nu}(z):\partial_w X^{\rho}(w)) $$

which can obviously be further expanded to

$$ ... = \eta_{\mu\nu}\langle \partial_z X^{\mu}(z)\partial_w X^{\rho}(w)\rangle \partial_z X^{\nu}(z) + \eta_{\mu\nu}\langle \partial_z X^{\nu}(z)\partial_w X^{\rho}(w)\rangle \partial_z X^{\mu}(z) $$

It is exactly this last step I dont understand. If this initial stumbling block is removed, I can understand the rest of the derivation, so can somebody help me remove it?

To generalize a bit, it seems I do not yet properly understand how such expressions involving normal and radial (time) ordered products are generally evaluated. So if somebody could give me a more general hint about this, I would probably be able to see how the last expression in my particular example is obtained.