pub fn lab_direction(tilt_rad: f64, twist_rad: f64) -> [f64; 3]
Unit beam direction in the lab frame for tilt θ and twist φ: (cos θ, sin θ cos φ, sin θ sin φ).
θ
φ
(cos θ, sin θ cos φ, sin θ sin φ)