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)
p |
|
---|---|
distribution |
|
numeric
value of integrand at p = c(p_perp, p_par)