This folder contains ASCII files for the (x,y,error) data plotted in Figs. 6A, 6B, 6C, 6D, 6E, 6F, 7A, 7B, 9B, 9C, 9D, 10, 11A, 11B, 13C, 13F
and (x,y) fits/calculations plotted in Figs. 5B, 6A, 6B, 6C, 6D, 6E, 6F, 7A, 7B, 9F, 9G, 9H, 9J, 9K, 9L, 11A, 11B, 12, 13A, 13B, 13D, 13E, 13G, 13H
and the intensity colour maps in Figs. 3A, 3B, 3C, 3D, 3E, 3F, 3G, 3H, 3A, 4B, 4C, 4D, 4E, 4F, 4G, 4H, 5A, 5B, 
	9A, 9B, 9C, 9D, 9E, 9F, 9G, 9H, 9I, 9J, 9K, 9L, 12, 13A, 13B, 13D, 13E, 13G, 13H

The folder also contains (x,y,e) data for the 2m1 continuum boundary in 2m1.xye, which is used in Figs. 6A, 6B, 7A and 7B,
and three RGB colour maps in colmap1.dat, colmap2.dat and colmap3.dat.
colmap1.dat is used for the single orientation slices in Figs. 3 and 4, colmap2.dat is used for the slices through
Horace scans in Figs. 5, 9 and 13, and colmap3.dat is used in Fig. 12.

The supplied MATLAB code plot_data.m plots 50 figures when run in the same directory as the unzipped data files.
This takes about a minute to run. Each figure panel in the paper is displayed in a separate matlab figure window, 
with the paper figure panel number, i.e. 3A to 13H, displayed in the figure window top banner. 
One can also run separate sections of the code for each individual figure.

(x,y) data : fits and lines/calculations
++++++++++++++++++++++++++++++++++++++++++++++
Energies are given in meV for modes ordered in increasing energy, and (h,k,l) in rlu.
Intensities are given in arbitrary units. 
Magnetic fields are given in T.
fig12.xy was extracted by eye from Fig. 2c of Ref. [33] so no error bars are provided.
In fig11A.xy and fig11B.xy, NaN values are used where the modes are not defined.

fig5B.xy			 6 columns  [l Energy(m1 primary) Energy(m1 shadow) Energy(m2 primary) Energy(m2 shadow) Energy(2m1 boundary)]
fig6A.xy		 	 3 columns  [Field Energy(eye guide) Intensity(eye guide)]
fig6B.xy		 	 3 columns  [Field Energy(eye guide) Intensity(eye guide)]
fig6C.xy		 	 4 columns  [Energy Intensity(fit to all peaks) Intensity(fit, m3 peak only) Intensity(fit, m4 peak only)]
						First row gives background intensity. Second row gives position of 2m1 continuum boundary at l=0.5.
fig6D.xy		 	 4 columns  [Energy Intensity(fit to all peaks) Intensity(fit, m3 peak only) Intensity(fit, m4 peak only)]
						First row gives background intensity. Second row gives position of 2m1 continuum boundary at l=0.5.
fig6E.xy		 	 3 columns  [Energy Intensity(fit to all peaks) Intensity(fit, m3 peak only)]
						First row gives background intensity. Second row gives position of 2m1 continuum boundary at l=0.5.
fig6F.xy		 	 3 columns  [Energy Intensity(fit to all peaks) Intensity(fit, m3 peak only)]
						First row gives background intensity. Second row gives position of 2m1 continuum boundary at l=0.5.						
