Overview flowchart for single-band modeling of a strongly lensed quasar system with our new automated code. The modeling consists of two main steps: the lens mass distribution modeling with source or image positions, and the reconstruction of the source surface brightness (SSB; which is the unlensed light distribution of the quasar host galaxy). Each step is divided into smaller modeling steps. The details are described in flowcharts for each step and in the corresponding subsection, as referenced in the flowchart.