def _create_plc(gamma_grid_size,log_values,normalize,skip_build,emitters_type):
n_e_plc = EmittersDistribution(name='plc',
spectral_type='plc',
normalize=normalize,
emitters_type=emitters_type,
log_values=log_values,
skip_build=skip_build,
gamma_grid_size=gamma_grid_size)
a_t, b_t = n_e_plc.set_bounds(1, 1E9, log_val=n_e_plc._log_values)
n_e_plc.add_par('gamma_cut', par_type='turn-over-energy', val=1E4, vmin=a_t, vmax=b_t,
unit='lorentz-factor',log=log_values)
n_e_plc.add_par('p', par_type='LE_spectral_slope', val=2.0, vmin=-10., vmax=10, unit='')
n_e_plc.set_distr_func(distr_func_plc)
return n_e_plc
def _create_pl(gamma_grid_size, log_values, normalize, skip_build, emitters_type):
n_e_pl = EmittersDistribution(name='pl',
spectral_type='pl',
normalize=normalize,
emitters_type=emitters_type,
log_values=log_values,
skip_build=skip_build,
gamma_grid_size=gamma_grid_size)
a_t, b_t = n_e_pl.set_bounds(1, 1E9, log_val=n_e_pl._log_values)
n_e_pl.add_par('p', par_type='LE_spectral_slope', val=2.0, vmin=-10., vmax=10, unit='')
n_e_pl.set_distr_func(distr_func_pl)
return n_e_pl
def _create_lpep(gamma_grid_size, log_values, normalize, skip_build, emitters_type):
n_lep = EmittersDistribution(name='lpep',
spectral_type='lp',
normalize=normalize,
emitters_type=emitters_type,
log_values=log_values,
skip_build=skip_build,
gamma_grid_size=gamma_grid_size)
a_t, b_t = n_lep.set_bounds(1, 1E9, log_val=n_lep._log_values)
n_lep.add_par('gamma_p', par_type='turn-over-energy', val=1E4, vmin=a_t, vmax=b_t,
unit='lorentz-factor',log=log_values)
n_lep.add_par('r', par_type='spectral_curvature', val=1.0, vmin=-15., vmax=15., unit='')
n_lep.set_distr_func(distr_func_lep)
return n_lep