def phi(self):
        """Polar angle / inclination of this vector in radians

        Based on sperical coordinate space
        returns angle between this vector and the positive z-azis
        range: (0 <= phi <= pi)
        """
        return np.arctan2(np.sqrt(self.x**2 + self.y**2), self.z)