[430] | 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> |
---|