Animate Clock Incomplete Circle

from manimlib.imports import *

class ClockOrganization(VGroup):
  CONFIG = {
    "numbers" : 4,
    "radius" : 3.1,
    "color" : WHITE
  }

  def __init__(self, **kwargs):
    digest_config(self, kwargs, locals())
    self.generate_nodes()
    VGroup.__init__(self, *self.node_list,**kwargs)

  def generate_nodes(self):
    self.node_list = []
    for i in range(self.numbers):
      mobject = VMobject()
      number = TexMobject(str(i+1))
      circle = Circle(radius=0.4,color=self.color)
      mobject.add(number)
      mobject.add(circle)
      mobject.move_to(
        self.radius * np.cos((-TAU / self.numbers) * i + 17*TAU / 84) * RIGHT
        + self.radius * np.sin((-TAU / self.numbers) * i + 17*TAU / 84) * UP
      )
      self.node_list.append(mobject)

  def select_node(self, node):
    selected_node = self.node_list[node]
    selected_node.scale(1.2)
    selected_node.set_color(RED)

  def deselect_node(self, selected_node):
    node = self.node_list[selected_node]
    node.scale(0.8)
    node.set_color(self.color)

class Testing4(Scene):
  def construct(self):
    test = ClockOrganization(numbers=21)
    self.play(Write(test), run_time=1.5)
    animation_steps=[]
    #animation_steps.append(test)
    num_circ=15
    for i in range(num_circ):
      thing = test.deepcopy()
      thing.select_node((1+i)%test.numbers-1)
      animation_steps.append(thing)
    test_normal=test.copy()
    test.save_state()
    self.play(Transform(test, animation_steps[0]))
    self.wait(2)
    anims=[Restore(test)]
    theta=180*DEGREES/num_circ
    lag_constant=1.25
    for i in range(1,num_circ+1):
        test.node_list[(1+i)%test.numbers-1].generate_target()
        test.node_list[(1+i)%test.numbers-1].target.scale(1.2)
        test.node_list[(1+i)%test.numbers-1].target.set_color(RED)
        if i in range(1,8):
          stop_smooth=(lag_constant*sigmoid(-10/7.5*(i-3.46))-sigmoid(-5))/(1-2*sigmoid(-5))
        else:
          stop_smooth=(lag_constant*sigmoid(10/7.5*(i-7.5-0.812-3.46))-sigmoid(-5))/(1-2*sigmoid(-5))
        if i in range(1,num_circ):
          anims.append(MoveToTarget(test.node_list[(1+i)%test.numbers-1],rate_func=there_and_back))
        else:
          anims.append(MoveToTarget(test.node_list[(1+i)%test.numbers-1]))
        anims.append(Animation(Mobject(),run_time=stop_smooth))
    self.play(
        AnimationGroup(*anims,lag_ratio=0.05)
        )
    self.wait()