skip to main content


Title: COMAP Early Science. III. CO Data Processing
Abstract

We describe the first-season CO Mapping Array Project (COMAP) analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and mapmaking. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High-efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, includingχ2and multiscale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a data set with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers.

 
more » « less
Award ID(s):
1910999
NSF-PAR ID:
10369007
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ;   « less
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
933
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 184
Size(s):
["Article No. 184"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed–Feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing, and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scalesk= 0.051–0.62 Mpc−1, we estimatePCO(k) = −2. 7 ± 1.7 × 104μK2Mpc3, the first direct 3D constraint on the clustering component of the CO(1–0) power spectrum in the literature.

     
    more » « less
  2. Abstract

    We present estimates of line-of-sight distortion fields derived from the 95 and 150 GHz data taken by BICEP2, BICEP3, and the Keck Array up to the 2018 observing season, leading to cosmological constraints and a study of instrumental and astrophysical systematics. Cosmological constraints are derived from three of the distortion fields concerning gravitational lensing from large-scale structure, polarization rotation from magnetic fields or an axion-like field, and the screening effect of patchy reionization. We measure an amplitude of the lensing power spectrumALϕϕ=0.95±0.20. We constrain polarization rotation, expressed as the coupling constant of a Chern–Simons electromagnetic termgaγ≤ 2.6 × 10−2/HI, whereHIis the inflationary Hubble parameter, and an amplitude of primordial magnetic fields smoothed over 1 MpcB1Mpc≤ 6.6 nG at 95 GHz. We constrain the rms of optical depth fluctuations in a simple “crinkly surface” model of patchy reionization, findingAτ< 0.19 (2σ) for the coherence scale ofLc= 100. We show that all of the distortion fields of the 95 and 150 GHz polarization maps are consistent with simulations including lensed ΛCDM, dust, and noise, with no evidence for instrumental systematics. In some cases, theEBandTBquadratic estimators presented here are more sensitive than our previous map-based null tests at identifying and rejecting spuriousB-modes that might arise from instrumental effects. Finally, we verify that the standard deprojection filtering in the BICEP/Keck data processing is effective at removing temperature to polarization leakage.

     
    more » « less
  3. Abstract

    The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background over 75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. This paper describes the CLASS data pipeline and maps for 40 GHz observations conducted from 2016 August to 2022 May. We demonstrate how well the CLASS survey strategy, with rapid (∼10 Hz) front-end modulation, recovers the large-scale Galactic polarization signal from the ground: the mapping transfer function recovers ∼67% (85%) ofEEandBB(VV) power at= 20 and ∼35% (47%) at= 10. We present linear and circular polarization maps over 75% of the sky. Simulations based on the data imply the maps have a white noise level of110μKarcminand correlated noise component rising at low-as−2.4. The transfer-function-corrected low-component is comparable to the white noise at the angular knee frequencies of≈ 18 (linear polarization) and≈ 12 (circular polarization). Finally, we present simulations of the level at which expected sources of systematic error bias the measurements, finding subpercent bias for the Λ cold dark matterEEpower spectra. Bias fromE-to-Bleakage due to the data reduction pipeline and polarization angle uncertainty approaches the expected level for anr= 0.01BBpower spectrum. Improvements to the instrument calibration and the data pipeline will decrease this bias.

     
    more » « less
  4. ABSTRACT

    The 1D power spectrum P1D of the Ly α forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of P1D with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of 54 600 quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with 2D image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly α emission line. In a companion paper, we present a similar analysis based on the Fast Fourier Transform estimate of the power spectrum. We conclude with a comparison of these two approaches and discuss the key sources of systematic error that we need to address with the upcoming DESI Year 1 analysis.

     
    more » « less
  5. Abstract

    We present a catalog of 536 fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst (CHIME/FRB) Project between 400 and 800 MHz from 2018 July 25 to 2019 July 1, including 62 bursts from 18 previously reported repeating sources. The catalog represents the first large sample, including bursts from repeaters and nonrepeaters, observed in a single survey with uniform selection effects. This facilitates comparative and absolute studies of the FRB population. We show that repeaters and apparent nonrepeaters have sky locations and dispersion measures (DMs) that are consistent with being drawn from the same distribution. However, bursts from repeating sources differ from apparent nonrepeaters in intrinsic temporal width and spectral bandwidth. Through injection of simulated events into our detection pipeline, we perform an absolute calibration of selection effects to account for systematic biases. We find evidence for a population of FRBs—composing a large fraction of the overall population—with a scattering time at 600 MHz in excess of 10 ms, of which only a small fraction are observed by CHIME/FRB. We infer a power-law index for the cumulative fluence distribution ofα=1.40±0.11(stat.)0.09+0.06(sys.), consistent with the −3/2 expectation for a nonevolving population in Euclidean space. We find thatαis steeper for high-DM events and shallower for low-DM events, which is what would be expected when DM is correlated with distance. We infer a sky rate of[820±60(stat.)200+220(sys.)]/sky/dayabove a fluence of 5 Jy ms at 600 MHz, with a scattering time at 600 MHz under 10 ms and DM above 100 pc cm−3.

     
    more » « less