We classify the probability distribution of partial sums of spins for a n-species mean field model in the thermodinamic limit. We show under which conditions this distribution converges to a multivariate normal distribution and we give the covariance matrix that is tied to the suscetivity matrix of the model.