sysrap/tests/sblackbody_test.py

0001 #!/usr/bin/env python
0002
0003 import numpy as np
0004 from opticks.ana.fold import Fold
0005 import matplotlib.pyplot as mp
0006 SIZE=np.array([1280, 720])
0007
0008 if __name__ == '__main__':
0009     f = Fold.Load(symbol="f")
0010     psr = f.psr
0011     psr_names = f.psr_names
0012
0013     print(psr)
0014     print(psr_names)
0015
0016     title = "sblackbody_test.sh : %s " % f.base
0017     fig, ax = mp.subplots(figsize=SIZE/100.)
0018     fig.suptitle(title)
0019     for i in range(len(psr)):
0020         nm = psr[i,:,0]
0021         bb = psr[i,:,1]
0022         bb /= bb.sum()
0023         ax.plot( nm, bb, label=psr_names[i] )
0024     pass
0025     ax.legend()
0026     fig.show()
0027
0028