def METHOD_NAME(self):
    self.zlens, self.zsource = 0.5, 1.5
    epl_kwargs = {
        "theta_E": 0.8,
        "center_x": 0.1,
        "center_y": 0.0,
        "e1": -0.2,
        "e2": 0.1,
        "gamma": 2.05,
    }
    shear_kwargs = {"gamma1": 0.09, "gamma2": -0.02}
    kwargs_macro = [epl_kwargs, shear_kwargs]
    self.x_image = np.array([0.65043538, -0.31109505, 0.78906059, -0.86222271])
    self.y_image = np.array([-0.89067493, 0.94851787, 0.52882605, -0.25403778])
    halo_list = ["SIS", "SIS", "SIS"]
    halo_z = [self.zlens - 0.1, self.zlens, self.zlens + 0.4]
    halo_kwargs = [
        {"theta_E": 0.05, "center_x": 0.3, "center_y": -0.9},
        {"theta_E": 0.01, "center_x": 1.3, "center_y": -0.5},
        {"theta_E": 0.02, "center_x": -0.4, "center_y": -0.4},
    ]
    self.kwargs_epl = kwargs_macro + halo_kwargs
    self.zlist_epl = [self.zlens, self.zlens] + halo_z
    self.lens_model_list_epl = ["EPL", "SHEAR"] + halo_list
    self.lensModel = LensModel(
        self.lens_model_list_epl,
        self.zlens,
        self.zsource,
        self.zlist_epl,
        multi_plane=True,
    )
    self.param_class = PowerLawFreeShear(self.kwargs_epl)