Difference between revisions of "Planning a Horace scan"

From Horace
Jump to: navigation, search
Line 32: Line 32:
  
 
<pre>
 
<pre>
en=[0:10:190];%energy bins (coarse for simulation)
+
en=[0:10:190]; % energy bins (coarse for simulation)
par_file='/usr/local/mprogs/Libisis/InstrumentFiles/maps/4to1_124.par';%detector parameter file
+
par_file='/usr/local/mprogs/Libisis/InstrumentFiles/maps/4to1_124.par'; % detector parameter file
sqw_file_fake='/home/maps/maps_users/Gruenwald/fake_200meV.sqw';%name of output file
+
sqw_file_fake='/home/maps/maps_users/Gruenwald/fake_200meV.sqw';         % name of output file
  
efix=200;%incident energy
+
efix=200; % incident energy
emode=1;%select 1 for direct geometry spectrometer
+
emode=1;   % select 1 for direct geometry spectrometer
  
 
alatt=[5.7,5.7,5.7];
 
alatt=[5.7,5.7,5.7];
Line 44: Line 44:
 
u=[1,0,0]; v=[0,1,1];
 
u=[1,0,0]; v=[0,1,1];
  
psi=[-90:0];%scan range in degrees
+
psi=[-90:0]; % scan range in degrees
  
 
omega=0; dpsi=0; gl=0; gs=0;
 
omega=0; dpsi=0; gl=0; gs=0;
Line 55: Line 55:
  
 
<pre>
 
<pre>
proj.u=[1,0,0]; proj.v=[0,1,0]; proj.type='rrr';%viewing axes for plots
+
proj.u=[1,0,0]; proj.v=[0,1,0]; proj.type='rrr'; % viewing axes for plots
  
my_l=[-2:2];%loop over a set of values of L
+
my_l=[-2:2]; % loop over a set of values of L
  
 
for i=1:numel(my_l)
 
for i=1:numel(my_l)

Revision as of 16:10, 15 January 2019

Simple Scan Planner

Horace has a useful little tool which allows you to work out the coverage of reciprocal space that will be achieved for a given instrument / incident energy / sample angle range.

You start it up by typing

horace_planner

The graphical user interface shown below should open up

2019-01-15 11 59 56-Horace Planner.png]

  • You first select the par file associated with your instrument (which specifies detector positions), by browsing to it, and then clicking Load.
  • Next you enter the vector u and v as comma-separated lists. u is the Q-space direction parallel to the incident beam when the crystal rotation angle, psi, is equal to 0, and v is a second vector lying in the equatorial plane of the detectors (horizontal for most instruments).
  • Enter the lattice parameters (in Angstroms) as a comma-separated list, and the lattice angles.
  • Enter the incident energy, Ei, and the energy transfer you are interested in, then the minimum and maximum values of psi you are considering.
  • Now hit Calculate and the coverage is calculated. The axes of the three plots are the Q directions parallel to u and v, and the direction out-of-plane, respectively. The axes are labelled with inverse Angstroms, and circles are shown at integer (h,k,l) positions.
  • The three views shown in the image above correspond to looking at the Q-coverage volume from three different perspectives. On the top left we see the view from above (i.e. down an axis perpendicular to u and v), on the top right we see the view from the right-hand side (i.e. along u from the positive side), and on the bottom left we see a view from the bottom (i.e. along v from the negative side)

More Detailed Scan Planner

The scan planner shown above is quick and simple to use. But sometimes you want to inspect very carefully exactly what reciprocal space coverage you will get for a given scan. To do this you can generate a fake dataset, where the 'signal' is represented as the angle of of the contributing scan at each point. Note that this process is rather slow, so is better for checking your decisions about what scan range and incident energy from the simple scan planner detailed above.

An example script for generating fake data is given below.

en=[0:10:190];  % energy bins (coarse for simulation)
par_file='/usr/local/mprogs/Libisis/InstrumentFiles/maps/4to1_124.par';  % detector parameter file
sqw_file_fake='/home/maps/maps_users/Gruenwald/fake_200meV.sqw';         % name of output file

efix=200;  % incident energy
emode=1;   % select 1 for direct geometry spectrometer

alatt=[5.7,5.7,5.7];
angdeg=[90,90,90];

u=[1,0,0]; v=[0,1,1];

psi=[-90:0];  % scan range in degrees

omega=0; dpsi=0; gl=0; gs=0;

fake_sqw (en, par_file, sqw_file_fake, efix, emode, alatt, angdeg,...
                     u, v, psi, omega, dpsi, gl, gs);

Once you have created this fake dataset you can take cuts and slices out of it in exactly the same way as you would a real dataset. For example:

proj.u=[1,0,0]; proj.v=[0,1,0]; proj.type='rrr';  % viewing axes for plots

my_l=[-2:2];  % loop over a set of values of L

for i=1:numel(my_l)
    my_slice(i)=cut_sqw(sqw_file_fake,proj,[-3,0.03,3],[-3,0.03,3],[my_l(i)-0.1,my_l(i)+0.1],[80,90],'-nopix');
    plot(compact(my_slice(i))); keep_figure;
end

Sample plot of fake dataset

Scan figure.png