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. """