Download Ver. 2.76 User`s Manual

4 SECTIONS FORMAT
104
4.15 [ M a g n e t i c F i e l d ] section
4.15.1 Charged particle
You can set a magnetic field in the PH I TS calculation. Cell number (reg), magnetic field type (typ), half
distance of magnets (gap)(cm), magnetic field intensity (mgf), transformation (trcl) and time dependence (time)
should be defined as
[ Magnetic
reg
1
2
3
4
( 150 <
( 150 <
( 150 <
( 150 <
( 150 <
( 150 <
...
...
Field ]
typ
4
4
2
2
61 )
4
62 )
4
63 )
4
64 )
4
65 )
4
66 )
4
...
...
gap
10.00000
10.00000
10.00000
10.00000
13.00000
13.00000
13.00000
13.00000
13.00000
13.00000
........
........
mgf
-5.956540
6.416140
-7.611980
3.516000
7.880140
-7.440800
9.441010
-8.295220
3.694830
-2.099350
........
........
trcl
3
1
0
0
2
2
2
2
2
2
...
...
time
non
non
0.0
pi/2
non
non
non
non
non
non
...
...
The column of trcl is omissible. The zero for trcl means no transformation. The time is a parameter of user
defined time dependent magnetic field. The column of time is also omissible. The “non” for time means no time
dependence. Two subroutines, usrmgt1.f and usrmgt2.f are included in the source as user defined subroutines for
the time dependent magnetic field. The former is for Wobbler magnet, and the latter is for pulse magnet for neutron
optics. You can choose these two subroutine by usrmgt=1, 2 in the parameter section. For the Wobbler magnet,
“time” means phase of the magnet, starting time for pulse magnet, respectively.
In the above expression, reg is region number, typ can take 2 or 4 for dipole electromagnet, or quadrupole
electromagnet respectively.
mgf denotes the strength of the magnetic field (kG), and trcl is the coordinate
transformation number defined in [transform] section.
You can use the format ( { 2 - 5 } 8 9 ), and you can use the lattice and universe style as
( 6 < 10[1 0 0] < u=3 ). But you need to close a value by (
) if it is not a single numeric value.
By using this format, you can set different magnetic field for each lattice. If a cell is re-defined, the value,
which is defined at first, is used.
In the case of dipole magnet, the distances gap make no sense, but set any numeric. The magnetic field is
available not only in the void region, but also in the material where the normal reaction can be occurred.
z-axis is assumed to be the center of the magnetic field. The direction of the magnetic field is positive direction
of y-axis for dipole, i.e. the positive charge particle is bent to positive direction of x-axis when it goes to positive
direction of z-axis. For quadrupole, the positive particle is converged in x-axis, diverged in y-axis when it goes to
positive direction of z-axis. You need the coordinate transformation by trcl for different geometrical situation.
When you specify charge number of the projectile particle with izst in [source] section, the motion of the
particle with the number in the magnetic field is described. Using izst, PH I TS can simulate the motion of the
particle with charge states. The charge number defined with izst doesn’t change while the particle moves. It
should be noted that particles produced from nuclear reactions are not affected by the value of izst; the charge of
the produced particle is given as its atomic number.