sysrap/tests/RINDEX.py

0001 #!/usr/bin/env python
0002
0003 import os, numpy as np
0004 import matplotlib.pyplot as mp
0005 SIZE=np.array([1280, 720])
0006
0007 if __name__ == '__main__':
0008     path = os.path.expandvars("$HOME/.opticks/GEOM/$GEOM/CSGFoundry/SSim/stree/material/$material/RINDEX.npy")
0009     a = np.load(path)
0010     print(path)
0011     print(a)
0012
0013     ri = a[:,1]
0014     eV = a[:,0]*1e6
0015     hc_eVnm = 1239.84198
0016     nm = hc_eVnm/eV    ## 330,331,...,599,600  very close to integer nm with 1239.84198
0017
0018
0019     fig, ax = mp.subplots(figsize=SIZE/100.)
0020     fig.suptitle(path)
0021
0022     ax.plot( nm, ri )
0023     ax.scatter( nm, ri )
0024
0025     ax.axvline(200, dashes=[2,1])
0026     ax.axvline(120, dashes=[1,1])
0027     ax.axvline(80,  dashes=[3,1])
0028
0029     fig.show()
0030
0031
0032 pass
0033