... | ... | @@ -55,8 +55,8 @@ Common Attributes for all Elements |
|
|
The following attributes are allowed on all elements:
|
|
|
|
|
|
TYPE::
|
|
|
A string value see Section [astring]. It specifies an ``engineering
|
|
|
type'' and can be used for element selection.
|
|
|
A string value see Section [astring]. It specifies an "engineering
|
|
|
type" and can be used for element selection.
|
|
|
APERTURE::
|
|
|
A string value see Section [astring] which describes the element
|
|
|
aperture. All but the last attribute of the aperture have units of
|
... | ... | @@ -239,7 +239,7 @@ DESIGNENERGY:: |
|
|
FMAPFN::
|
|
|
Name of the field map for the magnet. Currently maps of type
|
|
|
`1DProfile1` can be used see Section [1DProfile1]. The default option
|
|
|
for this attribute is `FMAPN` = ```1DPROFILE1-DEFAULT`''
|
|
|
for this attribute is `FMAPN` = `1DPROFILE1-DEFAULT`
|
|
|
see Section [benddefaultfieldmapopalt]. The field map is used to
|
|
|
describe the fringe fields of the magnet see Section [1DProfile1].
|
|
|
|
... | ... | @@ -286,7 +286,7 @@ DESIGNENERGY:: |
|
|
FMAPFN::
|
|
|
Name of the field map for the magnet. Currently maps of type
|
|
|
`1DProfile1` can be used see Section [1DProfile1]. The default option
|
|
|
for this attribute is `FMAPN` = ```1DPROFILE1-DEFAULT`''
|
|
|
for this attribute is `FMAPN` = `1DPROFILE1-DEFAULT`
|
|
|
see Section [benddefaultfieldmapopalt]. The field map is used to
|
|
|
describe the fringe fields of the magnet see Section [1DProfile1].
|
|
|
|
... | ... | @@ -305,7 +305,7 @@ edge angle. |
|
|
|
|
|
L::
|
|
|
Chord length of the bend reference arc in meters (see Figure [sbend]),
|
|
|
given by: latexmath:[\[L = 2 R \sin\left(\frac{\alpha}{2}\right)\]]
|
|
|
given by: latexmath:[L = 2 R \sin\left(\frac{\alpha}{2}\right)]
|
|
|
GAP::
|
|
|
Full vertical gap of the magnet (meters).
|
|
|
HAPERT::
|
... | ... | @@ -343,7 +343,7 @@ DESIGNENERGY:: |
|
|
FMAPFN::
|
|
|
Name of the field map for the magnet. Currently maps of type
|
|
|
`1DProfile1` can be used see Section [1DProfile1]. The default option
|
|
|
for this attribute is `FMAPN` = ```1DPROFILE1-DEFAULT`''
|
|
|
for this attribute is `FMAPN` = `1DPROFILE1-DEFAULT`
|
|
|
see Section [benddefaultfieldmapopalt]. The field map is used to
|
|
|
describe the fringe fields of the magnet see Section [1DProfile1].
|
|
|
|
... | ... | @@ -374,7 +374,7 @@ will be bent through a different angle if its mean kinetic energy |
|
|
doesn’t correspond to the `DESIGNENERGY`.
|
|
|
4. Internally the bend geometry is setup based on the ideal reference
|
|
|
trajectory, as shown in Figure [rbend,sbend].
|
|
|
5. If the default field map, ```1DPROFILE-DEFAULT`''
|
|
|
5. If the default field map, `1DPROFILE-DEFAULT`
|
|
|
see Section [benddefaultfieldmapopalt], is used, the fringe fields will
|
|
|
be adjusted so that the effective length of the real, soft edge magnet
|
|
|
matches the ideal, hard edge bend that is defined by the reference
|
... | ... | @@ -400,7 +400,7 @@ positive direction 30 degrees (towards the negative x axis as if |
|
|
Figure [rbend]). It has a design energy of 10MeV, a length of 0.5m, a
|
|
|
vertical gap of 2cm and a 0latexmath:[^{\circ}] entrance edge angle.
|
|
|
(Therefore the exit edge angle is 30latexmath:[^{\circ}].) We are
|
|
|
using the default, internal field map ``1DPROFILE1-DEFAULT''
|
|
|
using the default, internal field map "1DPROFILE1-DEFAULT"
|
|
|
see Section [benddefaultfieldmapopalt] which describes the magnet fringe
|
|
|
fields see Section [1DProfile1]. When _OPAL_ is run, you will get the
|
|
|
following output (assuming an electron beam) for this `RBEND`
|
... | ... | @@ -510,7 +510,8 @@ The output is effectively identical, to within a small numerical error. |
|
|
Now, let us modify this first example so that we bend instead in the
|
|
|
negative x direction. There are several ways to do this:
|
|
|
|
|
|
1. ....
|
|
|
1.
|
|
|
....
|
|
|
Bend: RBend, ANGLE = -30.0 * Pi / 180.0,
|
|
|
FMAPFN = "1DPROFILE1-DEFAULT",
|
|
|
ELEMEDGE = 0.25,
|
... | ... | @@ -518,16 +519,18 @@ Bend: RBend, ANGLE = -30.0 * Pi / 180.0, |
|
|
L = 0.5,
|
|
|
GAP = 0.02;
|
|
|
....
|
|
|
2. ....
|
|
|
2.
|
|
|
....
|
|
|
Bend: RBend, ANGLE = 30.0 * Pi / 180.0,
|
|
|
FMAPFN = "1DPROFILE1-DEFAULT",
|
|
|
ELEMEDGE = 0.25,
|
|
|
DESIGNENERGY = 10.0E6,
|
|
|
L = 0.5,
|
|
|
GAP = 0.02,
|
|
|
ROTATION = Pi;
|
|
|
ROTATION = Pi;
|
|
|
....
|
|
|
3.
|
|
|
....
|
|
|
3. ....
|
|
|
Bend: RBend, K0 = 0.0350195,
|
|
|
FMAPFN = "1DPROFILE1-DEFAULT",
|
|
|
ELEMEDGE = 0.25,
|
... | ... | @@ -535,14 +538,15 @@ Bend: RBend, K0 = 0.0350195, |
|
|
L = 0.5,
|
|
|
GAP = 0.02;
|
|
|
....
|
|
|
4. ....
|
|
|
4.
|
|
|
....
|
|
|
Bend: RBend, K0 = -0.0350195,
|
|
|
FMAPFN = "1DPROFILE1-DEFAULT",
|
|
|
ELEMEDGE = 0.25,
|
|
|
DESIGNENERGY = 10.0E6,
|
|
|
L = 0.5,
|
|
|
GAP = 0.02,
|
|
|
ROTATION = Pi;
|
|
|
ROTATION = Pi;
|
|
|
....
|
|
|
|
|
|
In each of these cases, we get the following output for the bend (to
|
... | ... | @@ -621,7 +625,7 @@ Bend4: RBend, ANGLE = 20.0 * Pi / 180.0, |
|
|
....
|
|
|
|
|
|
Up to now, we have only given examples of `RBEND` definitions. If we
|
|
|
replaced ``RBend'' in the above examples with ``SBend'', we would still
|
|
|
replaced "RBend" in the above examples with "SBend", we would still
|
|
|
be defining valid _OPAL-t_ bends. In fact, by adjusting the `L`
|
|
|
attribute according to Section [RBend,SBend], and by adding the
|
|
|
appropriate definitions of the `E2` attribute, we could even get
|
... | ... | @@ -908,8 +912,11 @@ To model a particular bend magnet, one must fit the field profile along |
|
|
the mid-plane of the magnet perpendicular to its face for the entrance
|
|
|
and exit fringe fields to the Enge function:
|
|
|
|
|
|
latexmath:[\[\label{eq:enge_func}
|
|
|
F(z) = \frac{1}{1 + e^{\sum\limits_{n=0}^{N_{order}} c_{n} (z/D)^{n}}}\]]
|
|
|
[#eq:enge_func]
|
|
|
[latexmath]
|
|
|
++++
|
|
|
F(z) = \frac{1}{1 + e^{\sum\limits_{n=0}^{N_{order}} c_{n} (z/D)^{n}}}
|
|
|
++++
|
|
|
|
|
|
where latexmath:[D] is the full gap of the magnet,
|
|
|
latexmath:[N_{order}] is the Enge function order and latexmath:[z]
|
... | ... | @@ -926,17 +933,20 @@ element as illustrated in Figure [rbend,sbend]. (As already stated, any |
|
|
bend can be described by a rotated positive bend.) _OPAL-t_ then has the
|
|
|
following information:
|
|
|
|
|
|
latexmath:[\[\begin{aligned}
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
B_0 &= \text{Field amplitude (T)} \\
|
|
|
R &= \text{Bend radius (m)} \\
|
|
|
n &= -\frac{R}{B_{y}}\frac{\partial B_y}{\partial x} \text{ (Field index, set using the parameter \texttt{K1})} \\
|
|
|
n &= -\frac{R}{B_{y}}\frac{\partial B_y}{\partial x} \text{ (Field index, set using the parameter } \texttt{K1} \text{)} \\
|
|
|
F(z) &= \left\{
|
|
|
\begin{array}{lll}
|
|
|
& F_{entrance}(z_{entrance}) \\
|
|
|
& F_{center}(z_{center}) = 1 \\
|
|
|
& F_{exit}(z_{exit})
|
|
|
\end{array}
|
|
|
\right.\end{aligned}\]]
|
|
|
\right.\end{aligned}
|
|
|
++++
|
|
|
|
|
|
Here, we have defined an overall Enge function, latexmath:[F(z)], with
|
|
|
three parts: entrance, center and exit. The exit and entrance fringe
|
... | ... | @@ -944,7 +954,9 @@ field regions have the form of Equation [enge_func] with parameters |
|
|
defined by the `1DProfile1` field map file given by the element
|
|
|
parameter `FMAPFN`. Defining the coordinates:
|
|
|
|
|
|
latexmath:[\[\begin{aligned}
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
y &\equiv \text{Vertical distance from magnet mid-plane} \\
|
|
|
\Delta_x &\equiv \text{Perpendicular distance to reference trajectory (see Figure~\ref{rbend,sbend})} \\
|
|
|
\Delta_z &\equiv \left\{
|
... | ... | @@ -953,41 +965,59 @@ y &\equiv \text{Vertical distance from magnet mid-plane} \\ |
|
|
& \text{Not defined, Enge function is always 1 in this region.} \\
|
|
|
& \text{Distance from exit Enge function origin perpendicular to magnet exit face.}
|
|
|
\end{array}
|
|
|
\right.\end{aligned}\]]
|
|
|
\right.\end{aligned}
|
|
|
++++
|
|
|
|
|
|
using the conditions
|
|
|
|
|
|
latexmath:[\[\begin{aligned}
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
\nabla \cdot \overrightarrow{B} &= 0 \\
|
|
|
\nabla \times \overrightarrow{B} &= 0\end{aligned}\]]
|
|
|
\nabla \times \overrightarrow{B} &= 0\end{aligned}
|
|
|
++++
|
|
|
|
|
|
and making the definitions:
|
|
|
|
|
|
latexmath:[\[\begin{aligned}
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
F'(z) &\equiv \frac{\mathrm{d} F(z)}{\mathrm{d} z} \\
|
|
|
F''(z) &\equiv \frac{\mathrm{d^{2}} F(z)}{\mathrm{d} z^{2}} \\
|
|
|
F'''(z) &\equiv \frac{\mathrm{d^{3}} F(z)}{\mathrm{d} z^{3}}\end{aligned}\]]
|
|
|
F'''(z) &\equiv \frac{\mathrm{d^{3}} F(z)}{\mathrm{d} z^{3}}\end{aligned}
|
|
|
++++
|
|
|
|
|
|
we can expand the field off axis, with the result:
|
|
|
|
|
|
latexmath:[\[\begin{aligned}
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
B_x(\Delta_x, y, \Delta_z) &= -\frac{B_0 \frac{n}{R}}{\sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z}}} e^{-\frac{n}{R} \Delta_x} \sin \left[ \left( \sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}} \right) y \right] F(\Delta_z) \\
|
|
|
B_y(\Delta_x, y, \Delta_z) &= B_0 e^{-\frac{n}{R} \Delta_x} \cos \left[ \left( \sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}} \right) y \right] F(\Delta_z) \\
|
|
|
B_z(\Delta_x, y, \Delta_z) &= B_0 e^{-\frac{n}{R} \Delta_x} \left\{\frac{F'(\Delta_z)}{\sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}}} \sin \left[ \left( \sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}} \right) y \right] \right. \\
|
|
|
&- \frac{1}{2 \sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}}} \left(F'''(\Delta_z) - \frac{F'(\Delta_z) F''(\Delta_z)}{F(\Delta_z)} \right) \left[ \frac{\sin \left[ \left( \sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}} \right) y \right]}{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}} \right. \\
|
|
|
&- \left. \left. y \frac{\cos \left[ \left( \sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}} \right) y \right]}{\sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}}} \right] \right\}\end{aligned}\]]
|
|
|
&- \left. \left. y \frac{\cos \left[ \left( \sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}} \right) y \right]}{\sqrt{\frac{n^2}{R^2} + \frac{F''(\Delta_z)}{F(\Delta_z)}}} \right] \right\}\end{aligned}
|
|
|
++++
|
|
|
|
|
|
These expression are not well suited for numerical calculation, so, we
|
|
|
expand them about latexmath:[y] to latexmath:[O(y^2)] to obtain:
|
|
|
|
|
|
* In fringe field regions: latexmath:[\[\begin{aligned}
|
|
|
* In fringe field regions:
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
B_x(\Delta_x, y, \Delta_z) &\approx -B_0 \frac{n}{R} e^{-\frac{n}{R} \Delta_x} y \\
|
|
|
B_y(\Delta_x, y, \Delta_z) &\approx B_0 e^{-\frac{n}{R} \Delta_x} \left[ F(\Delta_z) - \left( \frac{n^2}{R^2} F(\Delta_z) + F''(\Delta_z) \right) \frac{y^2}{2} \right] \\
|
|
|
B_z(\Delta_x, y, \Delta_z) &\approx B_0 e^{-\frac{n}{R} \Delta_x} y F'(\Delta_z)\end{aligned}\]]
|
|
|
* In central region: latexmath:[\[\begin{aligned}
|
|
|
B_z(\Delta_x, y, \Delta_z) &\approx B_0 e^{-\frac{n}{R} \Delta_x} y F'(\Delta_z)\end{aligned}
|
|
|
++++
|
|
|
* In central region:
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
B_x(\Delta_x, y, \Delta_z) &\approx -B_0 \frac{n}{R} e^{-\frac{n}{R} \Delta_x} y \\
|
|
|
B_y(\Delta_x, y, \Delta_z) &\approx B_0 e^{-\frac{n}{R} \Delta_x} \left[ 1 - \frac{n^2}{R^2} \frac{y^2}{2} \right] \\
|
|
|
B_z(\Delta_x, y, \Delta_z) &\approx 0\end{aligned}\]]
|
|
|
B_z(\Delta_x, y, \Delta_z) &\approx 0\end{aligned}
|
|
|
++++
|
|
|
|
|
|
These are the expressions _OPAL-t_ uses to calculate the field inside an
|
|
|
`RBEND` or `SBEND`. First, a particle’s position inside the bend is
|
... | ... | @@ -1003,15 +1033,18 @@ Default Field Map (_OPAL-t_) |
|
|
Rather than force users to calculate the field of a dipole and then fit
|
|
|
that field to find Enge coefficients for the dipoles in their
|
|
|
simulation, we have a default set of values we use from [enge] that are
|
|
|
set when the default field map, ```1DPROFILE1-DEFAULT`'' is used:
|
|
|
set when the default field map, `1DPROFILE1-DEFAULT` is used:
|
|
|
|
|
|
latexmath:[\[\begin{aligned}
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
c_{0} &= 0.478959 \\
|
|
|
c_{1} &= 1.911289 \\
|
|
|
c_{2} &= -1.185953 \\
|
|
|
c_{3} &= 1.630554 \\
|
|
|
c_{4} &= -1.082657 \\
|
|
|
c_{5} &= 0.318111\end{aligned}\]]
|
|
|
c_{5} &= 0.318111\end{aligned}
|
|
|
++++
|
|
|
|
|
|
The same values are used for both the entrance and exit regions of the
|
|
|
magnet. In general they will give good results. (Of course, at some
|
... | ... | @@ -1041,7 +1074,7 @@ map format) map: |
|
|
0.318111
|
|
|
....
|
|
|
|
|
|
As one can see, the default magnet gap for ```1DPROFILE1-DEFAULT’`'' is
|
|
|
As one can see, the default magnet gap for `1DPROFILE1-DEFAULT` is
|
|
|
set to 2.0cm. This value can be overridden by the `GAP` attribute of the
|
|
|
magnet (see Section [RBend,SBend]).
|
|
|
|
... | ... | @@ -1480,7 +1513,7 @@ locates in multiple field regions, all the field maps are superposed. if |
|
|
SUPERPOSE is set to false, then only one field map, which has highest
|
|
|
priority, is used to do interpolation for the particle tracking. The
|
|
|
priority ranking is decided by their sequence in the list of RFMAPFN
|
|
|
argument, i.e., ``e1.h5hart'' has the highest priority and ``e4.h5hart''
|
|
|
argument, i.e., "e1.h5hart" has the highest priority and "e4.h5hart"
|
|
|
has the lowest priority.
|
|
|
|
|
|
Another method to model an RF cavity is to read the RF voltage profile
|
... | ... | @@ -1604,7 +1637,7 @@ L:: |
|
|
The length of the cavity (default: 0 m)
|
|
|
VOLT::
|
|
|
The peak RF voltage (default: 0 MV). The effect of the cavity is
|
|
|
latexmath:[\delta E=\text{\texttt{VOLT}}\cdot\sin(2\pi(\text{\texttt{LAG}}-\text{\texttt{HARMON}}\cdot f_0 t))].
|
|
|
latexmath:[\delta E=\texttt{VOLT}\cdot\sin(2\pi(\texttt{LAG}-\texttt{HARMON}\cdot f_0 t))].
|
|
|
LAG::
|
|
|
The phase lag [rad] (default: 0). In _OPAL-t_ this phase is in general
|
|
|
relative to the phase at which the reference particle gains the most
|
... | ... | @@ -1724,7 +1757,10 @@ L:: |
|
|
Full length of the cavity [mm].
|
|
|
|
|
|
The field inside the cavity is given by
|
|
|
latexmath:[\[\mathbf{E} = \big(0, 0, E_0(t)\sin[2\pi f(t) t+\phi(t)]\big)\]]
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\mathbf{E} = \big(0, 0, E_0(t)\sin[2\pi f(t) t+\phi(t)]\big)
|
|
|
++++
|
|
|
with no field outside the cavity boundary. There is no magnetic field or
|
|
|
transverse dependence on electric field.
|
|
|
|
... | ... | @@ -1733,7 +1769,7 @@ Time Dependence |
|
|
^^^^^^^^^^^^^^^
|
|
|
|
|
|
The `POLYNOMIAL_TIME_DEPENDENCE` element is used to define time
|
|
|
dependent parameters in RF cavities in terms of a fourth order
|
|
|
dependent parameters in RF cavities in terms of a third order
|
|
|
polynomial.
|
|
|
|
|
|
P0::
|
... | ... | @@ -1746,8 +1782,10 @@ P3:: |
|
|
Third order term in the polynomial expansion [nslatexmath:[^{-3}]].
|
|
|
|
|
|
The polynomial is evaluated as
|
|
|
latexmath:[\[g(t) = p_0 + p_1 t + p_2 t^2 + p_3 t^3 %% + p_4 t^4
|
|
|
.\]]
|
|
|
[latexmath]
|
|
|
++++
|
|
|
g(t) = p_0 + p_1 t + p_2 t^2 + p_3 t^3.
|
|
|
++++
|
|
|
|
|
|
An example of a Variable Frequency RF cavity of cyclotron with
|
|
|
polynomial time dependence of parameters is given below:
|
... | ... | @@ -1772,7 +1810,7 @@ RF_CAVITY: VARIABLE_RF_CAVITY, PHASE_MODEL="RF_PHASE", AMPLITUDE_MODEL="RF_AMPLI |
|
|
Traveling Wave Structure
|
|
|
~~~~~~~~~~~~~~~~~~~~~~~~
|
|
|
|
|
|
. Figure 6:The on-axis field of an S-band (2997.924 MHz) `TRAVELINGWAVE` structure. The field of a single cavity is shown between its entrance and exit fringe fields. The fringe fields extend one half wavelength (latexmath:[\lambda/2]) to either side.
|
|
|
.Figure 6:The on-axis field of an S-band (2997.924 MHz) `TRAVELINGWAVE` structure. The field of a single cavity is shown between its entrance and exit fringe fields. The fringe fields extend one half wavelength (latexmath:[\lambda/2]) to either side.
|
|
|
image:./figures/traveling-wave-structure/FINSB-RAC-field.png[]
|
|
|
|
|
|
An example of a 1D `TRAVELINGWAVE` structure field map is shown in
|
... | ... | @@ -1784,14 +1822,20 @@ _OPAL-t_ reads in this field map and constructs the total field of the |
|
|
structure fields and the exit fringe field.
|
|
|
|
|
|
The fringe fields are treated as standing wave structures and are given
|
|
|
by: latexmath:[\[\begin{aligned}
|
|
|
\mathbf{E_{entrance}}(\mathbf{r}, t) &= \mathbf{E_{from-map}}(\mathbf{r}) \cdot \text{\texttt{VOLT}} \cdot \cos \left( 2\pi \cdot \text{\texttt{FREQ}} \cdot t
|
|
|
by:
|
|
|
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{aligned}
|
|
|
\mathbf{E_{entrance}}(\mathbf{r}, t) &= \mathbf{E_{from-map}}(\mathbf{r}) \cdot \texttt{VOLT} \cdot \cos \left( 2\pi \cdot \texttt{FREQ} \cdot t
|
|
|
+ \phi_{entrance} \right) \\
|
|
|
\mathbf{E_{exit}}(\mathbf{r}, t) &= \mathbf{E_{from-map}}(\mathbf{r}) \cdot \text{\texttt{VOLT}} \cdot \cos \left( 2\pi \cdot \text{\texttt{FREQ}} \cdot t
|
|
|
\mathbf{E_{exit}}(\mathbf{r}, t) &= \mathbf{E_{from-map}}(\mathbf{r}) \cdot \texttt{VOLT} \cdot \cos \left( 2\pi \cdot \texttt{FREQ} \cdot t
|
|
|
+ \phi_{exit} \right)
|
|
|
\end{aligned}\]] where VOLT and FREQ are the field magnitude and
|
|
|
\end{aligned}
|
|
|
++++
|
|
|
where VOLT and FREQ are the field magnitude and
|
|
|
frequency attributes (see below).
|
|
|
latexmath:[ \phi_{entrance}= \text{\texttt{LAG}}], the phase attribute
|
|
|
latexmath:[ \phi_{entrance}= \texttt{LAG}], the phase attribute
|
|
|
of the element (see below). latexmath:[ \phi_{exit} ] is dependent
|
|
|
upon both LAG and the NUMCELLS attribute (see below) and is calculated
|
|
|
internally by _OPAL-t_.
|
... | ... | @@ -1800,16 +1844,19 @@ The field of the main accelerating structure is reconstructed from the |
|
|
center section of the standing wave solution shown in
|
|
|
Figure [FINSB-RAC-field] using
|
|
|
|
|
|
latexmath:[\[\begin{split}
|
|
|
\mathbf{E} ( \mathbf{r},t ) &= \frac{\text{\texttt{VOLT}}}{\sin (2 \pi \cdot \text{\texttt{MODE}})} \\
|
|
|
[latexmath]
|
|
|
++++
|
|
|
\begin{split}
|
|
|
\mathbf{E} ( \mathbf{r},t ) &= \frac{\texttt{VOLT}}{\sin (2 \pi \cdot \texttt{MODE})} \\
|
|
|
&\phantom{=}
|
|
|
\times \Biggl\{ \mathbf{E_{from-map}} (x,y,z) \cdot \cos \biggl( 2 \pi \cdot \text{\texttt{FREQ}} \cdot t + \text{\texttt{LAG}}+ \frac{\pi}{2} \cdot \text{\texttt{MODE}} \Bigr) \\
|
|
|
\times \Biggl\{ \mathbf{E_{from-map}} (x,y,z) \cdot \cos \biggl( 2 \pi \cdot \texttt{FREQ} \cdot t + \texttt{LAG}+ \frac{\pi}{2} \cdot \texttt{MODE} \Bigr) \\
|
|
|
&\phantom{= \times \Biggl\{}
|
|
|
+ \mathbf{E_{from-map}}(x,y,z+d) \cdot \cos \biggl( 2 \pi \cdot \text{\texttt{FREQ}} \cdot t + \text{\texttt{LAG}} + \frac{3 \pi}{2} \cdot \text{\texttt{MODE}} \Bigr) \Biggr\}
|
|
|
\end{split}\]]
|
|
|
+ \mathbf{E_{from-map}}(x,y,z+d) \cdot \cos \biggl( 2 \pi \cdot \texttt{FREQ} \cdot t + \texttt{LAG} + \frac{3 \pi}{2} \cdot \texttt{MODE} \Bigr) \Biggr\}
|
|
|
\end{split}
|
|
|
++++
|
|
|
|
|
|
where d is the cell length and is defined as
|
|
|
latexmath:[\text{d} = \lambda \cdot \text{\texttt{MODE}} ]. MODE is an
|
|
|
latexmath:[\text{d} = \lambda \cdot \texttt{MODE} ]. MODE is an
|
|
|
attribute of the element (see below). When calculating the field from
|
|
|
the map (latexmath:[\mathbf{E_{from-map}}(x,y,z)]), the longitudinal
|
|
|
position is referenced to the start of the cavity fields at
|
... | ... | @@ -1835,7 +1882,7 @@ L:: |
|
|
cells.
|
|
|
VOLT::
|
|
|
The peak RF voltage (default: 0 MV). The effect of the cavity is
|
|
|
latexmath:[\delta E=\text{\texttt{VOLT}}\cdot\sin(\text{\texttt{LAG}}- 2\pi\cdot\text{\texttt{FREQ}}\cdot t)].
|
|
|
latexmath:[\delta E=\texttt{VOLT}\cdot\sin(\texttt{LAG}- 2\pi\cdot\texttt{FREQ}\cdot t)].
|
|
|
LAG::
|
|
|
The phase lag [rad] (default: 0). In _OPAL-t_ this phase is in general
|
|
|
relative to the phase at which the reference particle gains the most
|
... | ... | @@ -1952,7 +1999,7 @@ The reference system for a collimator is a Cartesian coordinate system |
|
|
_OPAL-t_ mode
|
|
|
^^^^^^^^^^^^^
|
|
|
|
|
|
The `CCOLLIMATOR` isn’t supported. `ECOLLIMATOR`s and `RCOLLIMATOR`s
|
|
|
The `CCOLLIMATOR` isn’t supported. ``ECOLLIMATOR``s and ``RCOLLIMATOR``s
|
|
|
detect all particles which are outside the aperture defined by `XSIZE`
|
|
|
and `YSIZE`. Lost particles are saved in an H5hut file defined by
|
|
|
`OUTFN`. The `ELEMEDGE` defines the location of the collimator and `L`
|
... | ... | @@ -2002,8 +2049,8 @@ PARTICLEMATTERINTERACTION:: |
|
|
`PARTICLEMATTERINTERACTION` is an attribute of the element. Collimator
|
|
|
physics is only a kind of particlematterinteraction. It can be applied
|
|
|
to any element. If the type of `PARTICLEMATTERINTERACTION` is
|
|
|
`COLLIMATOR`, the material is defined here. The material ``Cu'',
|
|
|
``Be'', ``Graphite'' and ``Mo'' are defined until now. If this is not
|
|
|
`COLLIMATOR`, the material is defined here. The material "Cu",
|
|
|
"Be", "Graphite" and "Mo" are defined until now. If this is not
|
|
|
set, the particle matter interaction module will not be activated. The
|
|
|
particle hitting collimator will be recorded and directly deleted from
|
|
|
the simulation.
|
... | ... | @@ -2091,7 +2138,7 @@ ystart=-1226.85, yend=-1241.3; |
|
|
|
|
|
The particles probed on the PROBE are recorded in the ASCII file
|
|
|
_<inputfilename>.loss_. Please note that these particles are not deleted
|
|
|
in the simulation, however, they are recorded in the ``loss'' file.
|
|
|
in the simulation, however, they are recorded in the "loss" file.
|
|
|
|
|
|
[[stripper-opal-cycl]]
|
|
|
Stripper (_OPAL-cycl_)
|
... | ... | @@ -2120,7 +2167,7 @@ OPCHARGE:: |
|
|
Charge number of the out-coming particle. Negative value represents
|
|
|
negative charge.
|
|
|
OPMASS::
|
|
|
Mass of the out-coming particles. [/c]
|
|
|
Mass of the out-coming particles. [latexmath:[\mathrm{GeV/c^2}]]
|
|
|
OPYIELD::
|
|
|
Yield of the out-coming particle (the outcome particle number per
|
|
|
income particle) , the default value is 1.
|
... | ... | |