Warning, /eic-opticks/ana/source_debug.rst is written in an unsupported language. File is not indexed.
0001 Wavelength Distribution Debugging
0002 ====================================
0003
0004 Compares simulated photon wavelengths against blackbody expectation.
0005
0006
0007 * still a hint of "ringing steps" from 200:400nm, but seems acceptable
0008 (TODO: try increasing icdf bins from 1024 to identify)
0009
0010
0011 [ISSUE] wp last bin elevated
0012 -----------------------------
0013
0014 ::
0015
0016 In [69]: plt.close();plt.hist(wp, bins=200)
0017
0018 , 2215., 2158., 2046., 2017., 2052., 2111., 2565.]),
0019
0020
0021 [FIXED] Bug with w0 sel.recwavelength(0)
0022 -----------------------------------------
0023
0024 Without selection sel.recwavelength(0) from ggv-newton:
0025
0026 * length of 500000
0027
0028 * three bin spike at lower bound around 60nm, comprising about 7000 photons
0029 (not present in the uncompressed wp)
0030
0031 **FIXED WHEN AVOID WAVELENGTH DOMAIN DISCREPANCY BETWEEN SOURCES AND COMPRESSION**
0032
0033 * plateau from 60~190 nm
0034
0035 **MADE MUCH LESS OBJECTIONABLE BY INCREASING ICDF BINS FROM 256 TO 1024**
0036
0037 * normal service resumes above 190nm with good
0038 match to Planck black body curve
0039
0040 * 256 unique linspaced values, a result of the compression::
0041
0042 In [36]: np.allclose(np.linspace(60,820,256),np.unique(w)) # upper changed 810 to 820 by the fix
0043 Out[36]: True
0044
0045
