Dear YT people,
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.
yt.ProjectionPlot(ds, 'z', 'density, center= 'm')
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')