LoadFullprofFile v1

../_images/LoadFullprofFile-v1_dlg.png

LoadFullprofFile dialog.

Summary

Load file generated by Fullprof.

Properties

Name Direction Type Default Description
Filename Input string Mandatory Name of [http://www.ill.eu/sites/fullprof/ Fullprof] .hkl or .prf file. Allowed extensions: [‘.hkl’, ‘.prf’, ‘.dat’]
PeakParameterWorkspace Output TableWorkspace Mandatory Name of table workspace containing peak parameters from .hkl file.
OutputWorkspace Output MatrixWorkspace Mandatory Name of data workspace containing the diffraction pattern in .prf file.

Description

This algorithm is to import Fullprof .irf file (peak parameters) and .hkl file (reflections) and record the information to TableWorkspaces, which serve as the inputs for algorithm LeBailFit.

Format of Instrument parameter TableWorkspace

Instrument parameter TableWorkspace contains all the peak profile parameters imported from Fullprof .irf file.

Presently these are the peak profiles supported

* Thermal neutron back to back exponential convoluted with pseudo-voigt (profile No. 10 in Fullprof)

Each row in TableWorkspace corresponds to one profile parameter.

Columns include Name, Value, FitOrTie, Min, Max and StepSize.

Format of reflection TableWorkspace

Each row of this workspace corresponds to one diffraction peak. The information contains the peak’s Miller index and (local) peak profile parameters of this peak. For instance of a back-to-back exponential convoluted with Gaussian peak, the peak profile parameters include Alpha, Beta, Sigma, centre and height.

How to use algorithm with other algorithms

This algorithm is designed to work with other algorithms to do Le Bail fit. The introduction can be found in the wiki page of LeBailFit v1.

Usage

Example - load a Fullprof .prf file:

LoadFullprofFile(Filename=r'LaB6_1bank3_C.prf',
    PeakParameterWorkspace='LaB6_InfoTable',
    OutputWorkspace='PG3_LaB6_Bank3')

infotablews = mtd["LaB6_InfoTable"]
dataws = mtd["PG3_LaB6_Bank3"]

print "LaB6: A = B = C = %.5f, Alpha = Beta = Gamma = %.5f" % (infotablews.cell(0, 1), infotablews.cell(5, 1))

maxy = max(dataws.readY(1))
print "Maximum peak value (calculated) = %.5f" % (maxy)

Output:

Data set counter =  5431
Data Size =  5431
LaB6: A = B = C = 4.15689, Alpha = Beta = Gamma = 90.00000
Maximum peak value (calculated) = 13.38550

Example - load a Fullprof .irf file:

LoadFullprofFile(Filename=r'LB4854b3.hkl',
    PeakParameterWorkspace='LaB6_Ref_Table',
    OutputWorkspace='Fake')

fakedataws = mtd["Fake"]
reftablews = mtd["LaB6_Ref_Table"]

print "Reflection table imported %d peaks.  Faked data workspace contains %d data points." % (
    reftablews.rowCount(), len(fakedataws.readX(0)))

index = 0
print "Peak %d of (%d, %d, %d): Alpha = %.5f, Beta = %.5f, FWHM = %.5f" % (index, reftablews.cell(index, 0),
    reftablews.cell(index, 1), reftablews.cell(index, 2), reftablews.cell(index, 3), reftablews.cell(index, 4), reftablews.cell(index, 7))
index = 75
print "Peak %d of (%d, %d, %d): Alpha = %.5f, Beta = %.5f, FWHM = %.5f" % (index, reftablews.cell(index, 0),
    reftablews.cell(index, 1), reftablews.cell(index, 2), reftablews.cell(index, 3), reftablews.cell(index, 4), reftablews.cell(index, 7))

Output:

Reflection table imported 76 peaks.  Faked data workspace contains 1 data points.
Peak 0 of (1, 1, 0): Alpha = 0.01963, Beta = 0.01545, FWHM = 289.07450
Peak 75 of (9, 3, 0): Alpha = 0.25569, Beta = 0.13821, FWHM = 14.67480

Categories: Algorithms | Diffraction\DataHandling