PFS-JIRA

PFS-JIRA https://pfspipe.ipmu.jp/jira This file is an XML representation of an issue en-us

8.3.4

803005

13-09-2019

[SIM2D-150] Provide model spectra of flux standards for 2D simulator https://pfspipe.ipmu.jp/jira/browse/SIM2D-150 DRP 2-D Simulator <p>I will provide the set of model spectra of flux standard stars for the 2D simulator. These spectra will be read in the simulator. The spectra are in the format of pfsSimObject. </p> SIM2D-150

Provide model spectra of flux standards for 2D simulator

Task Normal Done Done Takuji Yamashita Takuji Yamashita flux-calibration model-templates Mon, 15 Aug 2022 08:57:59 +0000 Fri, 16 Jun 2023 23:56:12 +0000 Fri, 16 Jun 2023 23:56:12 +0000 0 3 <p>Update:<br/> I will provide 200 spectra. 100 of them are uniformly distributed across stellar parameters, Teff = 6300-7300K, log(g/(cm/s^2)) = 3.5-4.5, <span class="error">[Fe/H]</span> = -1 - 0.4. The remaining 100 spectra are also uniformly distributed, but the range is wider: Teff = 6000-7500K, log(g/(cm/s^2)) = 1-5, <span class="error">[Fe/H]</span> = -5 - 1. In both samples, <span class="error">[alpha/Fe]</span> depends on <span class="error">[Fe/H]</span>.</p> <p>Then, the spectra are doppler shifted with a radial velocity randomly selected from -200 to 200 km/s. Galactic extinction is also applied to the spectra with `F99ExtinctionCurve` at a fixed ra/dec=0deg/0deg.</p> <p>200 model spectra of F stars were made in the pfsSimObject format. The relevant catalog contains catId, objId, tract/patch, magnitude, stellar parameters, etc. I use catId = 9 but I can change it if you have any recommendations. type is set to be FLUXSTD. I did not calculate J band magnitudes.<br/> <br/> model spectra and catalog:<br/> <a href="https://hscdata.mtk.nao.ac.jp/hsc_bin_dist/temp/deleted_on_20221101/model_fstars_2022_08_17.tar.gz" class="external-link" rel="nofollow">https://hscdata.mtk.nao.ac.jp/hsc_bin_dist/temp/deleted_on_20221101/model_fstars_2022_08_17.tar.gz</a><br/> <br/> The objId = 0 – 99 are stars of the sample with wider parameter ranges, which includes stars scarcely observed (the latter in the above description). On the other hand, the objId = 100-199 are of the sample composed of stars with likely stellar parameters (the former). For the integration test and weekly, it would be better to select stars of objId = 100-199.<br/> <br/> Paul, could you have a look at these spectra and catalog? <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=price" class="user-hover" rel="price">price</a><br/> </p> <p>The catalog file includes columns <tt>m_g</tt>, <tt>m_r</tt>, <tt>m_i</tt>, <tt>m_z</tt> and <tt>m_Y</tt>. What photometric system is this? I would have guessed HSC, except you use <tt>Y</tt> instead of <tt>y</tt>. We need to be able to map these to the bandpasses in <a href="https://github.com/Subaru-PFS/obs_pfs/tree/master/pfs/fluxCal/bandpass" class="external-link" rel="nofollow">obs_pfs</a>. <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=sogo.mineo" class="user-hover" rel="sogo.mineo">sogo.mineo</a>, is there a lookup-table somewhere that matches the filter name in the <tt>pfsConfig</tt> to the bandpass file?</p> <p>I think that we need to add noise to these photometric measurements. Do you have a prescription that you'd like to use, or should I invent something?</p> <p>Everything has r = 18 mag. Perhaps it would be better to have a variety of brightnesses, so that you can investigate performance as a function of brightness? (I could artificially dim spectra in the simulator if you'd like; but I need from you a prescription for selecting the magnitudes.)</p> <p>I've verified that the simulator can read the spectra, and that the flux scaling is reasonable.</p> <p>Is it OK to randomly select stars from your catalog for the simulated exposures? We currently have three modes for the simulator:</p> <ul> <li>The "integration test" has fifteen fibers, three of which are flux standards. The purpose of this dataset is to allow a relatively quick check that changes have not broken the pipeline. Due to the small number of flux standards, I would probably lean towards using brighter sources so that operations downstream from flux calibration are not affected by lower signal-to-noise. There are two designs (the second is just a shuffled version of the first), so we could potentially use up to 6 different spectra.</li> <li>The "weekly" has a single spectrograph with full fiber sampling, with 60 flux standards. The purpose of this dataset is to check that the pipeline is performing well on data that looks something like we'll get from the real instrument. There are two designs (the second is just a shuffled version of the first), so we could potentially use up to 120 different spectra.</li> <li>The "science" set has a single spectrograph with full fiber sampling, with 59 flux standards. The purpose of this dataset is to produce science-like outputs for the 1D team. There are two designs (one stars, one extragalactic) so we could use a total of 118 different spectra.</li> </ul> <blockquote><p>is there a lookup-table somewhere that matches the filter name in the pfsConfig to the bandpass file?</p></blockquote> <p>Here is it: <a href="https://github.com/Subaru-PFS/drp_stella/blob/13d4aef8fd33b22f49bd0d959b6334e68f762733/python/pfs/drp/stella/fitReference.py#L111" class="external-link" rel="nofollow">https://github.com/Subaru-PFS/drp_stella/blob/13d4aef8fd33b22f49bd0d959b6334e68f762733/python/pfs/drp/stella/fitReference.py#L111</a></p> <p>I am not sure whether <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=takuji.yamashita" class="user-hover" rel="takuji.yamashita">Takuji Yamashita</a> uses this mapping or not. In particular, this mapping maps <tt>"i"</tt> to <tt>"HSC-I"</tt> filter as opposed to <tt>"HSC-I2"</tt>. If he were to use it, <tt>m_i</tt> would not be <tt>HSC-I2</tt> but <tt>HSC-I</tt>.</p> <p>I am sorry for this late reply... <br/> The photometric system is HSC. I used g, r2, i2, z, and y of HSC in FilterCurve that Mineo-san mentioned above. </p> <p>For magnitude and error, I am sorry, I misunderstood. I meant 18 mag is just a default value. I was thinking of changing magnitudes by using --fluxStdMag option or by editing the catalog on my side for each run. I was also thinking of adding magnitude errors by editing pfsConfig. I will send an updated catalog which includes magnitude errors and a variety of magnitudes.</p> <p>For the selection of spectra, I will select 3 spectra for integration test, 60 for weekly, and 59 for science.</p> <p>This is the updated catalog. <br/> <a href="https://hscdata.mtk.nao.ac.jp/pfs_tempfiles/deleted_on_20221231/model_fstars_2022_10_06.tar.gz" class="external-link" rel="nofollow">https://hscdata.mtk.nao.ac.jp/pfs_tempfiles/deleted_on_20221231/model_fstars_2022_10_06.tar.gz</a></p> <p>I added columns of magnitude errors (e_m_x), fluxes, and flux errors. I use uniform flux errors for all stars, which corresponds to S/N = 100 at 20 mag (it is slightly worse than one of HSC-SSP). The magnitudes are distributed from g = 14 to 20 mag. This magnitude range is close to one of our flux standard F stars.</p> <p>I recommend objId = 103, 170, and 182 for the integration test. Because the parameters in the catalog do not depend on objId, we can use 60 spectra with objId = 100 - 159 for weekly. They are the main sample. Similarly, for the science mode, we can use 59 spectra with objId = 141 - 199.</p> <p>As part of this, we will need to check that we're using the correct filter names, following <a href="https://pfspipe.ipmu.jp/jira/browse/PIPE2D-1110" title="Make match filterName between pfsConfig and flux calibration" class="issue-link" data-issue-key="PIPE2D-1110"><del>PIPE2D-1110</del></a>.</p> <p>The filters I used in the above model spectra are g/r2/i2/z/y of HSC. But I have updated some minor points. Could you use the following new spectra and catalog?<br/> <br/> new: <br/> <a href="https://hscdata.mtk.nao.ac.jp/pfs_tempfiles/deleted_on_20230331/model_fstars_2022_12_22.tar.gz" class="external-link" rel="nofollow">https://hscdata.mtk.nao.ac.jp/pfs_tempfiles/deleted_on_20230331/model_fstars_2022_12_22.tar.gz</a><br/> </p> <ul class="alternate" type="square"> <li>Filters were changed from HSC to PS1. The filter curves, from g_ps1 to y_ps1, in fitReference.py were used.</li> <li>A flux error was changed to 350nJy, which corresponds to a sensitivity of PS1 g-band in 3pi survey.</li> </ul> <p>These spectra have been incorporated into a new integration test and weekly data set (<a href="https://pfspipe.ipmu.jp/jira/browse/SIM2D-151" title="Incorporate simulated flux standard spectra" class="issue-link" data-issue-key="SIM2D-151"><del>SIM2D-151</del></a>).</p> Blocks SIM2D-151 PIPE2D-992 Development Epic Link PIPE2D-303 Rank 0|02qpjd:00r20060i2004