SIMBAD references

With high-quality data from programs like the Hubble Frontier Fields, cluster lensing has reached the point that models are dominated by systematic rather than statistical uncertainties. We introduce a Bayesian framework to quantify systematic effects by determining how different lens modelling choices affect the results. Our framework includes a new two-sample test for quantifying the difference between posterior probability distributions that are sampled by methods like Monte Carlo Markov chains. We use the framework to examine choices related to the selection and treatment of cluster member galaxies in two of the Frontier Field clusters: Abell 2744 and MACS J0416.1-2403. When selecting member galaxies, choices about depth and area affect the models; we find that model results are robust for an I-band magnitude limit of m_lim >= 22.5 mag and a radial cut of r_lim >= 90 arcsec (from the centre of the field), although the radial limit likely depends on the spatial extent of lensed images. Mass is typically assigned to galaxies using luminosity/mass scaling relations. We find that the slopes of the scaling relations can have significant effects on lens model parameters but only modest effects on lensing magnifications. Interestingly, scatter in the scaling relations affects the two fields differently. This analysis illustrates how our framework can be used to analyse lens modelling choices and guide future cluster lensing programs.