Co-moving volume density of our HLIRGs in different redshift ranges compared with global stellar mass function measurements from the literature based on observations (red, green, yellow, orange, and purple lines as well as black pentagons) or simulations (black solid line). Error bars combine the Poisson error and the standard deviation among the three deep fields. Empty markers represent a subsample with good quality in the photo-z estimates and SED fits.