def find_peaks(x,num=8):
    d = x;
    #figure();
    #plot(d);

    peaks = r_[0.:num];
    i = 0;
    mmm = max(d)
    for p in range(0,int(num)):
        while (d[i] < mmm*.2):
            i += 1;
        peaks[p] = 0.;
        for n in range(-1,10):
            peaks[p] += r_[i+n]*d[i+n]
        peaks[p] = peaks[p]/sum(d[i-1:i+9])
        i += 13
    return peaks;


def fit_peaks(p):
    l=lambda a,x:a[0]*x+a[1]
    print p; 
    plotfit(fittodata_linear(r_[0.:size(p)],p,[25.,1.],err=.5*ones(size(p)),ret_all=True));

def plot_mcacalib():
    figure();
    plot(fmca_calib[0:450]);
    txt = text(0.22,.75,"delta-T/",
            horizontalalignment='center',
            verticalalignment='center',
            transform = ax.transAxes)