Galaxy stellar mass function up to z = 4 for the full sample. Each colour corresponds to different redshift bins of variable step size. Fits are shown in the mass range covered by our dataset. The filled areas correspond to the 68% confidence level regions, after accounting for Poissonian errors, the cosmic variance and the uncertainties created during the template fitting procedure. The open triangles and squares correspond to the local estimates by Moustakas et al. (2013) and Baldry et al. (2012), respectively.