idx = np.argwhere(np.diff(np.sign(z - alphaLine))).flatten()