function to calculate the integrand in an integral of a distribution over all momentum space. The integral is carried out in cylindrical coordinates, so the volume element in the integral is p_perp * dp_perp * dp_par. Therefore the distribution is here multiplied by p_perp. Also note the input p must be of order 1 and is scaled in the function by distribution[["p_scale"]].

integrand_distribution(p, distribution)

Arguments

p

numeric momentum vector: p = c(p_perp, p_par). Of order 1.

distribution

list with the momentum distribution.

Value

numeric value of integrand at p = c(p_perp, p_par)