Procedural Blender geometry study (Based on David Mignot work)

This project was made completely procedurally with python code.

This grid was made with 30x30 units.
There's a 20% chance for the cube to be a "tower", which can be from 2 up to 10 stacks.
Tower's cubes get rotated by 0.785398 radians, or 45°, in Z.
Towers get an emission material, normal cubes get a metallic (.780) and Roughness (.366)