csiborgtools/scripts/plot_concentration.ipynb
Richard Stiskalek 7f58b1f78c
JAX and fix conc (#6)
* change to log10 initlogRs

* add uncertainty

* add check if hessian negative

* update TODO

* update TODO

* output the error too

* save e_logRs

* add concentration calculation

* calcul concentration

* missed comma in a list

* Update TODO

* fix bug

* add box units and pretty status

* make uncertainty optional

* add BIC function

* remove BIC again

* add new overdensity calculation

* change defualt values to max dist and mass

* change to r200

* new halo find

* speed up the calculation

* add commented fucn

* update TODO

* add check whether too close to outside boundary

* Update TODO

* extract velocity as well

* calculate m200 and m500

* update nb

* update TODO
2022-11-05 21:17:05 +00:00

340 KiB

In [2]:
import numpy as np
import healpy as hp
%matplotlib notebook
import matplotlib.pyplot as plt
# Local imports
try:
    import csiborgtools
except ModuleNotFoundError:
    import sys
    sys.path.append("../")
    import csiborgtools
import utils

%load_ext autoreload
%autoreload 2
In [110]:
obs = utils.load_2mpp()


cols = ["ra", "dec", "mass_cl", "dist"]

mmains = utils.load_mmains(1)
mmains = csiborgtools.utils.list_to_ndarray(mmains, cols)
sim = mmains[0, ...]
sim = csiborgtools.utils.array_to_structured(sim ,cols)

planck = utils.load_planck2015()
100%|█████████████████████████████████████████████| 1/1 [00:00<00:00,  3.30it/s]
In [119]:
dx = 20
dmin = 125
dmax = dmin + dx

mask_obs = (dmin < obs["CDIST_CMB"]) & (obs["CDIST_CMB"] < dmax)
mask_sim = (dmin < sim["dist"]) & (sim["dist"] < dmax) & (sim["mass_cl"] > 1e12)

width = 6.4
plt.figure(figsize=(width, width*0.75))
plt.scatter(obs["RA"][mask_obs], obs["DEC"][mask_obs], s=1.5, label="2M++")
plt.scatter(sim["ra"][mask_sim] , sim["dec"][mask_sim], s=1.5, label="CSiBORG")
plt.scatter(planck["RA"], planck["DEC"], label="Planck SZ clusters < 214 MPc", c="red")


plt.legend()
plt.xlabel("RA")
plt.ylabel("dec")
plt.tight_layout()
# plt.savefig("../plots/2mpp_overlap.png", dpi=450)
plt.show()
No description has been provided for this image
In [ ]:
fpath = "/mnt/zfsusers/rstiskalek/csiborgtools/data/HFI_PCCS_SZ-union_R2.08.fits"
obs = csiborgtools.io.read_planck2015(fpath, cosmo, 214)


c1 = SkyCoord(ra=obs["RA"]*u.degree, dec=obs["DEC"]*u.degree, distance=obs["COMDIST"]*u.Mpc)
c2 = SkyCoord(ra=sim["ra"]*u.degree, dec=sim["dec"]*u.degree, distance=sim["dist"]*u.Mpc)

sep, indxs = csiborgtools.match.brute_spatial_separation(c1, c2, N=100, angular=False, verbose=True, )
In [ ]:

In [ ]:

In [ ]:

In [ ]:

In [ ]: