From 8fda8839ac79bdb1f54435e209a449bfd3016e20 Mon Sep 17 00:00:00 2001
From: Johannes Schumann <johannes.schumann@fau.de>
Date: Wed, 11 May 2022 17:35:25 +0200
Subject: [PATCH] Fix calculation
---
km3buu/output.py | 9 +++++++--
1 file changed, 7 insertions(+), 2 deletions(-)
diff --git a/km3buu/output.py b/km3buu/output.py
index acdde18..21395df 100644
--- a/km3buu/output.py
+++ b/km3buu/output.py
@@ -498,14 +498,19 @@ class GiBUUOutput:
def bjorken_y(roottuple_data):
"""
Calculate Bjorken y scaling variable for the GiBUU events
- (Lab. frame)
+
+ Definition: y = pq/pk
Parameters
----------
roottuple_data: awkward.highlevel.Array
"""
d = roottuple_data
- y = 1 - np.divide(np.array(d.lepOut_E), np.array(d.lepIn_E))
+ q = GiBUUOutput._q(d)
+ pq = q[0, :] * d.nuc_E - q[1, :] * d.nuc_Px - q[2, :] * d.nuc_Py - q[
+ 3, :] * d.nuc_Pz
+ pk = d.lepIn_E * d.nuc_E - d.lepIn_Px * d.nuc_Px - d.lepIn_Py * d.nuc_Py - d.lepIn_Pz * d.nuc_Pz
+ y = pq / pk
return y
@property
--
GitLab