  We discuss the finite size behaviour in the canonical ensemble of the balls
in boxes model. We compare theoretical predictions and numerical results for
the finite size scaling of cumulants of the energy distribution in the
canonical ensemble and perform a detailed analysis of the first and third order
phase transitions which appear for different parameter values in the model.