fig7A.xy		 	 5 columns  [Field Energy(m1) Energy(m2) Energy(m3) Energy(m4)]
fig7B.xy		 	 8 columns  [Field Energy(m1) Energy(m2) Energy(m3) Energy(m4) Energy(m5) Energy(m6) Energy(m7)]
fig9F-H.xy		     5 columns  [l Energy(m1 primary) Energy(m1 shadow) Energy(m2 primary) Energy(m2 shadow)]
fig9J-L.xy		 	 3 columns  [l Energy(\omega_1) Energy(\omega_2)]
fig11A.xy		 	 8 columns  [Field Energy(m1) Energy(m2) Energy(m3) Energy(m4) Energy(m5) Energy(m6) Energy(m7)]
fig11B.xy		 	 4 columns  [Field Energy(m1) Energy(m2) Energy(m3)]
fig12.xy		 	 2 columns  [Field Energy]		
fig13A.xy		 	 3 columns  [h Energy(m1 primary) Energy(m1 shadow)]
fig13B.xy		 	 3 columns  [k Energy(m1 primary) Energy(m1 shadow)]
fig13D.xy			 3 columns  [h Energy(\omega_1) Energy(\omega_3)]
fig13E.xy			 3 columns  [k Energy(\omega_1) Energy(\omega_3)]
fig13G.xy			 2 columns  [h Energy(\omega_2)]
fig13H.xy			 2 columns  [k Energy(\omega_2)]

(x,y,error) data : data points with errors
++++++++++++++++++++++++++++++++++++++++++++++
Energies are given in meV for modes ordered in increasing energy, and (h,k,l) in rlu.
Intensities are given in arbitrary units. 
Magnetic fields are given in T.
In fig7A.xye, fig7B.xye, fig11A.xye and fig11B.xye, NaN values are used when the mode cannot be resolved at the given field.

fig6A.xye			6 columns [Field Energy Error(Energy) FWHM Intensity Error(Intensity)]
fig6B.xye		 	6 columns [Field Energy Error(Energy) FWHM Intensity Error(Intensity)]
fig6C.xye		 	3 columns [Energy Intensity Error(Intensity)]
fig6D.xye		 	3 columns [Energy Intenstiy Error(Intensity)]
fig6E.xye		 	3 columns [Energy Intensity Error(Intensity)]
fig6F.xye		 	3 columns [Energy Intenstiy Error(Intensity)]
fig7A.xye		   15 columns [Field Energy(m1) Error(Energy(m1)) Energy(m2) Error(Energy(m2)) Energy(m3) Error(Energy(m3)) ...
								Energy(m4) Error(Energy(m4)) Energy(m5) Error(Energy(m5)) Energy(m6) Error(Energy(m6)) Energy(m7) Error(Energy(m7))]
fig7B.xye		   15 columns [Field Energy(m1) Error(Energy(m1)) Energy(m2) Error(Energy(m2)) Energy(m3) Error(Energy(m3)) ...
								Energy(m4) Error(Energy(m4)) Energy(m5) Error(Energy(m5)) Energy(m6) Error(Energy(m6)) Energy(m7) Error(Energy(m7))]
