[yt-users] halo mass function question

[Robert Thompson]

I'm running through different cookbook tutorials and I've ran into a snag
when trying to form the halo mass function. (
http://yt.enzotools.org/doc/howto/hmf_howto.html#hmf-howto)

I can run the halo finder, but the halo profiler is where I run into the
issue.  I've been unsuccessful with both N-body only and N-body+gas runs of
different sizes.  It seems that the filter the tutorial employs filters out
all of my halos and I get this for every halo:

yt.lagos   ERROR      2010-08-04 12:34:35,198 Skipping halo with r_max /
r_min = 0.153458.
yt.lagos   ERROR      2010-08-04 12:34:35,198 Skipping halo with r_max /
r_min = 0.167703.
yt.lagos   ERROR      2010-08-04 12:34:35,199 Skipping halo with r_max /
r_min = 0.209952.
yt.lagos   ERROR      2010-08-04 12:34:35,199 Skipping halo with r_max /
r_min = 0.184501.
yt.lagos   ERROR      2010-08-04 12:34:35,199 No halos in filtered list.

Now if I change any of the parameters such as must_be_virialized=False, or
reduce the overdensity_field, or reduce the total mass filter to a lower
mass I get the same result.  By simply removing the filter:

hp.add_halo_filter(HP.VirialFilter,must_be_virialized=True,
                   overdensity_field='ActualOverdensity',
                   virial_overdensity=200,
                   virial_filters=[['TotalMassMsun','>=','1e8']],
                   virial_quantities=['TotalMassMsun','RadiusMpc'])

 it goes through and profiles all of the halos, but doesn't output the
TotalMassMsun that I'm looking for.  I've tried simply removing the filter
and adding:

hp.add_profile('TotalMassMsun', weight_field=None, accumulation=False)

just before hp.make_profile() and it just returns yt.lagos ERROR no halos in
filtered list once again.  Am I missing something simple?  Thanks~

-Robert Thompson
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20100804/b4a137bb/attachment.htm>