code: def multi_gaussian_lens_light(self, kwargs_lens_light, model_bool_list=None, e1=0, e2=0, n_comp=20, deltaPix=None, numPix=None):
          
          if 'center_x' in kwargs_lens_light[0]:
              center_x = kwargs_lens_light[0]['center_x']
              center_y = kwargs_lens_light[0]['center_y']
          else:
              center_x, center_y = 0, 0
          r_h = self.half_light_radius_lens(kwargs_lens_light, center_x=center_x, center_y=center_y,
                                            model_bool_list=model_bool_list, deltaPix=deltaPix, numPix=numPix)
          r_array = np.logspace(-3, 2, 200) * r_h * 2
          x_coords, y_coords = param_util.transform_e1e2(r_array, np.zeros_like(r_array), e1=-e1, e2=-e2)
          x_coords += center_x
          y_coords += center_y
          #r_array = np.logspace(-2, 1, 50) * r_h
          flux_r = self._lens_light_internal(x_coords, y_coords, kwargs_lens_light,
                                             model_bool_list=model_bool_list)
          amplitudes, sigmas, norm = mge.mge_1d(r_array, flux_r, N=n_comp)
          return amplitudes, sigmas, center_x, center_y
docs: multi-gaussian decomposition of the lens light profile (in 1-dimension)
  
          :param kwargs_lens_light:
          :param n_comp:
          :return: