def test__tag(self):
hyper = al.SetupHyper(
hyper_galaxies_lens=True,
hyper_background_noise=al.hyper_data.HyperBackgroundNoise,
hyper_image_sky=al.hyper_data.HyperImageSky,
)
setup_mass = al.SetupMassLightDark(align_bulge_dark_centre=True)
setup = al.SetupPipeline(setup_hyper=hyper, setup_mass=setup_mass)
assert (
setup.tag == "setup__"
"hyper[galaxies_lens__bg_sky__bg_noise]__"
"mass[light_dark__bulge_sersic__disk_exp__mlr_free__dark_nfw_ludlow__with_shear__align_bulge_dark_centre]"
)
setup_source = al.SetupSourceInversion(
pixelization_prior_model=al.pix.Rectangular,
regularization_prior_model=al.reg.Constant,
)
setup_light = al.SetupLightParametric(light_centre=(1.0, 2.0))
setup_mass = al.SetupMassLightDark(mass_centre=(3.0, 4.0),
with_shear=False)
setup = al.SetupPipeline(setup_source=setup_source,
setup_light=setup_light,
setup_mass=setup_mass)
assert (
setup.tag == "setup__"
"light[parametric__bulge_sersic__disk_exp__align_bulge_disk_centre__centre_(1.00,2.00)]__"
"mass[light_dark__bulge_sersic__disk_exp__mlr_free__dark_nfw_ludlow__no_shear__centre_(3.00,4.00)]__"
"source[inversion__pix_rect__reg_const]")
setup_mass = al.SetupMassLightDark(align_bulge_dark_centre=True)
setup = al.SetupPipeline(setup_mass=setup_mass)
assert (
setup.tag == "setup__"
"mass[light_dark__bulge_sersic__disk_exp__mlr_free__dark_nfw_ludlow__with_shear__align_bulge_dark_centre]"
)
smbh = al.SetupSMBH(smbh_centre_fixed=True)
setup_subhalo = al.SetupSubhalo(
subhalo_instance=al.mp.SphericalNFWMCRLudlow(centre=(1.0, 2.0),
mass_at_200=1e8))
setup = al.SetupPipeline(setup_smbh=smbh, setup_subhalo=setup_subhalo)
assert (
setup.tag == "setup__"
"smbh[point_mass__centre_fixed]__"
"subhalo[nfw_sph_ludlow__mass_is_model__source_is_model__grid_5__centre_(1.00,2.00)__mass_1.0e+08]"
)
def test__mass_tag(self):
slam = al.SLaM(
pipeline_source_parametric=al.SLaMPipelineSourceParametric(),
pipeline_light_parametric=al.SLaMPipelineLightParametric(),
pipeline_mass=al.SLaMPipelineMass(),
)
assert (
slam.mass_tag == f"mass__"
f"light[parametric__bulge_sersic__disk_exp__align_bulge_disk_centre]__"
f"mass[total__power_law__with_shear]__"
f"source[parametric__bulge_sersic]")
pipeline_source_parametric = al.SLaMPipelineSourceParametric(
setup_light=al.SetupLightParametric(
bulge_prior_model=al.lp.SphericalExponential,
disk_prior_model=None,
align_bulge_disk_centre=False,
),
setup_mass=al.SetupMassTotal(
mass_prior_model=al.mp.EllipticalPowerLaw,
mass_centre=(0.0, 0.0)),
setup_source=al.SetupSourceParametric(
bulge_prior_model=al.lp.SphericalDevVaucouleurs),
)
pipeline_source_inversion = al.SLaMPipelineSourceInversion(
setup_source=al.SetupSourceInversion(
pixelization_prior_model=al.pix.VoronoiMagnification,
regularization_prior_model=al.reg.AdaptiveBrightness,
))
pipeline_light_parametric = al.SLaMPipelineLightParametric(
setup_light=al.SetupLightParametric(
bulge_prior_model=al.lp.SphericalDevVaucouleurs,
disk_prior_model=al.lp.SphericalExponential,
align_bulge_disk_centre=False,
light_centre=(0.0, 0.0),
))
pipeline_mass = al.SLaMPipelineMass(setup_mass=al.SetupMassLightDark(
bulge_prior_model=al.lmp.EllipticalSersicRadialGradient))
slam = al.SLaM(
pipeline_source_parametric=pipeline_source_parametric,
pipeline_source_inversion=pipeline_source_inversion,
pipeline_light_parametric=pipeline_light_parametric,
pipeline_mass=pipeline_mass,
)
assert (
slam.mass_tag == f"mass__"
f"light[parametric__bulge_dev_sph__disk_exp_sph__centre_(0.00,0.00)]__"
f"mass[light_dark__bulge_sersic_grad__disk_exp__mlr_free__dark_nfw_ludlow__with_shear]__"
f"source[inversion__pix_voro_mag__reg_adapt_bright]")
def test__source_inversion_tag(self):
slam = al.SLaM(
pipeline_source_parametric=al.SLaMPipelineSourceParametric(),
pipeline_source_inversion=al.SLaMPipelineSourceInversion(),
pipeline_mass=al.SLaMPipelineMass(),
)
assert (slam.source_inversion_tag == f"source__"
f"mass[total__sie__with_shear]__"
f"source[inversion__pix_rect__reg_const]")
assert slam.source_inversion_tag == slam.source_tag
pipeline_source_parametric = al.SLaMPipelineSourceParametric(
setup_light=al.SetupLightParametric(
bulge_prior_model=al.lp.SphericalExponential,
disk_prior_model=None),
setup_mass=al.SetupMassTotal(
mass_prior_model=al.mp.EllipticalPowerLaw,
mass_centre=(0.0, 0.0)),
setup_source=al.SetupSourceParametric(
bulge_prior_model=al.lp.SphericalDevVaucouleurs),
)
pipeline_source_inversion = al.SLaMPipelineSourceInversion(
setup_source=al.SetupSourceInversion(
pixelization_prior_model=al.pix.VoronoiMagnification,
regularization_prior_model=al.reg.AdaptiveBrightness,
))
slam = al.SLaM(
pipeline_source_parametric=pipeline_source_parametric,
pipeline_source_inversion=pipeline_source_inversion,
pipeline_light_parametric=al.SLaMPipelineLightParametric(),
pipeline_mass=al.SLaMPipelineMass(),
)
assert (slam.source_inversion_tag == f"source__"
f"light[parametric__bulge_exp_sph]__"
f"mass[total__power_law__with_shear__centre_(0.00,0.00)]__"
f"source[inversion__pix_voro_mag__reg_adapt_bright]")
assert slam.source_inversion_tag == slam.source_tag
By default, this again assumes _EllipticalIsothermal_ profile for the lens galaxy's mass model.
For this runner the _SLaMPipelineSourceInversion_ customizes:
- The _Pixelization_ used by the _Inversion_ of this pipeline.
- The _Regularization_ scheme used by of this pipeline.
The _SLaMPipelineSourceInversion_ use's the _SetupMass_ of the _SLaMPipelineSourceParametric_.
The _SLaMPipelineSourceInversion_ determines the source model used in the _SLaMPipelineLight_ and _SLaMPipelineMass_ pipelines, which in this
example therefore both use an _Inversion_.
"""
setup_source = al.SetupSourceInversion(
pixelization_prior_model=al.pix.VoronoiBrightnessImage,
regularization_prior_model=al.reg.AdaptiveBrightness)
pipeline_source_inversion = al.SLaMPipelineSourceInversion(
setup_source=setup_source)
# %%
"""
__SLaMPipelineMassTotal__
The _SLaMPipelineMassTotal_ pipeline fits the model for the lens galaxy's total mass distribution.
A full description of all options can be found ? and ?.
The model used to represent the lens galaxy's mass is input into _SLaMPipelineMassTotal_ and this runner uses the default of an
_EllipticalPowerLaw_ in this example.
"""
# %%
setup_mass = al.SetupMassTotal(with_shear=True)
# %%
"""
We also use the `SetupSourceInversion` object to customize the `Inversion` used for the source, specifically:
- The `Pixelization` used by the `Inversion` of this pipeline.
- The `Regularization` scheme used by the `Inversion` of this pipeline.
"""
# %%
setup_source = al.SetupSourceInversion(
pixelization_prior_model=al.pix.VoronoiMagnification,
regularization_prior_model=al.reg.Constant,
)
setup = al.SetupPipeline(
path_prefix=path.join("howtolens", "c4_t8_inversion"),
setup_mass=setup_mass,
setup_source=setup_source,
)
# %%
"""
__Pipeline Creation__
To create a pipeline we import it from the pipelines folder and run its `make_pipeline` function, inputting the
`Setup` and `SettingsPhase` above.
"""
