Vertical distribution of disk stars as a function of [Fe/H] and [Mg/Fe]. Colors according to the GMM groups as in Fig 2. The curve+points correspond to the median values of stars separated in bins equally populated. Error bands at the 1σ and 2σ level are computed by Monte Carlo resamplings to take into account the effect of errors in metallicity and abundances on the sample classification.