def METHOD_NAME(self, kT):
    self.system.thermostat.set_brownian(
        kT=kT,
        gamma=[
            self.gamma_ortho,
            self.gamma_ortho,
            self.gamma_parallel],
        seed=self.seed)
    self.system.integrator.set_brownian_dynamics()