docs: def get_smpar(self, accuracy='integrate', scale_sm=91.1876):
          """Compute the SM MS-bar parameters at the electroweak scale.
  
          This method can be used to validate the accuracy of the iterative
          extraction of SM parameters. If successful, the values returned by this
          method should agree with the values in the dictionary
          `wilson.run.smeft.smpar.p`."""
          if accuracy == 'integrate':
              C_out = self._rgevolve(scale_sm)
          elif accuracy == 'leadinglog':
              C_out = self._rgevolve_leadinglog(scale_sm)
          else:
              raise ValueError("'{}' is not a valid value of 'accuracy' (must be either 'integrate' or 'leadinglog').".format(accuracy))
          return smpar.smpar(C_out)