2024ApJ...966..158P

Blazars are characterized by largely aperiodic variability on timescales ranging from minutes to decades across the electromagnetic spectrum. The Transiting Exoplanet Survey Satellite (TESS) mission provides continuous sampling of blazar variability on timescales ranging from tens of minutes to 27 days for a single sector observation. Proper removal of the background, thermal ramping, and onboard systematic effects are crucial to the extraction of a reliable blazar light curve. Multiple publicly available procedures have been created to correct for these effects. Using ground-based observations from the Zwicky Transit Facility (ZTF) and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as “ground truth” observations, we compare six different methods (regression, cotrending basis vectors, pixel-level decorrelation, eleanor, quaver, and simple differential photometry) to each other, and to our “ground truth” observations, to identify which methods properly correct light curves of a sample of 11 bright blazars, including the suspected neutrino source TXS 0506+056. In addition to comparing the resulting light curves, we compare the slopes of the power spectral densities, perform least-square fitting to simultaneous ZTF data, and quantify other statistical qualities generated from the light curves of each method. We find that only three of the six methods compared (simple differential photometry, eleanor, and quaver) produce TESS light curves consistent with the ground-based ZTF and ATLAS observations.