[Yt-dev] EnzoFOF merger tree

[John Wise]

Hi all,

I just extended Matt's merger tree that works on enzo's inline FOF halo finder to work on multiple outputs, and build a complete merger tree.  It produces DOT (Graphviz) files.  I tried to make it simple to use and put in the necessary docstrings.  Here's an example to make a merger tree from z = 15->7 only including halos with >500 particles.

"""
from yt.mods import *
import yt.analysis_modules.halo_merger_tree.api as hm

tree = hm.EnzoFOFMergerTree((7.0,15.0), load_saved=True)
#tree.build_tree(0, max_children=4)
tree.build_tree(0, min_particles=500)
tree.write_dot()
"""

The load_saved argument will load the pickled results from previous calculations.  It will save it by default if it's false.

Here's the result from the most massive halo in a 256^3 sim.  Careful, it's large!  1.9MB, 4325x7461 pix
http://www.astro.princeton.edu/~jwise/pics/merger-tree.png

If you try it out, please let me know of any problems!

Cheers,
John