From 34c0408c16fda725f07f6da026c57f918b398714 Mon Sep 17 00:00:00 2001
From: Johannes Schumann <johannes.schumann@fau.de>
Date: Sat, 1 Apr 2023 09:09:27 +0200
Subject: [PATCH] Check if initial vertex position is in CAN in order to decide
for status of hadr. secondaries
---
km3buu/output.py | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/km3buu/output.py b/km3buu/output.py
index 1090986..df74a19 100644
--- a/km3buu/output.py
+++ b/km3buu/output.py
@@ -747,8 +747,12 @@ def write_detector_file(gibuu_output,
if prop_particles is not None:
lep_out_trk.status = km3io.definitions.trkmembers[
"TRK_ST_PROPLEPTON"]
- generator_particle_state = km3io.definitions.trkmembers[
- "TRK_ST_UNDEFINED"]
+ if geometry.in_can(vtx_pos):
+ generator_particle_state = km3io.definitions.trkmembers[
+ "TRK_ST_FINALSTATE"]
+ else:
+ generator_particle_state = km3io.definitions.trkmembers[
+ "TRK_ST_UNDEFINED"]
else:
lep_out_trk.status = km3io.definitions.trkmembers[
"TRK_ST_FINALSTATE"]
--
GitLab