Warning, /eic-opticks/u4/tests/U4SimulateTest_ph.rst is written in an unsupported language. File is not indexed.
0001 U4SimulateTest_ph.rst
0002 =======================
0003
0004 Looking for records that cross between multiple PMTs by looking at ReplicaNumber at each step point::
0005
0006 In [33]: np.c_[np.arange(len(i_rp)),i_rp,n_rp,u_rp]
0007 Out[33]:
0008 array([[ 0, 88, 3, -1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0009 [ 1, 116, 2, -1, -1, -1, 1000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0010 [ 2, 57, 10, -1, -1, -1, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0011 [ 3, 80, 1, -1, -1, -1, 4000, 4000, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0012 [ 4, 49, 1, -1, -1, -1, 4000, 4000, -1, 3000, 3000, -1, 2000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0013 [ 5, 77, 5, -1, -1, -1, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0014 [ 6, 938, 1, -1, -1, -1, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0015 [ 7, 936, 1, -1, -1, -1, 4000, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0016 [ 8, 921, 1, -1, -1, -1, 4000, 4000, 4000, 4000, 4000, 4000, -1, -1, 2000, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0017 [ 9, 111, 9, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0018 [ 10, 0, 70, -1, -1, 1000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0019 [ 11, 2, 1, -1, -1, 1000, 1000, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0020 [ 12, 34, 3, -1, -1, 1000, 1000, -1, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0021 [ 13, 926, 60, -1, -1, 3000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0022 [ 14, 934, 4, -1, -1, 3000, 3000, -1, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0023 [ 15, 942, 1, -1, -1, 3000, 3000, -1, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0024 [ 16, 81, 560, -1, -1, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0025 [ 17, 87, 4, -1, -1, 4000, 4000, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0026 [ 18, 132, 5, -1, -1, 4000, 4000, -1, 1000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0027 [ 19, 221, 1, -1, -1, 4000, 4000, -1, 1000, 1000, -1, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0028 [ 20, 199, 1, -1, -1, 4000, 4000, -1, 1000, 1000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0029 [ 21, 144, 1, -1, -1, 4000, 4000, -1, 1000, 1000, 1000, 1000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0030 [ 22, 152, 1, -1, -1, 4000, 4000, -1, 1000, 1000, 1000, 1000, 1000, 1000, -1, 4000, 0, 0, 0, 0, 0, 0, 0, 0],
0031 [ 23, 574, 1, -1, -1, 4000, 4000, -1, 2000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0032 [ 24, 447, 1, -1, -1, 4000, 4000, -1, 2000, 2000, -1, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000],
0033 [ 25, 444, 1, -1, -1, 4000, 4000, -1, 2000, 2000, 2000, 2000, 2000, 2000, -1, 4000, 0, 0, 0, 0, 0, 0, 0, 0],
0034 [ 26, 686, 7, -1, -1, 4000, 4000, -1, 3000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0035 [ 27, 848, 1, -1, -1, 4000, 4000, -1, 3000, 3000, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0036 [ 28, 770, 1, -1, -1, 4000, 4000, -1, 3000, 3000, 3000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0037 [ 29, 707, 1, -1, -1, 4000, 4000, -1, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000, 0, 0],
0038 [ 30, 167, 124, -1, -1, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0039 [ 31, 91, 44, -1, -1, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0040 [ 32, 107, 32, -1, -1, 4000, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0041 [ 33, 345, 1, -1, -1, 4000, 4000, 4000, 4000, 4000, -1, 2000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0042
0043 ::
0044
0045 In [34]: np.where( v_rp == 24 )
0046 Out[34]: (array([447]),)
0047
0048 In [35]: np.where( v_rp == 39 )
0049 Out[35]: (array([189, 194, 285, 776, 808]),)
0050
0051 In [36]: rp[v_rp == 39]
0052 Out[36]:
0053 array([[ -1, -1, 4000, 4000, 4000, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0054 [ -1, -1, 4000, 4000, 4000, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0055 [ -1, -1, 4000, 4000, 4000, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0056 [ -1, -1, 4000, 4000, 4000, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
0057 [ -1, -1, 4000, 4000, 4000, 4000, 4000, 4000, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]], dtype=int32)
0058
0059 In [37]: rp[v_rp == 24]
0060 Out[37]: array([[ -1, -1, 4000, 4000, -1, 2000, 2000, -1, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000]], dtype=int32)
0061
0062 In [38]: np.where(v_rp == 24)
0063 Out[38]: (array([447]),)
0064
0065 In [39]: seqhis_(t.seq[447,0])
0066 Out[39]: ['TO BT BR BT BT BR BT BT BT SR SR SR SR SR SR SR', 'SR SR SR SR SA']
0067
0068
0069
