[yt-users] HaloProfiler vs sum_mass_in_sphere.py

[Agarwal, Shankar]

Hi Matt,

Actually, that earlier post was a different question. There I was comparing HopFinder and HaloFinder. Stephen clarified that for me that these two calculate mass differently.

Here I am comparing HaolProfiler to itself. Let me know if my question is confusing.

shankar

-----Original Message-----
From: yt-users-bounces at lists.spacepope.org on behalf of Matthew Turk
Sent: Mon 11/30/2009 11:54 AM
To: Discussion of the yt analysis package
Subject: Re: [yt-users] HaloProfiler vs sum_mass_in_sphere.py
 
Hi Shankar,

You wrote to the list a week ago asking an almost identical question:

http://lists.spacepope.org/pipermail/yt-users-spacepope.org/2009-November/000252.html

and Stephen answered you:

http://lists.spacepope.org/pipermail/yt-users-spacepope.org/2009-November/000253.html

The same answer applies here.

-Matt

On Mon, Nov 30, 2009 at 9:48 AM, Agarwal, Shankar <sagarwal at ku.edu> wrote:
> Hi, I ran HaloProfiler on HaloAnalysis.out and here is the data for the first halo...
>
> # id    center[0]       center[1]       center[2]       RadiusMpc       TotalMassMsun
> 0000    0.169329105     0.144388856     0.8246570556    3.1986507886    1.07530961369e+15
>
>
>
>
> Then I ran sum_mass_in_sphere.py with the above center and radius ...
>
> from yt.mods import *
> fn = "RedshiftOutput0002"
> pf = load(fn)
> sp = pf.h.sphere([0.169329105,0.144388856,0.8246570556],3.1986507886/pf["mpc"])
> baryon_mass, particle_mass = sp.quantities["TotalQuantity"](
>        ["CellMassMsun", "ParticleMassMsun"], lazy_reader=True)
> print "Total mass in sphere is %0.5e (gas = %0.5e / particles = %0.5e)" % \
>            (baryon_mass + particle_mass, baryon_mass, particle_mass)
>
>
> And got...
>
> Total mass in sphere is 1.27226e+15 (gas = 2.06890e+14 / particles = 1.06537e+15)
>
>
> Let me know if the TotalMassMsun should match. Because they don't.
>
> regards
> shankar
> KU Cosmology
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

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