The superstructures of two cyclodextrin derivatives chemisorbed on a gold (111) surface were modeled using a novel random planting algorithm combined with the standard molecular dynamics simulation technique. Surface density and thickness of the modeled cyclodextrin layers are found to be in good agreement with corresponding experimental results. In addition, a number of structural quantities characterizing the lateral order in the films as well as the influence of solvent are investigated.