def METHOD_NAME(self):
    """Test for dist thermalization."""
    N = 100
    N_half = int(N / 2)
    self.system.part.clear()
    self.system.time_step = 0.02
    self.system.periodicity = [True, True, True]
    m1 = 1.0
    m2 = 10.0
    # Place particles
    partcls_m1 = self.system.part.add(
        pos=np.random.random(
            (N_half, 3)), mass=N_half * [m1])
    partcls_m2 = self.system.part.add(
        pos=np.random.random(
            (N_half, 3)), mass=N_half * [m2])
    t_dist = 2.0
    g_dist = 4.0
    t_com = 0.0
    g_com = 0.0
    thermalized_dist_bond = espressomd.interactions.ThermalizedBond(
        temp_com=t_com, gamma_com=g_com, temp_distance=t_dist,
        gamma_distance=g_dist, r_cut=9, seed=51)
    self.system.bonded_inter.add(thermalized_dist_bond)
    for p1, p2 in zip(partcls_m1, partcls_m2):
        p1.add_bond((thermalized_dist_bond, p2))
    # Warmup
    self.system.integrator.run(50)
    # Sampling
    loops = 150
    v_stored = np.zeros((N_half * loops, 3))
    for i in range(loops):
        self.system.integrator.run(5)
        v_dist = partcls_m2.v - partcls_m1.v
        v_stored[i * N_half:(i + 1) * N_half, :] = v_dist
    v_minmax = 5
    bins = 50
    error_tol = 0.015
    self.check_velocity_distribution(
        v_stored, v_minmax, bins, error_tol, t_dist)