mirror/sims_python_files
1
0
Fork 0
mirror of https://github.com/glatterf42/sims_python_files.git synced 2024-09-18 15:03:45 +02:00
Code Releases Activity

Added cumulative mass profile plotting.

Browse source
This commit is contained in:
glatterf42 2022-04-27 11:50:54 +02:00
parent c92465a95a
commit 2ba0971a3c
1 changed files with 8 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
get_mass_zoom_halo.py
View file

@ -14,9 +14,9 @@ import matplotlib.pyplot as plt
def V(r):
return 4 * np.pi * r**3 / 3
directory = Path(r"/home/ben/sims/swiftsim/examples/zoom_tests/")
directory = Path(r"/home/ben/sims/swiftsim/examples/zoom_tests/auriga6_halo_arj")
snap_number = 7
snap_number = 0 # 7 for our tests
fof_file = h5py.File(directory / f"fof_output_000{snap_number}.hdf5", "r")
file = h5py.File(directory / f'output_000{snap_number}.hdf5', 'r')
@ -49,7 +49,8 @@ number_of_members = groups['Sizes'][:]
table_width = 4
separate_unique_counter = 0
for i in range(len(groupids)-1):
# for i in range(len(groupids)-1):
for i in range(11):
if np.isin(groupids[i], unique_groups):
highres_members = particle_count[separate_unique_counter]
contamination = (1 - highres_members / number_of_members[i])
@ -96,7 +97,7 @@ Lbox = 100
softening = Lbox / Nres / 30
plt.axvline(4 * softening / group_radius, linestyle='--', color='grey')
plt.title('Density profile Auriga 6')
plt.title('Density profile Auriga 6 ARJ')
plt.xlabel(r'R / $\mathrm{R}_\mathrm{group}$')
plt.ylabel(r'ρ [$10^{10}\ \mathrm{M}_\odot\ /\ \mathrm{Mpc}^3$]')
@ -104,6 +105,9 @@ plt.loglog(plot_log_radial_bins, masses_in_radial_bins)
plt.show()
# kumulatives Dichteprofil -> sollte bei grossen Massen praktisch gleich sein, am besten direkt uebereinander plotten.
# plt.hist2D mit x und y, jeweils Zentrum der Gruppe +- 0.5 Mpc anschauen, alles andere wegwerfen, sollte direkt Struktur erkennen koennen, ob sie gleich ist.