How to manage this? Could change `peak_flux()` function to `max_flux_within_range()` which then potentially returns a list of values if the function is exactly periodic? Or: change entirely to focus on the 'is this detected?' problem