6.0

\(\renewcommand\AA{\unicode{x212B}}\)

MuonProcess v1

../_images/MuonProcess-v1_dlg.png

MuonProcess dialog.

Summary

Processes and analyses Muon workspace.

Properties

Name Direction Type Default Description
InputWorkspace Input Workspace Mandatory Input workspace loaded from file (e.g. by LoadMuonNexus)
Mode Input string Combined Mode to run in. CorrectAndGroup applies dead time correction and grouping; Analyse changes bin offset, crops, rebins and calculates asymmetry; Combined does all of the above. Allowed values: [‘Combined’, ‘Analyse’, ‘CorrectAndGroup’]
SummedPeriodSet Input int list   Comma-separated list of periods to be summed
SubtractedPeriodSet Input int list   Comma-separated list of periods to be subtracted from the SummedPeriodSet
ApplyDeadTimeCorrection Input boolean False Whether dead time correction should be applied to loaded workspace
DeadTimeTable Input TableWorkspace   Table with dead time information, e.g. from LoadMuonNexus.Must be specified if ApplyDeadTimeCorrection is set true.
DetectorGroupingTable Input TableWorkspace   Table with detector grouping information, e.g. from LoadMuonNexus.
TimeZero Input number Optional Value used for Time Zero correction
LoadedTimeZero Input number Mandatory Time Zero value loaded from file, e.g. from LoadMuonNexus.
RebinParams Input dbl list   Params used for rebinning. If empty - rebinning is not done.
Xmin Input number Optional Minimal X value to include
Xmax Input number Optional Maximal X value to include
OutputType Input string PairAsymmetry What kind of workspace required for analysis. Allowed values: [‘PairAsymmetry’, ‘GroupAsymmetry’, ‘GroupCounts’]
PairFirstIndex Input number Optional Workspace index of the first pair group
PairSecondIndex Input number Optional Workspace index of the second pair group
Alpha Input number 1 Alpha value of the pair
GroupIndex Input number Optional Workspace index of the group
OutputWorkspace Output Workspace Mandatory An output workspace.
CropWorkspace Input boolean True Determines if the input workspace should be cropped at Xmax, Xmin is still aplied.
WorkspaceName Input string   The name of the input workspace

Description

The algorithm replicates the sequence of actions undertaken by MuonAnalysis in order to produce a Muon workspace ready for fitting. It is a workflow algorithm used internally by this interface.

It acts on a workspace loaded from a Muon NeXus file, most commonly by LoadMuonNexus v2.

Specifically:

  1. Apply dead time correction
  2. Group the workspaces
  3. Offset, crop and rebin the workspace
  4. Perform counts or asymmetry calculation to get the resulting workspace.

MuonProcess can be applied in three different modes. * CorrectAndGroup applies the dead time correction and groups the workspaces only, returning the output of this step without performing the rest of the analysis. * Analyse performs the second half only, i.e. offset time zero, crop, rebin and calculate asymmetry. * Combined performs the whole sequence of actions above (CorrectAndGroup followed by Analyse).

Analysis

The asymmetry is calculated from either one or several provided data acquisition period workspaces (only the first one is mandatory). When more than one are supplied, the algorithm merges the counts and then calculates the asymmetry.

The way in which period data is merged before the asymmetry calculation is determined by SummedPeriodSet and SubtractedPeriodSet. For example, setting SummedPeriodSet to “1,2” and SubtractedPeriodSet to “3,4” would result in the period arithmetic \((1+2)-(3+4)\).

The algorithm supports three output types:

  • PairAsymmetry - asymmetry is calculated for a given pair of groups, using the alpha value provided. The pair to use is specified via PairFirstIndex and PairSecondIndex.
  • GroupAsymmetry - asymmetry between the given group (specified via GroupIndex) and the Muon exponential decay is calculated.
  • GroupCount - no asymmetry is calculated, pure counts of the specified group (via GroupIndex) are used.

