https://pfspipe.ipmu.jp/jira This file is an XML representation of an issue en-us 8.3.4 803005 13-09-2019 https://pfspipe.ipmu.jp/jira/browse/SIM2D-146 DRP 2-D Simulator <p>The dichroic throughputs from <tt>pfs&#95;thr&#95;20201231&#95;ext&#95;all&#95;&#42;.dat</tt> have a sum that exceeds unity, which is unphysical.</p> <p><span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/14534/14534_dichroic.png" style="border: 0px solid black" /></span></p> <div class="code panel" style="border-width: 1px;"><div class="codeContent panelContent"> <pre class="code-java"> <span class="code-keyword">import</span> numpy as np <span class="code-keyword">import</span> matplotlib.pyplot as plt blueData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20201231_ext_all_blu.dat"</span>) redData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20201231_ext_all_red.dat"</span>) wavelength = np.concatenate((blueData[:2600, 0], redData[:, 0])) <span class="code-keyword">assert</span>(np.all(blueData[2600:, 0] == redData[:401, 0])) blue = np.zeros_like(wavelength) blue[:3001] = blueData[:, 5] red = np.zeros_like(wavelength) red[2600:] = redData[:, 5] total = blue + red plt.plot(wavelength, blue, <span class="code-quote">"b-"</span>) plt.plot(wavelength, red, <span class="code-quote">"r-"</span>) plt.plot(wavelength, total, <span class="code-quote">"k-"</span>) plt.xlabel(<span class="code-quote">"Wavelength (nm)"</span>) plt.ylabel(<span class="code-quote">"Dichroic transmission"</span>) plt.suptitle(<span class="code-quote">"pfs_thr_20201231_ext_all_*.dat"</span>) plt.show() </pre> </div></div> SIM2D-146

Dichroic model is inconsistent between blue and red

Story Normal Done Done Kiyoto Yabe price Thu, 16 Dec 2021 19:41:19 +0000 Wed, 22 Dec 2021 00:17:49 +0000 Wed, 22 Dec 2021 00:17:49 +0000 0 3 <p>I just upload the new version (pfs_thr_20211220_ext_all_*.dat) to `drp_instdata`. The cause is that the reflectance of the blue dichroic used for the estimation is taken from the datasheet of the "witness sample". I could not find the proper estimation for the reflectance, so assume that reflectance = 1 - transmittance, but note that the absorption documented is 3-4% in average and might be severer in bluer part. Anyway, the updated version does not show any tip at the transition wavelength when blue and red are combined. </p> <p>It looks like that with the new curves, the dichroic cutoff in the blue arm has moved 15 nm to shorter wavelengths, and the ripples on the blue arm have disappeared.<br/> Comparing the ripples in the old curves with quartz spectra from PFI, I estimated the shift was about 10 nm.</p> <p><span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/15061/15061_sim2d-146.png" style="border: 0px solid black" /></span></p> <p>I confirmed that the NIR is unaffected.</p> <p><span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/15063/15063_sim2d-146-nir.png" style="border: 0px solid black" /></span></p> <div class="code panel" style="border-width: 1px;"><div class="codeContent panelContent"> <pre class="code-java"> <span class="code-keyword">import</span> numpy as np <span class="code-keyword">import</span> matplotlib.pyplot as plt blueData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20201231_ext_all_blu.dat"</span>) redData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20201231_ext_all_red.dat"</span>) nirData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20201231_ext_all_nir.dat"</span>) <span class="code-keyword">assert</span>(np.all(blueData[2600:, 0] == redData[:401, 0])) <span class="code-keyword">assert</span>(np.all(redData[3150:, 0] == nirData[:401, 0])) wavelength = np.concatenate((blueData[:2600, 0], redData[:, 0], nirData[401:, 0])) blue = np.zeros_like(wavelength) blue[:3001] = blueData[:, 9] red = np.zeros_like(wavelength) red[2600:6151] = redData[:, 9] nir = np.zeros_like(wavelength) nir[5750:] = nirData[:, 9] total = blue + red + nir plt.plot(wavelength, blue, <span class="code-quote">"b:"</span>) plt.plot(wavelength, red, <span class="code-quote">"g:"</span>) plt.plot(wavelength, nir, <span class="code-quote">"r:"</span>) plt.plot(wavelength, total, <span class="code-quote">"k:"</span>, label=<span class="code-quote">"Old"</span>) blueData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20211220_ext_all_blu.dat"</span>) redData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20211220_ext_all_red.dat"</span>) nirData = np.genfromtxt(<span class="code-quote">"../drp_instdata/data/throughput/pfs_thr_20211220_ext_all_nir.dat"</span>) <span class="code-keyword">assert</span>(np.all(blueData[2600:, 0] == redData[:401, 0])) <span class="code-keyword">assert</span>(np.all(redData[3150:, 0] == nirData[:401, 0])) wavelength = np.concatenate((blueData[:2600, 0], redData[:, 0], nirData[401:, 0])) blue = np.zeros_like(wavelength) blue[:3001] = blueData[:, 9] red = np.zeros_like(wavelength) red[2600:6151] = redData[:, 9] nir = np.zeros_like(wavelength) nir[5750:] = nirData[:, 9] total = blue + red + nir plt.plot(wavelength, blue, <span class="code-quote">"b-"</span>) plt.plot(wavelength, red, <span class="code-quote">"g-"</span>) plt.plot(wavelength, nir, <span class="code-quote">"r-"</span>) plt.plot(wavelength, total, <span class="code-quote">"k-"</span>, label=<span class="code-quote">"New"</span>) plt.xlabel(<span class="code-quote">"Wavelength (nm)"</span>) plt.ylabel(<span class="code-quote">"Transmission"</span>) plt.legend() plt.suptitle(<span class="code-quote">"Old vs New"</span>) plt.show() </pre> </div></div> <p>I've updated the simulator to use the new throughput curves provided by <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=kiyoto.yabe" class="user-hover" rel="kiyoto.yabe">Kiyoto Yabe</a>.</p> <p>No objections with suggested changes in pull requests.</p> <p>Thanks, <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=kiyoto.yabe" class="user-hover" rel="kiyoto.yabe">Kiyoto Yabe</a> and <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=hassan" class="user-hover" rel="hassan">hassan</a>.</p> <p>Merged to master.</p> Development Rank 0|02qpt6:508yl0chzzy0000q01r3609 Reviewers hassan