ResLibCal: compute and display the neutron TAS resolution function
To start the application and open the main GUI window, use
ResLibCal
To compute directly the resolution function, sending an EXP ResLib-like
configuration structure, use:
out = ResLibCal(EXP);
To convolve an iFunc model with a 4D resolution function, use:
out = ResLibCal(model);
To use ResLibCal from the command line, use:
out = ReslibCal(command, arguments...);
where 'command' is one of:
open open a configuration file (par, cfg, res, ini, m): ResLibCal('open','file')
save save the configuration in the Preference directory (ini format)
saveas save the configuration into a specified file/format: ResLibCal('saveas','file')
export dump the main ResLibCal window into a file
exit close all active views, and save current configuration
reset re-load the default configuration
print generate an HTML document to be printed
create open the main GUI (start interface), and read last saved configuration
compute only compute the matrix (no plotting/printing)
update compute, and then update open views, or send result to the console
view2 display the 2D view (resolution projections)
view3 display the 3D view (resolution)
tas display the spectrometer geometry
default same as create, but does not read last configuration (using reset configuration)
quit same as exit, but does not save the configuration
close close the 2D, 3D and TAS view windows
resol print-out the resolution matrix a la RESCAL
bragg print-out the Bragg widths a la RESCAL
list print-out the RESCAL parameter list
config return the current configuration (ResLib EXP)
hkle return the current HKLE location. Set it back with ResLibCal(hkle{:});
silent use silent computation (no plot/display) for further arguments
<PAR>=<VALUE> sets a parameter value, e.g. 'DM=3.355'
To compute the resolution at a given HKLW location, using the current settings
resolution = ResLibCal(QH,QK,QL,W)
where QH,QK,QL,W can be vectors, or empty to use current settings.
To only compute the resolution function, use:
ResLibCal('close'); out=ResLibCal('compute');
which will close the 2D,3D and TAS views, then compute the resolution
function, returned in out.resolution
The application contains a main interface with:
* Menu, Method, Scan and Instrument parameters (main)
* Resolution function plot (2D)
* Resolution function plot (3D)
* Instrument view
when changing any value in the main GUI:
* Method and Scan parameters, Instrument parameters
any opened view is updated after a re-computation of the resolution.
The 2D and 3D views can be closed without ending the application.
When the main window is closed, or Exit is selected all views are closed
Convolution in 4D for TAS ----------------------------------------------------
The 4D convolution syntax allow to simulate a TAS scan from a parametrised
dispersion. In the following example, we simulate a scan through a cubic crystal
dispersion, and show the ideal S(q,w) as well as the measured, broadened
measurement. Then, the simulated scan is inserted in a representation of the
S(q,w), to visualise the scan trajectory and simulated signal.
To use this tool, you need to input a 2D or 4D dispersion model (iFunc).
The dispersion is either a 2D model S(|q|,w) or a 4D model S(qh,qk,ql,w).
The axes of the dispersion are in the lattice reciprocal space, in r.l.u.
s=sqw_cubic_monoatomic; % create a 4D S(q,w) for a cubic pure material
t=ResLibCal(s); % convolute it with a TAS resolution, and open ResLibCal.
w=linspace(0.01,20,50); qh=0.3*ones(size(w)); qk=0*qh; ql=qk; % a scan
signal1=iData(t, [], qh,qk,ql,w);
signal0=iData(s, [], qh,qk,ql,w);
figure; plot(squeeze([signal1 signal0*100])); % plot the dispersion and simulated measurement
% now plot the 4D dispersion with the scan in it, for fun
qh=linspace(0.01,.5,50);qk=qh; ql=qh; w=linspace(0.01,10,51); % a 4D grid
f=iData(s,[],qh,qk,ql,w); % evaluate the model on the 4D grid
figure; surf(log(f(:,:,1,:)),'half'); hold on; % plot dispersion, and scan
scatter3(log(signal1(:,:,1,:)),'filled');
input: any combination of:
command: string among those listed above, which can be followed by any other
allowed parameter.
qh,qk,ql,w: 4 vectors or 4D matrices which delimit a region in the reciprocal
space where the TAS resolution function should be computed.
model: an iFunc model which is to be convolved with the TAS response.
EXP: a structure holding a ResLibCal, ResCal or ResLib configuration.
output:
a ResLibCal configuration with e.g. 'resolution' field, or a 4D convoluted model.
Version: Aug. 22, 2017
(c) E.Farhi, ILL. License: EUPL.
See also iFunc/conv