Short Manual for Acceptance Calculations using TRACK
  ----------------------------------------------------

  First : We do not calculate the acceptance with space charge, so in "track.dat" the parameter
         "current" should be set to 0. --> "current=0.".
          
          Why ? Because the particle distribution for the acceptance calculation is not realistic, 
               it is usually very large in one phase plane and very small in the other two planes.
               Calculating SC for this kind of distribution will not be correct and will result in 
               wrong acceptance calculation. 


  Second: We set the acceptance flag "iaccep" in "track.dat" to 1 --> "iaccep=1", the old flag
         "iflag_tac" is no longer used. 


  Third : We calculate the acceptance for every plane separately, so for the acceptance in the 
         horizontal plane (x,x'), we set a large emittance while we set small emittances for the 
         vertical plane (y,y') and the longitduinal plane (Dphi,DW/W). For example:

               epsnx = 10.0,  alfax = 0.0, betax = 20.0 
               epsny = .001,  alfay = 0.0, betay = 20.0  
               epsnz = .001,  alfaz = .01, betaz = 20.0

               phmax=.01, dwwmax=.01, dWWacc=.01

         the emittance value could be adjusted as needed. For the vertical plane (y,y') we use:  

               epsnx = .001,  alfax = 0.0, betax = 20.0 
               epsny = 10.0,  alfay = 0.0, betay = 20.0  
               epsnz = .001,  alfaz = .01, betaz = 20.0

               phmax=.01, dwwmax=.01, dWWacc=.01

          For the longitudinal plane (Dphi, DW/W) we use:

               epsnx = .001,  alfax = 0.0, betax = 20.0 
               epsny = .001,  alfay = 0.0, betay = 20.0  
               epsnz = 10.0,  alfaz = .01, betaz = 20.0

               phmax=180., dwwmax=100., dWWacc=100.

          phmax, dwwmax and dwwacc could be adjusted to smaller values as needed.

  Forth : We use the output file "accepall.out" to plot the calculated acceptance. This file 
         contains the input coordinates of all particles. Those accepted will have a loss flag
         value of 0, lost or stopped particles will have a non zero loss flag. There is also 
         the files "accepted.out" and "stopped.out" with accepted and stopped particles, which 
         could also be used. I'm actually thinking of removing those files to have all
         the information in one file "accepall.out".  

  That's all I think