Directory structure:

Directories are named according to the scheme L100_h030_m1000
where	L100 means size L=100 sites
	h030 means field h_x=0.30J
	m1000 means the number of kept states is up to 1000

Subdirectories:
SqwSum
SqwSumPF
SqwSumPFRes
Sxx
Syy
Szz


1. Sxx, Syy, and Szz contain the S^{xx}, S^{yy} and S^{zz} components of the
dynamical spin structure factor, respectively. The outSpectrum0.pgfplots
file is the raw data, and outSpectrum2.pgfplots has had the elastic peaks
removed by subtracting Lorentzians at each momentum value.
These folders also contain corresponding TeX files plotting. In some instances
(at least h_x=0, h_x=h_f) additional TeX files exist, which are used to draw
additional lines on top of the spectrum.

2. SqwSum contains outSpectrum3.pgfplots and associated TeX file.
- outSpectrum3.pgfplots is obtained by summing the outSpectrum2.pgfplots files
for S^{xx}, S^{yy}, and S^{zz} (with equal weights). It thus does not correspond
to the linear combination in the Cs2CoCl4 experiment, nor does it necessarily
respect sum rules. However, it can be appropriate for qualitative plots.
Enforcing sum rules on these files doesn't account for Bragg peak intensities,
and should only be used for an upper bound of QFI values.
- outSpectrum7.pgfplots was generated by summing the outSpectrum0.pgfplots files
for S^{xx}, S^{yy}, and S^{zz}, followed by enforcing the total moment sum rule,
and then subtracting Lorentzians to avoid elastic peaks. These spectra 
underestimate the real intensity, since the elastic features include both an
unphysical artifact due to finite size/boundary conditions, and Bragg peaks.
They may thus be used for a lower bound of QFI.

3. SqwSumPF contains spectra with experimental polarization factor applied
- outSpectrum4.pgfplots		is without artifact removal
- outSpectrum4_2.pgfplots	is with artifact removal

4. SqwSumPFRes contains spectra with experimental polarization factor and energy broadening applied.
- outSpectrum5.pgfplots		is without artifact removal
- outSpectrum5_2.pgfplots	is with artifact removal

The spectra in 3,4 were normalized in intensity relative to the experiment through the following procedure:
At zero field, remove elastic feature (which is an unphysical artifact) from DMRG data.
Scale DMRG intensities by sum rule.
Scale experimental zero-field data such that it agrees with DMRG (accounting for a 1.7 muB ordered moment)
This gives an absolute intensity scale to be used for finite field values too.




The special case of L100_h158114_m1000 has in S^{xx}, S^{yy},S^{zz} also files
outSpectrum8.pgfplots
outSpectrum9.pgfplots
obtained by summing the outSpectrum0.pgfplots files and enforcing the total
moment sum rule (outSpectrum8), then removing elastic features by subtracting
Lorentzians at each k (outSpectrum9).
