Table of Contents
This algorithm provides a method for making a rotating reference frame transformation, which provides a useful method of visualising high frequency signals. Since the transformation leaves a comparatively low frequency signal, the data can readily be rebinned to improve counting statistics - this is particularly useful during an experiment to provide an early check on the form of the measured signal. Data analysis (for example least squares fitting) to study the form of the signal relaxation can be completed in either the laboratory or rotating reference frame according to preference.
This algorithm is frequently run after making a phase quadrature transformation using the algorithm PhaseQuad v1. Both algorithms are fully described in the article by T.M. Riseman and J.H. Brewer [Hyp. Int., 65, (1990), 1107].
Name | Direction | Type | Default | Description |
---|---|---|---|---|
InputWorkspace | Input | MatrixWorkspace | Mandatory | Name of the input workspace containing the spectra in the lab frame |
OutputWorkspace | Output | MatrixWorkspace | Mandatory | Name of the output workspace containing the spectra in the RRF |
Frequency | Input | number | 0 | Frequency of the oscillations |
FrequencyUnits | Input | string | MHz | The frequency units. Allowed values: [‘MHz’, ‘Gauss’, ‘Mrad/s’] |
Phase | Input | number | 0 | Phase accounting for any misalignment of the counters |
Assuming InputWorkspace contains the real and imaginary parts of the asymmetry in the lab-fram, the algorithm returns the muon polarization in the Rotating Reference Frame (RRF), given the frequency of the oscillations and the phase of the detectors, as given by the formula:
where and are the real and imaginary part of the asymmetry in the lab-frame, is the input frequency, and the input phase.
Example - Computing asymmetry in RRF
import math
# Create an input workspace with two spectra
datax = [i/100. for i in range(1,300)]
datay1 = [ math.cos(2*3.14159*i/100.) for i in range(1,299) ]
datay2 = [ math.sin(2*3.14159*i/100.) for i in range(1,299) ]
datay = datay1 + datay2
input = CreateWorkspace(dataX=datax, dataY=datay,Nspec=2,UnitX="TOF")
# Compute polarization in RRF
output = RRFMuon(input,1.0,"MHz",0)
print("%.1f" % output.readY(0)[0])
print("%.1f" % output.readY(1)[0])
Output:
1.0
-0.0
Categories: Algorithms | Muon