tardis.energy_input.gamma_ray_transport module¶
- tardis.energy_input.gamma_ray_transport.calculate_average_energies(raw_isotope_abundance, gamma_ray_lines)[source]¶
Function to calculate average energies of positrons and gamma rays from a list of gamma ray lines from nndc.
- Parameters
- raw_isotope_abundancepd.DataFrame
isotope abundance in mass fractions
- gamma_ray_linespd.DataFrame
decay data
- Returns
- average_energies_listList
list of gamma ray energies
- average_positron_energies_listList
list of positron energies
- gamma_ray_line_array_listList
list of gamma ray lines
- tardis.energy_input.gamma_ray_transport.calculate_average_power_per_mass(energy_per_mass, time_delta)[source]¶
- tardis.energy_input.gamma_ray_transport.calculate_energy_per_mass(decay_energy, raw_isotope_abundance, cell_masses)[source]¶
Function to calculate decay energy per mass for each isotope chain.
- Parameters
- decay_energyDict
dictionary of decay chain energies for each isotope in each shell
- raw_isotope_abundancepd.DataFrame
isotope abundance in mass fractions.
- cell_massesnumpy.ndarray
shell masses in units of g
- Returns
- energy_per_masspd.DataFrame
decay energy per mass for each isotope chain For e.g Ni56 has 2 decay chains: Ni56 -> Co56 -> Fe56. It will calculate the decay energy per mass for each chain and store as a dataframe.
- tardis.energy_input.gamma_ray_transport.calculate_total_decays_old(inventories, time_delta)[source]¶
Function to create inventories of isotope
- Parameters
- modeltardis.Radial1DModel
The tardis model to calculate gamma ray propagation through
- time_endfloat
End time of simulation in days
- Returns
- Total decay listList
list of total decays for x g of isotope for time ‘t’
- tardis.energy_input.gamma_ray_transport.create_inventories_dict_old(isotope_dict)[source]¶
Function to create dictionary of inventories for each shell
- Parameters
- isotope_dictDict
dictionary of isotopes for each shell with their
masses.
- Returns
- invDict
dictionary of inventories for each shell For eg: {0: {‘Ni56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>,
‘Co56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>},
- {1: {‘Ni56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>,
‘Co56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>}} etc
- tardis.energy_input.gamma_ray_transport.create_isotope_dicts_old(raw_isotope_abundance, cell_masses)[source]¶
Function to create a dictionary of isotopes for each shell with their masses.
- Parameters
- raw_isotope_abundancepd.DataFrame
isotope abundance in mass fractions.
- cell_massesnumpy.ndarray
shell masses in units of g
- Returns
- isotope_dictsDict
dictionary of isotopes for each shell with their
masses. For eg: {0: {(28, 56): {‘Ni56’: 0.0001}, (27, 57): {‘Co56’: 0.0001}}{1: {(28, 56): {‘Ni56’: 0.0001}, (27, 57): {‘Co56’: 0.0001}}} etc
- tardis.energy_input.gamma_ray_transport.decay_chain_energies(average_energies, total_decays)[source]¶
Function to calculate decay chain energies.
- Parameters
- raw_isotope_abundancepd.DataFrame
isotope abundance in mass fractions
- average_energies_listList
list of gamma ray energies
- average_positron_energies_listList
list of positron energies
- gamma_ray_line_array_listList
list of gamma ray lines
- total_decaysDict
dictionary of total decays for each isotope in each shell
- Returns
- decay_energyDict
dictionary of decay chain energies for each isotope in each shell
- tardis.energy_input.gamma_ray_transport.distribute_packets(decay_energy, total_energy, num_packets)[source]¶
- tardis.energy_input.gamma_ray_transport.fractional_decay_energy(decay_energy)[source]¶
Function to calculate fractional decay energy Parameters ———- decay_energy : Dict
dictionary of decay chain energies for each isotope in each shell
Returns¶
- fractional_decay_energyDict
dictionary of fractional decay chain energies for each isotope in each shell
- tardis.energy_input.gamma_ray_transport.get_chain_decay_power_per_ejectamass(inventory, time, time_start, average_energies, average_positron_energies, tau_end)[source]¶
- tardis.energy_input.gamma_ray_transport.get_taus(raw_isotope_abundance)[source]¶
Function to calculate taus for each isotope
- Parameters
- raw_isotope_abundancepd.DataFrame
isotope abundance in mass fractions
- Returns
- tausDict
dictionary of taus for each isotope
- parentsDict
dictionary of parents for each isotope