hazma.rambo.generate_phase_space¶

hazma.rambo.generate_phase_space(num_ps_pts, masses, cme, mat_elem_sqrd=<function <lambda>>, num_cpus=None)¶

Generate a specified number of phase space points given a set of final state particles and a given center of mass energy. NOTE: weights are not normalized.

Parameters:

num_ps_pts (int) – Total number of phase space points to generate. masses (numpy.ndarray) – List of masses of the final state particles. cme (double) – Center-of-mass-energy of the process. mat_elem_sqrd ((double)(numpy.ndarray) {lambda klist: 1}) – Function for the matrix element squared. num_cpus (int {None}) – Number of cpus to use in parallel with rambo. If not specified, 75% of the cpus will be used.

Returns:

phase_space_points (numpy.ndarray) – List of phase space points. The phase space points are in the form {{ke11, kx11, ky11, kz11, …, keN1, kxN1, kyN1, kzN1, weight1}, . {ke1N, kx1N, ky1N, kz1N, …, keNN, kxNN, kyNN, kzNN, weightN}}

Examples

Generate 100000 phase space points for a 3 body final state.

>>> from hazma import rambo
>>> import numpy as np
>>> masses = np.array([100., 200., 0.0])
>>> cme = 1000.
>>> num_ps_pts = 100000
>>> num_fsp = len(masses)
>>>
>>> pts = rambo.generate_phase_space(num_ps_pts, masses, cme)