the temperature dependence of the EDCs at the X-point, stacked in Fig. \[fig:fig3\](c). At all temperatures, there is a large hump background corresponding to the large hole-like dispersion at high binding energy. At low temperatures, there is an additional peak around -0.05eV corresponding to the [$d_{xy}$]{} band bottom. We fit these EDCs with a Gaussian background for the large hump feature and a Lorentzian peak for the [$d_{xy}$]{} band. The integrated spectral weight for the fitted [$d_{xy}$]{} peak is plotted in Fig. \[fig:fig3\](d), which decreases toward zero with increasing temperature, seen as a non-trivial drop between 100K and 200K. As an independent check against trivial thermal effect, we choose small regions in the raw spectral image (marked in Fig. \[fig:fig3\](a)-(b)) dominated by [$d_{xy}$]{} (blue), [$d_{yz}$]{} (green), and mixed [$d_{xy}$]{} and [$d_{xz}$/$d_{yz}$]{} (magenta/cyan) characters and plot their integrated intensities as a function of temperature (Fig. \[fig:fig3\](e)). The spectral weight from [$d_{xy}$]{}-dominated region rapidly decreases, consistent with the fitted result in Fig. \[fig:fig3\](d), while that of [$d_{yz}$]{}-dominated region does not drop in a similar manner. The regions of mixed orbitals show a slower diminishing spectral weight compared to that for [$d_{xy}$]{}, reflecting the combined contributions from both [$d_{xz}$/$d_{yz}$]{} and [$d_{xy}$]{} orbitals. Although 