Current instruments and spectral analysis programs are now able to decompose the integrated spectrum of a stellar system into distributions of ages and metallicities. The reliability of these methods have rarely been tested on nearby systems with resolved stellar ages and metallicities. Here we derive the age-metallicity distribution of M54, the nucleus of the Sagittarius dwarf spheroidal galaxy, from its integrated MUSE spectrum. We find a dominant old (8-14 Gyr), metal-poor (-1.5 dex) and a young (1 Gyr), metal-rich (+0.25 dex) component - consistent with the complex stellar populations measured from individual stars in the same MUSE data set. There is excellent agreement between the (mass-weighted) average age and metallicity of the resolved and integrated analyses. Differences are only 3% in age and 0.2 dex metallicitiy. By co-adding individual stars to create M54's integrated spectrum, we show that the recovered age-metallicity distribution is insensitive to the magnitude limit of the stars or the contribution of blue horizontal branch stars - even when including additional blue wavelength coverage from the WAGGS survey. However, we find that the brightest stars can induce the spurious recovery of an old ($>8$ Gyr), metal-rich (+0.25 dex) stellar population, which is otherwise not expected from our understanding of chemical enrichment in M54. The overall derived stellar mass-to-light ratio of M54 is M/L$_{\mathrm{V}}=1.46$ with a scatter of 0.22 across the field-of-view, which we attribute to the stochastic contribution of a young, metal-rich component. These findings provide strong evidence that complex stellar population distributions can be reliably recovered from integrated spectra of extragalactic systems.

We present a new method of quantifying a galaxy's accretion history from its integrated spectrum alone. Using full spectral fitting and calibrated regularization techniques we show how we can accurately derive a galaxy's mass distribution in age-metallicity space and further separate this into stellar populations from different chemical enrichment histories. By exploiting the fact that accreted lower mass galaxies will exhibit an offset to lower metallicities at fixed age compared to the in-situ stellar population, we quantify the fraction of light that comes from past merger events, that are long since mixed in phase-space and otherwise indistinguishable. Empirical age-metallicity relations (AMRs) parameterized for different galaxy masses are used to identify the accreted stellar populations and link them back to the progenitor galaxy's stellar mass. This allows us to not only measure the host galaxy's total ex-situ mass fraction ($f_{acc}$), but also quantify the relative amount of accreted material deposited by satellite galaxies of different masses, i.e. the accreted satellite mass function in analogy to the subhalo mass function. Using mock spectra of simulated, present-day galaxies from the EAGLE suite we demonstrate that our method can recover the total accreted fraction to within $\approx 12 \%$, the stellar mass of the most massive accreted subhalo to within $\approx 26 \%$ and the slope of the accreted satellite mass function to within $\approx 16 \%$ of the true values from the EAGLE merger trees. Future application of this method to observations could potentially provide us accretion histories of hundreds of individual galaxies, for which deep integrated light spectroscopy is available.