mask[np.where(theta < 2*math.pi/3)] = 0