sysrap/tests/SEvtTest.py

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