Skip to content
Snippets Groups Projects

Draft: Resolve "Muon Propagation"

Closed Draft: Resolve "Muon Propagation"
8-muon-propagation into master
Closed Johannes Schumann requested to merge 8-muon-propagation into master
Compare
and
5 files
+ 248
− 109
Compare changes
  • Side-by-side
  • Inline
Files
5
km3buu/tests/test_output.py
+ 22
− 25
@@ -30,6 +30,7 @@ try:
libpath = environ.get("KM3NET_LIB")
if libpath is None:
libpath = Config().km3net_lib_path
print(libpath)
KM3NET_LIB_AVAILABLE = (ROOT.gSystem.Load(join(libpath,
"libKM3NeTROOT.so")) >= 0)
except ModuleNotFoundError:
@@ -112,21 +113,19 @@ class TestOfflineFile(unittest.TestCase):
11.90433897, 2.1818, 1.45689677, 0.49284856, 8.33975778, 0.28362369
])
np.testing.assert_array_almost_equal(evt.mc_tracks.dir_x, [
0.18255849, -0.2469, 0.48623089, 0.23767571, 0.24971059, 0.11284916
0.745621, 0.871178, 0.829481, 0.674661, 0.692116, -0.328162
])
np.testing.assert_array_almost_equal(evt.mc_tracks.dir_y, [
-0.80816248, -0.619212, -0.49241334, -0.84679953, -0.83055629,
-0.82624071
0.028316, -0.344457, 0.397546, 0.070903, 0.086895, -0.084494
])
np.testing.assert_array_almost_equal(evt.mc_tracks.dir_z, [
0.55995162, 0.745398, 0.72187854, 0.47585798, 0.4978161,
-0.55189796
-0.665769, -0.349854, -0.392325, -0.734715, -0.716537, -0.940835
])
# Test dataset is elec CC -> outgoing particles are placed at vertex pos
np.testing.assert_allclose(evt.mc_tracks.t, 8603022.62272017)
np.testing.assert_allclose(evt.mc_tracks.pos_x, -127.07940486)
np.testing.assert_allclose(evt.mc_tracks.pos_y, -122.54421157)
np.testing.assert_allclose(evt.mc_tracks.pos_z, 208.57726764)
np.testing.assert_allclose(evt.mc_tracks.t, 21670533.7327)
np.testing.assert_allclose(evt.mc_tracks.pos_x, 137.221704)
np.testing.assert_allclose(evt.mc_tracks.pos_y, 111.069631)
np.testing.assert_allclose(evt.mc_tracks.pos_z, 51.350148)
usr = evt.mc_tracks.usr[0]
# XSEC
np.testing.assert_almost_equal(evt.w2list[13], 40.62418521597373)
@@ -169,21 +168,19 @@ class TestMultiFileOutput(unittest.TestCase):
11.90433897, 2.1818, 1.45689677, 0.49284856, 8.33975778, 0.28362369
])
np.testing.assert_array_almost_equal(evt.mc_tracks.dir_x, [
0.18255849, -0.2469, 0.48623089, 0.23767571, 0.24971059, 0.11284916
0.745621, 0.871178, 0.829481, 0.674661, 0.692116, -0.328162
])
np.testing.assert_array_almost_equal(evt.mc_tracks.dir_y, [
-0.80816248, -0.619212, -0.49241334, -0.84679953, -0.83055629,
-0.82624071
0.028316, -0.344457, 0.397546, 0.070903, 0.086895, -0.084494
])
np.testing.assert_array_almost_equal(evt.mc_tracks.dir_z, [
0.55995162, 0.745398, 0.72187854, 0.47585798, 0.4978161,
-0.55189796
-0.665769, -0.349854, -0.392325, -0.734715, -0.716537, -0.940835
])
# Test dataset is elec CC -> outgoing particles are placed at vertex pos
np.testing.assert_allclose(evt.mc_tracks.t, 8603022.62272017)
np.testing.assert_allclose(evt.mc_tracks.pos_x, -127.07940486)
np.testing.assert_allclose(evt.mc_tracks.pos_y, -122.54421157)
np.testing.assert_allclose(evt.mc_tracks.pos_z, 208.57726764)
np.testing.assert_allclose(evt.mc_tracks.t, 21670533.7327)
np.testing.assert_allclose(evt.mc_tracks.pos_x, 137.221704)
np.testing.assert_allclose(evt.mc_tracks.pos_y, 111.069631)
np.testing.assert_allclose(evt.mc_tracks.pos_z, 51.350148)
usr = evt.mc_tracks.usr[0]
# XSEC
np.testing.assert_almost_equal(evt.w2list[13], 40.62418521597373)
@@ -205,16 +202,16 @@ class TestMultiFileOutput(unittest.TestCase):
np.testing.assert_array_almost_equal(
evt.mc_tracks.E, [7.043544, 3.274632, 4.429621, 0.21289])
np.testing.assert_array_almost_equal(
evt.mc_tracks.dir_x, [0.997604, 0.824817, 0.941969, 0.00302])
evt.mc_tracks.dir_x, [0.635764, 0.883134, 0.342816, -0.190106])
np.testing.assert_array_almost_equal(
evt.mc_tracks.dir_y, [-0.058292, -0.553647, 0.327013, -0.097914])
evt.mc_tracks.dir_y, [0.6503 , 0.181522, 0.896241, -0.334025])
np.testing.assert_array_almost_equal(
evt.mc_tracks.dir_z, [0.037271, 0.114683, -0.075871, 0.99519])
evt.mc_tracks.dir_z, [0.415829, 0.432579, 0.281477, 0.923194])
# Test dataset is elec CC -> outgoing particles are placed at vertex pos
np.testing.assert_allclose(evt.mc_tracks.t, 1951721.26185)
np.testing.assert_allclose(evt.mc_tracks.pos_x, -171.8025)
np.testing.assert_allclose(evt.mc_tracks.pos_y, -55.656482)
np.testing.assert_allclose(evt.mc_tracks.pos_z, 363.950535)
np.testing.assert_allclose(evt.mc_tracks.t, 14969858.655096)
np.testing.assert_allclose(evt.mc_tracks.pos_x, -223.792944)
np.testing.assert_allclose(evt.mc_tracks.pos_y, 240.693195)
np.testing.assert_allclose(evt.mc_tracks.pos_z, 244.288455)
usr = evt.mc_tracks.usr[0]
# XSEC
np.testing.assert_almost_equal(evt.w2list[13], 4.218262109165907)