eic-opticks/ana/lvn.py

0001 #!/usr/bin/env python
0002
0003 import numpy as np
0004 from opticks.ana.key import keydir
0005
0006 if __name__ == '__main__':
0007     kd = keydir()
0008     path = os.path.join(kd, "GNodeLib/all_volume_LVNames.txt")
0009
0010     lvn = np.loadtxt(path, dtype=np.object)
0011     print(" lvn.shape %s  path %s " % (str(lvn.shape), path) )
0012
0013     u_lvn, u_lvn_counts = np.unique(lvn, return_counts=True)
0014
0015     for i in range(len(u_lvn)):
0016         print(" %5d : %6d : %50s " % (i, u_lvn_counts[i], u_lvn[i] ))
0017     pass
0018
0019
0020
0021
0022
0023