Example Procedures
===================

Plot 2D Raman spectra
---------------------
Put histwv.m and all of the files labeled *2D*  into the working Matlab directory. Within Matlab, call the following:
Raman2D(N,bins,type,flag,nn,Jx,Jz)
2N the number of momentum points per dimension. bins the number of points on the Energy-axis (x-axis) of the resulting intensity-versus-energy plot. type is a string: either 'xx', 'xy', or 'xx-xy' for the three representative polarization combinations (see paper). flag is a boolean telling whether to estimate time for this run; typically choose 1 so that the estimate will be made. nn is the number of trials to average over, which is used to estimate error, reported at the end. Jx and Jz are the coupling constants of the lattice model with Jy=Jx assumed.
The function will create plots. The numerical output is described at the end of the script Raman2D.

Plot H-0 Raman spectra (with no magnetic field)
-----------------------------------------------
Put histwv.m, diag3D15, and all files labeled *14* into the working Matlab directory. Also unzip the contents of cell.zip into a folder called @cell in the same directory. Finally, call 
runner14_8.m
The assignment of variables is automated to make the plots in the paper. runner_H1_r and runner_H1_pi_r work similarly given the appropriate 'H1' and 'H1_pi' files in the working directory.

Plot H-0 Raman spectra with a magnetic field and resonant spectra
-----------------------------------------------------------------
Put histwv, diag3D15, and all *_9* files in the working directory with the folder cell unzipped into a folder @cell. Uncomment all lines starting in '% Ih' or '% I0'. Call 
runner_slab_max9_a1_6.m


Note: In the previous two cases the "runner" script sets a large value of N that will likely make the script take hours to run. These are chosen to make the plots smooth and accurate. Make N smaller to test the code much more quickly. In 3D N=15 may run in seconds and in 2D N=50 should be even faster.