aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/main.py
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2019-02-27 16:41:13 +0800
committerAaron LI <aly@aaronly.me>2019-02-27 16:41:13 +0800
commitd8143b2ba2de88bd76bfb4375fea2a8b77d3b31f (patch)
tree5927c616950a3a508aa686f36a74434f3798949b /fg21sim/extragalactic/clusters/main.py
parent5498d4224f5c386675ba0b0e017c1cabff51a7a2 (diff)
downloadfg21sim-d8143b2ba2de88bd76bfb4375fea2a8b77d3b31f.tar.bz2
clusters/halo: Update the spectral index calculation
Calculate the spectral index with a frequency separation of the width of 0.1 times the frequency value, instead of a fixed separation of 10 MHz.
Diffstat (limited to 'fg21sim/extragalactic/clusters/main.py')
-rw-r--r--fg21sim/extragalactic/clusters/main.py3
1 files changed, 1 insertions, 2 deletions
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py
index 994e0f3..e26b07f 100644
--- a/fg21sim/extragalactic/clusters/main.py
+++ b/fg21sim/extragalactic/clusters/main.py
@@ -375,8 +375,7 @@ class GalaxyClusters:
Tb_mean = haloem.calc_brightness_mean(
freq, flux=flux, pixelsize=self.sky.pixelsize)
- # Calculate spectral index with a narrow (10 MHz) band
- freq2 = freq + 10
+ freq2 = freq - freq*0.1
em2 = haloem.calc_emissivity(freq2)
index = -(np.log(em2)-np.log(em)) / (np.log(freq2)-np.log(freq))