def METHOD_NAME(peak_profile):
    x = numpy.arange(5000)
    peak_params = _PEAK_PARAMETERS[peak_profile]
    func = getattr(functions, peak_profile)
    y = func(x, *peak_params)
    estimated_peak_params = peaks.peak_search(y=y, fwhm=100, relevance_info=True)
    assert len(estimated_peak_params)==6, "Wrong number of peaks detected"
    for i, (peak_position, *_) in enumerate(estimated_peak_params):
        theoretical_peak_position = peak_params[i*(len(peak_params)//6) + 1]
        assert abs(peak_position - theoretical_peak_position) < 25