1 | <head> |
---|
2 | <title>Global Reference System</title> |
---|
3 | <!-- Changed by: Chris ISELIN, 17-Jul-1997 --> |
---|
4 | <!-- Changed by: Hans Grote, 10-Jun-2002 --> |
---|
5 | </head> |
---|
6 | |
---|
7 | <body bgcolor="#ffffff"> |
---|
8 | |
---|
9 | <center> |
---|
10 | EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH |
---|
11 | <IMG SRC="http://cern.ch/madx/icons/mx7_25.gif" align=right> |
---|
12 | <h2>Global Reference System</h2> |
---|
13 | </center> |
---|
14 | |
---|
15 | The <a href="#global">global reference orbit</a> of the accelerator |
---|
16 | is uniquely defined by the sequence of physical elements. |
---|
17 | The local reference system (<i>x</i>, <i>y</i>, <i>s</i>) |
---|
18 | may thus be referred to a global Cartesian coordinate system |
---|
19 | (<i>X</i>, <i>Y</i>, <i>Z</i>) (see <a href="#global">Figure 1</a>). |
---|
20 | The positions between beam elements are numbered 0,...,i,...n. |
---|
21 | The local reference system |
---|
22 | (<i>x<sub>i</sub>, y<sub>i</sub>, s<sub>i</sub></i>) |
---|
23 | at position <i>i</i>, |
---|
24 | i.e. the displacement and direction of the reference orbit |
---|
25 | with respect to the system (<i>X</i>, <i>Y</i>, <i>Z</i>) |
---|
26 | are defined by three displacements |
---|
27 | (<i>X<sub>i</sub></i>, <i>Y<sub>i</sub></i>, <i>Z<sub>i</sub></i>) |
---|
28 | and three angles |
---|
29 | (<i>Theta<sub>i</sub></i>, <i>Phi<sub>i</sub></i>, <i>Psi<sub>i</sub></i>) |
---|
30 | The above quantities are defined more precisely as follows: |
---|
31 | <ul> |
---|
32 | <li>X: |
---|
33 | Displacement of the local origin in <i>X</i>-direction. |
---|
34 | <li>Y: |
---|
35 | Displacement of the local origin in <i>Y</i>-direction. |
---|
36 | <li>Z: |
---|
37 | Displacement of the local origin in <i>Z</i>-direction. |
---|
38 | <li><a name=theta>THETA</a>: |
---|
39 | Angle of rotation (azimuth) about the global <i>Y</i>-axis, |
---|
40 | between the global <i>Z</i>-axis and the projection |
---|
41 | of the reference orbit onto the (<i>Z</i>, <i>X</i>)-plane. |
---|
42 | A positive angle THETA forms a right-hand screw with the <i>Y</i>-axis. |
---|
43 | <li><a name=phi>PHI</a>: |
---|
44 | Elevation angle, i.e. the angle between the reference orbit and its projection |
---|
45 | onto the (<i>Z</i>, <i>X</i>)-plane. |
---|
46 | A positive angle PHI correspond to increasing <i>Y</i>. |
---|
47 | If only horizontal bends are present, |
---|
48 | the reference orbit remains in the (<i>Z</i>, <i>X</i>)-plane. |
---|
49 | In this case PHI is always zero. |
---|
50 | <li><a name=psi>PSI</a>: |
---|
51 | Roll angle about the local <i>s</i>-axis, |
---|
52 | i.e. the angle between the intersection (<i>x</i>, <i>y</i>)- and |
---|
53 | (<i>Z</i>, <i>X</i>)-planes and the local <i>x</i>-axis. |
---|
54 | A positive angle PSI forms a right-hand screw with the <i>s</i>-axis. |
---|
55 | </ul> |
---|
56 | The angles (THETA, PHI, PSI) are <b>not</b> the Euler angles. |
---|
57 | The reference orbit starts at the origin and points by default |
---|
58 | in the direction of the positive <i>Z</i>-axis. |
---|
59 | The initial local axes (<i>x</i>, <i>y</i>, <i>s</i>) |
---|
60 | coincide with the global axes (<i>X</i>, <i>Y</i>, <i>Z</i>) in this order. |
---|
61 | The six quantities (X<sub>0</sub>, Y<sub>0</sub>, Z<sub>0</sub>, |
---|
62 | THETA<sub>0</sub>, PHI<sub>0</sub>, PSI<sub>0</sub>) |
---|
63 | thus all have zero initial values by default. |
---|
64 | The program user may however specify different initial conditions. |
---|
65 | <p> |
---|
66 | Internally the displacement is described by a vector <i>V</i> |
---|
67 | and the orientation by a unitary matrix <i>W</i>. |
---|
68 | The column vectors of <i>W</i> are the unit vectors spanning |
---|
69 | the local coordinate axes in the order (<i>x, y, s</i>). |
---|
70 | <i>V</i> and <i>W</i> have the values: |
---|
71 | <p> |
---|
72 | <img align=bottom src="../equations/VW.gif"> |
---|
73 | <p> |
---|
74 | where |
---|
75 | <p> |
---|
76 | <img align=bottom src="../equations/PhiThetaPsi.gif"> |
---|
77 | <p> |
---|
78 | The reference orbit should be closed and it should not be twisted. |
---|
79 | This means that the displacement of the local reference system |
---|
80 | must be periodic with the revolution frequency of the accelerator, |
---|
81 | while the position angles must be periodic modulo(2 pi) |
---|
82 | with the revolution frequency. |
---|
83 | If PSI is not periodic module(2 pi), coupling effects are introduced. |
---|
84 | When advancing through a beam element, |
---|
85 | MAD computes <i>V<sub>i</sub></i> and <i>W<sub>i</sub></i> |
---|
86 | by the recurrence relations |
---|
87 | <p> |
---|
88 | <i>V<sub>i</sub> = W<sub>i-1</sub>R<sub>i</sub> + V<sub>i-1</sub></i>, |
---|
89 | <i>W<sub>i</sub> = w<sub>i-1</sub>S<sub>i</sub></i>. |
---|
90 | <p> |
---|
91 | The vector <i>R<sub>i</sub></i> is the displacement and the matrix |
---|
92 | <i>S<sub>i</sub></i> is the rotation of the local reference system |
---|
93 | at the exit of the element <i>i</i> with respect to the entrance |
---|
94 | of the same element. |
---|
95 | The values of <i>R<sub>i</sub></i> and <i>S<sub>i</sub></i> |
---|
96 | are listed in the: |
---|
97 | |
---|
98 | <a href="local_system.html#straight">straight reference system</a> |
---|
99 | for each physical element type. |
---|
100 | <center> |
---|
101 | <a name=global><img align=bottom src="../figures/global.gif"></a> |
---|
102 | <p> |
---|
103 | <b>Figure 1:</b> Global Reference System |
---|
104 | </center> |
---|
105 | <p> |
---|
106 | <address> |
---|
107 | <a href="http://www.cern.ch/Hans.Grote/hansg_sign.html">hansg</a>, |
---|
108 | January 24, 1997 |
---|
109 | </address> |
---|
110 | |
---|
111 | </body> |
---|