Stars¶
Particle Stars¶
A particle.Starsobject is the core object in synthesizer for modelling stellar distributions described by discrete particles (e.g. in hydrodynamical simulations).
To initialise one you need equal length arrays describing the ages, metallicities and initial masses of your star particles. These can then be provided directly to the Stars object on initialisation. In the example below, we make some mock data (with units) and load these into a particle Stars object:
[1]:
import numpy as np
from unyt import Msun, Myr, Mpc
from synthesizer.particle import Stars
N = 100
ages = np.random.rand(N) * 100 * Myr
metallicities = np.random.rand(N) / 100
initial_masses = np.ones(N) * 1e6 * Msun
stars = Stars(
initial_masses=initial_masses,
ages=ages,
metallicities=metallicities,
)
The provided ages, metallicities and initial masses of the star particles describe the star formation and metal enrichment (SFZH) history. Some basic statistics of this object can be viewed by:
[2]:
print(stars)
+----------------------------------------------------------------------------------+
| PARTICLES |
+---------------------------+------------------------------------------------------+
| Attribute | Value |
+---------------------------+------------------------------------------------------+
| nparticles | 100 |
+---------------------------+------------------------------------------------------+
| metallicity_floor | 1.00e-05 |
+---------------------------+------------------------------------------------------+
| name | 'Stars' |
+---------------------------+------------------------------------------------------+
| component_type | 'Stars' |
+---------------------------+------------------------------------------------------+
| resampled | False |
+---------------------------+------------------------------------------------------+
| nstars | 100 |
+---------------------------+------------------------------------------------------+
| ages (100,) | 1.59e+05 yr -> 9.84e+07 yr (Mean: 4.88e+07 yr) |
+---------------------------+------------------------------------------------------+
| metallicities (100,) | 2.81e-04 -> 9.98e-03 (Mean: 5.20e-03) |
+---------------------------+------------------------------------------------------+
| initial_masses (100,) | 1.00e+06 Msun -> 1.00e+06 Msun (Mean: 1.00e+06 Msun) |
+---------------------------+------------------------------------------------------+
| log10ages (100,) | 5.20e+00 -> 7.99e+00 (Mean: 7.58e+00) |
+---------------------------+------------------------------------------------------+
| log10metallicities (100,) | -3.55e+00 -> -2.00e+00 (Mean: -2.40e+00) |
+---------------------------+------------------------------------------------------+
This can be visualised in synthesizer using a Grid object (see here), to assign these to SFZH bins.
[3]:
from synthesizer.grid import Grid
grid_dir = "../../../tests/test_grid"
grid_name = "test_grid"
grid = Grid(grid_name, grid_dir=grid_dir)
stars.plot_sfzh(grid)
[3]:
(<Figure size 350x350 with 3 Axes>,
<Axes: xlabel='$\\log_{10}(\\mathrm{age}/\\mathrm{yr})$', ylabel='$\\log_{10}Z$'>)
By providing a Grid to the Stars object one can generate the stellar emission. This can include the integrated SED (see Sed docs), line emission (see line docs), as well as the per–particle SED and line luminosity / equivalent widths.
Optional arguments¶
Any component can take arbitrary extra arguments. Some of these are pre-defined kwargs (e.g. coordinates, redshift, current_masses etc.), but any additional **kwargs will be unpacked and attached to the component instances (here a Stars instance). This is particularly useful when working with emission models that require component level attributes (e.g. particle specific optical depths).
[4]:
tau_v = np.random.rand(N)
coordinates = np.random.normal(50, 0.005, (N, 3)) * Mpc
centre = np.array([50, 50, 50]) * Mpc
redshift = 5
stars = Stars(
initial_masses=initial_masses,
ages=ages,
metallicities=metallicities,
# Pre defined Quantitys
coordinates=coordinates,
centre=centre,
# Optional argument
redshift=redshift,
current_masses=initial_masses - (np.random.rand(N) * 1e6 * Msun),
# Arbitrary kwarg for modelling
tau_v=tau_v,
)
print(stars)
+--------------------------------------------------------------------------------------+
| PARTICLES |
+-------------------------------+------------------------------------------------------+
| Attribute | Value |
+-------------------------------+------------------------------------------------------+
| redshift | 5 |
+-------------------------------+------------------------------------------------------+
| nparticles | 100 |
+-------------------------------+------------------------------------------------------+
| metallicity_floor | 1.00e-05 |
+-------------------------------+------------------------------------------------------+
| name | 'Stars' |
+-------------------------------+------------------------------------------------------+
| component_type | 'Stars' |
+-------------------------------+------------------------------------------------------+
| resampled | False |
+-------------------------------+------------------------------------------------------+
| nstars | 100 |
+-------------------------------+------------------------------------------------------+
| coordinates (100, 3) | 5.00e+01 Mpc -> 5.00e+01 Mpc (Mean: 5.00e+01 Mpc) |
+-------------------------------+------------------------------------------------------+
| masses (100,) | 7.41e+02 Msun -> 9.83e+05 Msun (Mean: 5.10e+05 Msun) |
+-------------------------------+------------------------------------------------------+
| centre | [50 50 50] Mpc |
+-------------------------------+------------------------------------------------------+
| tau_v (100,) | 9.71e-03 -> 9.75e-01 (Mean: 5.26e-01) |
+-------------------------------+------------------------------------------------------+
| ages (100,) | 1.59e+05 yr -> 9.84e+07 yr (Mean: 4.88e+07 yr) |
+-------------------------------+------------------------------------------------------+
| metallicities (100,) | 2.81e-04 -> 9.98e-03 (Mean: 5.20e-03) |
+-------------------------------+------------------------------------------------------+
| initial_masses (100,) | 1.00e+06 Msun -> 1.00e+06 Msun (Mean: 1.00e+06 Msun) |
+-------------------------------+------------------------------------------------------+
| current_masses (100,) | 7.41e+02 Msun -> 9.83e+05 Msun (Mean: 5.10e+05 Msun) |
+-------------------------------+------------------------------------------------------+
| centered_coordinates (100, 3) | -1.08e-02 Mpc -> 1.39e-02 Mpc (Mean: 1.99e-04 Mpc) |
+-------------------------------+------------------------------------------------------+
| log10ages (100,) | 5.20e+00 -> 7.99e+00 (Mean: 7.58e+00) |
+-------------------------------+------------------------------------------------------+
| log10metallicities (100,) | -3.55e+00 -> -2.00e+00 (Mean: -2.40e+00) |
+-------------------------------+------------------------------------------------------+
Computing characteristics radii¶
Characteristic radii can be computed for the stellar particle distribution. We can either calculate the half-mass radius.
[5]:
r = stars.get_half_mass_radius()
print(r)
0.007610819678567318 Mpc
Or use an aribtrary attribute for the radius weight. Here we’ll use the optical depth for a “half tau_V radius”.
[6]:
r = stars.get_attr_radius("tau_v")
print(r)
0.0074253486264112475 Mpc
We can also compute radii for different fractions by passing the fraction we want to calculate a radius for.
[7]:
r20 = stars.get_attr_radius("masses", frac=0.2)
r80 = stars.get_attr_radius("masses", frac=0.8)
print(r20, r80)
0.00501872865379495 Mpc 0.011212249232457767 Mpc
We can also calculate characteristic light radii from photometry. See the photometry documentation for more details.
Parametric Stars¶
The core element of the parametric modelling in synthesizer is the parametric.Stars object. This contains everything necessary to describe a stellar population including:
A 2D SFZH grid containing the stellar mass formed in stellar age and metallicity bins.
A
Morphologyobject describing it’s distribution (optional but required for imaging).The axes of the SFZH grid.
The total initial stellar mass.
And if calculated by the user:
There are a number of different methods to define a parametric stellar population (from here on a Stars object). In what follows we will focus on user defined values/arrays but note that each different method could instead use simulated data (e.g. from Semi Analytic Models) rather than user defined values.
The possible ways of defining a Stars object are (the age and metallicity axis of the SFZH must always be supplied):
Explictly providing the SFZH grid itself. Note that this will mean all other definitions below are ignored.
Providing singular age and metallicity values to define an “instantaneous” SFZH, i.e. one where all the stellar mass is in a single age and metallicity bin.
Providing arrays which contain the SFH and metallicity distribution in terms of stellar mass formed.
Providing SFH and metallcity distribution functions (from
parametric.sf_histandparametric.metal_dist). These will be used internally to compute the arrays from the previous point.Any combination of the above.
Below we simply import some packages and objects we’ll need.
[8]:
import numpy as np
from unyt import Myr
from synthesizer.parametric import SFH, Stars, ZDist
Defining an instantaneous SFZH¶
To define an instantaneous SFZH we need only to provide the age and metallicity of the SFZH bin, alongside the SFZH axis arrays and the total initial mass of the Stars object. Note that the log10ages and log10metallicites/metallicities will likely come from a Grid object in most workflows (see the Grid docs).
[9]:
# Create the SFZH axis arrays
log10ages = np.arange(6.0, 10.5, 0.1) # log10(age/yr)
log10metallicities = np.arange(-5.0, -1.5, 0.25)
metallicities = 10**log10metallicities
# Get the instantaneous Stars
inst_stars = Stars(
log10ages,
metallicities,
sf_hist=100 * Myr,
metal_dist=10**-3,
initial_mass=10**9 * Msun,
)
To see the created SFZH you can print some basic statistics by printing the Stars object.
[10]:
print(inst_stars)
+---------------------------------------------------------------------------+
| STARS |
+--------------------------+------------------------------------------------+
| Attribute | Value |
+--------------------------+------------------------------------------------+
| component_type | 'Stars' |
+--------------------------+------------------------------------------------+
| metallicity_grid_type | 'log10Z' |
+--------------------------+------------------------------------------------+
| log10ages_lims | [6.0, 10.399999999999984, ] |
+--------------------------+------------------------------------------------+
| metallicities_lims | [1e-05, 0.01778279410038923, ] |
+--------------------------+------------------------------------------------+
| log10metallicities_lims | [-5.0, -1.75, ] |
+--------------------------+------------------------------------------------+
| ages (45,) | 1.00e+06 yr -> 2.51e+10 yr (Mean: 2.71e+09 yr) |
+--------------------------+------------------------------------------------+
| metallicities (14,) | 1.00e-05 -> 1.78e-02 (Mean: 2.90e-03) |
+--------------------------+------------------------------------------------+
| sf_hist (45,) | 0.00e+00 -> 1.00e+09 (Mean: 2.22e+07) |
+--------------------------+------------------------------------------------+
| metal_dist (14,) | 0.00e+00 -> 1.00e+09 (Mean: 7.14e+07) |
+--------------------------+------------------------------------------------+
| initial_mass | 999999999.9999999 Msun |
+--------------------------+------------------------------------------------+
| sfzh (45, 14) | 0.00e+00 -> 1.00e+09 (Mean: 1.59e+06) |
+--------------------------+------------------------------------------------+
| log10ages (45,) | 6.00e+00 -> 1.04e+01 (Mean: 8.20e+00) |
+--------------------------+------------------------------------------------+
| log10metallicities (14,) | -5.00e+00 -> -1.75e+00 (Mean: -3.38e+00) |
+--------------------------+------------------------------------------------+
Or to probe the SFZH you can plot it using the provided helper function.
[11]:
fig, ax = inst_stars.plot_sfzh(show=True)
Defining a SFZH from arrays¶
Of course, most of the time you won’t want an instantaneous SFZH but instead a distribution. One way to achieve this is to pass arrays describing the SFH and metallicity distribution. This is particularly useful if you have data explictly describing the ages and metallicities of a stellar population. However, here we’ll demonstrate this with entirely unphysical arrays.
[12]:
# Create arrays for SFH and ZDist
sf_hist = np.zeros_like(log10ages)
metal_dist = np.zeros_like(metallicities)
sf_hist[10:15] = 1
sf_hist[15:20] = 0.5
metal_dist[5:6] = 1
metal_dist[6:8] = 0.5
# Get the Stars
arr_stars = Stars(
log10ages,
metallicities,
sf_hist=sf_hist,
metal_dist=metal_dist,
initial_mass=10**9 * Msun,
)
print(arr_stars)
fig, ax = arr_stars.plot_sfzh(show=True)
+---------------------------------------------------------------------------+
| STARS |
+--------------------------+------------------------------------------------+
| Attribute | Value |
+--------------------------+------------------------------------------------+
| component_type | 'Stars' |
+--------------------------+------------------------------------------------+
| metallicity_grid_type | 'log10Z' |
+--------------------------+------------------------------------------------+
| log10ages_lims | [6.0, 10.399999999999984, ] |
+--------------------------+------------------------------------------------+
| metallicities_lims | [1e-05, 0.01778279410038923, ] |
+--------------------------+------------------------------------------------+
| log10metallicities_lims | [-5.0, -1.75, ] |
+--------------------------+------------------------------------------------+
| ages (45,) | 1.00e+06 yr -> 2.51e+10 yr (Mean: 2.71e+09 yr) |
+--------------------------+------------------------------------------------+
| metallicities (14,) | 1.00e-05 -> 1.78e-02 (Mean: 2.90e-03) |
+--------------------------+------------------------------------------------+
| sf_hist (45,) | 0.00e+00 -> 1.00e+00 (Mean: 1.67e-01) |
+--------------------------+------------------------------------------------+
| metal_dist (14,) | 0.00e+00 -> 1.00e+00 (Mean: 1.43e-01) |
+--------------------------+------------------------------------------------+
| initial_mass | 999999999.9999999 Msun |
+--------------------------+------------------------------------------------+
| sfzh (45, 14) | 0.00e+00 -> 6.67e+07 (Mean: 1.59e+06) |
+--------------------------+------------------------------------------------+
| log10ages (45,) | 6.00e+00 -> 1.04e+01 (Mean: 8.20e+00) |
+--------------------------+------------------------------------------------+
| log10metallicities (14,) | -5.00e+00 -> -1.75e+00 (Mean: -3.38e+00) |
+--------------------------+------------------------------------------------+
Defining a SFZH using functions¶
You won’t always have explict arrays to describe the age and metallicity distributions. For these situations synthesizer provides a suite of functions to describe the SFH and metallicity distribution.
SFH parametrisations can be found in the SFH module, while metallicity distribution parametrisations can be found in the ZDist module. To see a full list of available parametrisations just print the parametrisations list from these modules`.
[13]:
print(SFH.parametrisations)
print(ZDist.parametrisations)
('Constant', 'Gaussian', 'Exponential', 'TruncatedExponential', 'DecliningExponential', 'DelayedExponentialLogNormal', 'DoublePowerLaw', 'DenseBasis')
('DeltaConstant', 'Normal')
To use one of these parametrisations you have to first define an instance of a parametrisation and then pass it to the Stars object at instantiation. Below we demonstrate an example using a delta function for the metallicity distribution, and a constant SFH. Note that you can combine any type of SFH and metallicity distribution inputs.
[14]:
# Define a delta function for metallicity
metal_dist = ZDist.DeltaConstant(log10metallicity=-2.5)
print(metal_dist)
# Define a constant SFH
sfh = SFH.Constant(100 * Myr)
print(sfh)
# Create the Stars object
const_stars = Stars(
log10ages,
metallicities,
sf_hist=sfh,
metal_dist=metal_dist,
initial_mass=10**9 * Msun,
)
print(const_stars)
fig, ax = const_stars.plot_sfzh(show=True)
----------
SUMMARY OF PARAMETERISED METAL ENRICHMENT HISTORY
<class 'synthesizer.parametric.metal_dist.DeltaConstant'>
metallicity: None
log10metallicity: -2.5
----------
----------
SUMMARY OF PARAMETERISED STAR FORMATION HISTORY
<class 'synthesizer.parametric.sf_hist.Constant'>
min_age: 0 yr
max_age: 100 Myr
median age: 50.00 Myr
mean age: 50.00 Myr
----------
+-------------------------------------------------------------------------------------------------------+
| STARS |
+--------------------------+----------------------------------------------------------------------------+
| Attribute | Value |
+--------------------------+----------------------------------------------------------------------------+
| component_type | 'Stars' |
+--------------------------+----------------------------------------------------------------------------+
| sf_hist_func | <synthesizer.parametric.sf_hist.Constant object at 0x7f5a67279ed0> |
+--------------------------+----------------------------------------------------------------------------+
| metal_dist_func | <synthesizer.parametric.metal_dist.DeltaConstant object at 0x7f5a6727a1a0> |
+--------------------------+----------------------------------------------------------------------------+
| metallicity_grid_type | 'log10Z' |
+--------------------------+----------------------------------------------------------------------------+
| log10ages_lims | [6.0, 10.399999999999984, ] |
+--------------------------+----------------------------------------------------------------------------+
| metallicities_lims | [1e-05, 0.01778279410038923, ] |
+--------------------------+----------------------------------------------------------------------------+
| log10metallicities_lims | [-5.0, -1.75, ] |
+--------------------------+----------------------------------------------------------------------------+
| ages (45,) | 1.00e+06 yr -> 2.51e+10 yr (Mean: 2.71e+09 yr) |
+--------------------------+----------------------------------------------------------------------------+
| metallicities (14,) | 1.00e-05 -> 1.78e-02 (Mean: 2.90e-03) |
+--------------------------+----------------------------------------------------------------------------+
| sf_hist (45,) | 0.00e+00 -> 1.85e+08 (Mean: 2.22e+07) |
+--------------------------+----------------------------------------------------------------------------+
| metal_dist (14,) | 0.00e+00 -> 1.00e+09 (Mean: 7.14e+07) |
+--------------------------+----------------------------------------------------------------------------+
| initial_mass | 999999999.9999999 Msun |
+--------------------------+----------------------------------------------------------------------------+
| sfzh (45, 14) | 0.00e+00 -> 1.85e+08 (Mean: 1.59e+06) |
+--------------------------+----------------------------------------------------------------------------+
| log10ages (45,) | 6.00e+00 -> 1.04e+01 (Mean: 8.20e+00) |
+--------------------------+----------------------------------------------------------------------------+
| log10metallicities (14,) | -5.00e+00 -> -1.75e+00 (Mean: -3.38e+00) |
+--------------------------+----------------------------------------------------------------------------+
Combining Stars¶
We can also combine individual Stars objects together to produce more complicated star formation and metal enrichment histories by simply adding them.
[15]:
# Create an exponential Stars object
exp_stars = Stars(
log10ages,
metallicities,
sf_hist=sfh,
metal_dist=metal_dist,
initial_mass=10**9.5 * Msun,
)
# Create a log-normal Stars object
logn_stars = Stars(
log10ages,
metallicities,
sf_hist=sfh,
metal_dist=0.005,
initial_mass=10**9 * Msun,
)
combined = arr_stars + const_stars + exp_stars + logn_stars
print(combined)
fig, ax = combined.plot_sfzh()
+---------------------------------------------------------------------------+
| STARS |
+--------------------------+------------------------------------------------+
| Attribute | Value |
+--------------------------+------------------------------------------------+
| component_type | 'Stars' |
+--------------------------+------------------------------------------------+
| metallicity_grid_type | 'log10Z' |
+--------------------------+------------------------------------------------+
| log10ages_lims | [6.0, 10.399999999999984, ] |
+--------------------------+------------------------------------------------+
| metallicities_lims | [1e-05, 0.01778279410038923, ] |
+--------------------------+------------------------------------------------+
| log10metallicities_lims | [-5.0, -1.75, ] |
+--------------------------+------------------------------------------------+
| ages (45,) | 1.00e+06 yr -> 2.51e+10 yr (Mean: 2.71e+09 yr) |
+--------------------------+------------------------------------------------+
| metallicities (14,) | 1.00e-05 -> 1.78e-02 (Mean: 2.90e-03) |
+--------------------------+------------------------------------------------+
| sf_hist (45,) | 0.00e+00 -> 1.02e+09 (Mean: 1.37e+08) |
+--------------------------+------------------------------------------------+
| metal_dist (14,) | 0.00e+00 -> 4.16e+09 (Mean: 4.40e+08) |
+--------------------------+------------------------------------------------+
| initial_mass | 0.9999999999999999 Msun |
+--------------------------+------------------------------------------------+
| sfzh (45, 14) | 0.00e+00 -> 7.68e+08 (Mean: 9.78e+06) |
+--------------------------+------------------------------------------------+
| log10ages (45,) | 6.00e+00 -> 1.04e+01 (Mean: 8.20e+00) |
+--------------------------+------------------------------------------------+
| log10metallicities (14,) | -5.00e+00 -> -1.75e+00 (Mean: -3.38e+00) |
+--------------------------+------------------------------------------------+
Defining your own parametrisations¶
If synthesizer doesn’t already include a parametrisation you need, you can add a custom parametrisation by defining a new class.
The important requirements of a custom function are:
They inherit the
Commonparent class.They instantiate the parent class (
SFH/ZDist.Common) with the model parameters and a meaningful name.They implement a private
_sfrmethod which takes a single age and returns a single SFR weight in theSFHcase, or they implement a_weightmethod which takes a single metallicity and returns a single metallicity weight in theZDistcase.The
SFHmethods should perform their calculations in years.
Below we demonstrate how to define a custom function using existing SFH and metallicity distribution functions for demonstration purposes.
[16]:
class Normal(ZDist.Common):
"""
A normally distributed metallicity distribution.
Attributes:
mean (float)
The mean of the normal distribution.
sigma (float)
The standard deviation of the normal distribution.
"""
def __init__(self, mean, sigma):
"""
Initialise the metallicity distribution and parent.
Either metallicity or log10metallicity must be provided.
Args:
mean (float)
The mean of the normal distribution.
sigma (float)
The standard deviation of the normal distribution.
"""
# Instantiate the parent
ZDist.Common.__init__(
self,
name="Normal",
mean=mean,
sigma=sigma,
)
# Define this models parameters
self.mean = mean
self.sigma = sigma
def _weight(self, metal):
"""
Return the distribution at a metallicity.
Args:
metal (float)
The (linear) metallicity at which to evaluate the distribution.
Returns:
float
The weight of the metallicity distribution at metal.
"""
norm = 1 / (self.sigma * np.sqrt(2 * np.pi))
exponent = ((metal - self.mean) / self.sigma) ** 2
return norm * np.exp(-0.5 * exponent)
class Constant(SFH.Common):
"""
A constant star formation history.
The SFR is defined such that:
sfr = 1; t<=duration
sfr = 0; t>duration
Attributes:
duration (float)
The duration of the period of constant star formation.
"""
def __init__(self, duration):
"""
Initialise the parent and this parametrisation of the SFH.
Args:
duration (unyt_quantity)
The duration of the period of constant star formation.
"""
# Initialise the parent
SFH.Common.__init__(self, name="Constant", duration=duration)
# Set the model parameters
self.duration = duration.to("yr").value
def _sfr(self, age):
"""
Get the amount SFR weight in a single age bin.
Args:
age (float)
The age (in years) at which to evaluate the SFR.
"""
# Set the SFR based on the duration.
if age <= self.duration:
return 1.0
return 0.0
# Define a scaled delta function for metallicity
metal_dist = Normal(mean=0.05, sigma=0.001)
print(metal_dist)
# Define an exponential SFH
sfh = Constant(duration=500 * Myr)
print(sfh)
# Create the Stars object
custom_stars = Stars(
log10ages,
metallicities,
sf_hist=sfh,
metal_dist=metal_dist,
initial_mass=10**9 * Msun,
)
print(custom_stars)
fig, ax = custom_stars.plot_sfzh(show=True)
----------
SUMMARY OF PARAMETERISED METAL ENRICHMENT HISTORY
<class '__main__.Normal'>
mean: 0.05
sigma: 0.001
----------
----------
SUMMARY OF PARAMETERISED STAR FORMATION HISTORY
<class '__main__.Constant'>
duration: 500 Myr
median age: 250.00 Myr
mean age: 250.00 Myr
----------
+---------------------------------------------------------------------------+
| STARS |
+--------------------------+------------------------------------------------+
| Attribute | Value |
+--------------------------+------------------------------------------------+
| component_type | 'Stars' |
+--------------------------+------------------------------------------------+
| sf_hist_func | <__main__.Constant object at 0x7f5a6731a230> |
+--------------------------+------------------------------------------------+
| metal_dist_func | <__main__.Normal object at 0x7f5a67577190> |
+--------------------------+------------------------------------------------+
| metallicity_grid_type | 'log10Z' |
+--------------------------+------------------------------------------------+
| log10ages_lims | [6.0, 10.399999999999984, ] |
+--------------------------+------------------------------------------------+
| metallicities_lims | [1e-05, 0.01778279410038923, ] |
+--------------------------+------------------------------------------------+
| log10metallicities_lims | [-5.0, -1.75, ] |
+--------------------------+------------------------------------------------+
| ages (45,) | 1.00e+06 yr -> 2.51e+10 yr (Mean: 2.71e+09 yr) |
+--------------------------+------------------------------------------------+
| metallicities (14,) | 1.00e-05 -> 1.78e-02 (Mean: 2.90e-03) |
+--------------------------+------------------------------------------------+
| sf_hist (45,) | 0.00e+00 -> 1.85e+08 (Mean: 2.22e+07) |
+--------------------------+------------------------------------------------+
| metal_dist (14,) | 0.00e+00 -> 1.00e+09 (Mean: 7.14e+07) |
+--------------------------+------------------------------------------------+
| initial_mass | 999999999.9999999 Msun |
+--------------------------+------------------------------------------------+
| sfzh (45, 14) | 0.00e+00 -> 1.85e+08 (Mean: 1.59e+06) |
+--------------------------+------------------------------------------------+
| log10ages (45,) | 6.00e+00 -> 1.04e+01 (Mean: 8.20e+00) |
+--------------------------+------------------------------------------------+
| log10metallicities (14,) | -5.00e+00 -> -1.75e+00 (Mean: -3.38e+00) |
+--------------------------+------------------------------------------------+
