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.

Sign up for GitHub

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

Jump to bottom

add cuts in quasar redshift #977

Merged
merged 13 commits into from
Mar 14, 2023
Merged

add cuts in quasar redshift #977

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
f0bc2c7
add cuts in quasar redshift
cramirezpe Mar 6, 2023
014767a
renamed files in test
iprafols Mar 13, 2023
02cacbc
typo in filename
cramirezpe Mar 13, 2023
9b1f229
remove zcut max
cramirezpe Mar 13, 2023
490be8d
Revert "remove zcut max"
cramirezpe Mar 13, 2023
6bc08a0
remove zcutmax
cramirezpe Mar 13, 2023
e5e57ad
remove zcutmin
cramirezpe Mar 13, 2023
2265da5
Revert "remove zcutmin"
cramirezpe Mar 13, 2023
6183509
Revert "remove zcutmax"
cramirezpe Mar 13, 2023
64cd599
rename test
cramirezpe Mar 14, 2023
deac3e4
remove print functions
cramirezpe Mar 14, 2023
f0381c3
undid rename test, added test reset to remove hidden dependencies
iprafols Mar 14, 2023
7057a4b
Merge branch 'add_qso_cut' of https://github.com/igmhub/picca into ad…
iprafols Mar 14, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
23 changes: 21 additions & 2 deletions bin/picca_cf.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -105,6 +105,20 @@ def main(cmdargs):
'of the last absorber redshift and the object '
'redshift smaller than z-cut-max'))

parser.add_argument('--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument('--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument('--lambda-abs',
type=str,
default='LYA',
Expand Down Expand Up @@ -249,6 +263,7 @@ def main(cmdargs):
t0 = time.time()

### Read data 1
print('z_min_sources', args.z_min_sources)
data, num_data, z_min, z_max = io.read_deltas(args.in_dir,
cf.nside,
cf.lambda_abs,
Expand All @@ -258,7 +273,9 @@ def main(cmdargs):
max_num_spec=args.nspec,
no_project=args.no_project,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources,)
del z_max
cf.data = data
cf.num_data = num_data
Expand Down Expand Up @@ -288,7 +305,9 @@ def main(cmdargs):
max_num_spec=args.nspec,
no_project=args.no_project,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources,)
del z_max2
cf.data2 = data2
cf.num_data2 = num_data2
Expand Down
22 changes: 20 additions & 2 deletions bin/picca_cf1d.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -71,6 +71,20 @@ def main(cmdargs):
required=False,
help='Upper limit on observed wavelength [Angstrom]')

parser.add_argument('--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument('--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument('--dll',
type=float,
default=3.e-4,
Expand Down Expand Up @@ -164,7 +178,9 @@ def main(cmdargs):
args.z_ref,
cosmo=None,
max_num_spec=args.nspec,
no_project=args.no_project)
no_project=args.no_project,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
cf.data = data
cf.num_data = num_data
del z_min, z_max
Expand All @@ -182,7 +198,9 @@ def main(cmdargs):
args.z_ref,
cosmo=None,
max_num_spec=args.nspec,
no_project=args.no_project)
no_project=args.no_project,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
cf.data2 = data2
cf.num_data2 = num_data2
del z_min2, z_max2
Expand Down
24 changes: 22 additions & 2 deletions bin/picca_cf_angl.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -107,6 +107,22 @@ def main(cmdargs):
'absorber redshift and the object redshift smaller than '
'z-cut-max'))

parser.add_argument(
'--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--lambda-abs',
type=str,
Expand Down Expand Up @@ -202,7 +218,9 @@ def main(cmdargs):
cf.z_ref,
cosmo=None,
max_num_spec=args.nspec,
no_project=args.no_project)
no_project=args.no_project,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
cf.data = data
cf.num_data = num_data
del z_min, z_max
Expand All @@ -228,7 +246,9 @@ def main(cmdargs):
cf.z_ref,
cosmo=None,
max_num_spec=args.nspec,
no_project=args.no_project)
no_project=args.no_project,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
cf.data2 = data2
cf.num_data2 = num_data2
del z_min2, z_max2
Expand Down
24 changes: 22 additions & 2 deletions bin/picca_dmat.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -118,6 +118,22 @@ def main(cmdargs):
'absorber redshift and the object redshift smaller than '
'z-cut-max'))

parser.add_argument(
'--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--lambda-abs',
type=str,
Expand Down Expand Up @@ -280,7 +296,9 @@ def main(cmdargs):
max_num_spec=args.nspec,
no_project=args.no_project,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
del z_max
cf.data = data
cf.num_data = num_data
Expand Down Expand Up @@ -310,7 +328,9 @@ def main(cmdargs):
max_num_spec=args.nspec,
no_project=args.no_project,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
del z_max2
cf.data2 = data2
cf.num_data2 = num_data2
Expand Down
24 changes: 22 additions & 2 deletions bin/picca_metal_dmat.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -129,6 +129,22 @@ def main(cmdargs):
'absorber redshift and the object redshift smaller than '
'z-cut-max'))

parser.add_argument(
'--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--lambda-abs',
type=str,
Expand Down Expand Up @@ -313,7 +329,9 @@ def main(cmdargs):
cf.cosmo,
max_num_spec=args.nspec,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
del z_max
cf.data = data
cf.num_data = num_data
Expand Down Expand Up @@ -345,7 +363,9 @@ def main(cmdargs):
cf.cosmo,
max_num_spec=args.nspec,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
del z_max2
cf.data2 = data2
cf.num_data2 = num_data2
Expand Down
20 changes: 19 additions & 1 deletion bin/picca_metal_xdmat.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -143,6 +143,22 @@ def main(cmdargs):
'absorber redshift and the object redshift smaller than '
'z-cut-max'))

parser.add_argument(
'--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--lambda-abs',
type=str,
Expand Down Expand Up @@ -300,7 +316,9 @@ def main(cmdargs):
cosmo,
max_num_spec=args.nspec,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)

xcf.data = data
xcf.num_data = num_data
Expand Down
24 changes: 22 additions & 2 deletions bin/picca_wick.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -109,6 +109,22 @@ def main(cmdargs):
'absorber redshift and the object redshift smaller than '
'z-cut-max'))

parser.add_argument(
'--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--lambda-abs',
type=str,
Expand Down Expand Up @@ -307,7 +323,9 @@ def main(cmdargs):
cosmo,
max_num_spec=args.nspec,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
for deltas in data.values():
for delta in deltas:
delta.fname = 'D1'
Expand Down Expand Up @@ -399,7 +417,9 @@ def main(cmdargs):
cosmo,
max_num_spec=args.nspec,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
for deltas in data.values():
for delta in deltas:
delta.fname = 'D2'
Expand Down
19 changes: 18 additions & 1 deletion bin/picca_xcf.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -131,6 +131,21 @@ def main(cmdargs):
'absorber redshift and the object redshift smaller than '
'z-cut-max'))

parser.add_argument(
'--z-min-sources',
type=float,
default=0.,
required=False,
help=('Min redshift of the background quasars'))

parser.add_argument(
'--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--lambda-abs',
type=str,
Expand Down Expand Up @@ -293,7 +308,9 @@ def main(cmdargs):
max_num_spec=args.nspec,
no_project=args.no_project,
nproc=args.nproc,
rebin_factor=args.rebin_factor)
rebin_factor=args.rebin_factor,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
xcf.data = data
xcf.num_data = num_data
userprint("")
Expand Down
20 changes: 19 additions & 1 deletion bin/picca_xcf1d.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -110,6 +110,22 @@ def main(cmdargs):
'absorber redshift and the object redshift smaller than '
'z-cut-max'))

parser.add_argument(
'--z-min-sources',
type=float,
default=0.,
required=False,
help=('Limit the minimum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--z-max-sources',
type=float,
default=10.,
required=False,
help=('Limit the maximum redshift of the quasars '
'used as sources for spectra'))

parser.add_argument(
'--lambda-abs',
type=str,
Expand Down Expand Up @@ -199,7 +215,9 @@ def main(cmdargs):
args.z_ref,
cosmo=None,
max_num_spec=args.nspec,
no_project=args.no_project)
no_project=args.no_project,
z_min_qso=args.z_min_sources,
z_max_qso=args.z_max_sources)
xcf.data = data
xcf.num_data = num_data
sys.stderr.write("\n")
Expand Down
Loading