This code propagates the half-wave plate (HWP) non-idealities through map-making, foreground cleaning (HILC), and angular power spectra estimation, to return the value of the observed tensor-to-scalar ratio. It was used to perform the analysis presented in https://arxiv.org/abs/2311.07999 and returns the same plots as those shown in the paper.
Dependencies (other than standard Python libraries): healpy and camb.
The avaliable input files are for an ideal HWP. Realistic input files can be provided upon request.