def atoms_inside(self, atoms):
"""Decides which atoms are inside the body (see Body class)."""
atoms_inside = Cylinder.atoms_inside(self, atoms)
atoms_inside *= self.periodicity.mask_unique(atoms - self.shift_vector,
atoms_inside)
return atoms_inside
