From 8a2be11e9fe7a8ff750f833e3d6b6294353066ad Mon Sep 17 00:00:00 2001
From: Johannes Schumann <johannes.schumann@fau.de>
Date: Tue, 20 Apr 2021 22:09:59 +0200
Subject: [PATCH] Neutrino energy parameter handling
---
km3buu/jobcard.py | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/km3buu/jobcard.py b/km3buu/jobcard.py
index 705fc2a..9b76a1a 100644
--- a/km3buu/jobcard.py
+++ b/km3buu/jobcard.py
@@ -111,8 +111,8 @@ def generate_neutrino_jobcard(events,
Interaction channel ["CC", "NC", "antiCC", "antiNC"]
flavour: str
Flavour ["electron", "muon", "tau"]
- energy: float
- Initial energy of the neutrino in GeV
+ energy: tuple
+ Initial energy range of the neutrino in GeV
target: (int, int)
(Z, A) describing the target nucleon
write_pert: boolean (default: True)
@@ -129,6 +129,7 @@ def generate_neutrino_jobcard(events,
The input path pointing to the GiBUU lookup data which should be used
"""
jc = read_jobcard(join(dirname(abspath(__file__)), "data/template.job"))
+ jc["input"]["path_to_input"] = input_path
# NEUTRINO
jc["neutrino_induced"]["process_ID"] = PROCESS_LOOKUP[process.lower()]
@@ -145,6 +146,9 @@ def generate_neutrino_jobcard(events,
runs = events // run_events
jc["input"]["numEnsembles"] = run_events
jc["input"]["num_runs_SameEnergy"] = runs
+ # ENERGY
+ jc["nl_neutrino_energyflux"]["eflux_min"] = energy[0]
+ jc["nl_neutrino_energyflux"]["eflux_max"] = energy[1]
# DECAY
if do_decay:
for i in DECAYED_HADRONS:
--
GitLab