[yt-users] r_min in HaloProfiler

[Agarwal, Shankar]

Sam,

I am trying to get the virial masses of the halos (with mass upto 1e+13 Msun) with this filter...

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


But when I look at the radial_profiles/Halo_xxxx_profile.dat files, I see that none of the halos are virialized. The peak halo densities are reaching about 20*mean_matter_density. I began my 200Mpc/h Box with 512^3 particles simulation at z=99. I have not looked at the literature in great detail but is there something fishy here ?


shankar



-----Original Message-----
From: yt-users-bounces at lists.spacepope.org on behalf of Sam Skillman
Sent: Tue 12/8/2009 3:23 PM
To: Discussion of the yt analysis package
Subject: Re: [yt-users] r_min in HaloProfiler
 
Shankar,

r_min is defined in your email:
r_min = 2 * self.pf.h.get_smallest_dx() * self.pf['mpc']

it's twice the smallest dx in units of mpc, dx being the smallest cell size.
 it is being rejected because your halo has something like 3 cells (radial)
in it, which is almost certainly not resolved.

sam

On Tue, Dec 8, 2009 at 2:11 PM, Agarwal, Shankar <sagarwal at ku.edu> wrote:

> Hi,
>
> While running HaloProfiler, I am seeing this...
>
> P001 yt.lagos   ERROR      2009-12-08 15:53:22,403 Skipping halo with r_max
> / r_min = 1.885682.
>
>
>
> I looked at yt/extensions/HaloProfiler.py...
>
>           r_min = 2 * self.pf.h.get_smallest_dx() * self.pf['mpc']
>            if (halo['r_max'] / r_min < PROFILE_RADIUS_THRESHOLD):
>                mylog.error("Skipping halo with r_max / r_min = %f." %
> (halo['r_max']/r_min))
>
>
> r_max would be the distance to the furthest particle in the halo as found
> by HaloFinder.
> But what is r_min ? And what is the basis for skipping this halo ?
>
>
> shankar
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>



-- 
Samuel W. Skillman
DOE Computational Science Graduate Fellow
Center for Astrophysics and Space Astronomy
University of Colorado at Boulder
samuel.skillman[at]colorado.edu

-------------- next part --------------
A non-text attachment was scrubbed...
Name: winmail.dat
Type: application/ms-tnef
Size: 3949 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20091208/9fb0c016/attachment-0001.bin>