McStas logo McStas - A neutron ray-trace simulation package ILL; PSI; ESS Niels Bohr Institute DTU Physics NEXMAP


About McStas
 Conditions of use
 Project funding

 Other Downloads (share)

Mailing list

Search web/mailinglist

 Wiki (GitHub)
 McStas manual


Report bugs (GitHub)

Code-repository (GitHub)

McStas: ESS_butterfly

[ Identification | Description | Input parameters | Links ]

The ESS_butterfly Component

ESS butterfly moderator, 2016 revision


  • Site:
  • Author: Peter Willendrup and Esben Klinkby
  • Origin: DTU
  • Date: August-September 2016


ESS butterfly moderator with automatic choice of coordinate system, with origin
placed at relevant "Moderator Focus Coordinate System" depending on sector location.

To select beamport N 5 simply use

COMPONENT Source = ESS_butterfly(sector="N",beamline=5,Lmin=0.1,Lmax=20,dist=2,
cold_frac=0.5, yheight=0.03,focus_xw=0.1, focus_yh=0.1)

The geometry corresponds correctly to the latest release of the butterfly moderator,
including changes warranted by the ESS CCB in July 2016, the so called BF1 type moderator.
A set of official release documents are available with this component, see the benchmarking
website mentioned below.

Brilliances, geometry adapted from earlier BF2 design
The geometry and brightness data implemented in the McStas ESS source component ESS_butterfly.comp,
are released as an updated component library for McStas 2.3, as well as a stand alone archive for
use with earlier versions of McStas.

The following features are worth highlighting:
  • The brightness data are still based on last years MCNP calculations, based on the Butterfly 2 geometry. As a result, the spatial variation of the brightness across the moderator face should be considered to have an uncertainty of the order of 10%. Detailed information on the reasoning behind the change to the Butterfly 1 geometry can be found in [1] and detailed information on horizontal spatial brightness variation can be found in [2]. The spectral shape has been checked and has not changed significantly.
  • A scaling factor has been introduced to in order to account for the decrease in brightness since 2015. To accommodate the influence of the changed geometry, this scaling factor has been applied independently for the cold and thermal contributions and is beamline dependent. It is adjusted to agree with the spectrally-integrated 6cm width data shown in [1],Figure 3.
  • To allow future user adjustments of brilliance, the scalar parameters c_performance and t_performance have been implemented. For now, we recommend to keep these at their default value of 1.0.
  • The geometry has been updated to correspond within about 2 mm to the geometry described in [1]. This has been done by ensuring that the position and apparent width of the moderators correspond to [1],Figure 2, which has been derived from current MCNP butterfly 1 model.
  • The beamport is now defined directly by its sector and number (e.g. 'W' and '5'), rather than giving the angle, as before. [1],Figure 5 shows the geometry of the moderator2, beamport insert and beamline axis for beamline W5. Since the underlying data is still from last years MCNP run, when the brightness was calculated at 10-degree intervals, this means that the spectral curve for the nearest beamport on the grid 5,15,25,35,45,55 degrees is used. The use of this grid has no effect on the accuracy of the geometry or brilliance because of the above- mentioned beamline-dependent adjustments to the brilliance and geometry. See the website [3] for details.
As before, the beamports all originate at the focal point of the sector. The beamline will in almost all cases be horizontally tilted in order to view the cold or thermal moderator, which should be done using an Arm component.

We expect to release an MCNP-event-based source model later in 2016, and possibly also new set of brilliance functions for ESS_butterfly.comp. These are expected to include more realistic brilliances in terms of variation across sectors and potentially also performance losses due to engineering reality. Engineering reality An ad-hoc method for future implementation of "engineering reality" is included, use the "c_performance/t_performance" parameters to down-scale performance uniformly across all wavelengths. References:

  1. Release document "Update to ESS Moderators, latest version"
  2. Release document "Description and performance of the new baseline ESS moderators, latest version"
  3. benchmarking website with comparative McStas-MCNP figures
  4. html-based, interactive 3D model of moderators and monolith, as seen from beamline N4.
  5. Source code for ESS_butterfly.comp at GitHub.

Input parameters

Parameters in boldface are required; the others are optional.
sectorstrDefines the 'sector' of your instrument position. Valid values are "N","S","E" and "W""N"
beamline1Defines the 'beamline number' of your instrument position. Valid values are 1..10 or 1..11 depending on sector1
yheightmDefines the moderator height. Valid values are 0.03 m and 0.06 m0.03
cold_frac1Defines the statistical fraction of events emitted from the cold part of the moderator0.5
target_index1Relative index of component to focus at, e.g. next is +1 this is used to compute 'dist' automatically.0
distmDistance from origin to focusing rectangle; at (0,0,dist) - alternatively use target_index0
focus_xwmWidth of focusing rectangle0
focus_yhmHeight of focusing rectangle0
c_performance1Cold brilliance scalar performance multiplicator c_performance > 01
t_performance1Thermal brilliance scalar performance multiplicator t_performance > 01
LminAAMinimum wavelength simulated
LmaxAAMaximum wavelength simulated
tmax_multiplier1Defined maximum emission time at moderator, tmax= tmax_multiplier * ESS_PULSE_DURATION.3
n_pulses1Number of pulses simulated. 0 and 1 creates one pulse.1
acc_powerMWAccelerator power in MW5
tfocus_distmPosition of time focusing window along z axis0
tfocus_timesTime position of time focusing window0
tfocus_widthsTime width of time focusing window0


[ Identification | Description | Input parameters | Links ]

Generated on 2023-09-19 20:01:36

Last Modified: Tuesday, 16-Apr-2024 20:05:21 CEST
Search website mailinglist archive GitHub repos