Results of three different collisional models: main belt alone which is discussed in Sect. 3 (left column), main belt and comets from Sect. 4 (middle column), main belt and disrupting comets from Sect. 7 (right column). 1st row: the initial and evolved SFDs of the main belt populations for 100 Boulder simulations; 2nd row: the resulting family production functions (in order to distinguish 100 lines we plot them using different colours ranging from black to yellow) and their comparison to the observations; 3rd row: the production function affected by comminution for a selected simulation; and 4th row: the distribution of synthetic families with D_PB ≥ 50 km in the (age, D_PB) plot for a selected simulation, without comminution. The positions of synthetic families in the 4th-row figures may differ significantly for a different Boulder simulation due to stochasticity and low-number statistics. Moreover, in the middle and right columns, many families were created during the LHB, so there are many overlapping crosses close to 4 Gyr.