def METHOD_NAME():
    data = np.ones((50, 50))
    aper = CircularAperture((25, 25), r=10.0)
    mask = aper.to_mask()
    cutout = mask.cutout(data, copy=True)
    data[25, 25] = 100.0
    assert cutout[10, 10] == 1.0
    # test quantity data
    data2 = np.ones((50, 50)) * u.adu
    cutout2 = mask.cutout(data2, copy=True)
    assert cutout2.unit == data2.unit
    data2[25, 25] = 100.0 * u.adu
    assert cutout2[10, 10].value == 1.0