qsonic-calib

usage: qsonic-calib [-h] --input-dir INPUT_DIR --outdir OUTDIR [--fbase FBASE]
                    [--catalog CATALOG] [--mock-analysis]
                    [--keep-surveys [KEEP_SURVEYS ...]] [--remove-bal-qsos]
                    [--remove-targetid-list REMOVE_TARGETID_LIST]
                    [--nvarbins NVARBINS] [--var-use-cov] [--nwbins NWBINS]
                    [--var1 VAR1] [--var2 VAR2] [--min-snr MIN_SNR]
                    [--max-snr MAX_SNR] [--wave1 WAVE1] [--wave2 WAVE2]
                    [--forest-w1 FOREST_W1] [--forest-w2 FOREST_W2]

Input/output parameters and selections

--input-dir, -i

Input directory.

--outdir, -o

Output directory to save files.

--fbase

Basename

Default: “qsonic-eta-fits”

--catalog

Catalog filename to enable catalog related removals.

--mock-analysis

Input folder is mock. Uses nside=16

Default: False

--keep-surveys

Surveys to keep. Empty keeps all.

--remove-bal-qsos

Removes BAL sightlines in the catalog option.

Default: False

--remove-targetid-list

Text file with TARGETIDs to exclude from analysis.

Variance fitting parameters

--nvarbins

Number of variance bins (logarithmically spaced).

Default: 100

--var-use-cov

Use covariance in varlss-eta fitting.

Default: False

--nwbins

Number of wavelength bins. None creates bins with 120 A spacing

--var1

Lower variance bin.

Default: 0.0001

--var2

Upper variance bin.

Default: 20.0

--min-snr

Minimum SNR of the forest.

Default: 0

--max-snr

Maximum SNR of the forest.

Default: 100

Wavelength analysis region

--wave1

First observed wavelength edge.

Default: 3600.0

--wave2

Last observed wavelength edge.

Default: 6000.0

--forest-w1

First forest wavelength edge.

Default: 1050.0

--forest-w2

Last forest wavelength edge.

Default: 1180.0

Note: Not every option is used in this script.

Functions

get_parser(add_help=True)

Constructs the parser needed for the script.

Parameters:

add_help (bool, default: True) – Add help to parser.

Returns:

parser

Return type:

argparse.ArgumentParser

mpi_set_targetid_list_to_remove(args, comm=None, mpi_rank=0)

Return a ndarray of TARGETIDs to remove from the sample.

Can be used without MPI by passing comm=None (which is the default.)

Parameters:

• args (argparse.Namespace) – Options passed to script.

• comm (MPI.COMM_WORLD or None, default: None) – Communication object broadcast data.

• mpi_rank (int, default: 0) – Rank of the MPI process.

Returns:

ids_to_remove – TARGETIDs to remove from the sample.

Return type:

ndarray

Raises:

QsonicException – If error occurs while reading args.remove_targetid_list or if --remove_bal_qsos is passed but the input catalog is missing BAL columns .

mpi_read_all_deltas(args, comm=None, mpi_rank=0, mpi_size=1)

mpi_stack_fluxes(args, comm, deltas_list)

mpi_run_all(comm, mpi_rank, mpi_size)

main()