Add a conversion dict from GiBUU to GENIE

km3buu/output.py

@@ -79,7 +79,7 @@ FLUX_INFORMATION_DTYPE = np.dtype([("energy", np.float64),
                                    ("flux", np.float64),
                                    ("events", np.float64)])
 
-EVENT_TYPE = {
+SCATTERING_TYPE = {
     1: "QE",
     2: "P33(1232)",
     3: "P11(1440)",
@@ -119,6 +119,46 @@ EVENT_TYPE = {
     37: "2pi background",
 }
 
+SCATTERING_TYPE_TO_GENIE = {
+    1: 1,  # QE        -> kScQuasiElastic
+    2: 4,  # P33(1232) -> kScResonant  
+    3: 4,  # P11(1440) -> kScResonant
+    4: 4,  # S11(1535) -> kScResonant
+    5: 4,  # S11(1650) -> kScResonant
+    6: 4,  # S11(2090) -> kScResonant
+    7: 4,  # D13(1520) -> kScResonant
+    8: 4,  # D13(1700) -> kScResonant
+    9: 4,  # D13(2080) -> kScResonant
+    10: 4,  # D15(1675) -> kScResonant
+    11: 4,  # G17(2190) -> kScResonant
+    12: 4,  # P11(1710) -> kScResonant
+    13: 4,  # P11(2100) -> kScResonant
+    14: 4,  # P13(1720) -> kScResonant
+    15: 4,  # P13(1900) -> kScResonant
+    16: 4,  # F15(1680) -> kScResonant
+    17: 4,  # F15(2000) -> kScResonant
+    18: 4,  # F17(1990) -> kScResonant
+    19: 4,  # S31(1620) -> kScResonant
+    20: 4,  # S31(1900) -> kScResonant
+    21: 4,  # D33(1700) -> kScResonant
+    22: 4,  # D33(1940) -> kScResonant
+    23: 4,  # D35(1930) -> kScResonant
+    24: 4,  # D35(2350) -> kScResonant
+    25: 4,  # P31(1750) -> kScResonant
+    26: 4,  # P31(1910) -> kScResonant
+    27: 4,  # P33(1600) -> kScResonant
+    28: 4,  # P33(1920) -> kScResonant
+    29: 4,  # F35(1750) -> kScResonant
+    30: 4,  # F35(1905) -> kScResonant
+    31: 4,  # F37(1950) -> kScResonant
+    32: 0,  # pi neutron-background -> kScNull
+    33: 0,  # pi proton-background -> kScNull
+    34: 3,  # DIS -> kScDeepInelastic
+    35: 0,  # 2p2h QE -> kScNull
+    36: 0,  # 2p2h Delta -> kScNull
+    37: 0,  # 2pi background -> kScNull
+}
+
 ROOTTUPLE_KEY = "RootTuple"
 
 EMPTY_KM3NET_HEADER_DICT = {