By default, this selects all atoms near gaps in the alignment for loop optimization (by calling model.loops()). You should redefine this routine if you do not have an alignment, or you wish to set a different region for loop optimization. See section 2.3 for an example.