From d82487346e87aedde6b9a476d43c7b1e1a821319 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Mon, 1 Jan 2018 14:51:55 +0800 Subject: Minor changes --- fg21sim/configs/20-extragalactic.conf.spec | 2 +- fg21sim/extragalactic/clusters/helper.py | 44 +++++++++++++++--------------- 2 files changed, 23 insertions(+), 23 deletions(-) (limited to 'fg21sim') 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. -- cgit v1.2.2