<div dir="ltr"><div><div><div><div>Dear YT people, <br></div><div><br></div><div>I need to do integration of density field along line of sight. I know the ProjectionPlot function plots the integration of a field along line of sight.  <br></div><div><br></div><div>yt.ProjectionPlot(ds, 'z', 'density, center= 'm')<br></div><br>However, I was wondering if it is possible to get the projection as an array (or new field) , so that I can make a ProfilePlot  of the integrated field along line of sight and another field ('Mass')<br><br></div>In other words, I want to do the following:<br></div><div><br></div><div>sp = ds.sphere("max", (200, "kpc"))<br>proj = yt.ProfilePlot(sp, 'Integrated_Density_along_Line_Of_Sight', 'Mass)<br><br></div>Thanks<br></div>Tazkera<br><div><div><br></div></div></div>