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
|
# 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 Galactic emission
# components, which currently includes the following components:
# - synchrotron
# - freefree
# - snr
#
# NOTE:
# - The input templates for simulations should be HEALPix full-sky maps.
# - The input catalog should be in CSV format.
[galactic]
# Synchrotron emission component (unpolarized)
[[synchrotron]]
# The template map for the simulation.
template = string(default=None)
# The frequency of the template map (same unit as in [frequency] section)
template_freq = float(default=None, min=0.0)
# The unit of the template map pixel
template_unit = string(default=None)
# Spectral index map
indexmap = string(default=None)
# Whether add fluctuations on the small scales according the angular
# power spectrum prediction?
add_smallscales = boolean(default=True)
# Range of multipole moments (l) of the angular power spectrum.
# The power spectrum will be cut off to a constant for multipole l < lmin.
# NOTE: Update the ``lmax`` accordingly w.r.t. ``sky/healpix/nside``.
# Generally, lmax = 3 * nside - 1
lmin = integer(min=0, default=10)
lmax = integer(min=1, default=3071)
# Filename prefix for this component
prefix = string(default="gsync")
# Whether save this component to disk
save = boolean(default=True)
# Output directory to save the simulated results
output_dir = string(default=None)
# Free-free bremsstrahlung emission component
[[freefree]]
# The H{\alpha} map used as the free-free emission template
halphamap = string(default=None)
# The unit of the H{\alpha} template (e.g., "Rayleigh")
halphamap_unit = string(default=None)
# The 100-{\mu}m dust map used for dust absorption correction
dustmap = string(default=None)
# The unit of the above dust map (e.g., "MJy/sr")
dustmap_unit = string(default=None)
# Effective dust fraction in the LoS actually absorbing Halpha
dust_fraction = float(default=0.33, min=0.1, max=1.0)
# Halpha absorption threshold:
# When the dust absorption goes rather large, the true Halpha
# absorption can not well determined. This configuration sets the
# threshold below which the dust absorption can be well determined,
# while the sky regions with higher absorption are masked out due
# to unreliable absorption correction.
# Unit: [mag]
halpha_abs_th = float(default=1.0)
# The electron temperature assumed for the ionized interstellar medium
# that generating H{\alpha} emission.
# Unit: [K]
electron_temperature = float(default=7000.0, min=1000)
# Filename prefix for this component
prefix = string(default="gfree")
# Whether save this component to disk
save = boolean(default=True)
# Output directory to save the simulated results
output_dir = string(default=None)
# Supernova remnants emission
[[snr]]
# The Galactic SNRs catalog data (CSV file)
catalog = string(default=None)
# Output the effective/inuse SNRs catalog data (CSV file)
catalog_outfile = string(default=None)
# Resolution for simulating each SNR template, which are finally
# mapped to the all-sky HEALPix map if used.
# Unit: [arcsec]
resolution = float(default=30.0, min=5.0)
# Filename prefix for this component
prefix = string(default="gsnr")
# Whether save this component to disk
save = boolean(default=True)
# Output directory to save the simulated results
output_dir = string(default=None)
|