def METHOD_NAME(x, y, center_x, center_y, r_in, r_out, phi0=0, delta_phi=2 * np.pi):
    """
    :param x:
    :param y:
    :param center_x:
    :param center_y:
    :param r_in:
    :param r_out:
    :param phi0:
    :param delta_phi:
    :return:
    :rtype: array of size of input grid with integers 0 or 1
    """
    x_shift = x - center_x
    y_shift = y - center_y
    R = np.sqrt(x_shift * x_shift + y_shift * y_shift)
    phi = np.arctan2(x_shift, y_shift)
    phi_min = phi0 - delta_phi / 2.0
    phi_max = phi0 + delta_phi / 2.0
    mask = np.zeros_like(x, dtype="int")
    if phi_max > phi_min:
        mask[(R < r_out) & (R > r_in) & (phi > phi_min) & (phi < phi_max)] = 1
    else:
        mask[(R < r_out) & (R > r_in) & (phi > phi_max)] = 1
        mask[(R < r_out) & (R > r_in) & (phi < phi_min)] = 1
    return mask