def ph_average(namef, epsilon):
    with open(namef, 'r') as f:
        lnew = []
        for prod in f:
            lprod = prod.strip().split(',')
            s = 0
            for measurement in lprod[1:]:
                s = s + float(measurement)
            m = s/(len(lprod)-1)
            if abs(m-7) < epsilon:
                lnew.append(lprod[0])
    return sorted(lnew)