diff options
| author | Aaron LI <aly@aaronly.me> | 2018-01-01 14:51:55 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2018-01-01 14:51:55 +0800 |
| commit | d82487346e87aedde6b9a476d43c7b1e1a821319 (patch) | |
| tree | 588d8c85d1026e5da95f74a37b66fb3bc304eb33 /fg21sim | |
| parent | ee80b0578613410dc4b66487ec4dacd474b025d5 (diff) | |
| download | fg21sim-d82487346e87aedde6b9a476d43c7b1e1a821319.tar.bz2 | |
Minor changes
Diffstat (limited to 'fg21sim')
| -rw-r--r-- | fg21sim/configs/20-extragalactic.conf.spec | 2 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/helper.py | 44 |
2 files changed, 23 insertions, 23 deletions
diff --git a/fg21sim/configs/20-extragalactic.conf.spec b/fg21sim/configs/20-extragalactic.conf.spec index 1992f3f..a75907c 100644 --- a/fg21sim/configs/20-extragalactic.conf.spec +++ b/fg21sim/configs/20-extragalactic.conf.spec @@ -100,7 +100,7 @@ # The minimum mass for clusters when to determine the galaxy clusters # total counts and their distributions. # Unit: [Msun] - mass_min = float(default=1e14, min=1e12) + mass_min = float(default=1e14, min=1e13) # Boost the number of expected cluster number within the sky coverage # by the specified times. diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py index 29134f7..88c566a 100644 --- a/fg21sim/extragalactic/clusters/helper.py +++ b/fg21sim/extragalactic/clusters/helper.py @@ -231,6 +231,28 @@ def density_energy_thermal(mass, z=0.0, configs=CONFIGS): return e_th +def density_energy_electron(spectrum, gamma): + """ + Calculate the energy density of relativistic electrons. + + Parameters + ---------- + spectrum : 1D float `~numpy.ndarray` + The number density of the electrons w.r.t. Lorentz factors + Unit: [cm^-3] + gamma : 1D float `~numpy.ndarray` + The Lorentz factors of electrons + + Returns + ------- + e_re : float + The energy density of the relativistic electrons. + Unit: [erg cm^-3] + """ + e_re = integrate.trapz(spectrum*gamma*AU.mec2, gamma) + return e_re + + def magnetic_field(mass, z=0.0, configs=CONFIGS): """ Calculate the mean magnetic field strength within the ICM, which is @@ -255,28 +277,6 @@ def magnetic_field(mass, z=0.0, configs=CONFIGS): return B -def density_energy_electron(spectrum, gamma): - """ - Calculate the energy density of relativistic electrons. - - Parameters - ---------- - spectrum : 1D float `~numpy.ndarray` - The number density of the electrons w.r.t. Lorentz factors - Unit: [cm^-3] - gamma : 1D float `~numpy.ndarray` - The Lorentz factors of electrons - - Returns - ------- - e_re : float - The energy density of the relativistic electrons. - Unit: [erg cm^-3] - """ - e_re = integrate.trapz(spectrum*gamma*AU.mec2, gamma) - return e_re - - def speed_sound(kT): """ The adiabatic sound speed in cluster ICM. |