diff --git a/km3buu/output.py b/km3buu/output.py index cbd51c2c9bdd61ca6ffd529495999bf74a40eda6..bc978f892702c2d4c647ee28839588636e007a3d 100644 --- a/km3buu/output.py +++ b/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 = {