Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometer
Monte Carlo ray tracing simulations, of direct geometry spectrometers, have been particularly useful in instrument design and characterization. However, these tools can also be useful for experiment planning and analysis. To this end, the McStas Monte Carlo ray tracing model of SEQUOIA, the fine res...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
EDP Sciences
2015-01-01
|
| Series: | EPJ Web of Conferences |
| Online Access: | http://dx.doi.org/10.1051/epjconf/20158303006 |
| _version_ | 1818882824851161088 |
|---|---|
| author | Granroth G.E. Hahn S.E. |
| author_facet | Granroth G.E. Hahn S.E. |
| author_sort | Granroth G.E. |
| collection | DOAJ |
| description | Monte Carlo ray tracing simulations, of direct geometry spectrometers, have been particularly useful in instrument design and characterization. However, these tools can also be useful for experiment planning and analysis. To this end, the McStas Monte Carlo ray tracing model of SEQUOIA, the fine resolution fermi chopper spectrometer at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (ORNL), has been modified to include the time of flight resolution sample and detector components. With these components, the resolution ellipsoid can be calculated for any detector pixel and energy bin of the instrument. The simulation is split in two pieces. First, the incident beamline up to the sample is simulated for 1 × 1011 neutron packets (4 days on 30 cores). This provides a virtual source for the backend that includes the resolution sample and monitor components. Next, a series of detector and energy pixels are computed in parallel. It takes on the order of 30 s to calculate a single resolution ellipsoid on a single core. Python scripts have been written to transform the ellipsoid into the space of an oriented single crystal, and to characterize the ellipsoid in various ways. Though this tool is under development as a planning tool, we have successfully used it to provide the resolution function for convolution with theoretical models. Specifically, theoretical calculations of the spin waves in YFeO3 were compared to measurements taken on SEQUOIA. Though the overall features of the spectra can be explained while neglecting resolution effects, the variation in intensity of the modes is well described once the resolution is included. As this was a single sharp mode, the simulated half intensity value of the resolution ellipsoid was used to provide the resolution width. A description of the simulation, its use, and paths forward for this technique will be discussed. |
| first_indexed | 2024-12-19T15:23:54Z |
| format | Article |
| id | doaj.art-c55ebb52408c4c75a8d786c6389c2eae |
| institution | Directory Open Access Journal |
| issn | 2100-014X |
| language | English |
| last_indexed | 2024-12-19T15:23:54Z |
| publishDate | 2015-01-01 |
| publisher | EDP Sciences |
| record_format | Article |
| series | EPJ Web of Conferences |
| spelling | doaj.art-c55ebb52408c4c75a8d786c6389c2eae2022-12-21T20:15:56ZengEDP SciencesEPJ Web of Conferences2100-014X2015-01-01830300610.1051/epjconf/20158303006epjconf-Qens-Wins2014_03006Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometerGranroth G.E.0Hahn S.E.1Neutron Data Analysis and Visualization DivisionQuantum Condensed Matter Division, Oak Ridge National LaboratoryMonte Carlo ray tracing simulations, of direct geometry spectrometers, have been particularly useful in instrument design and characterization. However, these tools can also be useful for experiment planning and analysis. To this end, the McStas Monte Carlo ray tracing model of SEQUOIA, the fine resolution fermi chopper spectrometer at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (ORNL), has been modified to include the time of flight resolution sample and detector components. With these components, the resolution ellipsoid can be calculated for any detector pixel and energy bin of the instrument. The simulation is split in two pieces. First, the incident beamline up to the sample is simulated for 1 × 1011 neutron packets (4 days on 30 cores). This provides a virtual source for the backend that includes the resolution sample and monitor components. Next, a series of detector and energy pixels are computed in parallel. It takes on the order of 30 s to calculate a single resolution ellipsoid on a single core. Python scripts have been written to transform the ellipsoid into the space of an oriented single crystal, and to characterize the ellipsoid in various ways. Though this tool is under development as a planning tool, we have successfully used it to provide the resolution function for convolution with theoretical models. Specifically, theoretical calculations of the spin waves in YFeO3 were compared to measurements taken on SEQUOIA. Though the overall features of the spectra can be explained while neglecting resolution effects, the variation in intensity of the modes is well described once the resolution is included. As this was a single sharp mode, the simulated half intensity value of the resolution ellipsoid was used to provide the resolution width. A description of the simulation, its use, and paths forward for this technique will be discussed.http://dx.doi.org/10.1051/epjconf/20158303006 |
| spellingShingle | Granroth G.E. Hahn S.E. Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometer EPJ Web of Conferences |
| title | Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometer |
| title_full | Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometer |
| title_fullStr | Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometer |
| title_full_unstemmed | Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometer |
| title_short | Monte Carlo simulation of the resolution volume for the SEQUOIA spectrometer |
| title_sort | monte carlo simulation of the resolution volume for the sequoia spectrometer |
| url | http://dx.doi.org/10.1051/epjconf/20158303006 |
| work_keys_str_mv | AT granrothge montecarlosimulationoftheresolutionvolumeforthesequoiaspectrometer AT hahnse montecarlosimulationoftheresolutionvolumeforthesequoiaspectrometer |