In the renormalization procedure, is writing things like
$$\varphi=\sqrt{Z_{\varphi}}\ \varphi_R\ ,\ \ m_0^2=Z_m\ m_R^2\ ,\ \ g_0=Z_g \mu^{\epsilon}\ g_R$$ and $$Z_i=1+\sum_{\nu=1}^\infty C_i^{(\nu)}(m_R,\mu,\Lambda\text{ or }\epsilon)·g_R^\nu\ , \ \ \ \ \ i=\varphi, m, g$$ really more than just an arbitrary ansatz?
I have no idea what principle people follow, when people have a Lagrangian, say for QED and then write down Lagrangians in the to-be-renormalized-stage. There seems to be a motivation to make them look similar to the old Lagrangian before introducing that coupling constrant expansion - and why in $g$, not other variables like $m$? Hence they write things like $m_{old}=c·m_{new}$, which seems faily conservative, because it doesn't introduce new terms, beyond maybe counter terms that look structurally list the old ones. But as far as I can see, the theory really just starts with the Lagrangian, which contains the to be found $Z$-expressions. You don't use the Lagrangian before that, do you? At least not beyond tree graphs. Therefore I think you could just begin with a buch of terms, with object that have to be fitted by renormalization. The theory effectively seems just to start with the non-bare object.
From all the possible 'unphysical numbers' in the expansion for the (finite number of) $Z$-terms, why does only the 'scale' $\mu$ survive? Do all scheme leave one number open, and if yes, why? I don't get the what this object '$\mu$' is, at all.