sysrap/tests/gaussQTables.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 cdat = os.path.expandvars("$CDAT")
0008 a = np.loadtxt(cdat, delimiter="f,      //"  )
0009
0010 tab0 = a[:250]
0011 tab1 = a[250:]
0012
0013
0014 fig, ax = mp.subplots(figsize=SIZE/100.)
0015
0016 ax.axvline( 2e-6 )
0017 ax.axvline( 5e-4 )
0018
0019
0020 ax.scatter( a[:,1], a[:,0], s=0.1 )
0021 #ax.scatter( tab0[:,1], tab0[:,0], s=0.1 )
0022 #ax.scatter( tab1[:,1], tab1[:,0], s=0.1 )
0023
0024
0025
0026 fig.show()
0027
0028
0029
0030
0031
0032
0033