aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic
Commit message (Collapse)AuthorAgeFilesLines
* clusters/halo: Add set_electron_spectrum() method; update doc a bitAaron LI2017-10-191-1/+26
|
* clusters/halo: Add and use attribute "self.halos_data_dumpfile"Aaron LI2017-10-191-2/+4
|
* clusters/halo: Rewrite halos catalog and raw data saveAaron LI2017-10-191-15/+39
| | | | | Replace option "halos_dumpfile" with "halos_catalog_outfile", and use option "dump_halos_data" to control whether to dump the raw data.
* clusters: "catalog_outfile" is requiredAaron LI2017-10-191-3/+1
|
* clusters/halo: Improve time step determination a bitAaron LI2017-10-192-3/+6
|
* clusters/halo: Increase the time step, but require a minimum of 20 stepsAaron LI2017-10-181-0/+7
|
* clusters: Update radio halo radius estimationAaron LI2017-10-182-23/+25
| | | | | | | | | The simulated radio halo is assumed to have a radius of the falling sub-cluster; while previously it is assumed to be 1/4 of the virial radius of the merged cluster. The new estimation can agree better with the currently observed radio halos, which generally have a angular diameter size of ~2-7 [arcmin].
* clusters: More regular shape for radio halosAaron LI2017-10-161-6/+5
|
* clusters: Update the documents/comments a b itAaron LI2017-10-142-11/+6
|
* clusters/solver.py: Add FokkerPlanckTests with 3 casesAaron LI2017-10-121-9/+158
| | | | | | | | The FokkerPlanckSolver is validated with all the 3 test cases! References: * Park & Petrosian 1996, ApJS, 103, 255 * Donnert & Brunetti 2014, MNRAS, 443, 3564
* clusters/solver.py: significantly improve fix_boundary()Aaron LI2017-10-121-10/+30
| | | | Also suggest that ``buffer_np`` be specified to 5%-10% of ``x_np``.
* clusters/solver.py: Small tweaksAaron LI2017-10-071-8/+6
|
* clusters/solver.py: Adjust parameter tc,uc ordersAaron LI2017-10-061-10/+7
|
* clusters/solver.py: Stop calc. at "tstop" by adding "tstep" parameterAaron LI2017-10-061-12/+19
|
* clusters/solver.py: Update comments/references/descriptions etc.Aaron LI2017-10-061-25/+31
|
* clusters/psformalism.py: Rewrite sample_z_m() methodAaron LI2017-10-042-38/+54
| | | | | Also update clusters/main.py to calculate the halo mass distributions before sampling the mass and redshifts for clusters.
* clusters: Use "hmf" to calculate halo mass functions/distributionsAaron LI2017-10-041-69/+136
| | | | | | | | | | | | * 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
* clusters/psformalism.py: Improve catalog commentsAaron LI2017-10-021-4/+7
|
* clusters/emission: Assume electron pitch angle be pi/2Aaron LI2017-09-201-23/+17
| | | | And some cleanup and small changes
* clusters/emission: Rewrite synchrotron kernel functionAaron LI2017-09-201-9/+38
| | | | | | | | | | | 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.
* clusters: Add option "use_output_catalog" to continue halo simulationAaron LI2017-09-031-4/+19
|
* clusters: Add config option "boost"Aaron LI2017-09-011-2/+5
| | | | | | | 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.
* Update "Fnu_to_Tb_fast()" to "Fnu_to_Tb()"Aaron LI2017-09-016-9/+9
|
* clusters/psformalism: Remove config "f_darkmatter"Aaron LI2017-08-291-3/+2
| | | | Use "1 - baryon_fraction = 1 - Ob0/Om0" to replace "f_darkmatter".
* clusters/halo.py: Update injection rate calculationAaron LI2017-08-261-10/+5
| | | | | The electrons are assumed to be injected throughout the whole cluster ICM/volume.
* Return the list of filepath to the written sky filesAaron LI2017-08-261-0/+8
|
* clusters/main.py: Update sky header; update comments etc.Aaron LI2017-08-261-5/+7
|
* configs: Remove option "save"; all enabled components will be savedAaron LI2017-08-261-1/+0
|
* clusters/halo.py: Rewrite electron acceleration coefficient calcAaron LI2017-08-142-26/+53
| | | | | | | | | | * 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>
* clusters/main.py: Use OrderedDict to easy keys manipulationAaron LI2017-08-141-29/+31
| | | | | | | Also avoid forgetting to add the newly added item for DataFrame conversion. Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Simplify two coefficients usageAaron LI2017-08-141-7/+5
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters: Rename property "injection_rate" and save into halos dataAaron LI2017-08-142-37/+38
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* Use "GalaxyClusters" from main.py, and remove clusters.pyAaron LI2017-08-142-734/+2
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* Fix several bugsAaron LI2017-08-131-11/+3
| | | | | | | | | * 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>
* clusters/main.py: Implement "simulate_frequency()" and "simulate()"Aaron LI2017-08-131-0/+46
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters: Halo templates are normalized to have mean of 1Aaron LI2017-08-132-4/+12
| | | | | | Also update the documents of "draw_halo()" function Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/main.py: Add method "_outfilepath()"Aaron LI2017-08-131-0/+21
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/main.py: Set sky header propertiesAaron LI2017-08-131-4/+5
| | | | | | | * Clean up the unused configurations * Add "FREQ" keyword to the sky FITS header Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/main.py: Implement "_draw_halos()" methodAaron LI2017-08-122-2/+89
| | | | | | New functions "halo_rprofile()" and "draw_halo()" added to helper.py Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/main.py: Inherent cluster catalog info to halos dataAaron LI2017-08-121-4/+9
| | | | | | Also rename "frequencies" to "frequency" Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/main.py: Add missing columns to DataFrame conversionAaron LI2017-08-121-2/+2
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters: Implement calc_{power,flux,brightness_mean} in helperAaron LI2017-08-112-40/+116
| | | | | | Help development & debug ... Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo: "volume" use unit "kpc^3" instead of "cm^3"Aaron LI2017-08-112-4/+7
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/main.py: Calculate and cache more halo dataAaron LI2017-08-061-4/+9
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* Fix a bug in clusters/halo.py; Add a comment to utils/io/_create_dir()Aaron LI2017-08-061-1/+1
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* Update against "Fnu_to_Tb_fast()" unit changesAaron LI2017-08-036-28/+36
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/main.py: Convert halos data to be a DataFrame and dumpAaron LI2017-08-031-0/+19
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/emission.py: Make "F()" be class method, fixing memory hole!Aaron LI2017-08-011-22/+35
| | | | | | | | | | | 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>
* Improve logging messagesAaron LI2017-08-012-4/+9
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/emission.py: Fix bug in "F()"Aaron LI2017-08-011-13/+14
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
Copyright © 2017-2018 Aaron LI. All Rights Reserved.