diff --git a/gammapkt.cc b/gammapkt.cc
index 2939ad8d2..379d21cc4 100644
--- a/gammapkt.cc
+++ b/gammapkt.cc
@@ -324,7 +324,7 @@ auto choose_f(const double xx, const double zrand) -> double
 auto thomson_angle() -> double {
   const double B_coeff = (8. * rng_uniform()) - 4.;
 
-  const double t_coeff = std::cbrt((std::sqrt((B_coeff * B_coeff) + 4) - B_coeff) / 2.2);
+  const double t_coeff = std::cbrt((std::sqrt((B_coeff * B_coeff) + 4) - B_coeff) / 2);
 
   const double mu = (1 / t_coeff) - t_coeff;