Tracking of particles

This operation tracks up to maxpar particles around the machine. Depending on the version maxpar may be set initially to 1000 or 10000. The initial values of the coordinates of the particles are lost in the process of tracking.They are replaced by the final values of the coordinates.

Input format 


TRACking of particles .....(up to 80 char) 


NPLOT NPRINT NPART NTURN 


particle data ( $x,x',y,y',l,\delta$ for all particles) 


MLOCAT  LIST  NGRAPH XMIN XMAX XPMIN XPMAX YMIN YMAX 


YPMIN YPMAX NCOL NLINE [ALMIN ALMAX DELMIN DELMAX]

Parameter definitions

NPLOT

0 $\parbox{10cm}{ plot action after every element.}$

-1 $\parbox{10cm}{ no plot action.}$

n $\parbox{10cm}{ action occurs after n turns (used in conjunction with MLOCAT and LIST) at MLOCAT locations specified by LIST elements.}$

NPRINT

0 $\parbox{10cm}{ print action after every element.}$

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

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

n $\parbox{10cm}{ action occurs after n turns (used in conjunction with MLOCAT and LIST) at MLOCAT locations specified by LIST elements.}$

NOTE: when nplot=-1 and nprint=-2 then mlocat and list do not appear.MLOCAT and LIST are the same for plot and print.

NPART $\parbox{10cm}{ number of particles traced.}$

$<$0 $\parbox{10cm}{ abs(npart) particles are added to particles already present from previous operation.}$

0 $\parbox{10cm}{ existing particles kept - none added.}$

$>$0 $\parbox{10cm}{ previously used particles deleted. NPART new particles introduced.}$

NTURN [t]$\parbox{10cm}{number of turns to be traced.}$

$x,x',y,y',l,\delta$$\parbox{8cm}{ particle coordinates for NPART particles.}$

MLOCAT $\parbox{10cm}{ indicates the number of intervals in which printing ... ...nting occurs at end of lattice only. When MLOCAT is 0, no number is in list.}$

LIST $\parbox{10cm}{ set of pairs of numbers, each of which indicates the... ...ting takes place. List may contain up to mxlist numbers (set at 40 initially)}$

NGRAPH

1 $\parbox{10cm}{ plot x,x' plane}$

2 $\parbox{10cm}{ plot y,y' plane}$

3 $\parbox{10cm}{ plot x,y plane}$

4 $\parbox{10cm}{ plot all planes}$

11,12,13,14 $\parbox{10cm}{ as above but the graphs are accumulated and the plot is printed at the end.}$

15,16 $\parbox{10cm}{ accumulates the $\delta,l$\ or the $\delta,\phi$ plo... ...uses the RF phase instead of path length differences and is always readable.}$

17 $\parbox{10cm}{ is equivalent to 14 as regards the xx',yy' and xy pl... ...$\delta,\phi$\ plot identical to that of produced by the ngraph value of 16.}$

XMIN,XMAX,XPMIN,XPMAX,YMIN,YMAX,YPMIN,YPMAX

$\parbox{10cm}{ limits for the plotting windows.}$

ALMIN,ALMAX,DELMIN,DELMAX

$\parbox{10cm}{limits for the plotting windows for the cases NGRAPH... ... for the other values of NGRAPH For NGRAPH=16 AL is to be interpreted as PHI.}$

NCOL,NLINE $\parbox{10cm}{ number of columns and lines to be used in plot matrix.}$

Examples

The first two examples are taken from demo3. The second one shows how to use the operation to track already existing particles.

The third example comes from demo13. It illustrates how to plot the longitudinal phase space using the RF phase advance and the relative momentum.

TRACKING OF PARTICLE TO MAKE TRIANGLE APPARENT
1 -1 3 500
.001 0.0 0.0003 0.0 0.0 0.00
.0015 0.0 0.0004 0.0 0.0 0.00
.0020 0.0 0.0005 0.0 0.0 0.00
0
11
-.005 0.005 -.001 0.001
-0.0006 0.0006 -.0003 0.0003
51 51,

MOVEMENT ANALYSIS TO FIND UNSTABLE CLOSED ORBIT NEAR .03 .002 0 0
1 3 3 -4 1 0 .00001 0.03 0.002 0.0 0.0 0.0 0.00 0,
tracking the 10 particles generated by the movemenet analysis just above
1 -2 0 1600
0
11
-.010 0.010 -.002 0.002
-.0006 0.0006 .0003 0.0003
51 51,

tracking
1 -2 3 1000
0 0 0 0 0.001 0
0 0 0 0 0.001 0.0005
0 0 0 0 0.2 0
0
16
-0.01 0.01 -0.01 0.01
-0.01 0.01 -0.01 0.01
51 31
0 12.56 -0.003 0.003,