Skip to content
Snippets Groups Projects
Commit 8a2be11e authored by Johannes Schumann's avatar Johannes Schumann
Browse files

Neutrino energy parameter handling

parent e6cec255
No related branches found
No related tags found
1 merge request!7Flux handling
......@@ -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:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment