aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim/extragalactic/clusters')
-rw-r--r--fg21sim/extragalactic/clusters/halo.py27
1 files changed, 20 insertions, 7 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 4d0e41b..c83e791 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -423,19 +423,20 @@ class RadioHalo:
n_inj = self.fp_injection(self.gamma)
n0_e = n_inj * (self.age_begin - self.time_init)
- logger.debug("Derive the initial electron spectrum ...")
+ logger.debug("Deriving the initial electron spectrum ...")
self._acceleration_disabled = True
- dt = self.fpsolver.tstep
tstart = self.age_begin
tstop = self.age_begin + self.time_init
- self.fpsolver.tstep = 3 * dt # Bigger step to save time
+ self.fpsolver.tstep = self.time_step * 3 # To save time
+
n_e = self.fpsolver.solve(u0=n0_e, tstart=tstart, tstop=tstop)
- self.fpsolver.tstep = dt
self._acceleration_disabled = False
+ self.fpsolver.tstep = self.time_step
return n_e
- def calc_electron_spectrum(self, tstart=None, tstop=None, n0_e=None):
+ def calc_electron_spectrum(self, tstart=None, tstop=None, n0_e=None,
+ fiducial=False):
"""
Calculate the relativistic electron spectrum by solving the
Fokker-Planck equation.
@@ -456,11 +457,15 @@ class RadioHalo:
The initial electron spectrum (number distribution).
Default: ``self.electron_spec_init``
Unit: [cm^-3]
+ fiducial : bool
+ Whether to disable the turbulent acceleration and derive the
+ fiducial electron spectrum?
+ Default: ``False``
Returns
-------
n_e : float 1D `~numpy.ndarray`
- The solved electron spectrum at ``tstop``.
+ The solved electron spectrum.
Unit: [cm^-3]
"""
if tstart is None:
@@ -469,7 +474,15 @@ class RadioHalo:
tstop = self.age_obs
if n0_e is None:
n0_e = self.electron_spec_init
- return self.fpsolver.solve(u0=n0_e, tstart=tstart, tstop=tstop)
+ if fiducial:
+ self._acceleration_disabled = True
+ self.fpsolver.tstep = self.time_step * 2 # To save time
+
+ n_e = self.fpsolver.solve(u0=n0_e, tstart=tstart, tstop=tstop)
+ self._acceleration_disabled = False
+ self.fpsolver.tstep = self.time_step
+
+ return n_e
def fp_injection(self, gamma, t=None):
"""