[yt-dev] Issue #717: tracking rays on top of cut region (yt_analysis/yt)

[yong zheng]

New issue 717: tracking rays on top of cut region
https://bitbucket.org/yt_analysis/yt/issue/717/tracking-rays-on-top-of-cut-region

yong zheng:

I am trying to use the .hierarchy.ray( ) to extract data along a particular
 light ray, but when I check the temperature field of that light ray data, I
find that the data points it extracted were not really what I want, the temperature range goes beyond my limit, here is part of my code:

>     sp = pf.h.sphere(center, radius)
>     warmhot = sp.cut_region(['grid["Temperature"] >= 1e5',
> 'grid["Temperature"] < 1e6'])
>     p0 = np.array([center[0], center[1]-250./pf['kpc'], center[2]])
>     p1 = np.array([center[0], center[1]+250./pf['kpc'], center[2]])
>     ray = warmhot.hierarchy.ray(p0, p1)
>     tempt = warmhot['Temperature']
>     ray_tempt = ray['Temperature']
>
when I examine the temperature range for warmnot, it is
      (100000.1171875, 999999.5625)
but for ray_tempt, it is:
    (6607.2802734375, 86239664.0)
It seems that the .hierarchy.ray() extracts data directly from the pf.h dataset instead of the cut warmhot region. Could anyone help?