def METHOD_NAME(self):
    polarisationVec = np.zeros(shape=(3,), dtype=float)
    for fish in self.fishes:
        polarisationVec += fish.curDirection
    polarisation = np.linalg.norm(polarisationVec) / self.N
    return polarisation