We present results from a comprehensive lensing analysis in HST data, of the complete Cluster Lensing And Supernova survey with Hubble (CLASH) cluster sample. We identify new multiple-images previously undiscovered, allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics -which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are ∼ 40% in the (dimensionless) mass density, κ, and ∼ 20% in the magnification, µ. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii (z s = 2) we find that all typically agree within 10% between the two models, and Einstein masses agree, typically, within ∼ 15%. At larger radii, the total projected, 2D integrated mass profiles of the two models, within r ∼ 2 , differ by ∼ 30%. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of d log(Σ)/d log(r) ∼ −0.64 ± 0.1, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities. All mass models and magnification maps are made publicly available for the community.
