Particles are randimized so they lay on a sphere in the phase space,
either 4D or 5D (4D+ delta momentum). The same emmittance is put in both planes.
At each observation point regular and normalized phase spaces are plotted.
The points are evenly distributed on a circle. To find its radius,
for each plane, normalized coordinates are histogrammed, and function of
the following form is fitted Y=sqrt(A^2 - B^2*x^2) (Ellipse equation).
Then, the radius R is B/A, and emittance is R^2.
STARTS ENDS
Element s[m] Element s[m]
Linac 1 00.0 1226 54.8 Starts at the end of solenoids
CT 1227 54.8 1317 83.55
DL 1318 83.55 1494 125.93
TL1 1495 125.94 1578 159.77
CR 1578 159.77 1876 244.53
Phase-Space evolution elong the line
Below is emittance evolution elong the line
Phase-Space evolution elong the line
Below is emittance evolution elong the line
Phase-Space evolution elong the line
Below is emittance evolution elong the line
tmp data