QSOnic
Lightining-fast continuum fitting
QSOnic is an MPI-parallelized, highly optimized quasar continuum fitting package for DESI built on the same algorithm as picca, but faster. It also provides an efficient API to read DESI quasar spectra. If you use this software, please cite the article in the Journal of Open Source Software (JOSS).
Picca to QSOnic translator
QSOnic comes with a picca ini file translator script called picca2qsonic. Usage: picca2qsonic <picca_ini_file> -o submit-qsonic.sh. This will translate the picca ini file to a SLURM script that can be submitted or run on an interactive node. It should be treated as a starting point for your own SLURM script, and you may need to modify it to suit your needs.
The key differences
You can use any desired input continuum from another continuum prediction method as input continuum model.
Coadding of spectrograph arms can be performed after continuum fitting or disabled entirely.
Continuum is multiplied by a fiducial mean flux when provided.
You can pass fiducial var_lss (column VAR_LSS) and mean flux (column MEANFLUX) for observed wavelength LAMBDA in STATS extention of a FITS file. Wavelength should be linearly and equally spaced. This is the same format as rawio output from picca, except VAR column in picca is the variance on flux not deltas. We break away from that convention by explicitly requiring variance on deltas in a new column.
If no fiducial is passed, we fit only for var_lss (no eta fitting by default). Eta fitting can be enabled by passing --var-fit-eta. A covariance matrix is calculated based on delete-one Jackknife over subsamples between var_pipe and var_obs data points, and used in var_lss, eta fitting if --var-use-cov passed (soft recommendation to always enable this).
Internal weights for continuum fitting and coadding are based on smoothed IVAR, and output WEIGHT is based on this smoothed ivar. This smoothing can be turned off.
Chi2 information as well as best fits are saved in continuum_chi2_catalog.fits. Chi2 is calculated using smooth ivar and var_lss, and does not subtract sum of ln(weights).
Similarities
Delta files are the same. CONT column is mean flux times continuum even when fiducial mean flux is passed.
MEANSNR in header file and chi2 catalog is average of flux times square root of positive ivar values. Header values are per arm, but catalog values are the average over all arms.
Eta fitting does not rescale IVAR output. Pipeline noise will be modified with explicit calibration option.
1 maintainer
Metadata
1 maintainer
Release files for qsonic 1.2.1
For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.
Source distribution (sdist)
| File | Size | Uploaded | |
|---|---|---|---|
| qsonic-1.2.1.tar.gz | 73.1 kB | Details |
Built distribution (wheel)
| File | Interpreter | ABI | Platform | Reset |
|---|---|---|---|---|
| qsonic-1.2.1-py3-none-any.whl | Python 3 | none | any | Details |
Total release size: 145.0 kB
Release files / qsonic-1.2.1.tar.gz
| Download URL | qsonic-1.2.1.tar.gz |
|---|---|
| Size | 73.1 kB |
| Tags | Source |
|
SHA-256 checksum How to use checksums |
99f107c78c6418018296ef43da044f3fa8c7edb1b08af05cc8ecddd0c07951e0
|
|
BLAKE2b-256 checksum How to use checksums |
c4244e11f026b5f53f8d8f2a60d9c56dc972c19250f6553784048a5c22f55736
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/7.0.0 CPython/3.14.8
|
Release files / qsonic-1.2.1-py3-none-any.whl
| Download URL | qsonic-1.2.1-py3-none-any.whl |
|---|---|
| Size | 71.9 kB |
| Tags | Python 3 |
|
SHA-256 checksum How to use checksums |
afcf0fd475ab7a71536462846669a9e013369c63d24ccf7eeb9bf2b4dbb5f9a2
|
|
BLAKE2b-256 checksum How to use checksums |
4d2691d8824b09a13b659fd210ac040e67a6aa2c427d3a223b5d83be2610e9c9
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/7.0.0 CPython/3.14.8
|
