Tutorial 1 for JetSeT¶

In [1]:

from IPython.core.display import display, HTML
display(HTML("<style>.container { width:95% !important; }</style>"))

import numpy as np

Basic setup and access to Jet class¶

In [2]:

import jetset

In [3]:

print(jetset.__version__)

1.2.0

See for more details:

https://jetset.readthedocs.io/en/latest/user_guide/jet_model_phys_SSC/Jet_example_phys_SSC.html

In [4]:

from jetset.jet_model import Jet
my_jet=Jet(electron_distribution='lppl')

In [5]:

Jet.available_electron_distributions()

lp: log-parabola
pl: powerlaw
lppl: log-parabola with low-energy powerlaw branch
lpep: log-parabola defined by peak energy
plc: powerlaw with cut-off
bkn: broken powerlaw
superexp: powerlaw with super-exp cut-off

In [6]:

my_jet.show_pars()

Table length=11

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
jet_leptonic	R	region_size	cm	5.000000e+15	1.000000e+03	1.000000e+30	False	False
jet_leptonic	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
jet_leptonic	B	magnetic_field	gauss	1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	beam_obj	beaming	lorentz-factor*	1.000000e+01	1.000000e-04	--	False	False
jet_leptonic	z_cosm	redshift		1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
jet_leptonic	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
jet_leptonic	N	emitters_density	1 / cm3	1.000000e+02	0.000000e+00	--	False	False
jet_leptonic	gamma0_log_parab	turn-over-energy	lorentz-factor*	1.000000e+04	1.000000e+00	1.000000e+09	False	False
jet_leptonic	s	LE_spectral_slope		2.000000e+00	-1.000000e+01	1.000000e+01	False	False
jet_leptonic	r	spectral_curvature		4.000000e-01	-1.500000e+01	1.500000e+01	False	False

In [7]:

my_jet.parameters.par_table

Out[7]:

Table length=11

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
str12	str16	str19	object	float64	float64	float64	bool	bool
jet_leptonic	R	region_size	cm	5.000000e+15	1.000000e+03	1.000000e+30	False	False
jet_leptonic	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
jet_leptonic	B	magnetic_field	gauss	1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	beam_obj	beaming	lorentz-factor*	1.000000e+01	1.000000e-04	--	False	False
jet_leptonic	z_cosm	redshift		1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
jet_leptonic	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
jet_leptonic	N	emitters_density	1 / cm3	1.000000e+02	0.000000e+00	--	False	False
jet_leptonic	gamma0_log_parab	turn-over-energy	lorentz-factor*	1.000000e+04	1.000000e+00	1.000000e+09	False	False
jet_leptonic	s	LE_spectral_slope		2.000000e+00	-1.000000e+01	1.000000e+01	False	False
jet_leptonic	r	spectral_curvature		4.000000e-01	-1.500000e+01	1.500000e+01	False	False

In [8]:

my_jet.show_model()

--------------------------------------------------------------------------------
jet model description
--------------------------------------------------------------------------------
name: jet_leptonic  

electrons distribution:
 type: lppl  
 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
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
jet_leptonic	R	region_size	cm	5.000000e+15	1.000000e+03	1.000000e+30	False	False
jet_leptonic	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
jet_leptonic	B	magnetic_field	gauss	1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	beam_obj	beaming	lorentz-factor*	1.000000e+01	1.000000e-04	--	False	False
jet_leptonic	z_cosm	redshift		1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
jet_leptonic	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
jet_leptonic	N	emitters_density	1 / cm3	1.000000e+02	0.000000e+00	--	False	False
jet_leptonic	gamma0_log_parab	turn-over-energy	lorentz-factor*	1.000000e+04	1.000000e+00	1.000000e+09	False	False
jet_leptonic	s	LE_spectral_slope		2.000000e+00	-1.000000e+01	1.000000e+01	False	False
jet_leptonic	r	spectral_curvature		4.000000e-01	-1.500000e+01	1.500000e+01	False	False

--------------------------------------------------------------------------------

In [9]:

my_jet.set_par('B',val=0.2)
my_jet.set_par('gamma0_log_parab',val=5E3)
my_jet.set_par('gmin',val=1E2)
my_jet.set_par('gmax',val=1E8)
my_jet.set_par('R',val=1E15)
my_jet.set_par('N',val=1E3)

