A diagrammatic method for the calculation of moments of operators and its application to nuclear spectroscopy