def make_list(x, m, modelcut, errorcut, CorD):
    """Makes list in coupling correction """
    if x == []:
        return []

    result_names_list = []
    count = 0
    for i in range(len(x.NAME)):
        bn = x.NAME[i].upper()
        if bn in m.indx:
            indx = x.indx[bn]         
            if ((abs(x.Cf1001r[indx]-m.Cf1001r[indx]) < modelcut)
                        and (abs(x.Cf1001i[indx]-m.Cf1001i[indx]) < modelcut)
                        and (x.Cf1001iERR[indx] < errorcut)
                        and (x.Cf1001rERR[indx] < errorcut)):
                result_names_list.append(x.NAME[i])
        else:
            print "Not in Response:", bn
            count=count+1

    if count > 0:
        print "Warning: ", count, "BPMs removed from data for not beeing in the model"

    return result_names_list

b = beat_input.computevectorEXP(X,Y,DX) - beat_input.zerovector
...
def computevectorEXP(self, x, y, dx, xbet=[], ybet=[]):
...
betx.append(xbet.BETX[xbet.indx[b]])

[deltas, varslist] = Python_Classes4MAD.GenMatrix.correctbeatEXP(phase_x, phase_y, dx, beat_inp, cut=_InputData.singular_value_cut, app=0, path=_InputData.output_path, xbet=beta_x, ybet=beta_y)

# In correctbeatEXP()
vector = beat_input.computevectorEXP(x, y, dx, xbet, ybet)

# correct.py
Python_Classes4MAD.BCORR.bNCorrNumeric(phase_x, phase_y, dx, beat_inp, cut=_InputData.singular_value_cut, ncorr=_InputData.num_of_correctors, app=0, path=_InputData.output_path, beta_x=beta_x, beta_y=beta_y)


# BCORR.py
def bNCorrNumeric(X,Y,DX, beat_input, cut=0.001, ncorr=3, app=0, tol=1e-9,path="./", beta_x=None, beta_y=None):
  ...
  b = beat_input.computevectorEXP(X,Y,DX, beta_x, beta_y)-beat_input.zerovector