Notes. The mass components are grouped into (HALOS), shears and foreground galaxies (S), and cluster galaxies following the scaling relations (G). Parameters with no errors are fixed. All cluster-scale halos are modeled as dPIEs. Halos making up the hot gas component are taken from B18. Sky x, y coordinates are the offsets in arcsec from the reference BCG positions (MACS 1206: RA = 12^h06^m1215, Dec = , MACS 0416: RA = 04^h16^m0915, Dec = and AS1063: RA = 22^h48^m4397, Dec = ). The ellipticity e is defined as , where a and b are the semi-major and semi-minor axis; when a shear term is present the ellipticity refers to γ_shear (see C17b). The position angle θ is computed from west to north; σ_LT is the LensTool fiducial velocity dispersion, r_core and r_cut are the core and truncation radius, respectively. For the galaxy component (G), the number of spherical sub-halos included in the model and the normalization F160W magnitude are quoted in the third column. For each optimized parameter, we quote the median, and the 16th and 84th percentiles as errors. The last column shows the total root-mean-square separation between the model-predicted and observed positions of the multiple images.