Ctrl+K

pynbody.analysis.decomp.decomp

Contents

pynbody.analysis.decomp.decomp#

pynbody.analysis.decomp.decomp(h, aligned=False, j_disk_min=0.8, j_disk_max=1.1, E_cut=None, j_circ_from_r=False, log_interp=False, angmom_size='3 kpc')[source]#

Creates an array ‘decomp’ for star particles in the simulation, with an integer specifying components.

The possible values of the components are:

1 – thin disk

2 – halo

3 – bulge

4 – thick disk

5 – pseudo bulge

This routine is based on an original IDL procedure by Chris Brook.

Parameters:

  • h (SimSnap) – The simulation snapshot to analyze

  • aligned (bool) – If True, the simulation is assumed to be already aligned so that the disk is in the xy plane. Otherwise, the simulation is recentered and aligned into the xy plane.

  • j_disk_min (float) – The minimum angular momentum as a proportion of the circular angular momentum which a particle must have to be part of the ‘disk’.

  • j_disk_max (float) – The maximum angular momentum as a proportion of the circular angular momentum which a particle can have to be part of the ‘disk’.

  • E_cut (float) – The energy boundary between bulge and spheroid. If None, this is taken to be the median energy of the stars.

  • j_circ_from_r (bool) – If True, the maximum angular momentum is determined as a function of radius, rather than as a function of orbital energy. Default False (determine as function of energy).

  • angmom_size (str) – The size of the disk to use for calculating the angular momentum vector. Default is “3 kpc”.

Contents