def get_integrated_intensity(ave_azav, peak_bin, delta_bin=3):
    """
    Get the average integrated intensity.  Sum the bin values from
    the peak bin to delta bin on each side.

    Parameters
    ----------
    ave_azav: ndarray
        radially binned average azimuthal intensity from used curves

    peak_bin: int
        peak radial bin

    delta_bin: int
        number of bins on each side of peak to include in sum

    Returns
    -------
    integrated_intensity: float
        the calculated integrated intensity
    """
    low = peak_bin - delta_bin
    high = peak_bin + delta_bin
    integrated_intensity = ave_azav[low: high].sum(axis=0)

    return integrated_intensity