Source code for ytree.frontends.lhalotree.fields

LHaloTreeArbor fields


# Copyright (c) ytree development team. All rights reserved.
# Distributed under the terms of the Modified BSD License.
# The full license is in the file COPYING.txt, distributed with this software.

import numpy as np

from ytree.data_structures.fields import \
from ytree.frontends.lhalotree.utils import \

m_unit = "Msun"
p_unit = "unitary"
r_unit = "kpc"
v_unit = "km/s"

# This is where we set aliases for LHaloTree field names to univeral
# field names.
# Syntax is: (<universal field name>, <field name on disk>, <units>)
# TODO: Use units from parameter file

id_type = np.int64

[docs]class LHaloTreeFieldInfo(FieldInfoContainer): alias_fields = ( # ("uid", "uid", None), # ("desc_uid", "desc_uid", None), # ("scale_factor", "scale_factor", None), ("mass", "Mvir", m_unit), ("virial_mass", "Mvir", m_unit), # ("virial_radius", "Rvir", r_unit), # ("scale_radius", "rs", r_unit), ("velocity_dispersion", "VelDisp", v_unit), ("position_x", "x", p_unit), ("position_y", "y", p_unit), ("position_z", "z", p_unit), ("velocity_x", "vx", v_unit), ("velocity_y", "vy", v_unit), ("velocity_z", "vz", v_unit), ("angular_momentum_x", "Jx", None), ("angular_momentum_y", "Jy", None), ("angular_momentum_z", "Jz", None), # ("spin_parameter", "Spin", None), ) data_types = tuple( [("uid", id_type), ("desc_uid", id_type)] + [(item[0], np.dtype(item[1])) for item in dtype_header_default])