-
Type:
Story
-
Status: Done (View Workflow)
-
Priority:
Normal
-
Resolution: Done
-
Affects Version/s: None
-
Fix Version/s: None
-
Component/s: None
-
Labels:None
The dichroic throughputs from pfs_thr_20201231_ext_all_*.dat have a sum that exceeds unity, which is unphysical.
import numpy as np import matplotlib.pyplot as plt blueData = np.genfromtxt("../drp_instdata/data/throughput/pfs_thr_20201231_ext_all_blu.dat") redData = np.genfromtxt("../drp_instdata/data/throughput/pfs_thr_20201231_ext_all_red.dat") wavelength = np.concatenate((blueData[:2600, 0], redData[:, 0])) assert(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, "b-") plt.plot(wavelength, red, "r-") plt.plot(wavelength, total, "k-") plt.xlabel("Wavelength (nm)") plt.ylabel("Dichroic transmission") plt.suptitle("pfs_thr_20201231_ext_all_*.dat") plt.show()