"""
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 \
FieldInfoContainer
from ytree.frontends.lhalotree.utils import \
dtype_header_default
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])