Post by Robin on Jul 25, 2017 16:00:17 GMT 8
Of course, projected 1D flux intensity plots can be very useful for quickly distinguishing the contributions from the bulge and disk component (albeit hiding much of the information in a 2D representation). profitEllipsePlot( ) is a good tool for deprojecting the image pixels into 1D surface brightness profiles for double component models; however, it does not allow for plotting of a single component model. If ever you had a reason to plot the 1D surface brightness profile of a single component fit, you can quickly get away with this by adding a pseudo second component with essentially zero flux:
This can then be fed into profitEllipsePlot( ) as follows:
Which yields something like the following:
The legend in the plot is, of course, misleading, but this is simply a diagnostic plot as it is probably more useful to investigate the 2D representation of the model residual.
add_pseudo_bulge = function(model) # Add a zero-point magnitude bulge
{
# <param: model [list]> - Single component modellist.
# <return: pseudo [list]> - A (pseudo) modellist containing a second null component.
pseudo = model # copy model list
for (key in names(model$sersic)){pseudo$sersic[[key]][2] = model$sersic[[key]][1]} # duplicate a second component
pseudo$sersic$mag[1] = ZERO_POINT # set magnitude of bulge (loc=1) to ZERO_POINT
return(pseudo)
}
This can then be fed into profitEllipsePlot( ) as follows:
profitEllipsePlot(Data=Data,modellist=add_pseudo_bulge(modellist))Which yields something like the following:
The legend in the plot is, of course, misleading, but this is simply a diagnostic plot as it is probably more useful to investigate the 2D representation of the model residual.
