def find_breakthrough(times, data):
    """return tuple of time,value for the start breakthrough and end breakthrough"""
    selected = data > 0.999 * max(data)
    selected_times = times[selected]
    return (selected_times[0], max(data)), (selected_times[-1], max(data))
