diff options
| author | Aaron LI <aly@aaronly.me> | 2019-02-27 16:41:13 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2019-02-27 16:41:13 +0800 |
| commit | d8143b2ba2de88bd76bfb4375fea2a8b77d3b31f (patch) | |
| tree | 5927c616950a3a508aa686f36a74434f3798949b /fg21sim | |
| parent | 5498d4224f5c386675ba0b0e017c1cabff51a7a2 (diff) | |
| download | fg21sim-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')
| -rw-r--r-- | fg21sim/extragalactic/clusters/main.py | 3 |
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)) |