Projected total mass density distribution of M0416 from the best-fit lens model presented in this work. Left: Projected total mass density distribution of M0416. Contour levels correspond to values of [0.70,1.85,3.00] × 10⁹ M_⊙ kpc⁻². Right: Difference between the cluster total mass density distributions inferred from the present reference lens model and that by B21. This map is obtained by considering 500 realizations of the lens models randomly extracting samples of free parameter values from the MCMC chains. The procedure adopted to generate this map is detailed in Sect. 4. Contour levels correspond to values of [−5,5] × 10⁷ M_⊙ kpc⁻² (i.e., the limits of the color bar). The observed positions of the 237 multiple images used as model constraints are shown as black crosses.