Function solid_angle_transform 

pub fn solid_angle_transform(
    d_omega_rest: f64,
    v: f64,
    c: f64,
    theta: f64,
) -> f64

Solid angle transformation under Lorentz boost. d_omega_obs = d_omega_rest / (gamma^2 * (1 - beta*cos(theta))^2)