fig10.xye  		 	5 columns [l Energy(empirical) Error(Energy(empirical)) Energy(ED) ModeNumber]
						Here, ModeNumber == 1 (8T//a m1), 2 (8T//a m2), 3 (zero field m1), 4 (zero field m2), 
						5 (zero field kinetic bound state) or 6 (9T//b magnon)
fig11A.xye		   17 columns [Field Energy(m1) Error(Energy(m1)) Energy(m2) Error(Energy(m2)) Energy(m3) Error(Energy(m3)) ...
								Energy(m4) Error(Energy(m4)) Energy(m5) Error(Energy(m5)) Energy(m6) Error(Energy(m6)) ...
								 Energy(m7) Error(Energy(m7)) Energy(2m1 boundary) Error(Energy(2m1 boundary))]
fig11B.xye		    7 columns [Field Energy(m1) Error(Energy(m1)) Energy(m2) Error(Energy(m2)) Energy(m3) Error(Energy(m3))]
fig13C.xye  		7 columns [h k l Energy(Observed) Error(Energy Observed) ModeNumber Energy(Calculated)]
						Here ModeNumber == 1 or 2 corresponding to E_- and E_+, respectively, as per (A2)
fig13F.xye  		7 columns [h k l Energy(Observed) Error(Energy Observed) ModeNumber Energy(Calculated)]
						Here ModeNumber == 1, 2 correspond to \omega_-(Q) and \omega_+(Q) respectively as per (A8), 3 is \omega_-(Q+a*)

Colour map figures
++++++++++++++++++++++++++++++++++++++++++++++
The first row of the file gives the dimensions of the plot arrays in the form [nRows, nCols].

The intensities are smoothed between neighboring bins by convolving a matrix with relative weights 
A=[0.1  0.2  0.1; 0.2  0.8  0.2; 0.1  0.2  0.1] to have a smoother colour variation between neighboring 
bins in the plotted intensity map. The plotted intensities with experimental data in Fig. 9A are 
produced by applying this smoothing matrix twice. For all other experimental data panels, 
the smoothing matrix is applied once. In Figs. 9B, 9C and 9D, and 12the calculation plots are unsmoothed.
In Figs. 5B, 9F, 9G, 9H, 9J, 9K, 9L, this smoothing matrix has been applied once to better mimic the experimental data.

To obtain the plots in the paper, use the matlab code plot_data.m, running this will plot all 50 figure panels. 
It is also possible to run separate sections of the code for each individual figure. The code uses an internal function to 
reshape the data into 2D matrices, convert from bin centers to edges, and then use pcolor with appropriate axes, 
colour limits and colour map.

plot_data.m contains the subroutine 'plot_color_fig' which plots the colour data specified in the first argument,
the second argument specifies whether the data contains a column for errors (i.e. for fig_1a.dat), and the third argument
specifies whether smoothed intensities (default, flag=1), raw intensities (flag=0), or errors (flag=-1) should be plotted.
The fourth argument is optionally used to widen bins at the edges of data sets to avoid gaps in plots in Fig. 9A, where the data
is made up of the data in fig9A_1.dat and fig9A_2.dat overlaid on each other.
For this panel, data files _1 and _2 correspond to the angle between the a-axis and the incident beam direction
psi = 155 and 90 respectively. The data in Fig. 9A was previously reported in Fig 3B of ref [5] R. Coldea et al, Science 327, 177 (2010).

Energies are given in meV and (h,k,l) in rlu.
Intensities and Errors are given in arbitrary units. 
Field is given in T.

fig3A.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig3B.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig3C.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig3D.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig3E.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig3F.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig3G.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig3H.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4A.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4B.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4C.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4D.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4E.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4F.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4G.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig4H.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig5A.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig5B.dat			3 columns [l Energy Intensity]
fig9A_1.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig9A_2.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig9B.dat			3 columns [l Energy Intensity]
fig9C.dat			3 columns [l Energy Intensity]
fig9D.dat			3 columns [l Energy Intensity]
fig9E_lower.dat		5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig9E_upper.dat		5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig9F_lower.dat		3 columns [l Energy Intensity]
fig9F_upper.dat		3 columns [l Energy Intensity]
fig9G_lower.dat		3 columns [l Energy Intensity]
fig9G_upper.dat		3 columns [l Energy Intensity]
fig9H_lower.dat		3 columns [l Energy Intensity]
fig9H_upper.dat		3 columns [l Energy Intensity]
fig9I_lower.dat		5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig9I_upper.dat		5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig9J_lower.dat		3 columns [l Energy Intensity]
fig9J_upper.dat		3 columns [l Energy Intensity]
fig9K_lower.dat		3 columns [l Energy Intensity]
fig9K_upper.dat		3 columns [l Energy Intensity]
fig9L_lower.dat		3 columns [l Energy Intensity]
fig9L_upper.dat		3 columns [l Energy Intensity]
fig12.dat			3 columns [Field Energy Intensity]
fig13A.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig13B.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig13D.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig13E.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig13G.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]
fig13H.dat			5 columns [l Energy Intensity Error(Intensity) Intenstity(Smoothed)]

