def METHOD_NAME(self):
    """Mark geometry for materials, boundary conditions, faults, etc.
    """
    materials = (
        MaterialGroup(tag=1, entities=[self.box]),
    )
    for material in materials:
        material.create_physical_group()
    vertex_groups = (
        VertexGroup(name="boundary_xneg", tag=10, dim=2, entities=[self.s_xneg]),
        VertexGroup(name="boundary_xpos", tag=11, dim=2, entities=[self.s_xpos]),
        VertexGroup(name="boundary_yneg", tag=12, dim=2, entities=[self.s_yneg]),
        VertexGroup(name="boundary_ypos", tag=13, dim=2, entities=[self.s_ypos]),
        VertexGroup(name="boundary_zneg", tag=14, dim=2, entities=[self.s_zneg]),
        VertexGroup(name="boundary_zpos", tag=15, dim=2, entities=[self.s_zpos]),
        VertexGroup(name="fault", tag=20, dim=2, entities=[self.s_fault]),
        VertexGroup(name="fault_edges", tag=21, dim=1, entities=[self.c_fault_yneg, self.c_faultbot, self.c_fault_ypos]),
    )
    for group in vertex_groups:
        group.create_physical_group()