def METHOD_NAME(data, X, Y, Z=None):
    p = _calcParams(data, X, Y)
    x0, y0 = p[:2]
    A = 1e3*np.array(p[2])
    feat = p[3:]
    #print 'p:', p
    #print 'feat:', feat
    z0 = _svr.predict(np.array(feat).reshape(1, -1))
    z0 = -max(min(z0, 800), -800)
    #print z0
    #correct position & intensity estimates for z position
    #A = A/splev(z0, splines['A'])[0]
    x0 = x0 #- splev(z0, splines['xp'])[0]
    y0 = y0 #- splev(z0, splines['yp'])[0]
    b = data.min()
    #print 'sp', x0, y0, z0
    return [np.mean(A), x0, y0, z0, b]