more modular, problems with assumption that mean(obs) + std(obs) cutoff makes sense.
debugging peak calling