Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Guang
x-cigale
Commits
d35b391b
Commit
d35b391b
authored
Jul 01, 2019
by
Guang
Committed by
Guang
Oct 23, 2019
Browse files
cut AGN disk SED to lam>10^1.5 nm; extend X-ray down to lam=10^1.2
parent
ec286843
Changes
2
Hide whitespace changes
Inline
Side-by-side
database_builder/__init__.py
View file @
d35b391b
...
@@ -794,7 +794,7 @@ def build_skirtor2016(base):
...
@@ -794,7 +794,7 @@ def build_skirtor2016(base):
disk
+=
scatt
disk
+=
scatt
# Apply wavelength cut to avoid X-ray wavelength
# Apply wavelength cut to avoid X-ray wavelength
lam_cut
=
5
.
lam_cut
=
10
**
1.
5
lam_idxs
=
wl_uncut
>=
lam_cut
lam_idxs
=
wl_uncut
>=
lam_cut
# Calculate the re-normalization factor to keep energy conservation
# Calculate the re-normalization factor to keep energy conservation
norm_fac
=
np
.
trapz
(
intrin_disk
,
x
=
wl_uncut
)
/
\
norm_fac
=
np
.
trapz
(
intrin_disk
,
x
=
wl_uncut
)
/
\
...
...
pcigale/sed_modules/xray.py
View file @
d35b391b
...
@@ -65,8 +65,8 @@ class Xray(SedModule):
...
@@ -65,8 +65,8 @@ class Xray(SedModule):
self
.
nu_2keV
=
c
/
lam_2keV
self
.
nu_2keV
=
c
/
lam_2keV
# We define the wavelength range for the non thermal emission
# We define the wavelength range for the non thermal emission
# corresponding to 0.
4
-1200 keV
# corresponding to 0.
08
-1200 keV
self
.
wave
=
np
.
logspace
(
-
3
,
0.5
,
1000
)
self
.
wave
=
np
.
logspace
(
-
3
,
1.2
,
1000
)
# X-ray emission from galaxies: 1.hot-gas & 2.X-ray binaries
# X-ray emission from galaxies: 1.hot-gas & 2.X-ray binaries
# 1.Hot-gas, assuming power-law index gamma=3
# 1.Hot-gas, assuming power-law index gamma=3
...
@@ -91,8 +91,8 @@ class Xray(SedModule):
...
@@ -91,8 +91,8 @@ class Xray(SedModule):
# Calculate total AGN corona X-ray luminosity
# Calculate total AGN corona X-ray luminosity
self
.
l_agn_xray_total
=
np
.
trapz
(
self
.
lumin_corona
,
x
=
self
.
wave
)
self
.
l_agn_xray_total
=
np
.
trapz
(
self
.
lumin_corona
,
x
=
self
.
wave
)
# Calculate 2-10 keV AGN corona X-ray luminosity
# Calculate 2-10 keV AGN corona X-ray luminosity
# 2-10 keV corredponds to self.wave[598:798]
lam_idxs
=
(
self
.
wave
<=
self
.
lam_2keV
)
&
(
self
.
wave
>=
self
.
lam_10keV
)
self
.
l_agn_xray_2to10keV
=
np
.
trapz
(
self
.
lumin_corona
[
598
:
798
],
x
=
self
.
wave
[
598
:
798
])
self
.
l_agn_xray_2to10keV
=
np
.
trapz
(
self
.
lumin_corona
[
lam_idxs
],
x
=
self
.
wave
[
lam_idxs
])
def
process
(
self
,
sed
):
def
process
(
self
,
sed
):
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment