[yt-users] How ot overplot star particles on density projections

[Latif]

Hi Matt and Nathan,
Thank you for your help.  I think Matt's idea is very good and bit easy.
Unfortunately, I didn't get success with it yet.  I am getting empty
panels.  Matt, Is it close to what you suggested?  It is most likely that I
am messing up some thing due to my poor understanding.  Here is my script.
Do you guys know what is going wrong here.

http://paste.yt-project.org/show/3708/

Cheers
Latif


On Fri, Jul 19, 2013 at 3:53 PM, Matthew Turk <matthewturk at gmail.com> wrote:

> Hi Latif,
>
> On Fri, Jul 19, 2013 at 9:48 AM, Latif <latifne at gmail.com> wrote:
> > Hi all,
> > Is there any way to annotate particles to the multi-plot
> porjections/slices?
> > I am using the following recipe from the webpage.
> >
> > http://yt-project.org/doc/cookbook/multi_plot_slice_and_proj.py
> >
> > thanks in advance,
> >
> > Cheers
> > Latif
>
> Unfortunately this is a lot harder, as the FRBs used there do not
> expose the same annotate_* methods that the sliceplot, projectionplot,
> etc do.
>
> However, you *may* be able to do something similar, although I have
> not tested it, by creating a ProjectionPlot or SlicePlot, then
> swapping out the .axes and .figure objects that resides on the plot
> object itself.  Then you can call ._setup_plots() on the
> ProjectionPlot or SlicePlot, and it should re-create all the necessary
> info.
>
> So it would look something like this, once you have the axes objects
> you're interested in from the recipe you linked to:
>
> p = ProjectionPlot( ... )
> p.plots["Density"].figure = fig
> p.plots["Density"].axes = dens_axes[0]
> p.annotate_whatever()
> p._setup_plots()
>
> Then you can call:
>
> fig.savefig("%s_3x2" % pf)
>
> This is all very rough, but I think it should get you there.
>
> -Matt
>
> >
> >
> > On Thu, Jul 4, 2013 at 9:43 PM, Latif <latifne at gmail.com> wrote:
> >>
> >> HI Matt,
> >> Thank you for a prompt and precise response.
> >> Cheers
> >> Latif
> >>
> >>
> >>
> >> On Thu, Jul 4, 2013 at 7:04 PM, Matthew Turk <matthewturk at gmail.com>
> >> wrote:
> >>>
> >>> Hi Latif,
> >>>
> >>> The callback you're looking for is "particles" and you can specify a
> >>> "ptype" argument for specifying the type of particles.  (If you're
> >>> using Enzo, this should be the number in the "particle_type" field you
> >>> are selecting.)  Note also that annotate_particles accepts a width
> >>> argument, in code units, which says how wide the selected region will
> >>> be around the center of the slice or the center of the box for
> >>> projections.  So if you are using a projection and you want the whole
> >>> box, you can do 1.0/pf['unitary'] to get the full domain.
> >>>
> >>> Here is an example:
> >>>
> >>> s = SlicePlot(pf, "x", "Density")
> >>> s.annotate_particles(1.0/pf['kpc'], p_size = 1.0, ptype = 1)
> >>>
> >>> which will choose particle_type = =1.
> >>>
> >>> -Matt
> >>>
> >>> On Thu, Jul 4, 2013 at 10:20 AM, Latif <latifne at gmail.com> wrote:
> >>> > Hi all,
> >>> > I want to overplot star particles on density projections/slices. It
> is
> >>> > probably a very simple question but could not figure out how to do
> it.
> >>> > Can
> >>> > I also get information about their position and velocities as well?
> >>> > thanks in advance,
> >>> > Cheers
> >>> > Latif
> >>> >
> >>> > _______________________________________________
> >>> > 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
> >>
> >>
> >
> >
> > _______________________________________________
> > 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 --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20130720/87a1646e/attachment.htm>