Add gibuu weight to an additional field in w2list

km3buu/output.py

@@ -201,6 +201,7 @@ W2LIST_LOOKUP = {
     "TARGETZ": 16,
     "VERINCAN": 17,
     "LEPINCAN": 18,
+    "GIBUU_WEIGHT": 19
 }
 
 W2LIST_LENGTH = len(W2LIST_LOOKUP)
@@ -646,6 +647,7 @@ def write_detector_file(gibuu_output,
         evt.w2list[W2LIST_LOOKUP["ICHAN"]] = event.evType
         evt.w2list[W2LIST_LOOKUP["VERINCAN"]] = 1
         evt.w2list[W2LIST_LOOKUP["LEPINCAN"]] = 1
+        evt.w2list[W2LIST_LOOKUP["GIBUU_WEIGHT"]] = event.weight
 
         # Vertex Position
         vtx_pos = np.array(geometry.random_pos())

km3buu/tests/test_output.py

@@ -130,3 +130,5 @@ class TestAANET(unittest.TestCase):
         np.testing.assert_equal(evt.w2list[9], 34)
         # CC/NC
         np.testing.assert_equal(evt.w2list[10], 2)
+        # GiBUU weight
+        np.testing.assert_almost_equal(evt.w2list[19], 0.8167222969153614)