def METHOD_NAME(self, basis: VibrationalBasis) -> VibrationalStructureProblem:  # type: ignore[override]
    # pylint: disable=arguments-differ
    if self._logfile is not None:
        glr = GaussianLogResult(self._logfile)
    else:
        glr = GaussianLogDriver(jcf=self._jcf).METHOD_NAME()
    watson = glr.get_watson_hamiltonian(normalize=self._normalize)
    problem = watson_to_problem(watson, basis=basis)
    return problem