for iterator in QGIS3 Graphical Modeler

I´d like to integrate the for loop in my model in ModelBuilder in QGIS3 (similar as the for interator in ArcGIS ModelBuilder).

data:

testdata: point layer with different classes (assignment to different classes)

hexagons: polygon layer (hexagrids)

aim:

For each class of the point layer testdata calculate the number of points per hexagrid and assign the count to the centroid of the grid where the points are situated.
The result should be again 1 point layer (hexagrid centroids) with the number of points per hexagrid for each class as fields in the attribute table. If a hexagrid contains points from different classes, for each class a separate feature (separate row in the attribute table) should be available.

Example:

I suppose the for-loop (iteration) is not integrated as a tool in the QGIS ModelBuilder so a PyQGIS script is necessary. I´ve tried something but it didn´t worked properly.

My idea (without using PyQGIS) was the following one:

1. select hexagrids with at least 1 point and create a new layer from the selection (extract by location)
2. create centroids for the hexagrids
3. join the grid ID to the points (join attributes by location)
4. join the points to the centroids based on the grid ID (join attributes by field)
5. split the point layer testdata by attribute (class)
6. run the tool "Count points in polygons" in the batch mode for each splitted point layer
7. merge the created polygons together (1 layer)
8. join the merged hexagrid layer with the number of points per grid to the centroid layer (grid ID as the join field)

(alternatively to 6. and 7., join the hexagrid layers with the point count with the centroid layer in the batch mode and then merge them together into 1 layer).

This workflow functions with a limited number of layers (point classes) but for thousands of classes this is impossible so a PyQGIS script seems as the only appropriate sollution.

Could you help me how the python script should look like or do you have any other suggestions?

---

Here is a picture for better understanding:

Each centroid (red) of the hexagrid represents the number of points within the same hexagrid. The assignment to the centroid is based on the hexagrid ID (join attributes by location). If a hexagrid contains points from different classes, for each class a new centroid is created (see hexagrids in circles).