Change 'theta_E' keyword for Chameleon profile to reflect its actual definition

Effective Einstein radius is calculated fine for SPEMD, but not for Chameleon profile.

import numpy as np
import matplotlib.pyplot as plt

from lenstronomy.LensModel.lens_model import LensModel
from lenstronomy.LensModel.lens_model_extensions import LensModelExtensions


theta_es = np.linspace(0.4, 2, 30)
e1 = 0.0
e2 = 0.0

chams_theta_e = []
spemds_theta_e = []


for theta_e in theta_es:
    kwargs_chameleon = {'theta_E': theta_e, 'w_c': 1.13, 'w_t': 0.03, 'e1': e1, 'e2': e2,
                             'center_x': 0, 'center_y': 0}

    kwargs_spemd = {'theta_E': theta_e, 'gamma': 1.8, 'e1': e1, 'e2': e2, 'center_x': 0, 'center_y': 0}

    lens_model_ext_cham = LensModelExtensions(LensModel(lens_model_list=['CHAMELEON']))
    lens_model_ext_spemd = LensModelExtensions(LensModel(lens_model_list=['SPEMD']))

    spacing = 1000

    chams_theta_e.append(lens_model_ext_cham.effective_einstein_radius([kwargs_chameleon], spacing=spacing))
    spemds_theta_e.append(lens_model_ext_spemd.effective_einstein_radius([kwargs_spemd], spacing=spacing))

plt.plot(theta_es, spemds_theta_e, label='SPEMD')
plt.plot(theta_es, chams_theta_e, label='CHAMELEON')
plt.xlabel('True')
plt.ylabel('Calculated')
plt.legend()
plt.show()

[sibirrer]

Yes, the definition was not well documented. I changed the parameter to ‘theta_1’ now and added some better documentation. Thanks! I close the issue for now. Please reopen if things are not set. With my latest push, I also changed a bit the output format of a fitting sequence and added a plotting routine that can better deal with it. Let me know if you have questions.

[ajshajib]

That’s great! Thanks for making the changes.