aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/halo.py
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2018-01-03 23:32:20 +0800
committerAaron LI <aly@aaronly.me>2018-01-03 23:32:20 +0800
commit68bc85f371d9ba258b8d9e970e1efaeda87d0122 (patch)
treef770e865a252566200afa440b4cc51d4cadc3f29 /fg21sim/extragalactic/clusters/halo.py
parent07a8673ef96964706ec0a5314c22cad16ad65b1e (diff)
downloadfg21sim-68bc85f371d9ba258b8d9e970e1efaeda87d0122.tar.bz2
clusters/halo: minor cleanups
Diffstat (limited to 'fg21sim/extragalactic/clusters/halo.py')
-rw-r--r--fg21sim/extragalactic/clusters/halo.py29
1 files changed, 15 insertions, 14 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 73924c6..5c7eb5b 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -279,24 +279,25 @@ class RadioHalo:
@property
@lru_cache()
- def kT_main(self):
+ def kT_obs(self):
"""
- The mean temperature of the main cluster ICM at ``z_merger``
- when the merger begins.
-
+ The ICM mean temperature of the cluster at ``z_obs``.
Unit: [keV]
"""
- return helper.kT_cluster(mass=self.M_main, z=self.z_merger,
+ return helper.kT_cluster(self.M_obs, z=self.z_obs,
configs=self.configs)
@property
@lru_cache()
- def kT_obs(self):
+ def kT_main(self):
"""
- The "current" cluster ICM mean temperature at ``z_obs``.
+ The mean temperature of the main cluster ICM at ``z_merger``
+ when the merger begins.
+
+ Unit: [keV]
"""
- return helper.kT_cluster(self.M_obs, z=self.z_obs,
- configs=self.configs) # [keV]
+ return helper.kT_cluster(mass=self.M_main, z=self.z_merger,
+ configs=self.configs)
@property
@lru_cache()
@@ -356,18 +357,18 @@ class RadioHalo:
ICM/volume, i.e., do not restricted inside the halo volume.
Qe(γ) = Ke * γ^(-s),
- int[ Qe(γ) γ me c^2 ]dγ * t_cluster = eta_e * e_th
+ int[ Qe(γ) γ me c^2 ]dγ * t_cluster = η_e * e_th
=>
- Ke = [(s-2) * eta_e * e_th * γ_min^(s-2) / (me * c^2 * t_cluster)]
+ Ke = [(s-2) * η_e * e_th * γ_min^(s-2) / (me * c^2 * t_cluster)]
References
----------
Ref.[cassano2005],Eqs.(31,32,33)
"""
s = self.injection_index
- e_thermal = helper.density_energy_thermal(self.M_obs, self.z_obs,
- configs=self.configs)
- term1 = (s-2) * self.eta_e * e_thermal # [erg cm^-3]
+ e_th = helper.density_energy_thermal(self.M_obs, self.z_obs,
+ configs=self.configs)
+ term1 = (s-2) * self.eta_e * e_th # [erg cm^-3]
term2 = self.gamma_min**(s-2)
term3 = AU.mec2 * self.age_obs # [erg Gyr]
Ke = term1 * term2 / term3 # [cm^-3 Gyr^-1]