14. 28. 79. 197. 3.438  0.		|Zp, Ap, Ztar, Atar, yp, ytar
92.9371698 .6833521 1.4747732  		|T, vT/c, et0
5.9650246 8.1543486 8.0032876 		|muB/T, tau, R
-0.8556867 1.3552316 1.			|aT, ycm, fs
0.647803  1.  1.  1.			|lampi, lamK, normpi, normK 
1  1  1  1  1  1  1  1  0  1  0  0  0	|vary T,vT,et0,muB/T, ...
1  1  1  1  0  0			|inv dist: pi-,pi+,K+,K-,pro,pbar
.15  .15  .15  .2  .15  .15		|sys. error in pi-,pi+,K+,K-,pro,pbar
0  1  1  0				|corr: pi-,pi+,K+,K-
0. 0. 0. 0.				|sys. error in pi-,pi+,K+,K- corr
1.25  0.				|Ymeas, rweak(cm)
0  1  0  2	    			|option
3  5  11  8  4  3			|n1 - 1-part dist. integration points
5  5  11  8  4  3			|n2 - 2-part dist. integration points

EXPLANATION OF INPUTS
1)  Zp (Ztar), Ap (Atar), and yp (ytar) are the charge, atomic number,
    and lab rapidity of the incoming projectile (target) nucleus.
2)  The temperature T should be given in MeV.  vT is the maximum radial
    expansion velocity, th(et0) is the maximum longitudinal exp. vel.
3)  muB is the baryon chemical potential -- actual parameter is muB/T.
    tau is the proper time of freezeout along the beam axis (in fm/c).
    R is the maximum transverse radius of the source (in fm).
4)  aT must be in the range -1 < aT < 1, it controls radial freezeout
    aT < 0 means outside to inside, aT > 0 means inside to outside
    ycm is the CM rapidity of the source (known for symmetric collisions).
    fs is a fudge factor to model incomplete strangeness chem equilibrium
      For complete chemical equilibrium, use a fixed value of fs=1.
5)  lampi and lamK are the incoherence parameters for pions and kaons.
    normpi and normK are arbitrary normalizations for pion and kaon
      one-particle distributions.  Default is to fix them both at 1.
6)  1 means allow parameter to vary, 0 means hold it fixed.
7)  1 (0) means read (do not read) inv. dist. data from the files
      pim.dat, pip.dat, kp.dat, km.dat, pro.dat, pbar.dat
8)  Percent systematic errors associated with each of the above data sets.
9)  1 (0) means read (do not read) 2-part. corr. data from the files
      pim.cor, pip.cor, kp.cor, km.cor
10) Percent systematic errors associated with each of the above data sets.
11) Ymeas identifies the measurement frame for qlong 
      (only when option(4) = +/- 2) 
    rweak is a measure of how much weakly-decaying strange baryons contribute
      to the pi- distribution.  For no contribution from these weak decays,
      set rweak=0.  To only include weak decays which occur within 5cm 
      transversally of the beam line, set rweak=5.  For 10cm, rweak=10.  
      Setting rweak>100. effectively includes pi- from all weak decays.
      rweak = -1. suppresses all resonances, considers only direct pi, K, p.
12) option(1) determines the minimization to use 
      = 0 does not minimize, calculates chi^2, errors, etc.
      = 1 does not use derivatives, is slow, but durable.
      = 2 uses derivatives, almost always the best choice.
      = 3 first uses no derivatives, then uses derivatives.
    option(2) determines units used for inv. dist. *.dat files
      = 0 MeV
      = 1 GeV
      momenta for correlation data are always assumed to be in MeV.
    option(3) determines what is meant by columns of *.dat files
       = 0 2nd column is mT-m, 3rd column in Ed^3N/dp^3
       = 1 2nd column is pT, 3rd column is d^2N/(dydpT)
    option(4) determines what is meant by 3rd column of *.cor files
       = 1 3rd column is y1-y2
       = 2 3rd column is qlong in the measurement frame defined by Ymeas
       = 3 3rd column is qlong in the LCMS frame
       negatives of these mean that the relative sign of ql and qout
       has been averaged over in the data (so should be in the model).
13) Number of 1-part dist. integration points for 
    y, x, eta, mTres, yres, Minv
    mTres and yres are mT and rapidity of resonances decaying into pi,K,p.
    Minv is the inv. mass of other 2 products in 3-body decay.
14) Analagous number of 2-part dist. integration points. 
    
My experience has shown that for AGS energies
2  4  9  7  3  2			|n1 - 1-part dist. integration points
4  4  9  7  3  2			|n2 - 2-part dist. integration points
   is quick and fairly accurate
3  5  11  8  4  3			|n1 - 1-part dist. integration points
5  5  11  8  4  3			|n2 - 2-part dist. integration points
   is a good accurate workhorse
4  6  12  9  5  4			|n1 - 1-part dist. integration points
6  6  12  9  5  4			|n2 - 2-part dist. integration points
   is slow but very accurate

while for SPS energies more eta integration points are needed, namely
2  4  14  7  3  2			|n1 - 1-part dist. integration points
4  4  14  7  3  2			|n2 - 2-part dist. integration points
   is quick and fairly accurate
3  5  18  8  4  3			|n1 - 1-part dist. integration points
5  5  18  8  4  3			|n2 - 2-part dist. integration points
   is a good accurate workhorse
4  6  22  9  5  4			|n1 - 1-part dist. integration points
6  6  22  9  5  4			|n2 - 2-part dist. integration points
   is slow but very accurate

For RHIC energies, probably even more points will be needed for the
eta integration.  One check as to whether enough integration points are
being used is to look at the proton dn/dy plotted from plotfrz.  An
insufficient number of data points can lead to unphysical bumps and 
wiggles.