code: def rho2sigma(self, rho0, Ra, Rs):
          
          return np.pi * rho0 * Ra * Rs / (Rs + Ra)
docs: converts 3d density into 2d projected density parameter
          :param rho0:
          :param Ra:
          :param Rs:
          :return: