igmhub / picca

set of tools for continuum fitting, correlation function calculation, cosmological fits...
GNU General Public License v3.0
30 stars 22 forks source link

new option per_r_par for covariance smoothing #1030

Closed julienguy closed 11 months ago

julienguy commented 11 months ago

This is useful for some correlation functions where the default scheme is not accurate and can lead to non positive matrices. For instance SIV and CIV side bands auto-correlations, or the lyb x lya cross-correlation.

Example:

# default (non positive output matrix)
picca_export.py --data /global/cfs/cdirs/desicollab/science/lya/y1-kp6/iron-tests/correlations/correlation-lyalyb-1-0-0/cf_lya_x_lyb.fits --out cf_lya_x_lyb_exp.fits --blind-corr-type lyaxlyb

# with new option (positive output matrix)
picca_export.py --data /global/cfs/cdirs/desicollab/science/lya/y1-kp6/iron-tests/correlations/correlation-lyalyb-1-0-0/cf_lya_x_lyb.fits --out cf_lya_x_lyb_exp_bis.fits --blind-corr-type lyaxlyb --smooth-per-r-par
julienguy commented 11 months ago

done