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

Remodel Galaxy WL object structure #91

Open
caioolivv opened this issue May 3, 2023 · 1 comment
Open

Remodel Galaxy WL object structure #91

caioolivv opened this issue May 3, 2023 · 1 comment

Comments

@caioolivv
Copy link
Member

Currently WL analysis is done by creating GalaxyWL objects containing GalaxyWLDist and GalaxyRedshift objects. Since shape and position information are grouped in the same object, separated from redshift data, binning processes are not straight forward. It also presents a difficulty if we were to hierarchically sample information for kernel density estimation and forward modelling. These problems would be avoided if we rewrote the WL likelihood the objects in the following structure

  • NcGalaxyWL object which would contain the following objects:
    • NcGalaxySampleDistributionPosition with information on image position and redshift;
    • NcGalaxySampleDistributionShape with information on intrinsic shape, observed shape and the shear mapping (this might also be capable of holding information on measuring error);
    • NcGalaxySampleDistributionZProxy with information on proxy redshift;

Mathematically, the above objects would represent the likelihood

$$\mathcal{P} (\epsilon, r, z_p) = \int dz dr \mathcal{P}(\epsilon | z, r, \theta) \mathcal{P}(z_p | z, \theta ) \mathcal{P} (r | \theta) \mathcal{P} (z | \theta),$$

with the following map

  • NcGalaxyWL -> $\mathcal{P} (\epsilon, r, z_p)$;
  • NcGalaxySampleDistributionPosition -> $\mathcal{P} (r | \theta) \mathcal{P} (z | \theta)$;
  • NcGalaxySampleDistributionShape -> $\mathcal{P}(\epsilon | z, r, \theta)$;
  • NcGalaxySampleDistributionZProxy -> $\mathcal{P}(z_p | z, \theta )$;

Also, we might also use GSD instead of GalaxySampleDistribution.

@vitenti
Copy link
Contributor

vitenti commented May 17, 2023

The object structure is fine. The name can be GalaxySDPosition that is instead of GalaxySampleDistribution or GSD we can use GalaxySD. What do you think?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
NumCosmo 1.0
  
To do
Development

No branches or pull requests

2 participants