  With mock strong gravitational lensing images, we investigate the performance
of the broken power-law (BPL) model proposed by \citet{2020ApJ...892...62D} on
the mass reconstruction of galaxy-scale lenses. An end-to-end test is carried
out, including the creation of mock strong lensing images, the subtraction of
lens light, and the reconstruction of lensed images, where the lenses are
selected from the galaxies in the Illustris-1 simulation. We notice that,
regardless of the adopted mass models (the BPL model or its special cases), the
Einstein radius can be robustly determined from imaging data alone, and the
median bias is typically less than $1\%$. Away from the Einstein radius, the
lens mass distribution tends to be harder to measure, especially at radii where
there are no lensed images detected. We find that, with rigid priors, the BPL
model can clearly outperform the single power-law models by achieving $<5\%$
median bias on the radial convergence profile within the Einstein radius. As
for the source light reconstructions, they are found to be sensitive to both
lens light contamination and lens mass models, where the BPL model with rigid
priors still performs best when there is no lens light contamination. We show
that, by correcting for the projection effect, the BPL model can estimate the
aperture and luminosity weighted line-of-sight velocity dispersions to an
accuracy of $\sim6\%$ scatter. These results highlight the great potential of
the BPL model in strong lensing related studies.
