Quantitative constraint on the contribution of resolved gamma-ray sources to the sub-PeV Galactic diffuse gamma-ray flux measured by the Tibet AS{\gamma} experiment

Motivated by the difference between the fluxes of sub-PeV Galactic diffuse gamma-ray emission (GDE) measured by the Tibet AS$\gamma$ experiment and the LHAASO collaboration, our study constrains the contribution to the GDE flux measured by Tibet AS$\gamma$ from the sub-PeV gamma-ray sources in the first LHAASO catalog plus the Cygnus Cocoon. After removing the gamma-ray emission of the sources masked in the observation by Tibet AS$\gamma$, the contribution of the sources to the Tibet diffuse flux is found to be subdominant; in the sky region of $25^{\circ} < l < 100^{\circ}$ and $|b| < 5^{\circ}$, it is less than 26.9% $\pm$ 9.9%, 34.8% $\pm$ 14.0%, and ${13.5%}^{+6.3%}_{-7.7%}$ at 121 TeV, 220 TeV, and 534 TeV, respectively. In the sky region of $50^{\circ} < l < 200^{\circ}$ and $|b| < 5^{\circ}$, the fraction is less than 24.1% $\pm$ 9.5%, 27.1% $\pm$ 11.1% and ${13.5%}^{+6.2%}_{-7.6%}$. After subtracting the source contribution, the hadronic diffusive nature of the Tibet diffuse flux is the most natural interpretation, although some contributions from very faint unresolved hadronic gamma-ray sources cannot be ruled out. Different source-masking schemes adopted by Tibet AS$\gamma$ and LHAASO for their diffuse analyses result in different effective galactic latitudinal ranges of the sky regions observed by the two experiments. Our study concludes that the effect of the different source-masking schemes leads to the observed difference between the Tibet diffuse flux measured in $25^{\circ} < l < 100^{\circ}$ and $|b| < 5^{\circ}$ and LHAASO diffuse flux in $15^{\circ} < l < 125^{\circ}$ and $|b| < 5^{\circ}$.