.. _galactic_guide:
Galactic object models
======================
Tow models have been added to handle galactic objects:
1) :class:`.GalacticBeamed` to hadle galatcit objects that are not beamed
2) :class:`.GalacticUnbeamed` to hadle galatcit objects that are beamed
The redshift will be set 0, and not used, and a new parameter, ``DL_cm``
will be used in place.
The only difference, is that the ``GalacticBeamed`` has beaming set as
default to 10 and set free, whilst ``GalacticUnbeamed`` has beaming set
as default to 1 and frozen
.. code:: ipython3
from jetset.jet_model import GalacticUnbeamed,GalacticBeamed
from astropy import constants as const
from astropy import units as u
import matplotlib.pyplot as plt
import numpy as np
.. code:: ipython3
import jetset
print('tested on jetset',jetset.__version__)
.. parsed-literal::
tested on jetset 1.2.2
Galactic Unbeamed leptonic example
----------------------------------
This is an example for a PWN.
**Please notice, that for a realistic model you should build a composite
model with Radio and Wind region, and define a custom emitting region. A
proper thermal external field will be provided in the next version. The
``CMB`` is the same used for jets**
.. code:: ipython3
pwn=GalacticUnbeamed(emitters_distribution='bkn',verbose=False,emitters_type='electrons',distance=2*u.kpc,name='pwn')
pwn.add_EC_component(['EC_CMB'])
pwn.parameters.R.val=(2*u.pc).to('cm').value
pwn.parameters.B.val=1E-4
pwn.parameters.gamma_break.val=1E8
pwn.parameters.gmax.val=8E9
pwn.parameters.gmin.val=5E5
pwn.parameters.p.val=3.2
pwn.parameters.p_1.val=3.8
pwn.parameters.N.val=1E-8
pwn.show_model()
.. parsed-literal::
--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticUnbeamed
name: pwn
electrons distribution:
type: bkn
gamma energy grid size: 201
gmin grid : 2.000000e+00
gmax grid : 1.000000e+06
normalization: True
log-values: False
ratio of cold protons to relativistic electrons: 1.000000e-01
radiative fields:
seed photons grid size: 100
IC emission grid size: 100
source emissivity lower bound : 1.000000e-120
spectral components:
name:Sum, state: on
name:Sync, state: self-abs
name:SSC, state: on
name:EC_CMB, state: on
external fields transformation method: disk
SED info:
nu grid size jetkernel: 1000
nu size: 500
nu mix (Hz): 1.000000e+06
nu max (Hz): 1.000000e+30
flux plot lower bound : 1.000000e-30
--------------------------------------------------------------------------------
.. raw:: html
Table length=10
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
pwn | R | region_size | cm | 6.171355e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
pwn | B | magnetic_field | gauss | 1.000000e-04 | 0.000000e+00 | -- | False | False |
pwn | NH_cold_to_rel_e | cold_p_to_rel_e_ratio | | 1.000000e-01 | 0.000000e+00 | -- | False | True |
pwn | gmin | low-energy-cut-off | lorentz-factor* | 5.000000e+05 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | gmax | high-energy-cut-off | lorentz-factor* | 8.000000e+09 | 1.000000e+00 | 1.000000e+15 | False | False |
pwn | N | emitters_density | 1 / cm3 | 1.000000e-08 | 0.000000e+00 | -- | False | False |
pwn | gamma_break | turn-over-energy | lorentz-factor* | 1.000000e+08 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | p | LE_spectral_slope | | 3.200000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | p_1 | HE_spectral_slope | | 3.800000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |
.. parsed-literal::
--------------------------------------------------------------------------------
.. code:: ipython3
pwn.eval()
p=pwn.plot_model(frame='src')
.. image:: galactic_files/galactic_11_0.png
.. code:: ipython3
pwn.save_model('pwn.pkl')
.. raw:: html
Table length=10
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
pwn | R | region_size | cm | 6.171355e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
pwn | B | magnetic_field | gauss | 1.000000e-04 | 0.000000e+00 | -- | False | False |
pwn | NH_cold_to_rel_e | cold_p_to_rel_e_ratio | | 1.000000e-01 | 0.000000e+00 | -- | False | True |
pwn | gmin | low-energy-cut-off | lorentz-factor* | 5.000000e+05 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | gmax | high-energy-cut-off | lorentz-factor* | 8.000000e+09 | 1.000000e+00 | 1.000000e+15 | False | False |
pwn | N | emitters_density | 1 / cm3 | 1.000000e-08 | 0.000000e+00 | -- | False | False |
pwn | gamma_break | turn-over-energy | lorentz-factor* | 1.000000e+08 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | p | LE_spectral_slope | | 3.200000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | p_1 | HE_spectral_slope | | 3.800000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |
.. code:: ipython3
GalacticUnbeamed.load_model("pwn.pkl")
.. parsed-literal::
==> par z_cosm is now depending on ['DL_cm'] according to expr:z_cosm =
def _dummy_z_par_func(self,DL_cm):
self.cosmo._DL_cm=DL_cm*u.cm
return 0
.. raw:: html
Table length=12
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
pwn | gmin | low-energy-cut-off | lorentz-factor* | 5.000000e+05 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | gmax | high-energy-cut-off | lorentz-factor* | 8.000000e+09 | 1.000000e+00 | 1.000000e+15 | False | False |
pwn | N | emitters_density | 1 / cm3 | 1.000000e-08 | 0.000000e+00 | -- | False | False |
pwn | gamma_break | turn-over-energy | lorentz-factor* | 1.000000e+08 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | p | LE_spectral_slope | | 3.200000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | p_1 | HE_spectral_slope | | 3.800000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | R | region_size | cm | 6.171355e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
pwn | R_H | region_position | cm | 1.000000e+17 | 0.000000e+00 | -- | False | True |
pwn | B | magnetic_field | gauss | 1.000000e-04 | 0.000000e+00 | -- | False | False |
pwn | NH_cold_to_rel_e | cold_p_to_rel_e_ratio | | 1.000000e-01 | 0.000000e+00 | -- | False | True |
pwn | beam_obj | beaming | lorentz-factor* | 1.000000e+00 | 1.000000e-04 | -- | False | False |
pwn | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |
.. parsed-literal::
--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticUnbeamed
name: pwn
electrons distribution:
type: bkn
gamma energy grid size: 201
gmin grid : 5.000000e+05
gmax grid : 8.000000e+09
normalization: True
log-values: False
ratio of cold protons to relativistic electrons: 1.000000e-01
radiative fields:
seed photons grid size: 100
IC emission grid size: 100
source emissivity lower bound : 1.000000e-120
spectral components:
name:Sum, state: on
name:Sync, state: self-abs
name:SSC, state: on
name:EC_CMB, state: on
external fields transformation method: disk
SED info:
nu grid size jetkernel: 1000
nu size: 500
nu mix (Hz): 1.000000e+06
nu max (Hz): 1.000000e+30
flux plot lower bound : 1.000000e-30
--------------------------------------------------------------------------------
.. raw:: html
Table length=12
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
pwn | gmin | low-energy-cut-off | lorentz-factor* | 5.000000e+05 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | gmax | high-energy-cut-off | lorentz-factor* | 8.000000e+09 | 1.000000e+00 | 1.000000e+15 | False | False |
pwn | N | emitters_density | 1 / cm3 | 1.000000e-08 | 0.000000e+00 | -- | False | False |
pwn | gamma_break | turn-over-energy | lorentz-factor* | 1.000000e+08 | 1.000000e+00 | 1.000000e+09 | False | False |
pwn | p | LE_spectral_slope | | 3.200000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | p_1 | HE_spectral_slope | | 3.800000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
pwn | R | region_size | cm | 6.171355e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
pwn | R_H | region_position | cm | 1.000000e+17 | 0.000000e+00 | -- | False | True |
pwn | B | magnetic_field | gauss | 1.000000e-04 | 0.000000e+00 | -- | False | False |
pwn | NH_cold_to_rel_e | cold_p_to_rel_e_ratio | | 1.000000e-01 | 0.000000e+00 | -- | False | True |
pwn | beam_obj | beaming | lorentz-factor* | 1.000000e+00 | 1.000000e-04 | -- | False | False |
pwn | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |
.. parsed-literal::
--------------------------------------------------------------------------------
.. parsed-literal::
None
.. code:: ipython3
pwn.energetic_report()
.. raw:: html
Table length=16
name | type | units | val |
BulkLorentzFactor | jet-bulk-factor | | 1.000000e+00 |
U_e | Energy dens. rest. frame | erg / cm3 | 7.480812e-09 |
U_p_cold | Energy dens. rest. frame | erg / cm3 | 1.503276e-12 |
U_B | Energy dens. rest. frame | erg / cm3 | 3.978874e-10 |
U_Synch | Energy dens. rest. frame | erg / cm3 | 4.538190e-11 |
U_Disk | Energy dens. rest. frame | erg / cm3 | 0.000000e+00 |
U_BLR | Energy dens. rest. frame | erg / cm3 | 0.000000e+00 |
U_DT | Energy dens. rest. frame | erg / cm3 | 0.000000e+00 |
U_CMB | Energy dens. rest. frame | erg / cm3 | 4.185297e-13 |
L_Sync | Lum. rest. frame. | erg / s | 2.170466e+38 |
L_SSC | Lum. rest. frame. | erg / s | 2.161006e+35 |
L_EC_Disk | Lum. rest. frame. | erg / s | 0.000000e+00 |
L_EC_BLR | Lum. rest. frame. | erg / s | 0.000000e+00 |
L_EC_DT | Lum. rest. frame. | erg / s | 0.000000e+00 |
L_EC_CMB | Lum. rest. frame. | erg / s | 1.785622e+35 |
NH_cold_to_rel_e | cold_p_to_rel_e_ratio | | 1.000000e-01 |
Galactic Unbeamed hadronic pp example
-------------------------------------
This is an example for a galactic unbeamed hadronic (pp) object
.. code:: ipython3
gal_hadronic=GalacticUnbeamed(emitters_distribution='plc',verbose=False,emitters_type='protons',distance=2*u.kpc,name='gal_hadronic_unbeamed')
gal_hadronic.parameters.R.val=1E18
gal_hadronic.parameters.N.val=1000
gal_hadronic.parameters.B.val=1E-3
gal_hadronic.show_model()
.. parsed-literal::
--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticUnbeamed
name: gal_hadronic_unbeamed
protons distribution:
type: plc
gamma energy grid size: 201
gmin grid : 2.000000e+00
gmax grid : 1.000000e+06
normalization: True
log-values: False
radiative fields:
seed photons grid size: 100
IC emission grid size: 100
source emissivity lower bound : 1.000000e-120
spectral components:
name:Sum, state: on
name:Sync, state: self-abs
name:SSC, state: on
name:PP_gamma, state: on
name:PP_neutrino_tot, state: on
name:PP_neutrino_mu, state: on
name:PP_neutrino_e, state: on
name:Bremss_ep, state: on
external fields transformation method: disk
SED info:
nu grid size jetkernel: 1000
nu size: 500
nu mix (Hz): 1.000000e+06
nu max (Hz): 1.000000e+30
flux plot lower bound : 1.000000e-30
--------------------------------------------------------------------------------
.. raw:: html
Table length=9
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
gal_hadronic_unbeamed | R | region_size | cm | 1.000000e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
gal_hadronic_unbeamed | B | magnetic_field | gauss | 1.000000e-03 | 0.000000e+00 | -- | False | False |
gal_hadronic_unbeamed | gmin | low-energy-cut-off | lorentz-factor* | 2.000000e+00 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_unbeamed | gmax | high-energy-cut-off | lorentz-factor* | 1.000000e+06 | 1.000000e+00 | 1.000000e+15 | False | False |
gal_hadronic_unbeamed | N | emitters_density | 1 / cm3 | 1.000000e+03 | 0.000000e+00 | -- | False | False |
gal_hadronic_unbeamed | NH_pp | target_density | 1 / cm3 | 1.000000e+00 | 0.000000e+00 | -- | False | False |
gal_hadronic_unbeamed | gamma_cut | turn-over-energy | lorentz-factor* | 1.000000e+04 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_unbeamed | p | LE_spectral_slope | | 2.000000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
gal_hadronic_unbeamed | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |
.. parsed-literal::
--------------------------------------------------------------------------------
.. code:: ipython3
gal_hadronic.eval()
p=gal_hadronic.plot_model(frame='src')
p.setlim(y_min=1E33)
.. image:: galactic_files/galactic_18_0.png
.. code:: ipython3
gal_hadronic.energetic_report()
.. raw:: html
Table length=17
name | type | units | val |
BulkLorentzFactor | jet-bulk-factor | | 1.000000e+00 |
U_e | Energy dens. rest. frame | erg / cm3 | 4.252875e-06 |
U_B | Energy dens. rest. frame | erg / cm3 | 3.978874e-08 |
U_p | Energy dens. rest. frame | erg / cm3 | 2.386056e+01 |
U_p_target | Energy dens. rest. frame | erg / cm3 | 1.503276e-03 |
U_Synch | Energy dens. rest. frame | erg / cm3 | 4.851771e-10 |
U_Disk | Energy dens. rest. frame | erg / cm3 | 0.000000e+00 |
U_BLR | Energy dens. rest. frame | erg / cm3 | 0.000000e+00 |
U_DT | Energy dens. rest. frame | erg / cm3 | 0.000000e+00 |
U_CMB | Energy dens. rest. frame | erg / cm3 | 0.000000e+00 |
L_Sync | Lum. rest. frame. | erg / s | 6.085926e+37 |
L_SSC | Lum. rest. frame. | erg / s | 2.273252e+35 |
L_EC_Disk | Lum. rest. frame. | erg / s | 0.000000e+00 |
L_EC_BLR | Lum. rest. frame. | erg / s | 0.000000e+00 |
L_EC_DT | Lum. rest. frame. | erg / s | 0.000000e+00 |
L_EC_CMB | Lum. rest. frame. | erg / s | 0.000000e+00 |
L_pp_gamma | Lum. rest. frame. | erg / s | 1.217784e+40 |
Galactic Beamed hadronic pp example
-----------------------------------
.. code:: ipython3
from jetset.jet_model import GalacticUnbeamed,GalacticBeamed
from astropy import constants as const
from astropy import units as u
import matplotlib.pyplot as plt
import numpy as np
gal_hadronic=GalacticBeamed(emitters_distribution='plc',verbose=False,emitters_type='protons',distance=2*u.kpc,name='gal_hadronic_beamed')
gal_hadronic.parameters.R.val=1E18
gal_hadronic.parameters.N.val=1000
gal_hadronic.parameters.B.val=1E-3
gal_hadronic.show_model()
.. parsed-literal::
--------------------------------------------------------------------------------
model description:
--------------------------------------------------------------------------------
type: GalacticBeamed
name: gal_hadronic_beamed
protons distribution:
type: plc
gamma energy grid size: 201
gmin grid : 2.000000e+00
gmax grid : 1.000000e+06
normalization: True
log-values: False
radiative fields:
seed photons grid size: 100
IC emission grid size: 100
source emissivity lower bound : 1.000000e-120
spectral components:
name:Sum, state: on
name:Sync, state: self-abs
name:SSC, state: on
name:PP_gamma, state: on
name:PP_neutrino_tot, state: on
name:PP_neutrino_mu, state: on
name:PP_neutrino_e, state: on
name:Bremss_ep, state: on
external fields transformation method: blob
SED info:
nu grid size jetkernel: 1000
nu size: 500
nu mix (Hz): 1.000000e+06
nu max (Hz): 1.000000e+30
flux plot lower bound : 1.000000e-30
--------------------------------------------------------------------------------
.. raw:: html
Table length=11
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
gal_hadronic_beamed | R | region_size | cm | 1.000000e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
gal_hadronic_beamed | R_H | region_position | cm | 1.000000e+17 | 0.000000e+00 | -- | False | True |
gal_hadronic_beamed | B | magnetic_field | gauss | 1.000000e-03 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | beam_obj | beaming | lorentz-factor* | 1.000000e+01 | 1.000000e-04 | -- | False | False |
gal_hadronic_beamed | gmin | low-energy-cut-off | lorentz-factor* | 2.000000e+00 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_beamed | gmax | high-energy-cut-off | lorentz-factor* | 1.000000e+06 | 1.000000e+00 | 1.000000e+15 | False | False |
gal_hadronic_beamed | N | emitters_density | 1 / cm3 | 1.000000e+03 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | NH_pp | target_density | 1 / cm3 | 1.000000e+00 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | gamma_cut | turn-over-energy | lorentz-factor* | 1.000000e+04 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_beamed | p | LE_spectral_slope | | 2.000000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
gal_hadronic_beamed | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |
.. parsed-literal::
--------------------------------------------------------------------------------
.. code:: ipython3
gal_hadronic.eval()
p=gal_hadronic.plot_model(frame='src')
p.setlim(y_min=1E36)
.. image:: galactic_files/galactic_22_0.png
.. code:: ipython3
gal_hadronic.energetic_report()
.. raw:: html
Table length=35
name | type | units | val |
BulkLorentzFactor | jet-bulk-factor | | 1.000000e+01 |
U_e | Energy dens. blob rest. frame | erg / cm3 | 4.252875e-06 |
U_B | Energy dens. blob rest. frame | erg / cm3 | 3.978874e-08 |
U_p | Energy dens. blob rest. frame | erg / cm3 | 2.386056e+01 |
U_p_target | Energy dens. blob rest. frame | erg / cm3 | 1.503276e-03 |
U_Synch | Energy dens. blob rest. frame | erg / cm3 | 4.851771e-10 |
U_Synch_DRF | Energy dens. disk rest. frame | erg / cm3 | 4.851771e-06 |
U_Disk | Energy dens. blob rest. frame | erg / cm3 | 0.000000e+00 |
U_BLR | Energy dens. blob rest. frame | erg / cm3 | 0.000000e+00 |
U_DT | Energy dens. blob rest. frame | erg / cm3 | 0.000000e+00 |
U_CMB | Energy dens. blob rest. frame | erg / cm3 | 0.000000e+00 |
U_Disk_DRF | Energy dens. disk rest. frame | erg / cm3 | 0.000000e+00 |
U_BLR_DRF | Energy dens. disk rest. frame | erg / cm3 | 0.000000e+00 |
U_DT_DRF | Energy dens. disk rest. frame | erg / cm3 | 0.000000e+00 |
U_CMB_DRF | Energy dens. disk rest. frame | erg / cm3 | 0.000000e+00 |
L_Sync_rf | Lum. blob rest. frame. | erg / s | 6.085926e+37 |
L_SSC_rf | Lum. blob rest. frame. | erg / s | 2.273252e+35 |
L_EC_Disk_rf | Lum. blob rest. frame. | erg / s | 0.000000e+00 |
L_EC_BLR_rf | Lum. blob rest. frame. | erg / s | 0.000000e+00 |
L_EC_DT_rf | Lum. blob rest. frame. | erg / s | 0.000000e+00 |
L_EC_CMB_rf | Lum. blob rest. frame. | erg / s | 0.000000e+00 |
L_pp_gamma_rf | Lum. blob rest. frame. | erg / s | 1.217784e+40 |
jet_L_Sync | jet Lum. | erg / s | 1.521482e+39 |
jet_L_SSC | jet Lum. | erg / s | 5.683129e+36 |
jet_L_EC_Disk | jet Lum. | erg / s | 0.000000e+00 |
jet_L_EC_BLR | jet Lum. | erg / s | 0.000000e+00 |
jet_L_EC_DT | jet Lum. | erg / s | 0.000000e+00 |
jet_L_EC_CMB | jet Lum. | erg / s | 0.000000e+00 |
jet_L_pp_gamma | jet Lum. | erg / s | 3.044459e+41 |
jet_L_rad | jet Lum. | erg / s | 3.059730e+41 |
jet_L_kin | jet Lum. | erg / s | 2.235985e+50 |
jet_L_tot | jet Lum. | erg / s | 2.235985e+50 |
jet_L_e | jet Lum. | erg / s | 3.985390e+43 |
jet_L_B | jet Lum. | erg / s | 3.728622e+41 |
jet_L_p | jet Lum. | erg / s | 2.235985e+50 |
.. code:: ipython3
gal_hadronic.save_model('gal_hadronic_beamed.pkl')
.. raw:: html
Table length=11
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
gal_hadronic_beamed | R | region_size | cm | 1.000000e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
gal_hadronic_beamed | R_H | region_position | cm | 1.000000e+17 | 0.000000e+00 | -- | False | True |
gal_hadronic_beamed | B | magnetic_field | gauss | 1.000000e-03 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | beam_obj | beaming | lorentz-factor* | 1.000000e+01 | 1.000000e-04 | -- | False | False |
gal_hadronic_beamed | gmin | low-energy-cut-off | lorentz-factor* | 2.000000e+00 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_beamed | gmax | high-energy-cut-off | lorentz-factor* | 1.000000e+06 | 1.000000e+00 | 1.000000e+15 | False | False |
gal_hadronic_beamed | N | emitters_density | 1 / cm3 | 1.000000e+03 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | NH_pp | target_density | 1 / cm3 | 1.000000e+00 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | gamma_cut | turn-over-energy | lorentz-factor* | 1.000000e+04 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_beamed | p | LE_spectral_slope | | 2.000000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
gal_hadronic_beamed | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |
.. code:: ipython3
gal_hadronic_new=GalacticBeamed.load_model('gal_hadronic_beamed.pkl')
.. parsed-literal::
==> par z_cosm is now depending on ['DL_cm'] according to expr:z_cosm =
def _dummy_z_par_func(self,DL_cm):
self.cosmo._DL_cm=DL_cm*u.cm
return 0
.. raw:: html
Table length=11
model name | name | par type | units | val | phys. bound. min | phys. bound. max | log | frozen |
gal_hadronic_beamed | gmin | low-energy-cut-off | lorentz-factor* | 2.000000e+00 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_beamed | gmax | high-energy-cut-off | lorentz-factor* | 1.000000e+06 | 1.000000e+00 | 1.000000e+15 | False | False |
gal_hadronic_beamed | N | emitters_density | 1 / cm3 | 1.000000e+03 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | NH_pp | target_density | 1 / cm3 | 1.000000e+00 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | gamma_cut | turn-over-energy | lorentz-factor* | 1.000000e+04 | 1.000000e+00 | 1.000000e+09 | False | False |
gal_hadronic_beamed | p | LE_spectral_slope | | 2.000000e+00 | -1.000000e+01 | 1.000000e+01 | False | False |
gal_hadronic_beamed | R | region_size | cm | 1.000000e+18 | 1.000000e+03 | 1.000000e+30 | False | False |
gal_hadronic_beamed | R_H | region_position | cm | 1.000000e+17 | 0.000000e+00 | -- | False | True |
gal_hadronic_beamed | B | magnetic_field | gauss | 1.000000e-03 | 0.000000e+00 | -- | False | False |
gal_hadronic_beamed | beam_obj | beaming | lorentz-factor* | 1.000000e+01 | 1.000000e-04 | -- | False | False |
gal_hadronic_beamed | DL_cm(M) | distance | cm | 6.171355e+21 | 0.000000e+00 | 3.085678e+24 | False | False |