Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Update redrock prior widths to be wider than QuasarNET "wiggles" in redshift #2293

Closed
dylanagreen opened this issue Jul 15, 2024 · 1 comment
Assignees

Comments

@dylanagreen
Copy link
Contributor

QuasarNET uses a grid based system (called boxes) to estimate line positions which causes unphysical "wiggles" in the N(z) distribution due to an unphysical centering of estimated line positions in these boxes, see the following plot near the peak of the distribution:

z_comp

Currently after running QuasarNET we rerun redrock with a redshift prior centered at the QuasarNET redshift and QSO templates only in order to recover these quasars with updated redshifts. The current redrock prior width is $\Delta z = 0.1$, which we have since discovered is narrower than the width of the wiggles in the QuasarNET N(z).

However, we can analytically calculate the required width of the redrock prior necessary to remove the wiggles in the final redshift catalog. We do this by calculating the width of each output QuasarNET box in redshift space and using this as the associated width of the redrock prior. This width is calculated assuming the QuasarNET has correctly estimated the "coarse" position of the line (it got the right box assignment) but is inaccurate in its "fine" estimation, possibly centering the line unphysically in the center of the box. Using the width of the box as a prior should allow redrock to correct the redshift and remove the wiggles.

Using a logarithmic output grid, which QuasarNET uses in Iron, creates a linear dependence of box width on estimated QuasarNET redshift:
box_width_log

Fitting the linear dependence provides the following equation:
$$\Delta z = 0.07842467 * z + 0.07842494$$

It is probably sufficient to round this to three or four decimal places.

@akremin
Copy link
Member

akremin commented Aug 7, 2024

Closed with PR #2298

@akremin akremin closed this as completed Aug 7, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants