Stellar orbit distribution of TNG50 galaxy. The orbits are represented as the probability density distribution in the phase space of circularity λ_z vs. radius r. On the left, we use the λ_z and r calculated from instantaneous positions and velocities of the stellar particles. On the right, we use the orbital λ_z and r obtained through a phase-space averaging approach, where particles with similar energy and angular momentum are grouped. Dashed lines are used to demarcate the boundaries separating the cold (λ_z > 0.8), warm (0.25 < λ_z < 0.8), hot (|λ_z| < 0.25), and counter-rotating (CR) (λ_z < −0.25) components based on circularity.