Line geometric aberrations :one turn computation

This operation traces npart particles, placed on ellipses with nominal emittances ${\epsilon_{x}}_{i}, {\epsilon_{x}}_{i}$ for one turn. It then fits an ellipse to the output points obtained. From this fitted ellipse it determines the average values for $\beta_{x},\ \alpha_{x},\ \beta_{y},\ \alpha_{y}, \ \nu_{x}, \ \nu_{y}, \ \epsilon_{x}, {\rm and}\ \epsilon_{y}$.

Input format 


LINEgeometric aberrations....(up to 80 characters) 

$\beta_{x},\alpha_{x},\beta_{y},\alpha_{y}$ 

$x_{co},\ x'_{co},\ y_{co},\ y'_{co},$ ener 


ncase,npart,ncoup 


nplot,nprint,mlocat,[list] 

${\epsilon_{x}}_{i} ,\ {\epsilon_{y}}_{i}$ for $i=1$ to ncase 

Parameter definitions

$\beta_{x},\alpha_{x},\beta_{y},\alpha_{y}$

$\parbox{10cm}{ input values of the twiss parameters at the entrance of the line.}$

$x_{co},\ x'_{co},\ y_{co},\ y'_{co},$ ener

$\parbox{10cm}{ coordinates and the momentum of the trajectory around which the aberrations are to be computed.}$

ncase $\parbox{10cm}{ number of cases analyzed(maximum 10).}$

npart $\parbox{10cm}{ number of particles to be traced(maximum 300).}$

ncoup $\parbox{10cm}{ not used presently, but a value must be inserted.}$

nplot

1 $\parbox{10cm}{ plot the resulting particles at the end of the job. It computes its own plotting windows.}$ -1 $\parbox{10cm}{ no plotting.}$

nprint

-2 $\parbox{10cm}{ no printing.}$

-1 $\parbox{10cm}{ printing at end of the line only.}$

0 $\parbox{10cm}{ printing after every element}$

mlocat

$\parbox{10cm}{ number of intervals in which printing is to occur; used in conjunction with list. This is not yet implemented.}$

list

$\parbox{10cm}{ intervals in which printing occurs. List may contain up to mxlist numbers (set at 40 initially)}$

${\epsilon_{x}}_{i} ,\ {\epsilon_{y}}_{i}$

$\parbox{10cm}{ ncase values for the chosen nominal emittances in x and y using the unit mm-mrad (E-06 m-rad)}$

Examples

This example is fictitious.

line geometric aberrations
0.19518E+01 -0.5862785371E-01 0.1816010291E+02  0.3831956152E-01
0 0 0 0 0
3 100 0
1 -1 0 40
38 38
75 75
150 150;