We have analyzed the as-implanted damage produced in silicon by B, Si and Ge ions using molecular dynamics (MD) simulations. Implantations were carried out at 50 K to avoid damage migration and annealing. In order to make a statistical study of the damage features, we have simulated hundreds of independent cascades for each ion for the same nuclear deposited energy. We have obtained that the average number of displaced atoms (DA) from perfect lattice positions and the size of defect clusters formed increases with ion mass. This dependence has not been obtained from equivalent binary collisions simulations. This indicates that multiple interactions play an important role in the generation of damage. Amorphous regions are directly formed during the collisional phase of the cascade of Ge and Si ions.