In [10]:

my_jet.parameters.B.val=0.2
my_jet.parameters.r.val=0.4

In [11]:

my_jet.show_electron_distribution()

--------------------------------------------------------------------------------
electrons distribution:
 type: lppl  
 gamma energy grid size:  201
 gmin grid : 2.000000e+00
 gmax grid : 1.000000e+06
 normalization  True
 log-values  False

Table length=11

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
jet_leptonic	B	magnetic_field	gauss	2.000000e-01	0.000000e+00	--	False	False
jet_leptonic	N	emitters_density	1 / cm3	1.000000e+03	0.000000e+00	--	False	False
jet_leptonic	R	region_size	cm	1.000000e+15	1.000000e+03	1.000000e+30	False	False
jet_leptonic	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
jet_leptonic	beam_obj	beaming	lorentz-factor*	1.000000e+01	1.000000e-04	--	False	False
jet_leptonic	gamma0_log_parab	turn-over-energy	lorentz-factor*	5.000000e+03	1.000000e+00	1.000000e+09	False	False
jet_leptonic	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+08	1.000000e+00	1.000000e+15	False	False
jet_leptonic	gmin	low-energy-cut-off	lorentz-factor*	1.000000e+02	1.000000e+00	1.000000e+09	False	False
jet_leptonic	r	spectral_curvature		4.000000e-01	-1.500000e+01	1.500000e+01	False	False
jet_leptonic	s	LE_spectral_slope		2.000000e+00	-1.000000e+01	1.000000e+01	False	False
jet_leptonic	z_cosm	redshift		1.000000e-01	0.000000e+00	--	False	False

In [12]:

p=my_jet.electron_distribution.plot()

In [13]:

p=my_jet.electron_distribution.plot(energy_unit='TeV')

In [14]:

p=my_jet.electron_distribution.plot3p()

In [15]:

my_jet.eval()

In [16]:

from jetset.plot_sedfit import PlotSED
my_plot=PlotSED()
my_plot=my_jet.plot_model(plot_obj=my_plot)

In [17]:

my_plot=my_jet.plot_model(frame='src')
my_plot.setlim(y_max=1E45,y_min=1E38,x_min=1E8)

In [18]:

my_plot=my_jet.plot_model(frame='src',density=True)

In [19]:

my_jet.list_spectral_components()

Sum
Sync
SSC

In [20]:

Sync=my_jet.spectral_components.Sync

In [21]:

Sync=my_jet.get_spectral_component_by_name('Sync')

In [22]:

nu_sync=Sync.SED.nu
nuFnu_sync=Sync.SED.nuFnu

In [23]:

nu_sync_src=Sync.SED.nu_src
nuLnu_sync_src=Sync.SED.nuLnu_src

In [24]:

my_jet.spectral_components.build_table(restframe='obs')
t_obs=my_jet.spectral_components.table

In [25]:

t_obs[::10]

Out[25]:

Table length=50

nu	Sum	Sync	SSC
Hz	erg / (cm2 s)	erg / (cm2 s)	erg / (cm2 s)
float64	float64	float64	float64
1000000.0	0.0	0.0	0.0
3026648.059395689	0.0	0.0	0.0
9160598.47544371	1.8603797396558676e-27	1.8603797396558676e-27	0.0
27725907.59860481	5.1595363417357764e-26	5.1595363417357764e-26	0.0
83916564.42830162	1.431379827478216e-24	1.431379827478216e-24	0.0
253985906.87807292	3.9735563564626905e-23	3.9735563564626905e-23	0.0
768725952.1663721	1.1046717952023323e-21	1.1046717952023323e-21	0.0
2326662911.331458	3.080805682170684e-20	3.080805682170684e-20	0.0
7041989785.449296	8.644485946297075e-19	8.644485946297075e-19	0.0
...	...	...	...
1.7317171337233599e+25	8.215170127233826e-14	0.0	8.215170127233826e-14
5.2412983022060615e+25	5.0064656764602724e-14	0.0	5.0064656764602724e-14
1.5863565335085865e+26	2.2032108833042557e-14	0.0	2.2032108833042557e-14
4.801342923653465e+26	6.823069626704561e-15	0.0	6.823069626704561e-15
1.4531975242368953e+27	1.4136805292714239e-15	0.0	1.4136805292714239e-15
4.3983174666502106e+27	1.9532706759307126e-16	0.0	1.9532706759307126e-16
1.3312159025043105e+28	1.7769601809721487e-17	0.0	1.7769601809721487e-17
4.029122027951344e+28	9.162225441406983e-19	0.0	9.162225441406983e-19
1.2194734366967333e+29	0.0	0.0	0.0
3.690916910662782e+29	0.0	0.0	0.0

