essdiffraction 26.5.2__tar.gz → 26.6.1__tar.gz

{essdiffraction-26.5.2/src/essdiffraction.egg-info → essdiffraction-26.6.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: essdiffraction
-Version: 26.5.2
+Version: 26.6.1
 Summary: Diffraction data reduction for the European Spallation Source
 Author: Scipp contributors
 License-Expression: BSD-3-Clause
@@ -21,39 +21,35 @@ Requires-Python: >=3.11
 Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: dask>=2022.1.0
-Requires-Dist: essreduce>=26.4.0
-Requires-Dist: graphviz
+Requires-Dist: graphviz>=0.20
+Requires-Dist: essreduce>=26.6.0
 Requires-Dist: numpy>=2
-Requires-Dist: plopp>=26.2.0
+Requires-Dist: plopp>=26.5.0
 Requires-Dist: pythreejs>=2.4.1
-Requires-Dist: sciline>=25.04.1
-Requires-Dist: scipp>=25.11.0
-Requires-Dist: scippneutron>=26.3.0
-Requires-Dist: scippnexus>=23.12.0
 Requires-Dist: tof>=25.12.0
 Requires-Dist: ncrystal[cif]>=4.1.0
-Requires-Dist: spglib!=2.7
+Requires-Dist: spglib!=2.7,>=2.0.0
 Provides-Extra: test
 Requires-Dist: pandas>=2.1.2; extra == "test"
 Requires-Dist: pooch>=1.5; extra == "test"
 Requires-Dist: pytest>=7.0; extra == "test"
 Requires-Dist: ipywidgets>=8.1.7; extra == "test"
 Provides-Extra: docs
-Requires-Dist: autodoc-pydantic; extra == "docs"
-Requires-Dist: ipykernel; extra == "docs"
-Requires-Dist: ipympl; extra == "docs"
-Requires-Dist: ipython!=8.7.0; extra == "docs"
-Requires-Dist: myst-parser; extra == "docs"
-Requires-Dist: nbsphinx; extra == "docs"
-Requires-Dist: pandas; extra == "docs"
-Requires-Dist: pooch; extra == "docs"
+Requires-Dist: autodoc-pydantic>=2.0.0; extra == "docs"
+Requires-Dist: ipykernel>=6.20.0; extra == "docs"
+Requires-Dist: ipympl>=0.9.0; extra == "docs"
+Requires-Dist: ipython!=8.7.0,>=8.8.0; extra == "docs"
+Requires-Dist: myst-parser>=2.0.0; extra == "docs"
+Requires-Dist: nbsphinx>=0.9.3; extra == "docs"
+Requires-Dist: pandas>=2.1.2; extra == "docs"
+Requires-Dist: pooch>=1.5; extra == "docs"
 Requires-Dist: pydata-sphinx-theme>=0.14; extra == "docs"
-Requires-Dist: sphinx; extra == "docs"
-Requires-Dist: sphinx-autodoc-typehints; extra == "docs"
-Requires-Dist: sphinx-copybutton; extra == "docs"
-Requires-Dist: sphinx-design; extra == "docs"
-Requires-Dist: sphinxcontrib-bibtex; extra == "docs"
-Requires-Dist: pyarrow; extra == "docs"
+Requires-Dist: sphinx>=7; extra == "docs"
+Requires-Dist: sphinx-autodoc-typehints>=1.24.0; extra == "docs"
+Requires-Dist: sphinx-copybutton>=0.5.2; extra == "docs"
+Requires-Dist: sphinx-design>=0.5.0; extra == "docs"
+Requires-Dist: sphinxcontrib-bibtex>=2.5.0; extra == "docs"
+Requires-Dist: pyarrow>=12.0.0; extra == "docs"
 Dynamic: license-file
 
 [![Contributor Covenant](https://img.shields.io/badge/Contributor%20Covenant-2.1-4baaaa.svg)](CODE_OF_CONDUCT.md)

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/docs/api-reference/index.md

@@ -1,8 +1,25 @@
 # API Reference
 
+## ESSdiffraction
+
+### Submodules
+
+```{eval-rst}
+.. currentmodule:: ess.diffraction
+
+.. autosummary::
+   :toctree: ../generated/modules
+   :template: module-template.rst
+   :recursive:
+
+   pdf
+   peaks
+```
+
 ## ESSpowder
 
 ### Module Attributes
+
 ```{eval-rst}
 .. currentmodule:: ess.powder
 
@@ -37,9 +54,7 @@
    grouping
    logging
    masking
-   peaks
    smoothing
-   transform
    types
 ```
 
@@ -47,7 +62,6 @@
 
 ### Workflows
 
-
 ```{eval-rst}
 .. currentmodule:: ess.dream
 
@@ -103,15 +117,20 @@
 
 ### Workflows
 
-
 ```{eval-rst}
 .. currentmodule:: ess.beer
 
 .. autosummary::
    :toctree: ../generated/functions
 
+   BeerMcStasWorkflowPulseShaping
+   BeerMcStasWorkflowPulseShapingAnalytical
    BeerModMcStasWorkflowKnownPeaks
    BeerModMcStasWorkflow
+   BeerPowderWorkflow
+   BeerPowderWorkflowAnalytical
+   BeerPowderMcStasWorkflow
+   BeerPowderMcStasWorkflowAnalytical
 ```
 
 ### Top-level functions
@@ -122,7 +141,8 @@
 .. autosummary::
    :toctree: ../generated/functions
 
-    load_beer_mcstas
+   dhkl_peaks_from_cif
+   load_beer_mcstas
 ```
 
 ### Submodules
@@ -137,6 +157,7 @@
    conversions
    data
    io
+   mcstas
    workflow
 ```
 
@@ -144,7 +165,6 @@
 
 ### Workflows
 
-
 ```{eval-rst}
 .. currentmodule:: ess.magic
 

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/docs/conf.py

@@ -249,8 +249,12 @@ nbsphinx_execute_arguments = [
 # build already tests if those plots can be made.
 # So we simply disable plots in doctests.
 doctest_global_setup = """
+# Private import because of https://github.com/mctools/ncrystal/issues/361
+from NCrystal._common import set_ncrystal_print_fct
 import numpy as np
 
+set_ncrystal_print_fct(lambda *args, **kwargs: None)
+
 try:
     import scipp as sc
 

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/docs/user-guide/dream/dream-make-wavelength-lookup-table.ipynb

@@ -18,8 +18,9 @@
    "outputs": [],
    "source": [
     "import scipp as sc\n",
+    "import scippnexus as snx\n",
     "from ess.reduce import unwrap\n",
-    "from ess.reduce.nexus.types import AnyRun\n",
+    "from ess.reduce.nexus.types import AnyRun, Position\n",
     "from ess.dream.beamline import InstrumentConfiguration, choppers"
    ]
   },
@@ -62,14 +63,13 @@
    "source": [
     "wf = unwrap.LookupTableWorkflow()\n",
     "\n",
-    "wf[unwrap.LtotalRange] = sc.scalar(5.0, unit=\"m\"), sc.scalar(80.0, unit=\"m\")\n",
+    "wf[unwrap.LtotalRange[AnyRun, snx.NXdetector]] = sc.scalar(5.0, unit=\"m\"), sc.scalar(80.0, unit=\"m\")\n",
     "wf[unwrap.NumberOfSimulatedNeutrons] = 200_000  # Increase this number for more reliable results\n",
-    "wf[unwrap.SourcePosition] = sc.vector([0, 0, 0], unit='m')\n",
+    "wf[Position[snx.NXsource, AnyRun]] = sc.vector([0, 0, 0], unit='m')\n",
     "wf[unwrap.DiskChoppers[AnyRun]] = disk_choppers\n",
     "wf[unwrap.DistanceResolution] = sc.scalar(0.1, unit=\"m\")\n",
     "wf[unwrap.TimeResolution] = sc.scalar(250.0, unit='us')\n",
     "wf[unwrap.PulsePeriod] = 1.0 / sc.scalar(14.0, unit=\"Hz\")\n",
-    "wf[unwrap.PulseStride] = 1\n",
     "wf[unwrap.PulseStrideOffset] = None"
    ]
   },
@@ -88,7 +88,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "table = wf.compute(unwrap.LookupTable)\n",
+    "table = wf.compute(unwrap.LookupTable[AnyRun, snx.NXdetector])\n",
     "table.array"
    ]
   },

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/pyproject.toml

@@ -30,18 +30,14 @@ requires-python = ">=3.11"
 # Make sure to list one dependency per line.
 dependencies = [
   "dask>=2022.1.0",
-  "essreduce>=26.4.0",
-  "graphviz",
+  "graphviz>=0.20",
+  "essreduce>=26.6.0",
   "numpy>=2",
-  "plopp>=26.2.0",
+  "plopp>=26.5.0",
   "pythreejs>=2.4.1",
-  "sciline>=25.04.1",
-  "scipp>=25.11.0",
-  "scippneutron>=26.3.0",
-  "scippnexus>=23.12.0",
   "tof>=25.12.0",
   "ncrystal[cif]>=4.1.0",
-  "spglib!=2.7",  # https://github.com/mctools/ncrystal/issues/320
+  "spglib>=2.0.0,!=2.7",  # https://github.com/mctools/ncrystal/issues/320
 ]
 
 dynamic = ["version"]
@@ -54,23 +50,23 @@ test = [
     "ipywidgets>=8.1.7",
 ]
 docs = [
-    "autodoc-pydantic",
-    "ipykernel",
-    "ipympl",
+    "autodoc-pydantic>=2.0.0",
+    "ipykernel>=6.20.0",
+    "ipympl>=0.9.0",
     # Breaks syntax highlighting in Jupyter code cells.
-    "ipython!=8.7.0",
-    "myst-parser",
-    "nbsphinx",
-    "pandas",
-    "pooch",
+    "ipython>=8.8.0,!=8.7.0",
+    "myst-parser>=2.0.0",
+    "nbsphinx>=0.9.3",
+    "pandas>=2.1.2",
+    "pooch>=1.5",
     "pydata-sphinx-theme>=0.14",
-    "sphinx",
-    "sphinx-autodoc-typehints",
-    "sphinx-copybutton",
-    "sphinx-design",
-    "sphinxcontrib-bibtex",
+    "sphinx>=7",
+    "sphinx-autodoc-typehints>=1.24.0",
+    "sphinx-copybutton>=0.5.2",
+    "sphinx-design>=0.5.0",
+    "sphinxcontrib-bibtex>=2.5.0",
     # needed by pandas < 3.0
-    "pyarrow",
+    "pyarrow>=12.0.0",
 ]
 
 [project.urls]

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/src/ess/beer/__init__.py

@@ -7,12 +7,17 @@ Components for BEER
 
 import importlib.metadata
 
-from .io import load_beer_mcstas
+from .mcstas import load_beer_mcstas
 from .peakfinding import dhkl_peaks_from_cif
 from .workflow import (
     BeerMcStasWorkflowPulseShaping,
+    BeerMcStasWorkflowPulseShapingAnalytical,
     BeerModMcStasWorkflow,
     BeerModMcStasWorkflowKnownPeaks,
+    BeerPowderMcStasWorkflow,
+    BeerPowderMcStasWorkflowAnalytical,
+    BeerPowderWorkflow,
+    BeerPowderWorkflowAnalytical,
     default_parameters,
 )
 
@@ -25,8 +30,13 @@ del importlib
 
 __all__ = [
     'BeerMcStasWorkflowPulseShaping',
+    'BeerMcStasWorkflowPulseShapingAnalytical',
     'BeerModMcStasWorkflow',
     'BeerModMcStasWorkflowKnownPeaks',
+    'BeerPowderMcStasWorkflow',
+    'BeerPowderMcStasWorkflowAnalytical',
+    'BeerPowderWorkflow',
+    'BeerPowderWorkflowAnalytical',
     '__version__',
     'default_parameters',
     'dhkl_peaks_from_cif',

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/src/ess/beer/clustering.py

@@ -2,7 +2,7 @@ import scipp as sc
 from ess.powder.types import RunType
 from scipy.signal import find_peaks, medfilt
 
-from .conversions import tof_from_t0_estimate_graph
+from .conversions import tof_from_nominal_time_at_chopper_graph
 from .types import (
     GeometryCoordTransformGraph,
     RawDetector,
@@ -13,7 +13,7 @@ from .types import (
 def cluster_events_by_streak(
     da: RawDetector[RunType], gg: GeometryCoordTransformGraph
 ) -> StreakClusteredData[RunType]:
-    graph = tof_from_t0_estimate_graph(da, gg)
+    graph = tof_from_nominal_time_at_chopper_graph(da, gg)
 
     da = da.transform_coords(['dspacing'], graph=graph)
     da.bins.coords['coarse_d'] = da.bins.coords.pop('dspacing').to(unit='angstrom')

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/src/ess/beer/conversions.py

@@ -19,7 +19,6 @@ def compute_wavelength_in_each_cluster(
     da: StreakClusteredData[RunType],
     chopper_delay: WavelengthDefinitionChopperDelay,
     mod_period: ModulationPeriod,
-    graph: GeometryCoordTransformGraph,
 ) -> WavelengthDetector[RunType]:
     """Fits a line through each cluster, the intercept of the line is t0.
     The line is fitted using linear regression with an outlier removal procedure.
@@ -38,7 +37,7 @@ def compute_wavelength_in_each_cluster(
     if isinstance(da, sc.DataGroup):
         return sc.DataGroup(
             {
-                k: compute_wavelength_in_each_cluster(v, mod_period)
+                k: compute_wavelength_in_each_cluster(v, chopper_delay, mod_period)
                 for k, v in da.items()
             }
         )
@@ -133,7 +132,7 @@ def _compute_d_given_list_of_peaks(
         / (scipp.constants.h / scipp.constants.m_n)
     ).to(unit=f'{time_of_arrival.unit}/angstrom')
     for dhkl in dhkl_list:
-        dt = sc.abs(t - dhkl * const)
+        dt = sc.abs(t - dhkl * const - pulse_length / 2)
         dt_in_range = dt < pulse_length / 2
         no_dt_found = sc.isnan(dtfound)
         dtfound = sc.where(dt_in_range, sc.where(no_dt_found, dt, dtfound), dtfound)
@@ -191,34 +190,16 @@ def _tof_from_dhkl(
     return out
 
 
-def t0_estimate(
-    wavelength_estimate: sc.Variable,
-    source_to_wavelength_definition_chopper_distance: sc.Variable,
-) -> sc.Variable:
-    """
-    Computes the time a neutron reaches a chopper at
-    ``source_to_wavelength_chopper_distance`` distance from the source
-    if it has wavelength ``wavelength_estimate``.
-    """
-    return (
-        sc.constants.m_n
-        / sc.constants.h
-        * wavelength_estimate
-        * source_to_wavelength_definition_chopper_distance.to(
-            unit=wavelength_estimate.unit
-        )
-    ).to(unit='s')
-
-
-def tof_from_t0_estimate_graph(
+def tof_from_nominal_time_at_chopper_graph(
     da: RawDetector[RunType],
     gg: GeometryCoordTransformGraph,
 ) -> ElasticCoordTransformGraph[RunType]:
     """Graph for computing ``wavelength`` in pulse shaping chopper modes."""
     return {
         **gg,
-        't0': t0_estimate,
-        'tof': lambda time_of_arrival, t0: time_of_arrival - t0,
+        'tof': lambda time_of_arrival, nominal_time_at_chopper: (
+            time_of_arrival - nominal_time_at_chopper
+        ),
         'time_of_arrival': time_of_arrival,
     }
 
@@ -287,7 +268,7 @@ convert_from_known_peaks_providers = (
 )
 convert_pulse_shaping = (
     geometry_graph,
-    tof_from_t0_estimate_graph,
+    tof_from_nominal_time_at_chopper_graph,
     wavelength_detector,
 )
 providers = (compute_wavelength_in_each_cluster, geometry_graph)

{essdiffraction-26.5.2 → essdiffraction-26.6.1}/src/ess/beer/data.py

@@ -31,7 +31,7 @@ _registry = make_registry(
         # Simulations from new model with corrected(?) moderator to detector
         # distance.
         # For correct reduction you need to use
-        # beer.io.mcstas_chopper_delay_from_mode_new_simulations
+        # beer.mcstas.mcstas_chopper_delay_from_mode_new_simulations
         # to obtain the correct WavelengthDefinitionChopperDelay for these files.
         "silicon-mode10-new-model.h5": "md5:98500830f27700fc719634e1acd49944",
         "silicon-mode16-new-model.h5": "md5:393f9287e7d3f97ceedbe64343918413",

essdiffraction-26.6.1/src/ess/beer/io.py

@@ -0,0 +1,21 @@
+# SPDX-License-Identifier: BSD-3-Clause
+# Copyright (c) 2025 Scipp contributors (https://github.com/scipp)
+
+"""Backward-compatible import location for BEER McStas helpers."""
+
+import warnings
+
+from . import mcstas as _mcstas
+from .mcstas import *  # noqa: F403
+from .mcstas import __all__  # noqa: F401
+
+warnings.warn(
+    "ess.beer.io is deprecated; use ess.beer.mcstas instead.",
+    DeprecationWarning,
+    stacklevel=2,
+)
+
+
+def __getattr__(name: str):
+    """Return attributes from the replacement :mod:`ess.beer.mcstas` module."""
+    return getattr(_mcstas, name)

essdiffraction-26.6.1/src/ess/beer/mcstas/__init__.py

@@ -0,0 +1,42 @@
+# SPDX-License-Identifier: BSD-3-Clause
+# Copyright (c) 2026 Scipp contributors (https://github.com/scipp)
+
+"""BEER McStas simulation helpers."""
+
+from .beamline import MCSTAS_T_OFFSET, PulseShapingMode, simulation_choppers
+from .load import (
+    MCSTAS_PULSE_SHAPING_MODES,
+    load_beer_mcstas,
+    load_beer_mcstas_monitor,
+    load_beer_mcstas_monitor_provider,
+    load_beer_mcstas_provider,
+    mcstas_chopper_delay_from_mode,
+    mcstas_chopper_delay_from_mode_new_simulations,
+    mcstas_detector_ltotal,
+    mcstas_modulation_period_from_mode,
+    mcstas_providers,
+    mcstas_pulse_shaping_choppers,
+    mcstas_source_position,
+    pulse_shaping_mcstas_providers,
+    pulse_shaping_mode_from_mcstas_mode,
+)
+
+__all__ = [
+    'MCSTAS_PULSE_SHAPING_MODES',
+    'MCSTAS_T_OFFSET',
+    'PulseShapingMode',
+    'load_beer_mcstas',
+    'load_beer_mcstas_monitor',
+    'load_beer_mcstas_monitor_provider',
+    'load_beer_mcstas_provider',
+    'mcstas_chopper_delay_from_mode',
+    'mcstas_chopper_delay_from_mode_new_simulations',
+    'mcstas_detector_ltotal',
+    'mcstas_modulation_period_from_mode',
+    'mcstas_providers',
+    'mcstas_pulse_shaping_choppers',
+    'mcstas_source_position',
+    'pulse_shaping_mcstas_providers',
+    'pulse_shaping_mode_from_mcstas_mode',
+    'simulation_choppers',
+]

essdiffraction-26.6.1/src/ess/beer/mcstas/beamline.py

@@ -0,0 +1,201 @@
+# SPDX-License-Identifier: BSD-3-Clause
+# Copyright (c) 2026 Scipp contributors (https://github.com/scipp)
+
+import enum
+from types import MappingProxyType
+from typing import TypedDict
+
+import scipp as sc
+from scippneutron.chopper import DiskChopper
+
+
+class PulseShapingMode(enum.Enum):
+    """Pulse-shaping chopper modes for BEER."""
+
+    ps0 = "PS0"
+    ps1 = "PS1"
+    ps2 = "PS2"
+    ps3 = "PS3"
+    ds1 = "DS1"
+
+
+class _ChopperParameters(TypedDict):
+    frequency: float
+    phase: float
+    distance: float
+    open: list[float]
+    close: list[float]
+
+
+Hz = sc.Unit("Hz")
+deg = sc.Unit("deg")
+meter = sc.Unit("m")
+
+MCSTAS_T_OFFSET = sc.scalar(1.6, unit="ms")
+"""Time offset applied by the BEER McStas simulation source model."""
+
+_HIGH_FLUX: dict[str, _ChopperParameters] = {
+    "PSC1": {
+        "frequency": 168.0,
+        "phase": 318.6929881679336,
+        "distance": 6.450,
+        "open": [0.0],
+        "close": [144.0],
+    },
+    "PSC3": {
+        "frequency": -168.0,
+        "phase": -318.6929881679336,
+        "distance": 7.375,
+        "open": [0.0],
+        "close": [144.0],
+    },
+    "FC1A": {
+        "frequency": -28.0,
+        "phase": -18.44878787209148,
+        "distance": 8.283,
+        "open": [0.0],
+        "close": [72.0],
+    },
+    "FC2A": {
+        "frequency": -14.0,
+        "phase": -134.52965314925247,
+        "distance": 79.975,
+        "open": [0.0],
+        "close": [175.0],
+    },
+}
+
+_parameters: dict[PulseShapingMode, dict[str, _ChopperParameters]] = {
+    PulseShapingMode.ps0: _HIGH_FLUX,
+    PulseShapingMode.ps1: _HIGH_FLUX,
+    PulseShapingMode.ps2: {
+        "PSC1": {
+            "frequency": 168.0,
+            "phase": 310.265456971683,
+            "distance": 6.450,
+            "open": [0.0],
+            "close": [144.0],
+        },
+        "PSC2": {
+            "frequency": -168.0,
+            "phase": -310.265456971683,
+            "distance": 6.850,
+            "open": [0.0],
+            "close": [144.0],
+        },
+        "FC1A": {
+            "frequency": -28.0,
+            "phase": -18.44878787209148,
+            "distance": 8.283,
+            "open": [0.0],
+            "close": [72.0],
+        },
+        "FC2A": {
+            "frequency": -14.0,
+            "phase": -134.52965314925247,
+            "distance": 79.975,
+            "open": [0.0],
+            "close": [175.0],
+        },
+    },
+    PulseShapingMode.ps3: {
+        "PSC1": {
+            "frequency": 168.0,
+            "phase": 307.05496889692084,
+            "distance": 6.450,
+            "open": [0.0],
+            "close": [144.0],
+        },
+        "PSC2": {
+            "frequency": -168.0,
+            "phase": -307.05496889692084,
+            "distance": 6.650,
+            "open": [0.0],
+            "close": [144.0],
+        },
+        "FC1A": {
+            "frequency": -28.0,
+            "phase": -18.44878787209148,
+            "distance": 8.283,
+            "open": [0.0],
+            "close": [72.0],
+        },
+        "FC2A": {
+            "frequency": -14.0,
+            "phase": -134.52965314925247,
+            "distance": 79.975,
+            "open": [0.0],
+            "close": [175.0],
+        },
+    },
+    PulseShapingMode.ds1: {
+        "PSC1": {
+            "frequency": 168.0,
+            "phase": 318.6929881679336,
+            "distance": 6.450,
+            "open": [0.0],
+            "close": [144.0],
+        },
+        "PSC3": {
+            "frequency": -168.0,
+            "phase": -318.6929881679336,
+            "distance": 7.375,
+            "open": [0.0],
+            "close": [144.0],
+        },
+        "FC1A": {
+            "frequency": -14.0,
+            "phase": -3.22439393604574,
+            "distance": 8.283,
+            "open": [0.0],
+            "close": [72.0],
+        },
+        "FC1B": {
+            "frequency": -63.0,
+            "phase": -46.41910994173803,
+            "distance": 8.317,
+            "open": [0.0],
+            "close": [180.0],
+        },
+        "FC2B": {
+            "frequency": -7.0,
+            "phase": -68.58171174285046,
+            "distance": 80.025,
+            "open": [0.0],
+            "close": [85.0],
+        },
+    },
+}
+
+
+def simulation_choppers(
+    mode: PulseShapingMode, source_position: sc.Variable
+) -> MappingProxyType[str, DiskChopper]:
+    """
+    Dict of ESS BEER McStas choppers for the selected pulse-shaping mode.
+
+    We make the chopper information available in this way as loading it directly from
+    the NeXus files is currently not available for these simulated BEER data.
+
+    Parameters
+    ----------
+    mode:
+        BEER pulse-shaping chopper mode.
+    source_position:
+        Position of the source in the coordinate system of the choppers.
+        The raw chopper positions are defined relative to the source position.
+    """
+    return MappingProxyType(
+        {
+            key: DiskChopper(
+                frequency=ch["frequency"] * Hz,
+                beam_position=sc.scalar(0.0, unit="deg"),
+                phase=ch["phase"] * deg,
+                axle_position=sc.vector(value=[0, 0, ch["distance"]], unit="m")
+                + source_position,
+                slit_begin=sc.array(dims=["cutout"], values=ch["open"], unit="deg"),
+                slit_end=sc.array(dims=["cutout"], values=ch["close"], unit="deg"),
+            )
+            for key, ch in _parameters[mode].items()
+        }
+    )