Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2026-04-09 07:49:55

0001 import numpy as np
0002 
0003 a = np.load("/tmp/fakeuser/opticks/GEOM/fakegeom/simg4ox/ALL0_no_opticks_event_name/A000/record.npy")
0004 b = np.load("/tmp/fakeuser/opticks/GEOM/fakegeom/simg4ox/ALL0_no_opticks_event_name/B000/f000/record.npy")
0005 
0006 print(a.shape)
0007 print(b.shape)
0008 
0009 assert a.shape == b.shape
0010 
0011 diff = [i for i, (a, b) in enumerate(zip(a[:, 1:], b[:, 0:-1])) if not np.allclose(a, b, rtol=0, atol=1e-5)]
0012 print(diff)
0013 
0014 assert diff == [14, 22, 32, 34, 40, 81, 85]