Skip to content

Commit 8c9dfde

Browse files
committed
docs: correct sigma_LT/sigma_0 attribution in dPIE docstrings
The sqrt(3/2) relation between Lenstool's fiducial dispersion (sigma_LT) and the central dispersion (sigma_0) was attributed to Elíasdóttir et al. (2007) as if it were that paper's physical definition. Per the contributed derivation note (H. Ding 2026, "On the definitions of b0 and velocity dispersion in Lenstool / dPIE"), it is Lenstool's parameter convention: Lenstool computes b0 = 6*pia_c2*sigma^2 (an E07-style coefficient) but evaluates deflections with the Kassiola & Kovner (1993) / Limousin et al. (2005) amplitude, so the input sigma must be read through the L05 parameterization. Docstrings now state this, and give the equivalent central-dispersion form b0 = 4*648000*(sigma_0/c)^2*(D_LS/D_S). No numerical behaviour changes — documentation only; the sigma_LT parameterization is retained deliberately for Lenstool results-table parity. Co-Authored-By: Claude Fable 5 <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_012Eacj9naBWSKrevsqnVJWR
1 parent 213278c commit 8c9dfde

1 file changed

Lines changed: 18 additions & 8 deletions

File tree

autogalaxy/profiles/mass/total/dual_pseudo_isothermal_mass.py

Lines changed: 18 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -571,9 +571,10 @@ def from_lenstool(
571571
572572
- ``sigma`` is Lenstool's **fiducial** velocity dispersion ``v_disp`` (sigma_LT),
573573
*not* the central velocity dispersion sigma_0 of the dPIE profile. They differ
574-
by sigma_0 = sqrt(3/2) * sigma_LT (Eliasdottir et al. 2007, App. A; Bergamini
575-
et al. 2019, Eq. 5) — quoting a measured central/aperture dispersion here
576-
overestimates the mass by 50%. The lens strength is
574+
by sigma_0 = sqrt(3/2) * sigma_LT — a Lenstool parameter convention (E07-style
575+
b0 coefficient paired with the K93/L05 deflection amplitude; see the dPIEMass
576+
class docstring), quoted by Bergamini et al. 2019, Eq. 5 — so a measured
577+
central/aperture dispersion entered here overestimates the mass by 50%. The lens strength is
577578
b0 = 6 * 648000 * (sigma_LT / c)^2 * (D_LS / D_S) arcsec, where Lenstool's
578579
stored ``b0 = 6 * pia_c2 * sigma^2`` (``set_potfile.c``) carries no distance
579580
ratio — Lenstool applies D_LS / D_S separately at deflection time
@@ -1095,18 +1096,27 @@ class dPIEMass(dPIEMassB0):
10951096
must be handled when ingesting real-data catalogues).
10961097
- ``sigma`` — Lenstool's **fiducial** velocity dispersion ``v_disp``
10971098
(sigma_LT, km/s), *not* the central dispersion:
1098-
sigma_0 = sqrt(3/2) * sigma_LT (Elíasdóttir et al. 2007, App. A;
1099-
Bergamini et al. 2019, Eq. 5). Quoting a measured stellar dispersion here
1100-
overestimates the mass by 50%.
1099+
sigma_0 = sqrt(3/2) * sigma_LT (the relation quoted by Bergamini et al.
1100+
2019, Eq. 5). The sqrt(3/2) is Lenstool's parameter convention, not a
1101+
physical definition from Elíasdóttir et al. (2007): it arises because
1102+
Lenstool computes ``b0 = 6 * pia_c2 * sigma^2`` (an E07-style
1103+
coefficient, ``set_potfile.c``) but evaluates deflections with the
1104+
Kassiola & Kovner (1993) / Limousin et al. (2005) amplitude, so the
1105+
input sigma must be read through the L05 parameterization (see the
1106+
contributed derivation note "On the definitions of b0 and velocity
1107+
dispersion in Lenstool / dPIE", H. Ding 2026). Quoting a measured
1108+
stellar dispersion here overestimates the mass by 50%.
11011109
- ``r_core`` / ``r_cut`` — Lenstool ``core_radius`` / ``cut_radius`` in
11021110
arcseconds (the internal ``ra`` / ``rs``). For ``.par`` files using the kpc
11031111
variants, pre-convert with
11041112
``r_core = r_core_kpc / cosmology.kpc_per_arcsec_from(redshift=redshift_object)``.
11051113
11061114
The lens strength is fully normalized internally:
11071115
b0 = 6 * 648000 * (sigma_LT / c)^2 * (D_LS / D_S) arcsec — equivalently
1108-
E_0 = 6 pi (D_LS / D_S) (sigma_LT / c)^2 in radians (Elíasdóttir et al. 2007,
1109-
Eq. A24) with the E_0-to-b0 prefactor folded in. Lenstool stores its ``b0``
1116+
E_0 = 6 pi (D_LS / D_S) (sigma_LT / c)^2 in radians (the coefficient of
1117+
Elíasdóttir et al. 2007, Eq. A24) with the E_0-to-b0 prefactor folded in;
1118+
written in the central dispersion this is
1119+
b0 = 4 * 648000 * (sigma_0 / c)^2 * (D_LS / D_S). Lenstool stores its ``b0``
11101120
without the distance ratio and applies D_LS / D_S at deflection time
11111121
(``e_grad.c``); the two conventions are verified equivalent against the
11121122
Lenstool C source and reference deflections

0 commit comments

Comments
 (0)