I'm running through different cookbook tutorials and I've ran into a snag when trying to form the halo mass function. (<a href="http://yt.enzotools.org/doc/howto/hmf_howto.html#hmf-howto">http://yt.enzotools.org/doc/howto/hmf_howto.html#hmf-howto</a>)<div>
<br></div><div>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:</div>
<div><br></div><div><div>yt.lagos ERROR 2010-08-04 12:34:35,198 Skipping halo with r_max / r_min = 0.153458.</div><div>yt.lagos ERROR 2010-08-04 12:34:35,198 Skipping halo with r_max / r_min = 0.167703.</div>
<div>yt.lagos ERROR 2010-08-04 12:34:35,199 Skipping halo with r_max / r_min = 0.209952.</div><div>yt.lagos ERROR 2010-08-04 12:34:35,199 Skipping halo with r_max / r_min = 0.184501.</div><div>yt.lagos ERROR 2010-08-04 12:34:35,199 No halos in filtered list.</div>
</div><div><br></div><div>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:</div>
<div><br></div><div><div>hp.add_halo_filter(HP.VirialFilter,must_be_virialized=True,</div><div> overdensity_field='ActualOverdensity',</div><div> virial_overdensity=200,</div><div>
virial_filters=[['TotalMassMsun','>=','1e8']],</div><div> virial_quantities=['TotalMassMsun','RadiusMpc'])</div></div><div><br></div><div> 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:</div>
<div><br></div><div>hp.add_profile('TotalMassMsun', weight_field=None, accumulation=False) </div><div><br></div><div>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~</div>
<div><br></div><div>-Robert Thompson</div>