def test_quantum_number_L_z(self, name, rho, expected_L_z):
        # Setup
        Job = InitJob("test_data/JobDef_pcase.json")

        # Fake
        Job.NAtom = 1
        Job.NOrb = [len(rho)/2]
        Job.Electron.NElectrons = sum(rho[ii][ii] for ii in range(len(rho)))
        Job.Electron.rho = np.matrix(rho, dtype='complex')

        # Action
        L_z = Job.Electron.quantum_number_L_z(Job)

        # Result
        assert L_z == expected_L_z