diff options
| author | Aaron LI <aly@aaronly.me> | 2017-07-30 10:19:18 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-07-30 10:19:18 +0800 |
| commit | 1abb6bfddb6858fd0140accaf7eece61b3764ae4 (patch) | |
| tree | 59a0df3353380f61ea5ca8c8630232199d1eecb9 | |
| parent | 4632691c206502cbbf5396515d181c6e5c34f026 (diff) | |
| download | fg21sim-1abb6bfddb6858fd0140accaf7eece61b3764ae4.tar.bz2 | |
clusters/halo.py: Add property "angular_radius" (unit: [arcsec])
Signed-off-by: Aaron LI <aly@aaronly.me>
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 10 | ||||
| -rw-r--r-- | fg21sim/utils/units.py | 2 |
2 files changed, 12 insertions, 0 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index b65cd08..adc2bc5 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -159,6 +159,16 @@ class RadioHalo: return r_halo @property + def angular_radius(self): + """ + The angular radius of the radio halo. + Unit: [arcsec] + """ + DA = COSMO.DA(self.z) * 1e3 # [Mpc] -> [kpc] + theta = self.radius / DA # [rad] + return theta * AUC.rad2arcsec + + @property def volume(self): """ The halo volume, calculated from the above radius. diff --git a/fg21sim/utils/units.py b/fg21sim/utils/units.py index e65f831..a6d7034 100644 --- a/fg21sim/utils/units.py +++ b/fg21sim/utils/units.py @@ -59,6 +59,8 @@ class UnitConversions: # Angle rad2deg = au.rad.to(au.deg) deg2rad = au.deg.to(au.rad) + rad2arcsec = au.rad.to(au.arcsec) + arcsec2rad = au.arcsec.to(au.rad) arcmin2deg = au.arcmin.to(au.deg) deg2arcmin = au.deg.to(au.arcmin) arcsec2deg = au.arcsec.to(au.deg) |