In [26]:

my_jet.spectral_components.build_table(restframe='src')
t_src=my_jet.spectral_components.table

In [27]:

t_obs['Sync'][::10].to('GeV/cm2 s')

Out[27]:

$[0,~0,~1.1611577 \times 10^{-24},~3.2203293 \times 10^{-23},~8.9339702 \times 10^{-22},~2.4800988 \times 10^{-20},~6.894819 \times 10^{-19},~1.9228877 \times 10^{-17},~5.3954637 \times 10^{-16},~1.439719 \times 10^{-14},~1.4724294 \times 10^{-13},~4.3081138 \times 10^{-13},~7.9304756 \times 10^{-13},~1.3786969 \times 10^{-12},~2.388583 \times 10^{-12},~4.1130803 \times 10^{-12},~6.9784192 \times 10^{-12},~1.1471237 \times 10^{-11},~1.7776511 \times 10^{-11},~2.5247891 \times 10^{-11},~3.2686208 \times 10^{-11},~3.8509885 \times 10^{-11},~4.1105436 \times 10^{-11},~3.9867003 \times 10^{-11},~3.5113825 \times 10^{-11},~2.7967606 \times 10^{-11},~2.0195405 \times 10^{-11},~1.3222962 \times 10^{-11},~7.8224764 \times 10^{-12},~4.1874255 \times 10^{-12},~2.0302231 \times 10^{-12},~8.8898466 \times 10^{-13},~3.5162291 \times 10^{-13},~1.256826 \times 10^{-13},~3.9983064 \times 10^{-14},~1.0001713 \times 10^{-14},~1.0060628 \times 10^{-15},~3.4708297 \times 10^{-18},~1.592186 \times 10^{-25},~0,~0,~0,~0,~0,~0,~0,~0,~0,~0,~0] \; \mathrm{\frac{GeV}{s\,cm^{2}}}$

In [28]:

t_src.write('test_SED.txt',format='ascii.ecsv',overwrite=True)

In [29]:

my_jet.energetic_report()

Table length=37

name	type	units	val
BulkLorentzFactor			1.000000e+01
U_e	Energy dens. blob rest. frame	erg / cm3	4.950261e-01
U_p_cold	Energy dens. blob rest. frame	erg / cm3	1.503276e-01
U_B	Energy dens. blob rest. frame	erg / cm3	1.591549e-03
U_p	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_p_target	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_Synch	Energy dens. blob rest. frame	erg / cm3	1.316122e-02
U_Synch_DRF	Energy dens. disk rest. frame	erg / cm3	1.316122e+02
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. frme.	erg / s	1.652654e+39
L_SSC_rf	Lum. blob rest. frme.	erg / s	2.158648e+39
L_EC_Disk_rf	Lum. blob rest. frme.	erg / s	0.000000e+00
L_EC_BLR_rf	Lum. blob rest. frme.	erg / s	0.000000e+00
L_EC_DT_rf	Lum. blob rest. frme.	erg / s	0.000000e+00
L_EC_CMB_rf	Lum. blob rest. frme.	erg / s	0.000000e+00
L_pp_gamma_rf	Lum. blob rest. frme.	erg / s	0.000000e+00
jet_L_Sync	jet Lum.	erg / s	4.131634e+40
jet_L_SSC	jet Lum.	erg / s	5.396619e+40
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	0.000000e+00
jet_L_rad	jet Lum.	erg / s	9.528254e+40
jet_L_kin	jet Lum.	erg / s	6.047641e+42
jet_L_tot	jet Lum.	erg / s	6.157838e+42
jet_L_e	jet Lum.	erg / s	4.638913e+42
jet_L_B	jet Lum.	erg / s	1.491449e+40
jet_L_p_cold	jet Lum.	erg / s	1.408728e+42
jet_L_p	jet Lum.	erg / s	0.000000e+00