Note that the first part of the algorithm will have grouped the spectra in the input workspaces, hence the term ‘group’ is used here instead of ‘spectrum’.

Usage

Note

To run these usage examples please first download the usage data, and add these to your path. In Mantid this is done using Manage User Directories.

Note

For examples of applying custom dead times, please refer to ApplyDeadTimeCorr v1 documentation.

For examples of applying custom grouping, please refer to MuonGroupDetectors v1 documentation.

Example - Integrated pair asymmetry for MUSR run (combined mode):

# Begin by loading data into a workspace to process
load_result = LoadMuonNexus(Filename = "MUSR0015189.nxs",
                            DetectorGroupingTable = "grouping",
                            OutputWorkspace = "loadedWS")
loaded_time_zero = load_result[2] # time zero loaded from file
grouping = mtd['grouping']        # detector grouping loaded from file
grouping = grouping.getItem(0)    # grouping here is a WorkspaceGroup - use table for first period

MuonProcess(InputWorkspace = "loadedWS",
            Mode = "Combined",
            DetectorGroupingTable = grouping,
            SummedPeriodSet = "1",
            LoadedTimeZero = loaded_time_zero,
            TimeZero = 0.55,
            Xmin = 0.11,
            Xmax = 12,
            OutputType = "PairAsymmetry",
            PairFirstIndex = 0,
            PairSecondIndex = 1,
            Alpha = 1.0,
            OutputWorkspace = "MuonProcess_output")

processed = mtd['MuonProcess_output']
output_int = Integration(processed)

print('Integrated asymmetry for the run: {0:.3f}'.format(output_int.readY(0)[0]))

Output:

Integrated asymmetry for the run: 1.701

Example - Pair asymmetry for a single period (analysis only):

# Create example workspace
y = [1,2,3] + [4,5,6]
x = [1,2,3] * 2
first_period = CreateWorkspace(x, y, NSpec=2)
input = GroupWorkspaces(first_period)

# Grouping
grouping = CreateEmptyTableWorkspace()
grouping.addColumn('vector_int', 'detectors')
grouping.addRow([range(1,2)])

output = MuonProcess(InputWorkspace = input,
                     Mode = "Analyse",
                     DetectorGroupingTable = grouping,
                     SummedPeriodSet = "1",
                     LoadedTimeZero = 0,
                     OutputType = 'PairAsymmetry',
                     PairFirstIndex = 1,
                     PairSecondIndex = 0,
                     Alpha = 0.5)

print('Output: {}'.format(output.readY(0)))

Output:

Output: [ 0.77777778  0.66666667  0.6       ]

Example - Group asymmetry for two periods (analysis only):

# Create example workspaces
y1 = [100,50,10]
y2 = [150,20,1]
x = [1,2,3]
CreateWorkspace(x, y1,OutputWorkspace="first")
first_period = mtd["first"]
first_period.mutableRun().addProperty('goodfrm',10,True)
CreateWorkspace(x, y2,OutputWorkspace="second")
second_period = mtd["second"]
second_period.mutableRun().addProperty('goodfrm',10,True)

input = GroupWorkspaces([first_period, second_period])

# Grouping
grouping = CreateEmptyTableWorkspace()
grouping.addColumn('vector_int', 'detectors')
grouping.addRow([range(1,1)])
grouping.addRow([range(2,2)])

output = MuonProcess(InputWorkspace = input,
                     Mode = "Analyse",
                     DetectorGroupingTable = grouping,
                     SummedPeriodSet = 1,
                     SubtractedPeriodSet = 2,
                     LoadedTimeZero = 0,
                     OutputType = 'GroupAsymmetry',
                     GroupIndex = 0,Xmin=0,Xmax=4)

print('Output: {}'.format(output.readY(0)))

Output:

Output: [-0.44618444  0.54537717  0.24908794]

Categories: AlgorithmIndex | Workflow\Muon

Source

C++ header: MuonProcess.h (last modified: 2020-03-25)

C++ source: MuonProcess.cpp (last modified: 2020-04-07)