Range and Resolution Expansion

[ICenturyI]

Hello,

Let me start by showing appreciation for this tool since it has been pretty helpful with my college work, even for someone who still isn't neither fond of or good at programming... yet... hopefully.

Thus, I've been following your 'example.ipynb' layout and has been pretty good, and does almost everything I require it to do. However, since my project is related with the upcoming MOONS' observations, I would need my program to not only fit in the NIR range as well as having MOONS' corresponding resolutions - it has 3 resolutions for HR spectra and 3 for LR spectra.

I've tried some straightforward/brutish methods like just adding the Paschen and other lines to 'line_priors' , adding the remaining range in the 'conti_windows' and tying each line after 'line_priors' with its correspondent MOONS' resolution. However, these didn't really produced the results I wanted.

So I wanted to know, firstly if it's possible to extend the fitting range and apply different resolutions to it, and if so, what is the best way to go about it.

I don't really know if this 'Issues' tab is supposed to be used to clarify these type of questions, but since I'm using the example you've provided, as a layout, I thought it was the best way to know relatively quickly the solution to my problem.

Thank you for your help in advance,

JM

[WenkeRen]

Hi JM,

I'm sorry for the wait. I appreciate your efforts and feedback. Also, I'm very excited to hear that our pipeline is tested on latest specstropic surveys.

Firstly, I believe our code should be able to handle NIR data, your problem is mostly a consequence of lack of coverage of PCA templates. Since the lack of NIR survey data, we haven't developed any work around for AGN-host decomposition covering rest-frame NIR. A quick test would be setting the decompose_host=False in the function q_mle.Fit(...). If I understand correct, by saying "these didn't really produced the results I wanted" you mean the code didn't throw out any errors but can only extend the fitting to about rest-frame 7450A. You can try this advice and return me your outcome if there are any errors. The could should at least work (well, maybe a bad fitting) as my expectation.

By following my advice, you are skipping the AGN-host decomposition section, meaning the code will deal the whole spectra as a pure AGN and starting the fitting directly from continuum (i.e., power-law, polynomial component, FeII psedo continuum et al.). You must be very cautious to this results especially when you are dealing with rest-frame NIR where a strong galaxy contribution is expected and the long-range fitting with pure AGN continuum probably will be very weird. Nevertheless, the measurement of emission lines, including the width, flux, wavelength center, but be careful about the properties related to the continuum, like EW.

Besides, one minor suggestion to your fitting is that you definitely should adjust the conti_windows parameter, extending it to NIR according to your situation. This probably will improve the continuum fitting but not guaranteed.

Let me know if you have further issues and I'm always gald to help.

Best,
Wenke