Workflow for the joint modeling using RXJ 1131 as an example. The input datasets consisted of photometric images and the spatial kinematics of the lens galaxy. The red contours in the middle right of the green panel represent the isocontours of both the light and mass density distributions of the lens galaxy, derived from the MGE method (see Sect. 3.2). The modeled I_light represents the light fitted from optical imaging, and I_light, IFU corresponds to the light near the spectral absorption lines in the IFU data. In the paper, I_light is equivalent to I_light, IFU. We employed an MCMC sampler to simultaneously sample the parameter space η_LD for the lensing and the dynamical modeling.