File Applications/ResLibCal/ResLibCal.m

  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: Nov. 26, 2018
  (c) E.Farhi, ILL. License: EUPL.
  See also iFunc/conv