tardis.montecarlo.estimators.radfield_mc_estimators module

class tardis.montecarlo.estimators.radfield_mc_estimators.RadiationFieldMCEstimators(*args, **kwargs)

Bases: RadiationFieldMCEstimators

class_type = jitclass.RadiationFieldMCEstimators#7ff7a9e23ed0<j_estimator:array(float64, 1d, A),nu_bar_estimator:array(float64, 1d, A),j_blue_estimator:array(float64, 2d, A),Edotlu_estimator:array(float64, 2d, A),photo_ion_estimator:array(float64, 2d, A),stim_recomb_estimator:array(float64, 2d, A),bf_heating_estimator:array(float64, 2d, A),stim_recomb_cooling_estimator:array(float64, 2d, A),photo_ion_estimator_statistics:array(int64, 2d, A)>

tardis.montecarlo.estimators.radfield_mc_estimators.initialize_estimator_statistics(tau_sobolev_shape, gamma_shape)

Initializes the estimators used in the Monte Carlo simulation.

Parameters

tau_sobolev_shapetuple

Shape of the array with the Sobolev optical depth.

gamma_shapetuple

Shape of the array with the photoionization rate coefficients.

Returns

Estimators

The initialized estimators.

Examples

>>> tau_sobolev_shape = (10, 20)
>>> gamma_shape = (5, 5)
>>> initialize_estimators(tau_sobolev_shape, gamma_shape)
<Estimators object at 0x...>

tardis.montecarlo.estimators.radfield_mc_estimators.update_base_estimators(r_packet, distance, estimator_state, comov_nu, comov_energy)

Updating the estimators

tardis.montecarlo.estimators.radfield_mc_estimators.update_bound_free_estimators(comov_nu, comov_energy, shell_id, distance, estimator_state, t_electron, x_sect_bfs, current_continua, bf_threshold_list_nu)

Update the estimators for bound-free processes.

Parameters

comov_nufloat

comov_energyfloat

shell_idint

distancefloat

numba_estimatortardis.montecarlo.montecarlo_numba.numba_interface.Estimators

t_electronfloat

Electron temperature in the current cell.

x_sect_bfsnumpy.ndarray, dtype float

Photoionization cross-sections of all bound-free continua for which absorption is possible for frequency comov_nu.

current_continuanumpy.ndarray, dtype int

Continuum ids for which absorption is possible for frequency comov_nu.

bf_threshold_list_nunumpy.ndarray, dtype float

Threshold frequencies for photoionization sorted by decreasing frequency.

tardis.montecarlo.estimators.radfield_mc_estimators.update_line_estimators(radfield_mc_estimators, r_packet, cur_line_id, distance_trace, time_explosion)

Function to update the line estimators

Parameters

estimatorstardis.montecarlo.montecarlo_numba.numba_interface.Estimators

r_packettardis.montecarlo.montecarlo_numba.r_packet.RPacket

cur_line_idint

distance_tracefloat

time_explosionfloat