Plotting

From Horace
Jump to: navigation, search

An exhaustive list of commands to do with plotting your data - one of the primary functions of Horace!

plot

Plot a dnd or sqw object. This does not work for 0-dimensional objects (single points), or 4-dimensional objects (we couldn't think of a way of displaying 4 dimensions plus intensity!). For 1-dimensional objects a series of markers with errorbars are connected by a line. For 2-dimensional objects a 2D colourmap is displayed. For a 3-dimensional object the sliceomatic program is used, where a series of 2D slices within a box are plotted.

plot(w)

[figureHandle_, axesHandle_, plotHandle_] = plot(w)

The second line of code (with outputs) gives outputs that are Matlab handles to the figure window, the axes, and the plot respectively. These are useful if, for example, you wish to resize the axes, change the font size of labels, etc.

smooth

wsmooth=smooth(w)

wsmooth=smooth(w,[wid_x,wid_y,..],shape)

smooth allows you to smooth the data for plotting. Can optionally specify the width and smoothing function ('hat' or 'gaussian'). The default is 'hat', and 3 bins either side in all directions.

IMPORTANT NOTE You can only apply smoothing to dnd objects, not to sqw objects, since with the latter you would be destroying the very pixel information that the object is designed to hold. To convert an sqw object (win) to a dnd, simply type wout=d1d(win), if win is 1-dimensional, wout=d2d(win) if it is 2-dimensional, and so on.

Altering plot characteristics and other useful commands

Color

To change the marker and line color in a 1-dimensional Horace plot, type

acolor <the color that you want>

e.g. for a red plot:

acolor red

Line style

To change the line style in a 1-dimensional plot, type e.g.

aline --

to make a dashed line. You can also change the line thickness with something like

aline 2

To change multiple characteristics, type e.g.

aline(2,'--')

Type help aline in Matlab for a full list of options.


Marker style

You can change the style of marker in a 1-dimensional plot in a similar way to the above, e.g.

amark o %gives a circular marker
amark 6 %sets marker size to 6
amark(6,'o') %sets a circular marker with size 6

Type help amark in Matlab for a full list of options.


Axes

To change the x, y, or z limits of a plot, type e.g.

lx 0 3
ly -6 1
lz -2 2

which sets the x-axis limits to be 0 to 3, the y-axis limits to be -6 to 1, and the z-axis limits to be -2 to 2. Note the lz is used to change the color scale on a colormap plot.

You can also change the axes scales to be linear or logarithmic using

linx
logx
liny
...


Cursor

In order to get a cursor on your Horace plots, type

xycursor

You then left click the mouse on a position in a figure, and the x and y values are printed in the Matlab window. You can do this multiple times. To turn off the cursor, hit the carriage return key.

Keeping

To store a figure in your current session (i.e. so that the next plot you make opens in a new window, with the current plot preserved), type

keep_figure

If you have multiple figures open and you wish to alter one of them (e.g. by appending a line or more data to it)that has been kept using the above command, click on it and then type

make_current

Note that both of these options are also available in drop-down menus in the figures windows themselves.


1-dimensional plots

pd

Overplot command for 1-dimensional objects only (in addition to plot). Designed to mimic the p command from Genie, plotting markers, errorbars, and connecting lines. This means that if you use it for two objects in succession without closing the first plot window, it will overplot the second object on the first, provided it is also a plot of a 1-dimensional object. If there is no current figure window then it plot a new one. If you use this command and the current figure window does not correspond to a 1-dimensional object, then a new figure window will also be created.

pd(w_1d)

[figureHandle_, axesHandle_, plotHandle_] =  pd(w_1d)

pl

Overplot line between points for a 1-dimensional object. No markers or errorbars displayed.

pl(w_1d)

[figureHandle_, axesHandle_, plotHandle_] = pl(w_1d)

pm

Overplot markers at points for a 1-dimensional object. No line or errorbars displayed.

pm(w_1d)

[figureHandle_, axesHandle_, plotHandle_] = pm(w_1d)

pp

Overplot markers and errorbars for a 1-dimensional object. No lines linking points are displayed.

pp(w_1d)

[figureHandle_, axesHandle_, plotHandle_] = pp(w_1d)


pe

Overplot errorbars at points for a 1-dimensional object. No linking lines or markers are displayed.

pe(w_1d)

[figureHandle_, axesHandle_, plotHandle_] = pd(w_1d)

ph

Overplot histogram of a 1-dimensional object.

ph(w_1d);

[figureHandle_, axesHandle_, plotHandle_] = ph(w_1d)

ploc

Overplot pl, overplotting current figure regardless of type (i.e. can plot a 1d curve on top of a 2d dataset, such as when plotting a dispersion relation over a 2d Q-E slice).

ploc(w_1d);

peoc

Overplot pe, overplotting current figure regardless of type.

peoc(w_1d);

phoc

Overplot ph, overplotting current figure regardless of type.

phoc(w_1d);

pmoc

Overplot pm, overplotting current figure regardless of type.

pmoc(w_1d);


dd

Does exactly the same as the command pd, except that it does not overplot.

dh

Does exactly the same as the command ph, except that it does not overplot.

dl

Does exactly the same as the command pl, except that it does not overplot.

dm

Does exactly the same as the command pm, except that it does not overplot.

dp

Does exactly the same as the command pp, except that it does not overplot.

de

Does exactly the same as the command pe, except that it does not overplot.


2-dimensional plots

da

Area plot for a single 2-dimensional object, with colour-scale signifying intensity. It is this that is called when plot is used for a 2-dimensional object.

da(w_2d);

[figureHandle_, axesHandle_, plotHandle_] = da(w_2d)

ds

Surface plot for a sinlge 2-dimensional object, with colour scale and contour signifying intensity.

ds(w_2d);

[figureHandle_, axesHandle_, plotHandle_] = ds(w_2d)

ps

Overplot a surface plot of a d2d dataset or array of datasets

ps(w)

Optionally return figure, axes and plot handles:

[fig_handle, axes_handle, plot_handle] = ps(w,...)

spaghetti_plot

Plots data in sqw-file or sqw-object along HKL directions.

wsp = spaghetti_plot([0 0 0; 0.5 0.5 0.5; 0.5 0.5 0],sqw_file,'labels',{'\Gamma','R','M'})

3-dimensional plots

sliceomatic

Sliceomatic plot of multiple area plots, for a 3-dimensional object. This function is called by the plot routine.

sliceomatic(w_3d);

sliceomatic_overview

As sliceomatic, but the default view is from above. In effect this means you see a 2d slice which can be animated/changed by the third slider bar. Useful for e.g. following a spin wave dispersion ring/cone as a function of energy.

sliceomatic_overview(w_3d);

ps2

(Note that this routine is especially useful for making surface plots of dispersion relations, when used in conjunction with dispersion

Make a surface plot of a 2d SQW object, with the error array (or another data source) providing the intensity

ps2(w)       % Use error bars to set colour scale

ps2(w,wc)    % Signal in wc sets colour scale (sqw or d2d object with same array size as w, or a numeric array)

Differs from ps in that the signal sets the z axis, and the colouring is set by the error bars, or another object. This enable a function of three variables to be plotted (e.g. dispersion relation where the 'signal' array hold the energy and the error array hold the spectral weight).

One can optionally return figure, axes and plot handles:

[fig_handle, axes_handle, plot_handle] = ps2(w,...) 


Miscellaneous functions

meta(fig) allows you to copy the figure into a metafile. On Windows, this function puts the file in the clipboard so that it can be pasted directly into Word, Powerpoint etc.