pub fn gamma(m1: f64, m2: f64) -> f64
Kinematic factor gamma = 4 M1 M2 / (M1 + M2)^2.
gamma = 4 M1 M2 / (M1 + M2)^2