[PIPE2D-703] Provide F-star standards covering the full PFS spectral range Created: 21/Jan/21 Updated: 21/Jul/21 Resolved: 21/Jul/21 |
|
| Status: | Done |
| Project: | DRP 2-D Pipeline |
| Component/s: | None |
| Affects Version/s: | None |
| Fix Version/s: | None |
| Type: | Story | Priority: | Normal |
| Reporter: | hassan | Assignee: | Takuji Yamashita |
| Resolution: | Done | Votes: | 0 |
| Labels: | None | ||
| Remaining Estimate: | Not Specified | ||
| Time Spent: | Not Specified | ||
| Original Estimate: | Not Specified | ||
| Attachments: |
|
| Story Points: | 2 |
| Sprint: | 2DDRP-2021 A 2, 2DDRP-2021 A 7 |
| Reviewers: | hassan |
| Description |
|
Following the 2D DRP tech telecon 2020-01-20, it was pointed out that the simulated data for F-star standards used by DRP do not cover the full wavelength range (they cut off at around 1.23 micron while the NIR limit is 1.26 micron). Please provide spectra for these standards over the full PFS range. For specific details of these F-star standards, contact Paul Price. The GA team may also need to be contacted. |
| Comments |
| Comment by price [ 21/Jan/21 ] |
|
The spectra currently used for flux standards in the simulator are from AMBRE, provided by Masayuki Tanaka. |
| Comment by Masayuki Tanaka [ 21/Jan/21 ] |
|
Takuji Yamashita, could you provide the extrapolated version of these templates?
|
| Comment by Takuji Yamashita [ 21/Jan/21 ] |
|
OK! Should I make a branch for this on /drp_instdata/data/objects/fluxCal/ and upload the extrapolated templates there? |
| Comment by price [ 22/Jan/21 ] |
|
That would be great, thanks. Could you please upload only the spectra corresponding to the ones that already exist (we're trying to keep the size of the package from becoming overwhelmingly large) and update the README.txt file appropriately? |
| Comment by Takuji Yamashita [ 25/Jan/21 ] |
|
I have placed the extrapolated spectra and an updated README.txt in a new branch. These spectra were generated in the ticket of PIPE2D-364. The data format was changed from the original version. Only flux density in nJy is stored, you need to use FITST headers for wavelength solution, and the pixel scale is 0.01 nm/pix now. If you have any questions, please feel free to contact me. |
| Comment by price [ 30/Jan/21 ] |
|
I got a warning from git: Encountered 6 file(s) that should have been pointers, but weren't: data/objects/fluxCal/p6500_g+4.0_m0.0_t01_z+0.00_a+0.00.AMBRE_Extp.fits data/objects/fluxCal/p6500_g+4.0_m0.0_t01_z-1.00_a+0.00.AMBRE_Extp.fits data/objects/fluxCal/p7000_g+4.0_m0.0_t01_z+0.00_a+0.00.AMBRE_Extp.fits data/objects/fluxCal/p7000_g+4.0_m0.0_t01_z-1.00_a+0.00.AMBRE_Extp.fits data/objects/fluxCal/p7500_g+4.0_m0.0_t01_z+0.00_a+0.00.AMBRE_Extp.fits data/objects/fluxCal/p7500_g+4.0_m0.0_t01_z-1.00_a+0.00.AMBRE_Extp.fits Could you please check that you're using git lfs? |
| Comment by price [ 30/Jan/21 ] |
|
Adding files with a different format will require coordinated changes in drp_instmodel. |
| Comment by price [ 30/Jan/21 ] |
|
The spectrum appears to be discontinuous at the red end: the continuum doesn't match. >>> import astropy.io.fits >>> ff = astropy.io.fits.open("data/objects/fluxCal/p6500_g+4.0_m0.0_t01_z+0.00_a+0.00.AMBRE_Extp.fits") >>> flux = ff[1].data.Flux >>> import numpy as np >>> wavelength = np.arange(len(flux))*0.01 + 300 >>> import matplotlib.pyplot as plt >>> plt.plot(wavelength, flux, 'k-') >>> plt.show()
|
| Comment by Takuji Yamashita [ 02/Feb/21 ] |
|
I did not use git lfs. I do again. According to DAMD-50, the discussion on the data format is not fixed yet. We need to decide that, Masayuki Tanaka and hassan?
I made those extrapolation for the spectra which were convolved with LSF (Gaussian), while the stored spectra at 1200nm are ones before the LSF convolution. Please use them after the LSF convolution on your side. I can send the LSF-convolved spectra but I thought the unconvolved data is better because you can apply your LSF. |
| Comment by Takuji Yamashita [ 02/Feb/21 ] |
|
The spectrum that was convolved by an LSF (a Gaussian with FWHM = 2.4 \AA) is overplotted. A small "jump" remains at 1200 nm. This jump is smaller than 1% which is our target accuracy of flux calibration. I believe this is not a serious problem for flux calibration, but it is nice if we can resolve this. So, please let me add this issue to my development item list after my first loop of development. A function fitting for extrapolation was done on a spectrum with small absorption lines and therefore the extrapolation function is an "average" including small absorption lines and is a bit lower than the continuum. This is why we find the jump there.
|
| Comment by Takuji Yamashita [ 10/Feb/21 ] |
|
The FITS files were uploaded to github again using git lfs. Please check it. |
| Comment by price [ 11/Feb/21 ] |
|
The FITS files are still present in git, from the previous commit. This makes the git repo larger than it should be (this is the whole reason for git lfs). I'll fix it, as it requires editing history. |
| Comment by price [ 11/Feb/21 ] |
|
I've fixed the commits on the branch. If you have any problems on your local branch, please git fetch && git reset --hard origin/tickets/ |
| Comment by price [ 11/Feb/21 ] |
|
I've updated the simulator code that reads these files, and also the pipeline code that uses them. We'll also want to regenerate the weekly raw dataset before this merges. hassan, are you aware of any other coming changes that will require regenerating the weekly, so we can coordinate them? |
| Comment by hassan [ 11/Feb/21 ] |
|
To my understanding price, there are no other upcoming changes that require regenerating the weekly. |
| Comment by rhl [ 11/Feb/21 ] |
|
When are we getting Sean's new object lists? And when are we going to use different target lists for brn and bmn, that's a structural change to the weekly. That'll require |
| Comment by hassan [ 11/Feb/21 ] |
|
I've agreed with Paul offline about 30 mins ago that we will regenerate the weekly after |
| Comment by price [ 17/Jul/21 ] |
|
I believe this is ready to merge, and we should include it in the new generation of sims. |
| Comment by hassan [ 20/Jul/21 ] |
|
Paul has reviewed the changes. |
| Comment by price [ 21/Jul/21 ] |
|
Merged. |