def classify_clusters_given_diploid(mus, diploidInd):
	"""
	Classifies clusters given that one cluster has been classified as diploid.

	Arguments:
		mus (list of lists of floats): List of cluster means
		diploidInd (int): Index of the cluster that has been classified as the diploid cluster.

	Returns:
		singleCopyParamInds (list of ints): List of indices of clusters that have been classified as
										single copy states.
		zeroCopyParamInds (list of ints): List of indices of clusters that have been classified as
										zero copy states.
		ampParamInds (list of ints): List of indices of clusters that have been classified as
									 amplifications.
	"""

	diploidMuX = mus[diploidInd][0]
	diploidMuY = mus[diploidInd][1]

	delParamInds = []
	ampParamInds = []
	for i in range(len(mus)):
		if i == diploidInd: continue

		#deletions have lower RDR than the diploid; amplifications have greater RDR than the diploid
		if mus[i][0] < diploidMuX:
			delParamInds.append(i)
		else:
			ampParamInds.append(i)

	singleCopyParamInds = []
	zeroCopyParamInds = []
	for i in delParamInds:
		muX = mus[i][0]
		muY = mus[i][1]

		#condition for deciding if a deletion is a zero copy state (note: needs to be revised)
		if muX < diploidMuX - 0.2 and muY < diploidMuY + 0.1:
			zeroCopyParamInds.append(i)
		else:
			singleCopyParamInds.append(i)

	return singleCopyParamInds, zeroCopyParamInds, ampParamInds