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

MuonProcess v1

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’.

Workflow

../_images/MuonProcess-v1_wkflw.svg

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

C++ source: MuonProcess.cpp