Galactic object models
Tow models have been added to handle galactic objects:
1) GalacticBeamed
to hadle galatcit objects that are not beamed
2) 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
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
import jetset
print('tested on jetset',jetset.__version__)
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
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()
--------------------------------------------------------------------------------
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
--------------------------------------------------------------------------------
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 |
--------------------------------------------------------------------------------
pwn.eval()
p=pwn.plot_model(frame='src')
pwn.save_model('pwn.pkl')
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 |
GalacticUnbeamed.load_model("pwn.pkl")
==> 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
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 |
--------------------------------------------------------------------------------
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
--------------------------------------------------------------------------------
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 |
--------------------------------------------------------------------------------
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
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()
--------------------------------------------------------------------------------
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
--------------------------------------------------------------------------------
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 |
--------------------------------------------------------------------------------
gal_hadronic.eval()
p=gal_hadronic.plot_model(frame='src')
p.setlim(y_min=1E33)
gal_hadronic.energetic_report()
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
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()
--------------------------------------------------------------------------------
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
--------------------------------------------------------------------------------
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 |
--------------------------------------------------------------------------------
gal_hadronic.eval()
p=gal_hadronic.plot_model(frame='src')
p.setlim(y_min=1E36)
gal_hadronic.energetic_report()
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 |
gal_hadronic.save_model('gal_hadronic_beamed.pkl')
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 |
gal_hadronic_new=GalacticBeamed.load_model('gal_hadronic_beamed.pkl')
==> 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
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 |