Non-linear perturbation theory extension of the Boltzmann code CLASS

We present a new open-source code that calculates one-loop power spectra and cross spectra for matter fields and biased tracers in real and redshift space. These spectra incorporate all ingredients required for a direct application to data: nonlinear bias and redshift-space distortions, infrared resummation, counterterms, and the Alcock-Paczynski effect. Our code is based on the Boltzmann solver class and inherits its advantageous properties: user friendliness, ease of modification, high speed, and simple interface with other software. We present detailed descriptions of the theoretical model, the code structure, approximations, and accuracy tests. A typical end-to-end run for one cosmology takes 0.3 seconds, which is sufficient for Markov chain Monte Carlo parameter extraction. As an example, we apply the code to the Baryon Oscillation Spectroscopic Survey (BOSS) data and infer cosmological parameters from the shape of the galaxy power spectrum.