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