def find_rt(peak_data):
    """Determine retention time for a peak, given Pandas dataframe of rt and intensity"""
    max_intensity = peak_data["intensity"].max()
    rt = round(peak_data[peak_data["intensity"] == max_intensity].iloc[0]["rt"], 4)

    return rt