models.ensemble.ensemble#

This module contains a calculation of the different combinations needed to calculate the voxel-averaged emission, as well as the calculations of the probabilities to observe each combination as a function of the internal velocity \(v_i\). The procedure cycles through the available ensembles and performs the calculations separately.

Note

This module owes itself largely to the work done by Silke Andree-Labsch and Christoph Bruckmann. A clear deficiency in the calculations in this module is that we use a list of the necessary probabilities, which are of different shapes depending on the number of combinations observed for a given velocity bin. Since this computational difficulty is common in (astro) particle physics, there should exist a python package to properly handle this in python. This is a possible development that still needs to be tested and implemented.

Functions

calculate

!!!! <<DEGENERATE>>

calculate_combinations

This function calculates all of the different combinations of clump masses that may be in a line-of-sight.

create_clump_combinations

This function removes all of the unnecessary degenerate looping during this calculation.

initialise

Initialise the combinations of clumps for each ensemble and calculate the probabilities.

print_ensembleparameters

This function prints the number of clumps in the ensemble, both standard and normalised.

reinitialise

Reinitialise all subpackage variables to the correct number of clump sets.