Skip to content
forked from igmhub/picca

set of tools for continuum fitting, correlation function calculation, cosmological fits...

License

Notifications You must be signed in to change notification settings

SamYoules/picca

 
 

Repository files navigation

picca

Package for Igm Cosmological-Correlations Analyses.

This package contains tools used for the analysis of the Lyman-alpha forest sample from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the Dark Energy Spectroscopic Instrument (DESI). Here you will find tools to

  • fit continua of forests
  • compute correlation functions (1D and 3D) and power-spectra (1D)
  • compute covariance matrices
  • fit models for the correlation functions

The current reference is du Mas des Bourboux et al. 2020 (https://arxiv.org/abs/2007.08995).

Installation

download

git clone https://github.com/igmhub/picca.git

add to your bashrc

export PICCA_BASE=<path to your picca>

then make sure you have all required modules by running

pip install -r requirements.txt --user

and finally run

python setup.py install --user

(assuming you run as user; for a system-wide install omit --user option).

Alternatively, you can just add picca/py/ to your PYTHONPATH.

Examples

example run over 1000 spectra (the DLA catalog is not required):

delta field

picca_deltas.py
--in-dir data/
--drq ../DR14Q_v1_1.fits
--dla-vac ../dlas/DLA_DR14_v1b.dat
--out-dir deltas/
--mode pix
  • --mode can be pix (Anze/Jose format), spec (spec- files) or corrected-spec (corrected-spec files)
  • --in-dir points to the directory containing the data
  • NOTE: reading the spec files is very slow

correlation function

picca_cf.py
--in-dir deltas/
--out cf.fits.gz
--nside 32
  • nside determines the healpixelization used for the subsamples. nside=32 gives ~3200 subsamples for DR12.

distortion matrix

picca_dmat.py
--in-dir deltas/
--out dmat.fits.gz
--rej 0.95
  • --rej is 1-fraction of pairs used for the calculation

wick covariance (optional)

Only T123 implemented

# first calculate cf_1d from data
picca_cf1d.py
--in-dir deltas/
--out cf1d.fits.gz

# then use it for wick
picca_wick.py
--in-dir deltas/
--out t123.fits.gz
--rej 0.999
--cf1d cf1d.fits.gz

## use the export script to export to picca fitter format
picca_export.py
--data cf.fits.gz
--dmat dmat.fits.gz
--out cf-exp.out.gz

Name of tags

The tags name follow the names of the king of France:
https://fr.wikipedia.org/wiki/Liste_des_monarques_de_France#Liste_des_monarques_de_France

About

set of tools for continuum fitting, correlation function calculation, cosmological fits...

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages

  • Python 56.0%
  • Jupyter Notebook 44.0%