def combineSpectraCompleteLinkage(masses1, masses2, epsilon=0.04):
combMasses = np.append(masses1, masses2)
combMasses = np.sort(combMasses)
clusters = DNS.getClustersCompleteLinkage(combMasses, epsilon)
combMasses = []
for cluster in clusters:
if len(cluster) > 1:
combMasses += [sum(cluster) / len(cluster)]
else:
combMasses += cluster
return np.sort(np.array(combMasses))
