def METHOD_NAME(self):
    a = ase.Atoms('CC', positions=[[0.5, 0.5, 0.5], [0.5, 0.5, 1.0]],
                  cell=[2, 2, 2], pbc=True)
    i, j, dr = neighbour_list("ijD", a, 1.1)
    hist = angle_distribution(i, j, dr, 20, 1.1)
    self.assertEqual(hist.sum(), 0)