From 9653032f0380d9a5980985e65f606240ee5f66d6 Mon Sep 17 00:00:00 2001
From: Johannes Schumann <johannes.schumann@fau.de>
Date: Thu, 29 Apr 2021 17:25:22 +0200
Subject: [PATCH] Update w2list w.r.t. xsection
---
km3buu/output.py | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/km3buu/output.py b/km3buu/output.py
index db2e3e2..32d644e 100644
--- a/km3buu/output.py
+++ b/km3buu/output.py
@@ -126,7 +126,7 @@ PARTICLE_MC_STATUS = {
W2LIST_LOOKUP = {
"PS": 0,
"EG": 1,
- "XSEC": 2,
+ "XSEC_MEAN": 2,
"COLUMN_DEPTH": 3,
"P_EARTH": 4,
"WATER_INT_LEN": 5,
@@ -137,6 +137,7 @@ W2LIST_LOOKUP = {
"CC": 10,
"DISTAMAX": 11,
"WATERXSEC": 12,
+ "XSEC": 13,
"TARGETA": 15,
"TARGETZ": 16,
"VERINCAN": 17,
@@ -507,6 +508,8 @@ def write_detector_file(gibuu_output,
head.set_line(k, v)
head.Write("Head")
+ mean_xsec_func = gibuu_output.mean_xsec
+
for mc_event_id, event in enumerate(event_data):
evt.clear()
evt.id = mc_event_id
@@ -517,6 +520,7 @@ def write_detector_file(gibuu_output,
evt.w.push_back(-1.0) #w3 (= w2*flux)
# Event Information (w2list)
evt.w2list.resize(W2LIST_LENGTH)
+ evt.w2list[W2LIST_LOOKUP["XSEC_MEAN"]] = mean_xsec_func(event.lepIn_E)
evt.w2list[W2LIST_LOOKUP["XSEC"]] = event.xsec
evt.w2list[W2LIST_LOOKUP["TARGETA"]] = gibuu_output.A
evt.w2list[W2LIST_LOOKUP["TARGETZ"]] = gibuu_output.Z
--
GitLab