Ability to fit foreground stars and AGN at the same time?

[rremigio]

Hello,

I have been using BADASS to fit IFU spectra of nearby Seyfert 1s. One problem I have encountered in my work is that there are some galaxies that have foreground stars that are nearby (e.g. Mrk 110) or in front of the host galaxy (e.g. see archival HST ACS images of RXJ 2044.0+2833) such that there are spaxels where the observed spectrum essentially has two different components at different redshifts (say, AGN emission and the foreground star). This is not as much of a problem if one component dominates over the other (e.g. very weak/non-existent AGN continuum and strong foreground star):

I do run into issues if the spectra has detectable features from both objects, like here (H\beta absorption and Mg I 5167, 5173, 5182 A triplet from the foreground star; BLR emission from H\beta, some weak Mg I absorption redshifted to ~5360 A):

My question is that would it be possible to add/implement a fitting option or some way for BADASS to deal with a foreground object?