$s^{-1} Mpc^{-1}$, $n_s = 0.96$,$ \sigma_8 = 0.76$, $\Omega_m = 0.268$ and $ \Omega_{\lambda} = 0.732$). The cosmological volume contains a 1 Mpc scale zoom-in region with three levels of refinement that lead to a minimum mass of the dark matter particle of 1.2 $\cdot$ $10^5 M_{\odot}$ and a minimum softening length of 124 pc ([@2014ApJ...784..161P]). In the zoom-in region an isolated MW-sized halo assembles, with mass consistent with the lower limit for the estimated halo mass of the Milky Way, which provides a good fit with current kinematical data on the Milky Way stellar halo ([@2013ApJ...773L..32R]). The companion hydrodynamical simulations, Eris and ErisLE, have shown to produce realistic spiral galaxies by $z=0$, which reproduce nearly all the main properties of the disk and bulge of the Milky Way ([@2011ApJ...742...76G]; [@2012ASPC..453..289M]; [@2013ApJ...772...36G]; [@2013ApJ...773...43B]).

Satellite galaxies falling into the main host after the last major merger, which takes place at $z \sim 3$, are considered for replacement with hi-resolution dwarf galaxy models including stars. We analize the entire available set of outputs of the original simulation with the use of the Amiga Halo Finder (AHF) ([@2004MNRAS.351..410G];[@2009ApJS..182..608K]) with a virial overdensity criterion that varies with redshift according to the cosmology ([@1998ApJ...495...80B]). In 