| Commit message (Collapse) | Author | Age | Files | Lines |
|
|
|
|
| |
The electrons are assumed to be injected throughout the whole cluster
ICM/volume.
|
| |
|
| |
|
| |
|
|
|
|
|
|
|
|
|
|
| |
* Adopt the electron acceleration coefficient formula from [cassano2005]
* Rename method "_tau_acceleration()" to "_chi_acceleration()", and
rewrite
* Add property "kT_merger"
* Also save "kT_merger" and "chi" into halos data
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
Also avoid forgetting to add the newly added item for DataFrame
conversion.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
| |
* Brought back "clobber" property
* Delete the wrong "self.halos=[]"
* Fix the wrong assignment of "cimax0"
* Assign default values for {c,r}i{min,max}1
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Also update the documents of "draw_halo()" function
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
* Clean up the unused configurations
* Add "FREQ" keyword to the sky FITS header
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
New functions "halo_rprofile()" and "draw_halo()" added to helper.py
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Also rename "frequencies" to "frequency"
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Help development & debug ...
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
|
|
| |
For unknown reason! Each instance of "SynchrotronEmission" will cost
~5 MiB memory when interpolating the kernel function cached as
"self._F_interp", which won't be collected/freed by GC, for UNKNOWN
reason!
Therefore, make the method "F()" be class-wide, to avoid the memory hole!
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
| |
|
| |
|
|
|
|
|
|
| |
Also update logging message for RadioHalo a bit.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Also updated configuration options.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
This method simulates the radio halo properties for each cluster with a
recent major merger.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
| |
* Methods "power()", "flux()" already moved to "halo.py";
* Method "brightness()" will also be moved over there;
* Clean up thus unused parameters "z" and "radius", as well as imports.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
Calculate the synchrotron power (a.k.a. specific luminosity) and flux
density from the emissivity for the halo.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
This method invokes ``SynchrotronEmission`` to calculate the synchrotron
emissivity at specified frequencies for the derived electron spectrum.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
|
|
| |
NOTE:
Since we just substitute ``p`` (electron momentum) for
``γ`` (electron Lorentz factor), i.e.,
N = int_p n(p) dp = int_γ n(γ) dγ ,
therefore the coefficient ``sqrt(3) e^3 B / (m c^2)`` does NOT need
modification (e.g., multiply an extra ``m c``)!
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
|
| |
Use 2D grid of (discrete) samples to optimize integration speed,
avoiding the more complicated integration w.r.t. functions.
WARNING:
Current calculation results seems wrong!
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|