Let $G$ be a finite group and $H_G$ be the set of harmonic polynomials. In the case of the symmetric group these polynomials are spanned by the Vandermonde determinant and all its partial derivatives. This result seems to generalize to all finite reflection groups. For a general finite group, it is clear that $H_G$ is a finite dimensional vector space, but is there a similar construction to obtain a basis?