Script for making cuts and plotting

From Horace
Revision as of 12:09, 31 October 2014 by Russell Ewings (Talk | contribs) (Created page with "<pre> %Make different kinds of cuts and slices, and make basic plots of them %=============================== %First define the viewing axes that you wish to use. These do not ...")

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

%Make different kinds of cuts and slices, and make basic plots of them

%===============================
%First define the viewing axes that you wish to use. These do not need to have any particular relation to the
%spectrometer axes
proj.u=[1,1,0]; proj.v=[-1,1,0]; proj.uoffset=[0,0,0,0]; proj.type='rrr';
proj2.u=[1,0,0]; proj2.v=[0,1,0]; proj2.uoffset=[0,0,0,0]; proj2.type='rrr';
proj3.u=[1,1,1]; proj3.v=[-1,1,0]; proj3.uoffset=[0,0,0,0]; proj3.type='rrr';

%The u and v fields of proj define the first two viewing axes, and the third Q dimension is the cross-product
%of them. The uoffset is if you wish to explicitly give an offset from the origin for your cuts / slices.
%The type field (the best choice is 'rrr') is to say that all 3 Q axes should be in reciprocal lattice units.
%If you prefer inverse Angstroms then use 'aaa'

%================================

%3d volume slice and plot
my_vol=cut_sqw(sqw_file,proj,[0,0.1,8],[2,0.05,6],[-2,-1],[0,10,1000],'-nopix');%Makes a Q,Q,E volume plot
%the -nopix option ensures that the output is a "d3d", which takes less memory on your computer

plot(my_vol);%default style plot of the 3d volume - opens a new sliceomatic window


%================================

%2d slice and plot
my_slice=cut_sqw(sqw_file,proj,[0,0.1,8],[4,5],[-2,-1],[0,10,10000]);%Make as Q,E slice
%this time we did not use the nopix option, so the output is an "sqw", in which all contributing detector pixel information is retained

plot(my_slice);%default style plot of 2d slice

%Customise the slice plot:

%Make axes tight
plot(compact(my_slice));

%Default style smoothing (notice conversion to d2d first, since smoothing of sqw objects is not allowed, due to their strong connection to the raw data)
plot(smooth(d2d(my_slice)));

%Smoothing options
plot(smooth(d2d(my_slice),[2,2],'gaussian'));

%Set colour scale and other axes scales in script:
lz 0 0.5; %colour scale
ly 50 250
lx -1.5 1.5

%Reset a limit
lx

%Make a plot, and retain the figure window so that the next plot appears in a new window, and does not replace this one
plot(my_slice);
keep_figure;


%================================

%1d cut

my_cut1=cut_sqw(sqw_file,proj,[0,0.1,8],[5,6],[-0.1,0.1],[130,150]);

plot(my_cut1);%the lx and ly axes adjustment can be used here too

%Specify the colour, marker and other characteristics of the plot. Note you must set these first, then make the plot
acolor red
amark s
plot(my_cut1);

%Overplot another cut on the same set of axes (in this example just the first cut plus a constant):
acolor black
amark o
pp(my_cut1+3); %the pp command overplots a "plot"

%you can overplot a line using "pl", just markers using "pm", just errorbars using "pe".


%================================

%Taking a cut from a cut (CAN GIVE BIG SPEED SAVINGS):

%If you have a volume or slice, you can take a lower dimensional cut from it. Because the data are already in memory this 
%is much faster than going back to the data file on disk

my_cut2=cut(my_slice,[],[130,150]);
%note you do not need to specify all 4 axes, because some (in this case 2) have already been integrated
%the empty array [] specifies that you want to keep the same binning as that used in my_slice