CSGOptiX/tests/CSGOptiXSimTest.py

0001 #!/usr/bin/env python
0002
0003 import os, numpy as np
0004 from opticks.ana.fold import Fold
0005 from opticks.ana.p import *
0006
0007 if __name__ == '__main__':
0008     t = Fold.Load()
0009     PIDX = int(os.environ.get("PIDX","-1"))
0010
0011     p = t.photon if hasattr(t, "photon") else None
0012     r = t.record if hasattr(t, "record") else None
0013     seq = t.seq if hasattr(t, "seq") else None
0014     nib = seqnib_(seq[:,0])  if not seq is None else None
0015
0016
0017     for i in range(len(p)):
0018         if not (PIDX == -1 or PIDX == i): continue
0019         if PIDX > -1: print("PIDX %d " % PIDX)
0020         print("r[%d,:,:3]" % i)
0021         print(r[i,:nib[i],:3])
0022         print("\n\nbflagdesc_(r[%d,j])" % i)
0023         for j in range(nib[i]):
0024             print(bflagdesc_(r[i,j]))
0025         pass
0026
0027         #print("ridiff_(r[%d])*1000." % i)
0028         #print(ridiff_(r[i])*1000.)
0029
0030         print("\n")
0031         print("p[%d]" % i)
0032         print(p[i])
0033         print("\n")
0034         print("bflagdesc_(p[%d])" % i)
0035         print(bflagdesc_(p[i]))
0036         print("\n")
0037         if not seq is None:
0038             print("seqhis_(seq[%d,0]) nib[%d]  " % (i,i) )
0039             print(" %s : %s "% ( seqhis_(seq[i,0]), nib[i] ))
0040             print("\n")
0041         pass
0042         print("\n\n")
0043     pass
0044
0045
0046
0047
0048