computes the sector of rotational axes for orientations within region
sR = axisSector(oR) sR = axisSector(oR,omgea)