unit = "SATUN06"
bpmlist = ("SATUN01-DBPM070","SATUN02-DBPM070", "SATUN03-DBPM070", "SATUN04-DBPM070", "SATUN05-DBPM410", "SATUN06-DBPM070", "SATUN07-DBPM070", "SATUN08-DBPM070", "SATUN09-DBPM070", "SATUN10-DBPM070", "SATUN11-DBPM070", "SATUN12-DBPM070", "SATUN13-DBPM070", "SATUN14-DBPM410", "SATUN15-DBPM070", "SATUN16-DBPM070", "SATUN17-DBPM070", "SATUN18-DBPM070", "SATUN19-DBPM070", "SATUN20-DBPM070", "SATUN21-DBPM070", "SATUN22-DBPM070", "SATMA02-DBPM030", "SATMA02-DBPM060", "SATBD01-DBPM020", "SATBD01-DBPM060", "SATBD01-DBPM100")
bpmlist2 = ("SATUN01-DBPM070","SATUN02-DBPM070", "SATUN03-DBPM070", "SATUN04-DBPM070", "SATUN05-DBPM410", "SATUN06-DBPM070", "SATUN07-DBPM070", "SATUN08-DBPM070", "SATUN09-DBPM070", "SATUN10-DBPM070")

def read(bpmlist):
    X, Y = [], []
    for bpm in bpmlist:
        X.append(caget(bpm + ":X2"))
        Y.append(caget(bpm + ":Y2"))
    return X, Y

def average(bpmlist):
    X,Y = [0 for i in range(27)], [0 for i in range(27)]
    for i in range(10):
        x, y = read(bpmlist)
        X = arradd(X,x)
        Y = arradd(Y,y)
        print(i)
    for val in X: val = val / 10
    for val in Y: val = val / 10
    return X, Y

X1, Y1 = average(bpmlist)

print(X1, Y1)
res = Lscan