def test_sis_density_deflection(self):
"""
tests whether the unit conversion between the lensing parameter 'sigma0' and the units in the density profile are ok
:return:
"""
from lenstronomy.LensModel.Profiles.sis import SIS as Model
lensModel = Model()
theta_E = 1.
rho0 = lensModel.theta2rho(theta_E)
kwargs_lens = {'theta_E': theta_E}
kwargs_density = {'rho0': rho0}
r = .5
mass_2d = lensModel.mass_2d(r, **kwargs_density)
alpha_mass = mass_2d / r
alpha_r, _ = lensModel.derivatives(r, 0, **kwargs_lens)
npt.assert_almost_equal(alpha_mass / np.pi, alpha_r, decimal=5)
lensModel.density_2d(1, 1, rho0=1)
class PJaffe(object):
"""
class for pseudo Jaffe lens light (2d projected light/mass distribution)
"""
param_names = ['amp', 'Ra', 'Rs', 'center_x', 'center_y']
lower_limit_default = {
'amp': 0,
'Ra': 0,
'Rs': 0,
'center_x': -100,
'center_y': -100
}
upper_limit_default = {
'amp': 100,
'Ra': 100,
'Rs': 100,
'center_x': 100,
'center_y': 100
}
def __init__(self):
from lenstronomy.LensModel.Profiles.p_jaffe import PJaffe as PJaffe_lens
self.lens = PJaffe_lens()
def function(self, x, y, amp, Ra, Rs, center_x=0, center_y=0):
"""
:param x:
:param y:
:param amp:
:param a:
:param s:
:param center_x:
:param center_y:
:return:
"""
rho0 = self.lens.sigma2rho(amp, Ra, Rs)
return self.lens.density_2d(x, y, rho0, Ra, Rs, center_x, center_y)
def light_3d(self, r, amp, Ra, Rs):
"""
:param y:
:param amp:
:param Rs:
:param center_x:
:param center_y:
:return:
"""
rho0 = self.lens.sigma2rho(amp, Ra, Rs)
return self.lens.density(r, rho0, Ra, Rs)
def assert_integrals(self, Model, kwargs):
lensModel = Model()
int_profile = ProfileIntegrals(lensModel)
r = 2.
density2d_num = int_profile.density_2d(r, kwargs)
density2d = lensModel.density_2d(r, 0, **kwargs)
npt.assert_almost_equal(density2d / density2d_num, 1., decimal=1)
mass_2d_num = int_profile.mass_enclosed_2d(r, kwargs)
mass_2d = lensModel.mass_2d(r, **kwargs)
npt.assert_almost_equal(mass_2d / mass_2d_num, 1, decimal=1)
mass_3d_num = int_profile.mass_enclosed_3d(r, kwargs)
mass_3d = lensModel.mass_3d(r, **kwargs)
npt.assert_almost_equal(mass_3d / mass_3d_num, 1, decimal=2)