Tutorial 1 for JetSeT¶

In [18]:

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

import numpy as np

Basic setup and access to Jet class¶

In [19]:

import jetset

In [20]:

print(jetset.__version__)

1.3.0rc7

See for more details:

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

Starting from version 1.3.0 the Jet class from the jet_model module, and all the derived classes, perform the C computation using threads. This increase the computational speed Each time you create a new Jet object, you will get a log noticing how many C threads have been created The number of threads is automatically determined according to the number of cores and threads of your computation You can revert back to single thread using the set_num_c_threads and passing the number of threads:

In [22]:

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

===> setting C threads to 12

In [23]:

my_jet.set_num_c_threads(1)

===> setting C threads to 1

In [25]:

my_jet.set_num_c_threads(12)

===> setting C threads to 12

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=12

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
str12	str16	str21	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	NH_cold_to_rel_e	cold_p_to_rel_e_ratio		1.000000e+00	0.000000e+00	--	False	True
jet_leptonic	beam_obj	beaming		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=12

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
str12	str16	str21	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	NH_cold_to_rel_e	cold_p_to_rel_e_ratio		1.000000e+00	0.000000e+00	--	False	True
jet_leptonic	beam_obj	beaming		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()

--------------------------------------------------------------------------------
model description: 
--------------------------------------------------------------------------------
type: Jet
name: jet_leptonic  
geometry: spherical  

electrons distribution:
 type: lppl  
 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+00

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:Sum, hidden: False
   name:Sync, state: self-abs
   name:Sync, hidden: False
   name:SSC, state: on
   name:SSC, hidden: False
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=12

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
str12	str16	str21	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	NH_cold_to_rel_e	cold_p_to_rel_e_ratio		1.000000e+00	0.000000e+00	--	False	True
jet_leptonic	beam_obj	beaming		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=12

model name	name	par type	units	val	phys. bound. min	phys. bound. max	log	frozen
str12	str16	str21	object	float64	float64	float64	bool	bool
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	NH_cold_to_rel_e	cold_p_to_rel_e_ratio		1.000000e+00	0.000000e+00	--	False	True
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		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 [ ]: