We investigate fine global properties of nonnegative, integrable solutions to the Cauchy problem for the Fast Diffusion Equation with weights (WFDE) $u_t=|x|^\gamma\mathrm{div}\left(|x|^{-\beta}\nabla u^m\right)$ posed on $(0,+\infty)\times\mathbb{R}^d$, with $d\ge 3$, in the so-called good fast diffusion range $m_c<m<1$, within the range of parameters $\gamma, \beta$, optimal for the validity of the so-called Caffarelli-Kohn-Nirenberg inequalities. It is a natural question to ask in which sense such solutions behave like the Barenblatt $\mathfrak{B}$ (fundamental solution): for instance, asymptotic convergence, i.e. $\|u(t)-\mathfrak{B}(t)\|_{{\rm L}^p(\mathbb{R}^d)}\xrightarrow[]{t\to\infty}0$, is well known for all $1\le p\le \infty$, while only few partial results tackle a finer analysis of the tail behaviour. We characterize the maximal set of data $\mathcal{X}\subset{\rm L}^1_+(\mathbb{R}^d)$ that produces solutions which are pointwise trapped between two Barenblatt (Global Harnack Principle), and uniformly converge in relative error (UREC), i.e. ${\rm d}_\infty(u(t))=\|u(t)/\mathcal{B}(t)-1\|_{{\rm L}^\infty(\mathbb{R}^d)}\xrightarrow[]{t\to\infty}0$. Such characterization is in terms of an integral condition on $u(t=0)$. To the best of our knowledge, analogous issues for the linear heat equation $m=1$, do not possess such clear answers. Our characterization is also new for the classical, non-weighted, FDE. We are able to provide minimal rates of convergence to $\mathcal{B}$ in different norms. Such rates are almost optimal in the non weighted case, and become optimal for radial solutions. To complete the panorama, we show that solutions with data in ${\rm L}^1_+(\mathbb{R}^d)\setminus\mathcal{X}$, preserve the same "fat" spatial tail for all times, hence UREC fails.