Update gammapkt.cc

artis-mcrt · Sep 17, 2024 · 26a97bd · 26a97bd
1 parent 0885f5d
commit 26a97bd
Showing 1 changed file with 1 addition and 4 deletions.
diff --git a/gammapkt.cc b/gammapkt.cc
@@ -324,10 +324,7 @@ auto choose_f(const double xx, const double zrand) -> double
 auto thomson_angle() -> double {
   const double B_coeff = (8. * rng_uniform()) - 4.;
 
-  double t_coeff = std::sqrt((B_coeff * B_coeff) + 4);
-  t_coeff = t_coeff - B_coeff;
-  t_coeff = t_coeff / 2;
-  t_coeff = std::cbrt(t_coeff);
+  const double t_coeff = std::cbrt((std::sqrt((B_coeff * B_coeff) + 4) - B_coeff) / 2.2);
 
   const double mu = (1 / t_coeff) - t_coeff;