In present and future magnetic confined fusion devices with metallic plasma-facing components (PFCs) such as JET-ILW and ITER, the calculation of the plasma composition must account for multiple impurities of a wide range of mass and charge, resolve their poloidal asymmetries and account for different central peakings for various elements. Single measurements of radiation and effective charge are not enough to characterize this complex system and a self-consistent analysis of data from multiple diagnostics is required. This contribution describes a method to calculate the plasma composition simultaneously accounting for contributions of up to two low-Z impurities, and two mid-/high-Z impurities. The analysis stems from methodologies explained in Sertoli et al. (Rev. Sci. Instrum., vol. 89 (11), 2018, 113501), expanded to include more impurities and to coherently analyse multiple diagnostics within the same framework. The example Ne-seeded JET-ILW hybrid discharge reported here shows that Be, Ne, Ni and W are necessary to simultaneously explain the observed soft X-ray emission, the W concentration measured by passive vacuum ultra-violet spectroscopy, the line-of-sight integrated measurement of the effective charge, the observed poloidal asymmetry of the soft X-ray (SXR) emission, the Ne density measured by charge-exchange-recombination spectroscopy and the line-of-sight integrals of the total radiation as measured by bolometry. This consistent picture of the elemental composition enables the calculation of the radial profiles of the effective charge, the dilution and total radiation. For the cases analysed up to now, these are often very different from the typical assumptions presently used when modelling JET-ILW discharges. This will affect, among others, the calculation of neutron rates, current density profile and heat transport. These considerations are of course valid for all present and future magnetic-controlled fusion devices which exhibit multi-material plasma-facing components, including ITER.