JetTimeEvol

JetTimeEvol#

class jetset.jet_timedep.JetTimeEvol(jet_rad, only_radiation=False, Q_inj=None, name='jet_time_ev', inplace=True, log_sampling=False, jet_gamma_grid_size=200, setup=True)[source]#

Bases: object

Parameters:

  • jet

  • Q_inj

  • name

  • inplace

  • jet_gamma_grid_size

  • log_sampling

The model parameters can be set after object creation, eg:

from jetset.jet_emitters_factory import EmittersFactory
q_inj=InjEmittersFactory().create_inj_emitters('pl',emitters_type='electrons',normalize=True)
q_inj.parameters.gmax.val=5


from jetset.jet_timedep import  JetTimeEvol
from jetset.jet_model import Jet

my_jet=Jet()

temp_ev_acc=JetTimeEvol(jet=my_jet_acc,Q_inj=q_inj,inplace=True)

temp_ev_acc.parameters.duration.val=1E4

model parameters#

duration: float (s)

duration of the evolution in seconds

gmin_grid: float

lower bound of the gamma grid for the temporal evolution

gmax_grid: float

upper bound of the gamma grid for the temporal evolution

gamma_grid_size: int

size of the gamma grid

TStart_Acc: float (s)

starting time of the acceleration process

TStop_Acc: float (s)

stop time of the acceleration process

TStart_Inj: float (s)

starting time of the injection process

TStop_Inj: float (s)

stop time of the injection process

T_esc: float (R/c)

escape time scale in (R/c)

Esc_Index: float

power-law index for the escapce term t_esc(gamma)=T_esc*gamma^(Esc_Index)

t_D0: float (s)

acceleration time scale for diffusive acceleration term

t_A0: float (s)

acceleration time scale for the systematic acceleration term

Diff_Index: float (s)

power-law index for the diffusive acceleration term

Acc_Index: float (s)

power-law index for the systematic acceleration term

Lambda_max_Turb: float (cm)

max scale of the turbulence

Lambda_choer_Turb_factor: float (cm)

max scale for the quasi-linear resonant regime

t_size: int

number of steps for the time grid

num_samples: int

number of the time steps saved in the output

L_inj: float (erg/s)

if not None, the inj function is rescaled to match L_inj

Attributes Summary

Adiabatic_cooling

Delta_R_acc

IC_cooling

Sync_cooling

custom_acc_profile

custom_q_jnj_profile

delta_t

blob delta t in s

log_sampling

logarithmically spaced bool

region_expansion

t_unit_acc

blob time in R/c

t_unit_rad

blob time in R/c

temp_ev

tempev_table

time_steps_array

Methods Summary

V_acc()

V_rad()

eval_L_tot_inj()

get_region(region)

init_TempEv()

load_model(file_name)

plot_pre_run_plot([figsize, dpi])

plot_tempev_emitters([region, figsize, dpi, ...])

plot_tempev_model([comp, region, frame, t1, ...])

plot_time_profile([figsize, dpi])

run([only_injection, do_injection, ...])

param only_injection:

save_model(file_name)

set_time([time_slice, time, frame])

show_model([getstring, names_list, sort_key])

show_pars([sort_key])

Attributes Documentation

Adiabatic_cooling#
Delta_R_acc#
IC_cooling#
Sync_cooling#
custom_acc_profile#
custom_q_jnj_profile#
delta_t#

blob delta t in s

log_sampling#

logarithmically spaced bool

region_expansion#
t_unit_acc#

blob time in R/c

t_unit_rad#

blob time in R/c

temp_ev#
tempev_table#
time_steps_array#

Methods Documentation

V_acc()[source]#
V_rad()[source]#
eval_L_tot_inj()[source]#
get_region(region)[source]#
init_TempEv()[source]#
classmethod load_model(file_name)[source]#
plot_pre_run_plot(figsize=(8, 6), dpi=120)[source]#
plot_tempev_emitters(region='rad', figsize=(8, 8), dpi=120, energy_unit='gamma', loglog=True, plot_Q_inj=True, pow=None)[source]#
plot_tempev_model(comp='Sum', region='rad', frame='obs', t1=None, t2=None, time_slice=None, time_slice_bin=None, time=None, time_bin=None, density=False, use_cached=False, sed_data=None, plot_obj=None, average=False)[source]#
plot_time_profile(figsize=(8, 8), dpi=120)[source]#
run(only_injection=True, do_injection=True, cache_SEDs_rad=True, cache_SEDs_acc=False)[source]#
Parameters:

  • only_injection (if TRUE radiative regions is considered empty and is filled by particles escaped from acceleration region,)

  • region. (otherwise particle escaped from acc regions are injected and mixed with pre-existing particles in the radiative)

  • object (in the radiative region are cooled. if Q_inj is not provided in the)

  • FALSE (then do_injection will be considered)

  • do_injection (if TRUE, particles are injected in the acceleration region, otherwise acc region is skipped, and only particles)

  • object

  • FALSE

  • cache_SEDs_rad (if TRUE, all the SEDs of the radiative region will be cached)

  • cache_SEDs_acc (if TRUE, all the SEDs of the acceleration region will be cached)

save_model(file_name)[source]#
set_time(time_slice=None, time=None, frame='blob')[source]#
show_model(getstring=False, names_list=None, sort_key=None)[source]#
show_pars(sort_key='par type')[source]#
Contents