From 314b6bf7756996fcb099a2a81c97f492ed18cac4 Mon Sep 17 00:00:00 2001
From: Johannes Schumann <johannes.schumann@fau.de>
Date: Wed, 10 Nov 2021 10:41:22 +0100
Subject: [PATCH] Clean up target density calculation
---
km3buu/output.py | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/km3buu/output.py b/km3buu/output.py
index da3a873..b6bebe8 100644
--- a/km3buu/output.py
+++ b/km3buu/output.py
@@ -615,12 +615,12 @@ def write_detector_file(gibuu_output,
tau_secondaries = propagate_lepton(event_data, np.sign(nu_type) * 15)
media = read_default_media_compositions()
- density = media["SeaWater"]["density"]
+ density = media["SeaWater"]["density"] # [g/cm^3]
element = mendeleev.element(gibuu_output.Z)
target = media["SeaWater"]["elements"][element.symbol]
- target_density = 1e3 * density * target[1]
- targets_per_volume = target_density * (1e3 * constants.Avogadro /
- target[0].atomic_weight)
+ target_density = 1e3 * density * target[1] # [kg/m^3]
+ targets_per_volume = target_density / target[
+ 0].atomic_weight / constants.atomic_mass
w2 = gibuu_output.w2weights(geometry.volume, targets_per_volume, 4 * np.pi)
global_generation_weight = gibuu_output.global_generation_weight(4 * np.pi)
--
GitLab