The standard model has 12 massive leptons and 2 massive bosons other than the Higgs. My understanding of the Higgs mechanism is at about the level of this article, which goes as follows. Start with the massless Klein-Gordon equation for two fields, H and Z: $\nabla^2 H=0$ and $\nabla^2 Z=0$. Then modify the two equations to let H (the Higgs) have a v.e.v., and by adding an interaction term to the Z field's wave equation, $\nabla^2 Z=kH^2Z$. Since $H$ has a v.e.v., this looks exactly like the Klein-Gordon equation for a particle of mass $m$, with $m^2=kH^2$.
So extrapolating naively, I'd guess we have 14 coupling constants, $k_1$ through $k_{14}$, that have to be put in arbitrarily by hand and that determine the 14 masses that we observe in the low-energy limit for the standard model's 14 massive, fundamental particles.
Is this naive extrapolation accurate, or does $m_j$ actually depend not just on $k_j$ but also on the $k_i$ with $i\ne j$? If there is an interdependence, is it something that can be expressed linearly? Does chiral symmetry breaking (which I know nothing about) change the whole picture qualitatively? If there is an interdependence, do we get any prospect of a natural explanation for the features of the mass spectrum we see (e.g., why the masses cover such a huge range, from meV to TeV)?
Regardless of whether there is interdependence, is there any naturalness argument to explain why we see real-valued masses, as opposed to, say, imaginary masses, which would seem perfectly natural if there were no reason to prefer $k>0$?