def test_quantum_number_L_z_dimers(self, job_file, rho_file, expected_L_z):
        # Setup
        Job = InitJob(job_file)

        # Fake
        rho = load_json_file(rho_file)
        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