In [30]:

my_jet.energetic_report_table

Out[30]:

Table length=37

name	type	units	val
str17	str29	object	float64
BulkLorentzFactor			1.000000e+01
U_e	Energy dens. blob rest. frame	erg / cm3	4.950261e-01
U_p_cold	Energy dens. blob rest. frame	erg / cm3	1.503276e-01
U_B	Energy dens. blob rest. frame	erg / cm3	1.591549e-03
U_p	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_p_target	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_Synch	Energy dens. blob rest. frame	erg / cm3	1.316122e-02
U_Synch_DRF	Energy dens. disk rest. frame	erg / cm3	1.316122e+02
U_Disk	Energy dens. blob rest. frame	erg / cm3	0.000000e+00
U_BLR	Energy dens. blob rest. frame	erg / cm3	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	0.000000e+00
jet_L_rad	jet Lum.	erg / s	9.528254e+40
jet_L_kin	jet Lum.	erg / s	6.047641e+42
jet_L_tot	jet Lum.	erg / s	6.157838e+42
jet_L_e	jet Lum.	erg / s	4.638913e+42
jet_L_B	jet Lum.	erg / s	1.491449e+40
jet_L_p_cold	jet Lum.	erg / s	1.408728e+42
jet_L_p	jet Lum.	erg / s	0.000000e+00

In [31]:

my_jet.save_model('test_model.pkl')

In [32]:

my_jet_new=Jet.load_model('test_model.pkl')

Table length=11

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
jet_leptonic	gmin	low-energy-cut-off	lorentz-factor*	1.000000e+02	1.000000e+00	1.000000e+09	False	False
jet_leptonic	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+08	1.000000e+00	1.000000e+15	False	False
jet_leptonic	N	emitters_density	1 / cm3	1.000000e+03	0.000000e+00	--	False	False
jet_leptonic	gamma0_log_parab	turn-over-energy	lorentz-factor*	5.000000e+03	1.000000e+00	1.000000e+09	False	False
jet_leptonic	s	LE_spectral_slope		2.000000e+00	-1.000000e+01	1.000000e+01	False	False
jet_leptonic	r	spectral_curvature		4.000000e-01	-1.500000e+01	1.500000e+01	False	False
jet_leptonic	R	region_size	cm	1.000000e+15	1.000000e+03	1.000000e+30	False	False
jet_leptonic	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
jet_leptonic	B	magnetic_field	gauss	2.000000e-01	0.000000e+00	--	False	False
jet_leptonic	beam_obj	beaming	lorentz-factor*	1.000000e+01	1.000000e-04	--	False	False
jet_leptonic	z_cosm	redshift		1.000000e-01	0.000000e+00	--	False	False

In [33]:

my_plot=my_jet_new.plot_model()
my_plot.setlim(y_max=1E-11,y_min=3E-17,x_min=1E8)

Define custom emiotters distribution¶

See for more details:

https://jetset.readthedocs.io/en/latest/user_guide/custom_emitters_distr/custom_emitters.html

In [34]:

from jetset.jet_emitters import EmittersDistribution

In [35]:

def distr_func_super_exp(gamma,gamma_cut,s,a):
    return np.power(gamma,-s)*np.exp(-(1/a)*(gamma/gamma_cut)**a)

In [ ]:

In [36]:

n_e_super_exp=EmittersDistribution('super_exp',spectral_type='plc',normalize=False)
n_e_super_exp.add_par('gamma_cut',par_type='turn-over-energy',val=50000.,vmin=1., vmax=None, unit='lorentz-factor')
n_e_super_exp.add_par('s',par_type='LE_spectral_slope',val=2.3,vmin=-10., vmax=10, unit='')
n_e_super_exp.add_par('a',par_type='spectral_curvature',val=1.8,vmin=0., vmax=100., unit='')

In [37]:

n_e_super_exp.set_distr_func(distr_func_super_exp)

In [38]:

n_e_super_exp.parameters.show_pars()

Table length=6

