u4/tests/U4SimulateTest_af.py

0001 #!/usr/bin/env python
0002 """
0003 U4SimulateTest_af.py
0004 ========================
0005
0006 """
0007
0008
0009 import os, numpy as np
0010 from opticks.ana.fold import Fold
0011 from opticks.ana.p import *
0012
0013 MODE = int(os.environ.get("MODE", 0))
0014 if MODE > 0:
0015     from opticks.ana.pvplt import *
0016 pass
0017
0018 PID = int(os.environ.get("PID", -1))
0019 if PID == -1: PID = int(os.environ.get("OPTICKS_G4STATE_RERUN", -1))
0020
0021
0022 if __name__ == '__main__':
0023
0024     print("PID : %d " % (PID))
0025     a = Fold.Load("$BASE/ALL0", symbol="a")
0026     print(repr(a))
0027
0028     b = Fold.Load("$BASE/ALL1", symbol="b")
0029     print(repr(b))
0030
0031     ## modern way to do this is with sevt.py sevt.Load
0032
0033     print("seqhis_(a.seq[PID,0] : %s " % seqhis_(a.seq[PID,0] ))
0034     print("seqhis_(b.seq[PID,0] : %s " % seqhis_(b.seq[PID,0] ))
0035
0036     ar = a.record[PID]
0037     br = b.record[PID]
0038
0039     print("\nar = a.record[PID] ; ar[:4] \n", ar[:4])
0040     print("\nbr = b.record[PID] ; br[:4] \n", br[:4] )
0041
0042