This folder contains ASCII files for the (x,y,error) data plotted in Figs. 4, 6, 7, 8, 10, 11, and 13,
and the intensity colormaps in Figs. 2, 3, 5, 6, 7, 9, 11, 12, and 13.

The supplied MATLAB code plot_data.m plots all the figures when run in the same directory as the unzipped data files.
Note that running plot_data can take a couple of minutes as it plots in total 39 figures.

Also included is MATLAB code that uses the empirical model of the dispersion described in Appendix B:
empirical_dispersion.m     - computes the renormalized dispersion relation
empirical_intensity.m	   - computes the associated intensities as parameterized by the Z pre-factors in Eqs. (A3--A6)
empirical_model_example.m  - example script demonstrating how to use the above functions

To obtain the renormalized dispersion at a given q = (h,k,l) (rlu) with the parameters as per Table I, the calling syntax is
>>omega=empirical_dispersion(h,k,l,flag,load('best_fit_pars.txt', '-ascii'));
with the parameter flag = 1, 2 or 3 for the mode w(q), w(q-Q) or w(q+Q), respectively.

To obtain the intensity at a given q = (h,k,l) (rlu) and E with the parameters as per Table I, the calling syntax is
>>int=empirical_intensity(h,k,l,E,flag,load('best_fit_pars.txt', '-ascii'),[100 1 1 0 0.1 1 1 1]);
with the parameter flag = 1, 2 or 3 for the mode w(q), w(q-Q) or w(q+Q), respectively.

More information is obtained from
>>help empirical_dispersion
and
>>help empirical_intensity


(x,y,error) data
+++++++++++++++++++++++
Energies are given in meV and (h,k,l) in rlu.
Intensities and Errors are given in arbitrary units. 
|q| is given in 1/Angstroms and is used as the x coordinate in plots along wave vector paths, which is defined
as the distance travelled along a wave vector path in reciprocal space from the starting wave vector.

fig_4.xy	10 columns [|q| h k l Energy(LSWT) Energy(best-fit) Energy(K/2) Energy(K/2,mixed) Energy(M) Energy(M,mixed)]
fig_6.xye	6 columns [|q| h k l Energy Error_Energy]
fig_6.xy	7 columns [|q| h k l Energy(q) Energy(q-Q) Energy(q+Q)]
fig_7.xy	7 columns [|q| h k l Energy(q) Energy(q-Q) Energy(q+Q)]
fig_8a.xy	9 columns [|q| h k l Energy(q) Energy(q-Q) Energy(q+Q) Min_Energy(2-magnon) Max_Energy(2-magnon)]
fig_8b.xy	9 columns [|q| h k l Energy(q) Energy(q-Q) Energy(q+Q) Min_Energy(2-magnon) Max_Energy(2-magnon)]
fig_8c.xy	9 columns [|q| h k l Energy(q) Energy(q-Q) Energy(q+Q) Min_Energy(2-magnon) Max_Energy(2-magnon)]
fig_8d.xy	9 columns [|q| h k l Energy(q) Energy(q-Q) Energy(q+Q) Min_Energy(2-magnon) Max_Energy(2-magnon)]
fig_10a.xye	18 columns [Energy(data,low-energy) Intensity(data,low-energy) Error_Intensity(data,low-energy) ...
			Energy(data,high-energy) Intensity(data,high-energy) Error_Intensity(data,high-energy) ...
			Energy(continuum) Intensity(continuum) ...
			Energy(bkg,low-energy) Intensity(bkg,low-energy) Error_Intensity(bkg,low-energy) ...
			Energy(bkg,high-energy) Intensity(bkg,high-energy) Error_Intensity(bkg,high-energy) ...
			Energy(sim,1-magnon) Intensity(sim,1-magnon) Energy(sim,2-magnon) Intensity(sim,2-magnon)]
fig_10b.xye	21 columns [Energy(data,low-energy) Intensity(data,low-energy) Error_Intensity(data,low-energy) ...
			Energy(data,high-energy) Intensity(data,high-energy) Error_Intensity(data,high-energy) ...
			Energy(continuum) Intensity(continuum) ...
			Energy(bkg,M2) Intensity(bkg,M2) Error_Intensity(bkg,M2) ...
			Energy(bkg,Gamma,high-energy) Intensity(bkg,Gamma,high-energy) Error_Intensity(bkg,Gamma,high-energy) ...
			Energy(bkg,Gamma,low-energy) Intensity(bkg,Gamma,low-energy) Error_Intensity(bkg,Gamma,low-energy) ...
			Energy(sim,1-magnon) Intensity(sim,1-magnon) Energy(sim,2-magnon) Intensity(sim,2-magnon)]