name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
N	emitters_density	1 / cm3	1.000000e+02	0.000000e+00	--	False	False
gamma_cut	turn-over-energy	lorentz-factor*	5.000000e+04	1.000000e+00	--	False	False
s	LE_spectral_slope		2.300000e+00	-1.000000e+01	1.000000e+01	False	False
a	spectral_curvature		1.800000e+00	0.000000e+00	1.000000e+02	False	False

In [39]:

p=n_e_super_exp.plot()

In [40]:

p=n_e_super_exp.plot(energy_unit='eV')

In [41]:

from jetset.jet_model import Jet
my_jet=Jet(electron_distribution=n_e_super_exp)

In [42]:

n_e_super_exp.normalize

Out[42]:

False

In [43]:

my_jet.electron_distribution.normalize

Out[43]:

False

In [44]:

my_jet.parameters.N.val=5E4
my_jet.show_model()
my_jet.IC_nu_size=100
my_jet.eval()

--------------------------------------------------------------------------------
jet model description
--------------------------------------------------------------------------------
name: jet_leptonic  

electrons distribution:
 type: super_exp  
 gamma energy grid size:  201
 gmin grid : 2.000000e+00
 gmax grid : 1.000000e+06
 normalization  False
 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
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
jet_leptonic	R	region_size	cm	5.000000e+15	1.000000e+03	1.000000e+30	False	False
jet_leptonic	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
jet_leptonic	B	magnetic_field	gauss	1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	beam_obj	beaming	lorentz-factor*	1.000000e+01	1.000000e-04	--	False	False
jet_leptonic	z_cosm	redshift		1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
jet_leptonic	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
jet_leptonic	N	emitters_density	1 / cm3	5.000000e+04	0.000000e+00	--	False	False
jet_leptonic	gamma_cut	turn-over-energy	lorentz-factor*	5.000000e+04	1.000000e+00	--	False	False
jet_leptonic	s	LE_spectral_slope		2.300000e+00	-1.000000e+01	1.000000e+01	False	False
jet_leptonic	a	spectral_curvature		1.800000e+00	0.000000e+00	1.000000e+02	False	False

--------------------------------------------------------------------------------

In [45]:

my_jet.eval()

In [46]:

p=my_jet.plot_model()
p.setlim(y_min=1E-16,y_max=1E-12)

In [47]:

my_jet.electron_distribution.normalize=True

In [48]:

my_jet.parameters.N.val=5E4
my_jet.show_model()
my_jet.IC_nu_size=100
my_jet.eval()

--------------------------------------------------------------------------------
jet model description
--------------------------------------------------------------------------------
name: jet_leptonic  

electrons distribution:
 type: super_exp  
 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
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
jet_leptonic	R	region_size	cm	5.000000e+15	1.000000e+03	1.000000e+30	False	False
jet_leptonic	R_H	region_position	cm	1.000000e+17	0.000000e+00	--	False	True
jet_leptonic	B	magnetic_field	gauss	1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	beam_obj	beaming	lorentz-factor*	1.000000e+01	1.000000e-04	--	False	False
jet_leptonic	z_cosm	redshift		1.000000e-01	0.000000e+00	--	False	False
jet_leptonic	gmin	low-energy-cut-off	lorentz-factor*	2.000000e+00	1.000000e+00	1.000000e+09	False	False
jet_leptonic	gmax	high-energy-cut-off	lorentz-factor*	1.000000e+06	1.000000e+00	1.000000e+15	False	False
jet_leptonic	N	emitters_density	1 / cm3	5.000000e+04	0.000000e+00	--	False	False
jet_leptonic	gamma_cut	turn-over-energy	lorentz-factor*	5.000000e+04	1.000000e+00	--	False	False
jet_leptonic	s	LE_spectral_slope		2.300000e+00	-1.000000e+01	1.000000e+01	False	False
jet_leptonic	a	spectral_curvature		1.800000e+00	0.000000e+00	1.000000e+02	False	False

--------------------------------------------------------------------------------

In [49]:

my_jet.plot_model(p,comp='Sum',label='Normalized distr')
p.setlim(y_min=1E-16,y_max=1E-12)

In [50]:

p.fig

Out[50]:

In [ ]: