We study the polarization properties of extragalactic sources at 95 and 150 GHz in the SPTpol 500 deg2 survey. We estimate the polarized power by stacking maps at known source positions, and correct for noise bias by subtracting the mean polarized power at random positions in the maps. We show that the method is unbiased using a set of simulated maps with similar noise properties to the real SPTpol maps. We find a flux-weighted mean-squared polarization fraction 〈p2〉= [8.9 ± 1.1] × 10−4 at 95 GHz and [6.9 ± 1.1] × 10−4 at 150 GHz for the full sample. This is consistent with the values obtained for a subsample of active galactic nuclei. For dusty sources, we find 95 per cent upper limits of 〈p2〉95 < 16.9 × 10−3 and 〈p2〉150 < 2.6 × 10−3. We find no evidence that the polarization fraction depends on the source flux or observing frequency. The 1σ upper limit on measured mean-squared polarization fraction at 150 GHz implies that extragalactic foregrounds will be subdominant to the CMB E and B mode polarization power spectra out to at least ℓ ≲ 5700 (ℓ ≲ 4700) and ℓ ≲ 5300 (ℓ ≲ 3600), respectively, at 95 (150) GHz.

Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3 × 2pt, had to discard a lot of signal to noise from small scales due to our inability to accurately model non-linearities and baryonic effects. Galaxy–galaxy lensing, or the position–shear correlation between lens and source galaxies, is one of the three two-point correlation functions that are included in such analyses, usually estimated with the mean tangential shear. However, tangential shear measurements at a given angular scale θ or physical scale R carry information from all scales below that, forcing the scale cuts applied in real data to be significantly larger than the scale at which theoretical uncertainties become problematic. Recently, there have been a few independent efforts that aim to mitigate the non-locality of the galaxy–galaxy lensing signal. Here, we perform a comparison of the different methods, including the Y-transformation, the point-mass marginalization methodology, and the annular differential surface density statistic. We do the comparison at the cosmological constraints level in a combined galaxy clustering and galaxy–galaxy lensing analysis. We find that all the estimators yield equivalent cosmological results assuming a simulated Rubin Observatory Legacy Survey of Space and Time (LSST) Year 1 like set-up and also when applied to DES Y3 data. With the LSST Y1 set-up, we find that the mitigation schemes yield ∼1.3 times more constraining S8 results than applying larger scale cuts without using any mitigation scheme.