Radial surface brightness profiles of disrupted satellite material in the final remnants of all the models. The plotted physical magnitudes are provided in simulation units. The letters identifying each panel are those used in Table 3. Circles: all the stars originally belonging to the satellite. Red dashed lines: satellite bulge stars. Blue dotted lines: satellite disc stars. Black solid straight lines: exponential fits to the radial density profiles considering all the satellite stars. The data points plotted in green have been excluded from the fits. In the panels corresponding to remnants that are more accurately fitted with two exponential discs, the two scale-lengths are shown in the frame. When just one exponential disc provides a better fit to the global satellite stellar profile than two nested ones, only one scale-length is provided. The root mean square of the global fit is also shown in each remnant. Orange light solid curved lines: total fit to the radial surface density profiles in those cases that are more accurately fitted by two exponential discs. Vertical dashed lines: extent of each IC characterized by each fitted exponential radial profile, as defined in the text. The radial extension of the ICs is indicated using a shaded background with a different colour for the different ICs identified in each model (green for the innermost ones and orange for the outer ones). Sub-panels: relative residuals of the fits as a function of the radial position. The models with small primary bulges obey a different scaling to the rest of models (see Sect. 2). Therefore, the simulation units in these two experiments differ from those in the others, as observed in Figs. 2–4.