We investigated binding of hydrogen atoms to small Polycyclic Aromatic Hydrocarbons (PAHs) -i.e. graphene dots with hydrogen-terminated edges -using density functional theory and correlated wavefunction techniques. We considered a number of PAHs with 3 to 7 hexagonal rings and computed binding energies for most of the symmetry unique sites, along with the minimum energy paths for significant cases. The chosen PAHs are small enough to not present radical character at their edges, yet show a clear preference for adsorption at the edge sites which can be attributed to electronic effects. We show how the results, as obtained at different level of theory, can be rationalized in detail with the help of few simple concepts derivable from a tight-binding model of the π electrons.