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

Fixes PROPOSAL interface

parent 83e57c4f
No related branches found
No related tags found
1 merge request!20Proposal 7
Pipeline #23572 failed
...@@ -18,6 +18,7 @@ import proposal as pp ...@@ -18,6 +18,7 @@ import proposal as pp
from particle import Particle from particle import Particle
import awkward as ak import awkward as ak
from collections import defaultdict from collections import defaultdict
import pathlib
from .config import Config from .config import Config
...@@ -40,8 +41,9 @@ PROPOSAL_TARGET = pp.medium.AntaresWater() ...@@ -40,8 +41,9 @@ PROPOSAL_TARGET = pp.medium.AntaresWater()
def _setup_propagator(max_distance, particle_definition): def _setup_propagator(max_distance, particle_definition):
pp.InterpolationSettings.tables_path = Config().proposal_itp_tables itp_table_path = Config().proposal_itp_tables
pathlib.Path(itp_table_path).mkdir(parents=True, exist_ok=True)
pp.InterpolationSettings.tables_path = itp_table_path
crosssection = pp.crosssection.make_std_crosssection( crosssection = pp.crosssection.make_std_crosssection(
particle_def=particle_definition, particle_def=particle_definition,
target=PROPOSAL_TARGET, target=PROPOSAL_TARGET,
...@@ -50,12 +52,12 @@ def _setup_propagator(max_distance, particle_definition): ...@@ -50,12 +52,12 @@ def _setup_propagator(max_distance, particle_definition):
collection = pp.PropagationUtilityCollection() collection = pp.PropagationUtilityCollection()
collection.displacement = pp.make_displacement(crosssection, True) collection.displacement = pp.make_displacement(crosssection, True)
collection.interaction = pp.make_interaction(crosssection, True) collection.interaction = pp.make_interaction(crosssection, True)
collection.decay = pp.make_decay(crosssection, True) collection.decay = pp.make_decay(crosssection, particle_definition, True)
collection.time = pp.make_time(crosssection, True) collection.time = pp.make_time(crosssection, particle_definition, True)
utility = pp.PropagationUtility(collection=collection) utility = pp.PropagationUtility(collection=collection)
geometry = pp.geometry.Sphere(pp.Vector3D(), max_distance, 0) geometry = pp.geometry.Sphere(pp.Cartesian3D(), max_distance)
density_distr = pp.density_distribution.density_homogeneous( density_distr = pp.density_distribution.density_homogeneous(
PROPOSAL_TARGET.mass_density) PROPOSAL_TARGET.mass_density)
......
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