def wrap_batch_gbtxwigglez_data_run(inifile, generate=False,
outdir="./plots/"):
r"""Wrapper to the GBT x WiggleZ calculation"""
params_init = {"gbt_mapkey": "cleaned GBT map",
"wigglez_deltakey": "WiggleZ overdensity map",
"wigglez_mockkey": "WiggleZ overdensities from mocks",
"wigglez_selectionkey": "WiggleZ selection function",
"mode_transfer_1d_ini": "ini file -> 1d trans. function",
"mode_transfer_2d_ini": "ini file -> 2d trans. function",
"beam_transfer_ini": "ini file -> 2d beam trans. function",
"spec_ini": "ini file for the spectral estimation",
"output_tag": "tag identifying the output somehow"}
prefix = "cwx_"
params = parse_ini.parse(inifile, params_init, prefix=prefix)
print params
output_tag = "%s_%s" % (params['gbt_mapkey'], params['output_tag'])
output_root = "%s/%s/" % (outdir, output_tag)
if generate:
output_tag = None
print output_root
print output_tag
file_tools.mkparents(output_root)
parse_ini.write_params(params, output_root + 'params.ini',
prefix=prefix)
datapath_db = data_paths.DataPath()
mode_transfer_1d = None
if params["mode_transfer_1d_ini"]:
mode_transfer_1d = cct.wrap_batch_crosspwr_transfer(
params["mode_transfer_1d_ini"],
generate=generate,
outdir=outdir)
batch_gbtxwigglez_data_run(params["gbt_mapkey"],
params["wigglez_deltakey"],
params["wigglez_mockkey"],
params["wigglez_selectionkey"],
inifile=params["spec_ini"],
datapath_db=datapath_db,
outdir=output_root,
output_tag=output_tag,
beam_transfer=None,
mode_transfer_1d=mode_transfer_1d,
mode_transfer_2d=None,
theory_curve=None)
def wrap_batch_gbtxwigglez_data_run(inifile,
generate=False,
outdir="./plots/"):
r"""Wrapper to the GBT x WiggleZ calculation"""
params_init = {
"gbt_mapkey": "cleaned GBT map",
"wigglez_deltakey": "WiggleZ overdensity map",
"wigglez_mockkey": "WiggleZ overdensities from mocks",
"wigglez_selectionkey": "WiggleZ selection function",
"mode_transfer_1d_ini": "ini file -> 1d trans. function",
"mode_transfer_2d_ini": "ini file -> 2d trans. function",
"beam_transfer_ini": "ini file -> 2d beam trans. function",
"spec_ini": "ini file for the spectral estimation",
"output_tag": "tag identifying the output somehow"
}
prefix = "cwx_"
params = parse_ini.parse(inifile, params_init, prefix=prefix)
print params
output_tag = "%s_%s" % (params['gbt_mapkey'], params['output_tag'])
output_root = "%s/%s/" % (outdir, output_tag)
if generate:
output_tag = None
print output_root
print output_tag
file_tools.mkparents(output_root)
parse_ini.write_params(params, output_root + 'params.ini', prefix=prefix)
datapath_db = data_paths.DataPath()
mode_transfer_1d = None
if params["mode_transfer_1d_ini"]:
mode_transfer_1d = cct.wrap_batch_crosspwr_transfer(
params["mode_transfer_1d_ini"], generate=generate, outdir=outdir)
batch_gbtxwigglez_data_run(params["gbt_mapkey"],
params["wigglez_deltakey"],
params["wigglez_mockkey"],
params["wigglez_selectionkey"],
inifile=params["spec_ini"],
datapath_db=datapath_db,
outdir=output_root,
output_tag=output_tag,
beam_transfer=None,
mode_transfer_1d=mode_transfer_1d,
mode_transfer_2d=None,
theory_curve=None)
def wrap_batch_gbtpwrspec_data_run(inifile, generate=False,
outdir="./plots/"):
r"""Wrapper to the GBT x GBT calculation"""
params_init = {"gbt_mapkey": "cleaned GBT map",
"mode_transfer_1d_ini": "ini file -> 1d trans. function",
"mode_transfer_2d_ini": "ini file -> 2d trans. function",
"beam_transfer_ini": "ini file -> 2d beam trans. function",
"square_1dmodetrans": False,
"spec_ini": "ini file for the spectral estimation",
"output_tag": "tag identifying the output somehow"}
prefix="cp_"
params = parse_ini.parse(inifile, params_init, prefix=prefix)
print params
output_tag = "%s_%s" % (params['gbt_mapkey'], params['output_tag'])
output_root = "%s/%s/" % (outdir, output_tag)
if generate:
output_tag = None
print output_root
print output_tag
file_tools.mkparents(output_root)
parse_ini.write_params(params, output_root + 'params.ini',
prefix=prefix)
datapath_db = data_paths.DataPath()
mode_transfer_1d=None
if params["mode_transfer_1d_ini"]:
mode_transfer_1d = cct.wrap_batch_crosspwr_transfer(
params["mode_transfer_1d_ini"],
generate=generate,
outdir=outdir)
return batch_gbtpwrspec_data_run(params["gbt_mapkey"],
inifile=params["spec_ini"],
datapath_db=datapath_db,
outdir=output_root,
output_tag=output_tag,
beam_transfer=None,
square_1dmodetrans = params["square_1dmodetrans"],
mode_transfer_1d=mode_transfer_1d,
mode_transfer_2d=None)
def wrap_batch_gbtpwrspec_data_run(inifile, generate=False, outdir="./plots/"):
r"""Wrapper to the GBT x GBT calculation"""
params_init = {
"gbt_mapkey": "cleaned GBT map",
"mode_transfer_1d_ini": "ini file -> 1d trans. function",
"mode_transfer_2d_ini": "ini file -> 2d trans. function",
"beam_transfer_ini": "ini file -> 2d beam trans. function",
"square_1dmodetrans": False,
"spec_ini": "ini file for the spectral estimation",
"output_tag": "tag identifying the output somehow"
}
prefix = "cp_"
params = parse_ini.parse(inifile, params_init, prefix=prefix)
print params
output_tag = "%s_%s" % (params['gbt_mapkey'], params['output_tag'])
output_root = "%s/%s/" % (outdir, output_tag)
if generate:
output_tag = None
print output_root
print output_tag
file_tools.mkparents(output_root)
parse_ini.write_params(params, output_root + 'params.ini', prefix=prefix)
datapath_db = data_paths.DataPath()
mode_transfer_1d = None
if params["mode_transfer_1d_ini"]:
mode_transfer_1d = cct.wrap_batch_crosspwr_transfer(
params["mode_transfer_1d_ini"], generate=generate, outdir=outdir)
return batch_gbtpwrspec_data_run(
params["gbt_mapkey"],
inifile=params["spec_ini"],
datapath_db=datapath_db,
outdir=output_root,
output_tag=output_tag,
beam_transfer=None,
square_1dmodetrans=params["square_1dmodetrans"],
mode_transfer_1d=mode_transfer_1d,
mode_transfer_2d=None)
