Actually yes. You could check this yourself by examining the last column in HopAnalysis out. If max_r < 2*(1/512)=3.9e-3 then yes, it will skip it. In one of my simulations, a 5e11 Msun halo has a max_r of 2.64e-03 corresponding to 168 kpc/h. The halo profiler was not built for profiling halos with < 2 cells in radius. In my mind you have at least four options:<div>
<br></div><div>1. Run with AMR to resolve the dark matter potential so that you can use the halo profiler as built. Moderate computational cost. I'm guessing this was the case for Stephen's simulation where a decent fraction of the halos in the HopAnalysis were able to be profiled.</div>
<div>2. Run with a gigantic unigrid such that you can resolve these halos with > 2 (and preferably many) cells. This is a very computationally expensive route as you'd need more than a 2048^3 unigrid in a 200 Mpc box. I do not recommend this option.<br>
3. Write a dark matter particle profiler that uses the particle positions (and not grid cells) to bin up the profile. Little computational cost but requires coding on your part.</div><div>4. Run the same size unigrid on a smaller physical volume to capture a different range of the mass function.</div>
<div><br></div><div>Cheers,</div><div>Sam</div><div><br></div><div><br></div><div><br><div class="gmail_quote">On Sun, Feb 28, 2010 at 12:07 PM, Agarwal, Shankar <span dir="ltr"><<a href="mailto:sagarwal@ku.edu" target="_blank">sagarwal@ku.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Sam,<br>
<br>
Here is how I am using it right now...<br>
<div><br>
hp.add_halo_filter(HP.VirialFilter,must_be_virialized=True,<br>
overdensity_field='ActualOverdensity',<br>
virial_overdensity=200,<br>
</div> virial_filters=[['TotalMassMsun','>=','5e11']],<br>
virial_quantities=['TotalMassMsun','RadiusMpc'])<br>
hp.make_profiles(filename="VirialHaloes.out")<br>
<br>
And i get 500.<br>
<br>
<br>
Do you think it has something to do with the grid spacing (for 512 ^3 grid in 200 Box, it is 390 Kpc)<br>
<br>
HaloProfiler may be skipping small halos completely since they are smaller than 390 Kpc.<br>
<font color="#888888"><br>
<br>
Shankar<br>
</font><div><div></div><div><br>
<br>
-----Original Message-----<br>
From: <a href="mailto:yt-users-bounces@lists.spacepope.org" target="_blank">yt-users-bounces@lists.spacepope.org</a> on behalf of Sam Skillman<br>
Sent: Sun 2/28/2010 12:34 PM<br>
To: Discussion of the yt analysis package<br>
Subject: Re: [yt-users] Dealing with HOP groups<br>
<br>
Hi Shankar,<br>
<br>
Do you have a line that looks like:<br>
hp.add_halo_filter(HP.VirialFilter,must_be_virialized=True,<br>
overdensity_field='ActualOverdensity',<br>
virial_overdensity=200,<br>
virial_filters=[['TotalMassMsun','>=','1e14']],<br>
virial_quantities=['TotalMassMsun','RadiusMpc'])<br>
<br>
If so, please post it. The number of halos > 10^14 Msun is fairly close to<br>
500 in a 200 Mpc box. If you want the halo profiler to use all halos, set<br>
virial_filters=[['TotalMassMsun','>=','1e14']] accordingly.<br>
<br>
Sam<br>
<br>
On Sun, Feb 28, 2010 at 11:08 AM, Agarwal, Shankar <<a href="mailto:sagarwal@ku.edu" target="_blank">sagarwal@ku.edu</a>> wrote:<br>
<br>
> Hi Stephen,<br>
><br>
> I read<br>
> <a href="http://lists.spacepope.org/pipermail/yt-users-spacepope.org/2009-March/000110.html" target="_blank">http://lists.spacepope.org/pipermail/yt-users-spacepope.org/2009-March/000110.html</a><br>
><br>
> Do you remember the cosmology parameters for your run: particles, grid,<br>
> box, refinement, AMR ON or OFF ?<br>
><br>
><br>
> I have run a 200 Mpc/h box, 512^3 particle, 512^3 grid, AMR-OFF simulation<br>
> and I get about 30,000 halos using HOP. When I try to finding virial halos,<br>
> I get only 500. I may have to increase the grid to 1024^3.<br>
><br>
> Let me know how you got 40,000 virial halos.<br>
><br>
> Thanks<br>
> Shankar<br>
> _______________________________________________<br>
> yt-users mailing list<br>
> <a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><br>
<br>
<br>
<br>
--<br>
Samuel W. Skillman<br>
DOE Computational Science Graduate Fellow<br>
Center for Astrophysics and Space Astronomy<br>
University of Colorado at Boulder<br>
samuel.skillman[at]<a href="http://colorado.edu" target="_blank">colorado.edu</a><br>
<br>
</div></div><br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Samuel W. Skillman<br>DOE Computational Science Graduate Fellow<br>Center for Astrophysics and Space Astronomy<br>University of Colorado at Boulder<br>samuel.skillman[at]<a href="http://colorado.edu" target="_blank">colorado.edu</a><br>
</div>