It took me a while to understand the formula. Write it as
  std::sqrt (global_rad..._integral / this->get_volume())