PhaseQuad v1

../_images/PhaseQuad-v1_dlg.png

PhaseQuad dialog.

Summary

Generates a quadrature phase signal.

Overview

This algorithm provides a method for combining signals from multiple detector segments to form two output signals in quadrature phase. It is of particular use when working with muon spin rotation signals measured using highly segmented detector arrays that are typically found at pulsed muon sources (ISIS has instruments with up to 600 detector elements). The method allows information from individual detectors to be combined in a way that makes full use of the dataset.

This algorithm is frequently run as a precursor to making a Rotating Reference Frame transformation of the dataset using the algorithm RRFMuon v1. Both algorithms are fully described in the article by T.M. Riseman and J.H. Brewer [Hyp. Int., 65, (1990), 1107].

Properties

Name Direction Type Default Description
InputWorkspace Input MatrixWorkspace Mandatory Name of the input workspace containing the spectra
PhaseTable Input TableWorkspace Mandatory Name of the table containing the detector phases and asymmetries
OutputWorkspace Output MatrixWorkspace Mandatory Name of the output workspace

Description

Assuming that the InputWorkspace contains measured counts as a function of time, and PhaseTable contains the detector phases and asymmetries, the algorithm returns a workspace containing two spectra (squashograms) as a function of the same time binning. PhaseTable is expected to have three columns, corresponding to the detector id, its asymmetry and its phase.

Usage

Note

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

Example - Computing squashograms:

from math import pi
Load(Filename='MUSR00022725.nxs', OutputWorkspace='MUSR00022725')
CropWorkspace(InputWorkspace='MUSR00022725', OutputWorkspace='MUSR00022725', XMin=0, EndWorkspaceIndex=63)

# Create a PhaseTable with some arbitrary detector information
tab = CreateEmptyTableWorkspace()
tab.addColumn('int', 'DetID')
tab.addColumn('double', 'Phase')
tab.addColumn('double', 'Asym')
for i in range(0,32):
    phi = 2*pi*i/32.
    tab.addRow([1, 0.2, phi])
for i in range(0,32):
    phi = 2*pi*i/32.
    tab.addRow([1, 0.2, phi])
ows = PhaseQuad(InputWorkspace='MUSR00022725', PhaseTable='tab')
print "Output workspace has", ows.getNumberHistograms(), "histograms"

Output:

Output workspace has 2 histograms

Categories: Algorithms | Muon

Source

C++ source: PhaseQuadMuon.cpp

C++ header: PhaseQuadMuon.h