1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
|
# Configurations for "fg21sim"
# -*- mode: conf -*-
#
# Syntax: `ConfigObj`, https://github.com/DiffSK/configobj
#
# There are various configuration options that specify the input data
# and/or templates required by the simulations, the properties of the input
# data, the output products, as well as some parameters affecting the
# simulation behaviors.
#
# This file contains the options corresponding the extragalactic emission
# components, which currently includes the following components:
# - clusters
# - halos
# - pointsources
#
# NOTE:
# - The input templates for simulations should be HEALPix full-sky maps.
# - The input catalog should be in CSV format.
[extragalactic]
# Emissions from the clusters of galaxies
[[clusters]]
# The clusters catalog derived from the Hubble Volume Project (CSV file)
catalog = string(default=None)
# Output the effective/inuse clusters catalog data (CSV file)
catalog_outfile = string(default=None)
# The fraction that a cluster hosts a radio halo
halo_fraction = float(default=None, min=0.0, max=1.0)
# Resolution (unit: arcmin) for simulating each cluster, which are finally
# mapped to the HEALPix map of Nside specified in "[common]" section.
resolution = float(default=0.5, min=0.0)
# Filename prefix for this component
prefix = string(default="egcluster")
# Whether save this component to disk
save = boolean(default=True)
# Output directory to save the simulated results
output_dir = string(default=None)
# Emission of giant radio halos from galaxy clusters
[[halos]]
# Maximum redshift until where to tracing the cluster merging history
# (e.g., when calculating the electron spectrum)
zmax = float(default=3.0, min=0.0)
# Mass threshold of the sub-cluster to be regarded as a significant
# merger. (unit: Msun)
merger_mass_th = float(default=1e13, min=1e12)
# Radius of the giant radio halo in clusters (unit: kpc)
# XXX: currently only support a constant radius of halos
radius_halo = float(default=500, min=100)
# Magnetic field assumed for the cluster (unit: uG)
# XXX: currently only support a constant magnetic field
magnetic_field = float(default=0.5, min=0.1, max=10)
# Fraction of the turbulence energy in the form of magneto-sonic waves.
eta_t = float(default=0.3, min=0.0, max=1.0)
# Ratio of the total energy injected in cosmic-ray electrons during the
# cluster life to the present-day total thermal energy of the cluster.
eta_e = float(default=0.003, min=0.0, max=0.1)
# Minimum and maximum Lorentz factor (i.e., energy) of the relativistic
# electron spectrum.
pmin = float(default=1.0)
pmax = float(default=1e5)
# Number of points for the grid used during solving the Fokker-Planck
# equation to calculate the electron spectrum.
pgrid_num = integer(default=100, min=10)
# Number of grid points used as the buffer region near the lower
# boundary, and the value within this buffer region will be fixed to
# avoid unphysical pile-up of low-energy electrons.
buffer_np = integer(default=5, min=0)
# Time step for solving the Fokker-Planck equation (unit: Gyr)
time_step = float(default=0.01, min=1e-5, max=1.0)
# Index of the power-law spectrum assumed for the injected electrons.
injection_index = float(default=2.5)
# Extragalactic point sources
[[pointsources]]
# Whether save this point source catelogue to disk
save = boolean(default=True)
# Output directory to save the simulated catelogues
output_dir = string(default="PS_tables")
# PS components to be simulated
pscomponents=string_list(default=list())
# Resolution [arcmin]
resolution=float(default=0.6, min=0.0)
# Number of each type of point source
# Star forming
[[[starforming]]]
# Number of samples
numps = integer(default=1000)
# Prefix
prefix = string(default="SF")
[[[starbursting]]]
# Number of samples
numps = integer(default=1000)
# Prefix
prefix = string(default="SB")
[[[radioquiet]]]
# Number of samples
numps = integer(default=1000)
# Prefix
prefix = string(default="RQ")
[[[FRI]]]
# Number of samples
numps = integer(default=1000)
# Prefix
prefix = string(default="FRI")
[[[FRII]]]
# Number of samples
numps = integer(default=1000)
# Prefix
prefix = string(default="FRII")
|