Commit f21cb44f authored by Thomas Planche's avatar Thomas Planche

initial commit

parents
# Compiled source #
###################
# outputs (not inputs) #
########################################
*.fai
*.plt
*.ps
# logs and backup files #
#########################
*.log
*~
# Other stuffs#
#########################
**HRS FIT DIPOLE
!!Simplest HRS (i.e. magnifying section quads OFF) but only up to the center of the first dipole;
!!the purpose here is to fit dipole position and strength to get the righ beam position and angle in the center of the dipole
'OBJET' 1
546.198 !60 keV U1+ = 546.198
2 !2: All the initial coordinates must be entered explicitly
1 1 !total number of particles; number of distinct momenta
0. 0. 0. 0. 0. 1. ’o’ !Y; T; Z; P; S; D; 'marker'
1 !1 or -9 (-9 disables the tracking of this particle)
'DRIFT' 2
80.0 !80 cm long drift
'CHANGREF' 3
YS -120.0 ZR 20. !shift the optical axis by 120 cm (the reference radius of the dipole); then rotate around the vertical axis by 20 deg (20=AT-W)
'DIPOLE' 4
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
65 120. !AT, RM
65 4.6 0. 0. 0. !ACENT;B0; N; B; GX
7.0 -1. !ENTRANCE FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.47 225.50 0.0 0.0 225.50 !W+; edge angle; R1; U1; U2; R2
7.0 -1. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.47 -225.50 0.0 0.0 -225.50 0. !W-; edge angle; R1; U1; U2; R2
0. 0. !face 3 (Not used here)
0 0. 0. 0. 0. 0. 0. 0.
0. 0. 0. 0. 0. 0. 0.
2 4. !IRD(=2, 25 or 4)
0.5 !Integration step size[cm]
2 0. 0. 0. 0. !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'FIT' 5
2 !Number of physical parameters to be varied
3 1 0. 0.1 !in element #3 (i.e. CHANGREF); vary parameter #1 (i.e YS); coupling switch OFF=0; relative range +/- 10%
4 5 0. 0.2 !in element #5 (i.e. DIPOLE); vary parameter #6 (i.e. B0); coupling switch OFF=0; relative range +/- 20%
2 1e-10 500 !Number of constraints; Convergence threshold; maximum number of iterations
3 1 2 4 120.0 1. 0 !IC (type of constraint. =3 is for a constraint on a particle coordinate); Particle number; Particle coordinate (=2 for Y);Where?:at the end of element #4(i.e. DIPOLE);Wanted value(=120 cm); Weigth ; 0:no additional parameters
3 1 3 4 0.0 1. 0 !IC (type of constraint. =3 is for a constraint on a particle coordinate); Particle number; Particle coordinate (=3 for T);Where?:at the end of element #4(i.e. DIPOLE);Wanted value(=0 deg.); Weigth ; 0:no additional parameters
'END' 6
!!The FIT should converge to YS=-120.2408355 and B0=4.552771183 after 105 iteration.
!!Yes I know, I don't need so many significant figures, but I keep them just to track possible changes with different versions of zgoubi
** HSR FIT edge angle
!!Simplest HRS (i.e. magnifying section quads OFF)
!!The purpose here is to find the right edge angle to have 180 deg. of horizontal phase advance through the system
'OBJET' 1
546.198 !60 keV U1+ = 546.198
5 !5: generate 11 particles used to calculate transfer matrix ('MATRIX')
0.01 0.01 0.01 0.01 .1 .001 !step size in Y; T; Z; P; S; D
0.0 0. 0. 0. 0. 1.
'DRIFT' 2
80.0 !80 cm long drift
'CHANGREF' 3
YS -120.24084 ZR 20. !shift the optical axis (the value comes from runing 1_FIT_DIPOLE.in); then rotate around the vertical axis by 20 deg (20=AT-W)
'DIPOLE' 4
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT, RM
65 4.552771 0. 0. 0. !ACENT;B0 (the value comes from runing 1_FIT_DIPOLE.in); N; B; GX
7.0 -1. !ENTRANCE FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.47 225.50 0.0 0.0 225.50 !W+; edge angle; R1; U1; U2; R2
7.0 -1. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.47 -225.50 0.0 0.0 -225.50 0. !W-; edge angle; R1; U1; U2; R2
0. 0. !face 3 (Not used here)
0 0. 0. 0. 0. 0. 0. 0.
0. 0. 0. 0. 0. 0. 0.
2 4. !IRD(=2, 25 or 4)
0.5 !Integration step size[cm]
2 0. 0. 0. 0. !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'CHANGREF' 5
ZR 20. YS 120.24084
'DRIFT' 6
80.0 !80 cm long drift
'DRIFT' 7
80.0
'CHANGREF' 8
YS -120.24084 ZR 20. !shift the optical axis (the value comes from runing 1_FIT_DIPOLE.in); then rotate around the vertical axis by 20 deg (20=AT-W)
'DIPOLE' 9
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT, RM
65 4.552771 0. 0. 0. !ACENT;B0 (the value comes from runing 1_FIT_DIPOLE.in); N; B; GX
7.0 -1. !ENTRANCE FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.47 225.50 0.0 0.0 225.50 !W+; edge angle; R1; U1; U2; R2
7.0 -1. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.47 -225.50 0.0 0.0 -225.50 0. !W-; edge angle; R1; U1; U2; R2
0. 0. !face 3 (Not used here)
0 0. 0. 0. 0. 0. 0. 0.
0. 0. 0. 0. 0. 0. 0.
2 4. !IRD(=2, 25 or 4)
0.5 !Integration step size[cm]
2 0. 0. 0. 0. !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'CHANGREF' 10
ZR 20. YS 120.24084
'DRIFT' 11
80.0
'MATRIX' 12
1 0 !order of the matrix/map; 0:means calculate the matrix here
'FIT' 13
3 !Number of physical parameters to be varied
4 36 -4.020 0.01 !element#4 (DIPOLE); parameter #36 (exit edge angle); inverse (minus sign) correlation with element#4 parameter#20 (entrance edge angle); relative range
9 20 4.020 0.01 !similar than before. Note: constraint w.r.t. the same element-parameter 4-20
9 36 -4.020 0.01 !similar than before. Note: constraint w.r.t. the same element-parameter 4-20
2 1e-12 100 !Number of constraints; Convergence threshold; maximum number of iterations
1 1 2 12 0.0 1. 0 !transfer matrix; matrix element 1 2; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
1 2 1 12 0.0 1. 0 !transfer matrix; matrix element 2 1; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
'END' 14
** HSR FIT edge curvature
!!Simplest HRS (i.e. magnifying section quads OFF)
!!The purpose here is to find the edge curvature that minimized the x-x'^2 (also called x-aa, or 1-22) aberation
'OBJET' 1
546.198 !60 keV U1+ = 546.198
6 !6: generate 61 particles used to calculate 2nd order map ('MATRIX')
0.01 0.01 0.01 0.01 .1 .001 !step size in Y; T; Z; P; S; D
0.0 0. 0. 0. 0. 1.
'DRIFT' 2
80.0 !80 cm long drift
'CHANGREF' 3
YS -120.24076 ZR 20. !shift the optical axis (the value comes from runing 1_FIT_DIPOLE.in); then rotate around the vertical axis by 20 deg (20=AT-W)
'DIPOLE' 4
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT, RM
65 4.552787 0. 0. 0. !ACENT;B0 (the value comes from runing 1_FIT_DIPOLE.in); N; B; GX
7.0 -1. !ENTRANCE FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.47452 225.58 0.0 0.0 225.58 !W+; edge angle; R1; U1; U2; R2
7.0 -1. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.47452 -225.58 0.0 0.0 -225.58 0. !W-; edge angle; R1; U1; U2; R2
0. 0. !face 3 (Not used here)
0 0. 0. 0. 0. 0. 0. 0.
0. 0. 0. 0. 0. 0. 0.
2 4. !IRD(=2, 25 or 4)
0.5 !Integration step size[cm]
2 0. 0. 0. 0. !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'CHANGREF' 5
ZR 20. YS 120.24076
'DRIFT' 6
80.0
'DRIFT' 7
80.0
'CHANGREF' 8
YS -120.24076 ZR 20. !shift the optical axis (the value comes from runing 1_FIT_DIPOLE.in); then rotate around the vertical axis by 20 deg (20=AT-W)
'DIPOLE' 9
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT, RM
65 4.552787 0. 0. 0. !ACENT;B0 (the value comes from runing 1_FIT_DIPOLE.in); N; B; GX
7.0 -1. !ENTRANCE FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.47452 225.58 0.0 0.0 225.58 !W+; edge angle; R1; U1; U2; R2
7.0 -1. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.47452 -225.58 0.0 0.0 -225.58 0. !W-; edge angle; R1; U1; U2; R2
0. 0. !face 3 (Not used here)
0 0. 0. 0. 0. 0. 0. 0.
0. 0. 0. 0. 0. 0. 0.
2 4. !IRD(=2, 25 or 4)
0.5 !Integration step size[cm]
2 0. 0. 0. 0. !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'CHANGREF' 10
ZR 20. YS 120.24076
'DRIFT' 11
80.0
'MATRIX' 12
2 0 !order of the map; 0:means calculate the matrix here
'FIT' 13
7 !Number of physical parameters to be varied
4 24 4.021 0.02
4 37 -4.021 0.02
4 40 -4.021 0.02
9 21 4.021 0.02
9 24 4.021 0.02
9 37 -4.021 0.02
9 40 -4.021 0.02
1 1e-15 100
2 1 22 12 0.0 1. 0 !2:2nd order map element; transfer map element 1-22 (i.e. x-x'^2); Desired value; 0:no additional parameters
'END' 14
! 'FIT' 13
! 10
! 4 36 -4.020 0.01
! 9 20 4.020 0.01
! 9 36 -4.020 0.01
! 4 24 4.021 0.01
! 4 37 -4.021 0.01
! 4 40 -4.021 0.01
! 9 21 4.021 0.01
! 9 24 4.021 0.01
! 9 37 -4.021 0.01
! 9 40 -4.021 0.01
! 3 1e-15 100
! 1 1 2 12 0.0 1. 0
! 1 2 1 12 0.0 1. 0
! 2 1 22 12 0.0 1. 0
** HSR no quad using DIPOLE
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1
1 121 1 1 1 1
0.0 1.0 0.0 0.0 0.0 1.
0.0 0.0 0.0 0.0 0.0 1.
'DRIFT' 2
80.5 !80 cm long drift
'CHANGREF' 3
YS -120.24076 ZR 20. !shift the optical axis (the value comes from runing 1_FIT_DIPOLE.in); then rotate around the vertical axis by 20 deg (20=AT-W)
'DIPOLE' 4
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT, RM
65 4.552774 0. 0. 0. !ACENT;B0 (the value comes from runing 1_FIT_DIPOLE.in); N; B; GX
7.0 -1. !ENTRANCE FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.56 223.79 0.0 0.0 223.79 !W+; edge angle; R1; U1; U2; R2
7.0 -1. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.56 -223.79 0.0 0.0 -223.79 0. !W-; edge angle; R1; U1; U2; R2
0. 0. !face 3 (Not used here)
0 0. 0. 0. 0. 0. 0. 0.
0. 0. 0. 0. 0. 0. 0.
2 4. !IRD(=2, 25 or 4)
0.5 !Integration step size[cm]
2 0. 0. 0. 0. !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'CHANGREF' 5
ZR 20. YS 120.24076
'DRIFT' 6
80.0
'DRIFT' 7
80.0
'CHANGREF' 8
YS -120.24076 ZR 20. !shift the optical axis (the value comes from runing 1_FIT_DIPOLE.in); then rotate around the vertical axis by 20 deg (20=AT-W)
'DIPOLE' 9
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT, RM
65 4.552774 0. 0. 0. !ACENT;B0 (the value comes from runing 1_FIT_DIPOLE.in); N; B; GX
7.0 -1. !ENTRANCE FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.56 223.79 0.0 0.0 223.79 !W+; edge angle; R1; U1; U2; R2
7.0 -1. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.56 -223.79 0.0 0.0 -223.79 0. !W-; edge angle; R1; U1; U2; R2
0. 0. !face 3 (Not used here)
0 0. 0. 0. 0. 0. 0. 0.
0. 0. 0. 0. 0. 0. 0.
2 4. !IRD(=2, 25 or 4)
0.5 !Integration step size[cm]
2 0. 0. 0. 0. !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'CHANGREF' 10
ZR 20. YS 120.24076
'DRIFT' 11
80.5
'FAISTORE' 12
zgoubi.fai
1
'END' 13
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "zgoubifai.eps"
set xlabel 'X [mm]'
set ylabel 'X` [mrad]'
set grid
titletext="Simplest HRS, multipole OFF"
labeltext="(c) `whoami`, "
today="`date +%Y/%b/%d`"
set label labeltext.today at screen .01, screen .02
set title titletext
set key top left
cm2mm=10.0
plot "thinBeam.fai" u ($10*cm2mm):11 with points title "edge angle 26.47 deg, curvature= 1/2.255m",\
"thinBeam_noCurvature.fai" u ($10*cm2mm):11 with points lc 3 title "edge angle 26.47 deg, curvature= 0",\
"zgoubi.fai" u ($10*cm2mm):11 with points lc 4 title "as build: 26.56 deg., 1/223.79m"
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "plots/BzXplot.eps"
set xlabel 'Theta [deg.]'
set ylabel 'Bz [T]'
set grid
titletext="HRS TOSCA reference trajectory"
labeltext="(c) `whoami`, "
today="`date +%Y/%b/%d`"
set label labeltext.today at screen .01, screen .02
set title titletext
set key top left
cm2mm=10.0
plot "zgoubi.plt" using ($22/pi*180.):($25/10.) with p pt 0
This diff is collapsed.
This source diff could not be displayed because it is too large. You can view the blob instead.
** HSR using TOSCA map
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1
1 121 1 1 1 1
0.0 1.0 0.0 0.0 0.0 1.
0.0 0.0 0.0 0.0 0.0 1.
! 'OBJET' 1
! 546.198 !60 keV U1+ = 546.198
! 5 !5: generate 11 particles used to calculate transfer matrix ('MATRIX')
! 0.01 0.01 0.01 0.01 .1 .001 !step size in Y; T; Z; P; S; D
! 0.0 0. 0. 0. 0. 1.
'PARTICUL' 2
22172.3 1.602176487E-19 0.0 0.0 0.0
'DRIFT' 3
25.500
'TOSCA' 4
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling; X coordinate scaling; Y; Z
HEADER_2 ! Title
281 65 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
Dipole2D.map !File name
0 0 0 0. ! ID (see manual);
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.1 !Integration step size[cm]
2 !KPOS, normally=2
132.5586547 -0.436332313 132.5586547 0.436332313 !RE has been fitted to be =RS; TE=-25deg(in rad) whic is the distance between the effective edge and the edge of th efield map; RS=RE; TS=-TE
'DRIFT' 5
2.0430 !24.0430 - 22.0 (which is the half length of the multipole field map)
'MAP2D-E' 6
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
0.000001 0.1 0.1!Normalization factor for E (to be put in MV/cm), X (cm), and Y(cm)
HEADER_0 !this actually only works if you have n=0 header lines. It crashes if n>0 and I could not figure out why
45 27 !Number of longitudinal and horizontal-transverse nodes of the mesh (the Z elevation is arbitrary)
Multipole2D.map !file name
0 0. 0. 0.
4
1.
1 0. 0. 0.
'DRIFT' 7
2.0430 !24.0430 - 22.0 (which is the half length of the multipole field map)
'TOSCA' 8
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling; X coordinate scaling; Y; Z
HEADER_2 ! Title
281 65 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
Dipole2D.map !File name
0 0 0 0. ! ID (see manual);
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.1 !Integration step size[cm]
2 !KPOS, normally=2
132.5586547 -0.436332313 132.5586547 0.436332313 !RE has been fitted to be =RS; TE=-25deg(in rad) whic is the distance between the effective edge and the edge of th efield map; RS=RE; TS=-TE
'DRIFT' 9
25.500
'FAISTORE' 10
zgoubi.fai
1
! 'MATRIX' 9
! 1 0 !order of the matrix/map; 0:means calculate the matrix here
! 'FIT' 10
! 1 !Number of physical parameters to be varied
! 2 1 7.001 0.5
! 2
! 1 1 2 9 0.0 1. 0 !transfer matrix;
! 1 2 1 9 0.0 1. 0 !transfer matrix;
'END' 11
THREED TYPE=SURFACE VECTOR=NO XORIGIN=0 YORIGIN=1200 ZORIGIN=106 ROTX=1.0E-04 ROTY=179.999 ROTZ=1.0E-04 XASPECT=1 YASPECT=1 ZASPECT=1 SIZE=704.188791082898 FACETANGLE=10 PERSPECTIVE=NO LINECOLOUR=YES OPTION=SETVIEW
/++ GRID X0=0 Y0=500 Z0=0 DXG=20 DYG=20 DZG=1 NXG=101 NYG=51 NZG=1 FILE='forZgoubi' BINARY=NO FORMAT=2 F1=Y UNIT1G=LENGU F2=Z UNIT2G=LENGU F3=X UNIT3G=LENGU F4=0 UNIT4G=FLUXU F5=Bz UNIT5G=FLUXU F6=0 UNIT6G=FLUXU UNIT7G=1 UNIT8G=1 UNIT9G=1 UNIT10G=1 UNIT11G=1 UNIT12G=1
/++ yes
//global WCS
SET XLOCAL=0 YLOCAL=0 ZLOCAL=0 TLOCAL=0,PLOCAL=0,SLOCAL=0 LOOK=ANYWHERE ABORT=YES
$const #mm2cm 0.1
$const #T2kG 10.
$const #pi ATAN(1)*4
$const #rmin 104. |//cm
$const #rmax 136. |//cm good field region : 120+/-16 cm
$const #Dr 1.0 |//cm
$const #Dtheta 0.5 |//deg
$const #Dz 0 |//cm
//open output file
$open 1 forZgoubi_MOD22.dat overwrite
$format 1 fixed 16 8
$format 2 character 16
/++ R0,DR,DTTA,DZ
$assign 1 1 1 1
$write 1 #rmin #Dr #Dtheta #Dz
$assign 2 2 2 2 2 2
$write 1 'Y[cm]' 'Z[cm]' 'X[cm]' 'Br[kG]' 'Bth[kG]' 'Bz[kG]'
$assign 1 1 1 1 1 1
/++ $do #r 900 1500 20
$do #r #rmin/#mm2cm #rmax/#mm2cm #Dr/#mm2cm
$do #th -25 115 #Dtheta
$const #x #r*SIND(-#th+45.0)
$const #y #r*COSD(-#th+45.0)
$const #z 0
POINT XP=#x YP=#y ZP=0 COMP=bz
$write 1 #y*#mm2cm #z*#mm2cm #x*#mm2cm 0.0 BZ*#T2kG 0.0
/++ $write 1 #r*#mm2cm #th*#pi/180. #z*#mm2cm 0.0 0.0 BZ*#T2kG
$end do
$end do
$close 1
../../opera/multipole/geneMultipole2DMap.comi
\ No newline at end of file
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "plots/thinBeam.eps"
set xlabel 'X [mm]'
set ylabel 'X` [mrad]'
set grid
titletext="HRS 'pure separetor'"
labeltext="`whoami`, "
today="`date +%Y/%b/%d`"
set label labeltext.today at screen .01, screen .02
set title titletext
set key top right
cm2mm=10.0
set xrange[-1:1]
plot 'ref.fai' u ($10*cm2mm):11 with points lc 0 pt 5 ps 0.5 title "Multipole OFF",\
'zgoubi.fai' u ($10*cm2mm):11 with points lc 1 pt 1 title "100 V on the rods A"
This diff is collapsed.
**HRS FIT TOSCA !!Simplest HRS (i.e. magnifying section quads OFF) using a TOSCA map
'OBJET' 1
546.198 !60 keV U1+ = 546.198
2 !2: All the initial coordinates must be entered explicitly
1 1 !total number of particles; number of distinct momenta
0. 0. 0. 0. 0. 1. ’o’ !Y; T; Z; P; S; D; 'marker'
1 !1 or -9 (-9 disables the tracking of this particle)
'TOSCA' 3
0 2 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling; X coordinate scaling; Y; Z
HEADER_2 ! Title
281 61 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
2D.map !File name
0 0 0 0. ! ID (see manual);
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.1 !Integration step size[cm]
2 !KPOS, normally=2
132.5586547 -0.436332313 132.5586547 0.436332313 !RE has been fitted to be =RS; TE=-25deg(in rad) whic is the distance between the effective edge and the edge of th efield map; RS=RE; TS=-TE
'FAISTORE' 3
zgoubi.fai
1
'FIT' 4
1 !number of variables
2 70 2.072 0.5 !element#2 (TOSCA); parameter #70 (RE); correlation set so that RE=RS (which is parameters#72 of element #2)
1 !number of constraints
3 1 2 2 0.0 1. 0 !3:particle coordinate; 1: particle #1; 2: Y(i.e. R); where?:at the end of element #2; wanted value=0.0; weight; 0:no additional parameters
'END' 6
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "plots/BzXplot.eps"
set xlabel 'Theta [deg.]'
set ylabel 'Bz [T]'
set grid
titletext="HRS TOSCA reference trajectory"
labeltext="(c) `whoami`, "
today="`date +%Y/%b/%d`"
set label labeltext.today at screen .01, screen .02
set title titletext
set key top left
cm2mm=10.0
plot "zgoubi.plt" using ($22/pi*180.):($25/10.) with p pt 0
** HSR using TOSCA map
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1
1 121 1 1 1 1
0.0 1.0 0.0 0.0 0.0 1.
0.0 0.0 0.0 0.0 0.0 1.
! 'OBJET' 1
! 546.198 !60 keV U1+ = 546.198
! 5 !5: generate 11 particles used to calculate transfer matrix ('MATRIX')
! 0.01 0.01 0.01 0.01 .1 .001 !step size in Y; T; Z; P; S; D
! 0.0 0. 0. 0. 0. 1.
'DRIFT' 2
25.500
'TOSCA' 3
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling; X coordinate scaling; Y; Z
HEADER_2 ! Title
281 61 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
2D.map !File name
0 0 0 0. ! ID (see manual);
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.1 !Integration step size[cm]
2 !KPOS, normally=2
132.5586547 -0.436332313 132.5586547 0.436332313 !RE has been fitted to be =RS; TE=-25deg(in rad) whic is the distance between the effective edge and the edge of th efield map; RS=RE; TS=-TE
'DRIFT' 4
24.0430
'DRIFT' 5
24.0430
'TOSCA' 6
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling; X coordinate scaling; Y; Z
HEADER_2 ! Title
281 61 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
2D.map !File name
0 0 0 0. ! ID (see manual);
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.1 !Integration step size[cm]