| Commit message (Collapse) | Author | Age | Files | Lines |
|
|
|
|
| |
Also update clusters/main.py to calculate the halo mass distributions before
sampling the mass and redshifts for clusters.
|
|
|
|
|
|
|
|
|
|
|
|
| |
* New dependency "hmf" (halo mass functions) module
* Calculate halo mass distributions/functions (dndlnm) with respect to masses
and redshifts, instead of use the previous data file ("ps_data")
* New section "[extragalactic][psformalism]" in configurations
* New functions to write and read the dndlnm data
TODO:
* update the method to sample (mass, redshift) for clusters from the dndlnm
data
|
| |
|
|
|
|
| |
And some cleanup and small changes
|
|
|
|
|
|
|
|
|
|
|
| |
Add asymptotic functions to calculate the values beyond the
interpolation bounds (e.g., <1e-3 and >10), otherwise, the calculated
synchrotron emissivity is overestimated at the higher frequencies.
By rewrite this synchrotron kernel function, the calculated results is
consistent with the theoretical/analytical results, e.g., the
synchrotron radiation of a population of electrons of power-law index n
is also a power-law with index (n-1)/2.
|
| |
|
|
|
|
|
|
|
| |
NOTE:
The cluster number (especially the bright ones) may be rather small
within a sky patch. This ``boost`` config option increase the expected
cluster number by a specified times, for better testing purpose.
|
| |
|
|
|
|
| |
Use "1 - baryon_fraction = 1 - Ob0/Om0" to replace "f_darkmatter".
|
|
|
|
|
| |
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>
|