Creates the function to define the exclusion zones

def _hammer_function(self, x,x0,r_x0, s_x0):
        '''
        Creates the function to define the exclusion zones
        '''
        return norm.logcdf((np.sqrt((np.square(np.atleast_2d(x)[:,None,:]-np.atleast_2d(x0)[None,:,:])).sum(-1))- r_x0)/s_x0)