eic-opticks/ana/bidx.py

0001 #!/usr/bin/env python
0002 """
0003 Hmm this duplicates ana/boundary_counts.py
0004 """
0005 import numpy as np, os
0006 from opticks.ana.key import keydir
0007 from opticks.ana.blib import BLib
0008
0009 KEYDIR=keydir()
0010 blib = BLib()
0011
0012 if __name__ == '__main__':
0013     avi = np.load(os.path.join(KEYDIR, "GNodeLib/all_volume_identity.npy"))
0014
0015     bidx = ( avi[:,2] >>  0)  & 0xffff
0016     midx = ( avi[:,2] >> 16)  & 0xffff
0017
0018     b,n = np.unique( bidx, return_counts=True)
0019
0020     for i in range(len(b)):
0021         print("%3d : %3d : %7d : %s " % (b[i],b[i]+1, n[i],blib.bname(b[i])))
0022     pass
0023
0024
0025
0026