The shell-model-like approach is implemented to treat the cranking many-body
Hamiltonian based on the covariant density functional theory including pairing
correlations with exact particle number conservation. The self-consistency is
achieved by iterating the single-particle occupation probabilities back to the
densities and currents. As an example, the rotational structures observed in
the neutron-rich nucleus 60Fe are investigated and analyzed. Without
introducing any \emph{ad hoc} parameters, the bandheads, the rotational
spectra, and the relations between the angular momentum and rotational
frequency for the positive parity band A, and negative parity bands B and C are
well reproduced. The essential role of the pairing correlations is revealed. It
is found that for band A, the bandcrossing is due to the change of the last two
occupied neutrons from the 1f5/2 signature partners to the 1g9/2
signature partners. For the two negative parity signature partner bands B and
C, the bandcrossings are due to the pseudo-crossing between the 1f7/2,5/2
and the 1f5/2,1/2 orbitals. Generally speaking, the deformation
parameters β for bands A, B, and C decrease with rotational frequency.
For band A, the deformation jumps from β∼0.19 to β∼0.29
around the bandcrossing. In comparison with its signature partner band C, band
B exhibits appreciable triaxial deformation