fig_11.xy	16 columns [|q| h k l Energy(q-Q) Energy(q+Q) ...
			|q|(Gamma) h(Gamma) k(Gamma) l(Gamma) Energy(Gamma) ...
			|q|(K') h(K') k(K') l(K') Energy(K')] 
fig_13a.xy	7 columns [|q| h k l E(q) E(q-Q) E(q+Q)]
fig_13b.xy	7 columns [|q| h k l E(q) E(q-Q) E(q+Q)]


Color 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 detectors and energy bins by a convolving matrix with weights 
[0.1  0.2  0.1; 0.2  0.8  0.2; 0.1  0.2  0.1] to have a smoother color variation between neighboring pixels in 
the plotted intensity map. All the plotted intensities are produced by applying this smoothing matrix twice. 
The energy plot in Fig. 5 and the density of states in Fig. 11 do not include any smoothing.

To obtain the plots in the paper, one could use the matlab code plot_data.m to reshape data into 2D matrices, 
convert from bin centers to edges, and then use pcolor with appropriate axes, color limits and colormap.
The RGB colormap in cmap_fig2.dat was used for Figs. 2, 6 and 7, the colormap cmap_fig3.dat was used for 
Figs. 3, 12 and 13, cmap_fig9.dat was used for Fig. 9, and cmap_fig11.dat was used for Fig. 11.

Energies are given in meV and (h,k,l) in rlu.
Intensities and Errors are given in arbitrary units. 
The density of states (DoS) in Fig. 11 is in units of 1/(meV Co).
|q| is given in 1/Angstroms and is used as the x coordinate in plots along wave vector paths, which is defined
as the distance travelled along a wave vector path in reciprocal space from the starting wave vector.

fig_2.dat	7 columns [|q| h k l Energy Intensity Error_Intensity]
fig_3a.dat	3 columns [(h,h) (-k,k) Intensity]
fig_3b.dat	3 columns [(h,h) (-k,k) Intensity]
fig_3c.dat	3 columns [(h,h) (-k,k) Intensity]
fig_3d.dat	3 columns [(h,h) (-k,k) Intensity]
fig_3e.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_3f.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_3g.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_3h.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_3i.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_3j.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_3k.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_3l.dat	4 columns [(h,h) (-k,k) Intensity Error_Intensity]
fig_5a.dat	3 columns [(h,h) (-k,k) Energy]
fig_5b.dat	3 columns [(h,h) (-k,k) Energy]
fig_6a.dat	7 columns [|q| h k l Energy Intensity Error_Intensity]
fig_6b.dat	6 columns [|q| h k l Energy Intensity]
fig_7a.dat	7 columns [|q| h k l Energy Intensity Error_Intensity]
fig_7b.dat	6 columns [|q| h k l Energy Intensity]
fig_9a.dat	4 columns [(h,h,0) Energy Intensity Error_Intensity]
fig_9b.dat	3 columns [(h,h,0) Energy Intensity]
fig_11.dat	6 columns [|q| h k l Energy DoS]
fig_12a.dat	4 columns [(h,h,0) Energy Intensity Error_Intensity]
fig_12b.dat	4 columns [(h,h,0) Energy Intensity Error_Intensity]
fig_12c.dat	4 columns [(h,h,0) Energy Intensity Error_Intensity]
fig_12d.dat	4 columns [(h,h,0) Energy Intensity Error_Intensity]
fig_12e.dat	3 columns [(h,h,0) Energy Intensity]
fig_12f.dat	3 columns [(h,h,0) Energy Intensity]
fig_12g.dat	3 columns [(h,h,0) Energy Intensity]
fig_12h.dat	3 columns [(h,h,0) Energy Intensity]
fig_13a.dat	3 columns [(h,h,0) Energy Intensity]
fig_13b.dat	3 columns [(h,h,0) Energy Intensity]	
