We study in detail the spectrum of heavy quarkonia with different orbital angular momentum along with their radial and gluonic excitations. Using an anisotropic formulation of Lattice QCD we achieved an unprecedented control over statistical errors and were able to study systematic errors such as lattice spacing artefacts, finite volume effects and relativistic corrections. First results on the spin structure in heavy hybrids are also presented.