Commit ab6bf61c authored by Thomas Planche's avatar Thomas Planche

got DIPOLE/edge-curvature to work OK

parent 3b9d12e4
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !BORO: Brho of 60 keV 238U+ = 544.12 kG.cm
543.687 !BORO: Brho of 60 keV 238U+ = 543.687 kG.cm
2 !KOBJ=2: 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'. Note: Brho=BORO*D
......@@ -10,7 +10,7 @@ zgoubi HRS tuto: DIPOLE
'DIPOLE'
2 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
65 120. !AT[deg.], RM[cm]
65 4.5343342 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
65 4.5307259 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
0 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.5 1e99 0.0 0.0 1e99 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
......@@ -29,12 +29,13 @@ atDIPOLEcenter.fai
'FIT'
2 !Number of physical parameters to be varied
2 5 0. 0.2 !element #2; vary parameter #64; coupling switch; relative range
2 64 0. 0.2 !element #2; vary parameter #64; coupling switch; relative range
2 1e-15 100 !Number of constraints; Convergence threshold; maximum number of iterations
3 1 2 2 120.0 1. 0 !IC (=3 for constraint on particle coordinate); Particle #; Particle coordinate (=2 for Y);element #;Wanted value; Weigth ; 0
3 1 3 2 0.0 1. 0 !IC (=3 for constraint on particle coordinate); Particle #; Particle coordinate (=3 for T);element #;Wanted value; Weigth ; 0
2 !Number of physical parameters to be varied
2 5 0. 0.2 !element #2; vary parameter #5 (BO); coupling switch; relative range
2 64 0. 0.2 !element #2; vary parameter #64 (RE); coupling switch; relative range
2 1e-8 100 !Number of constraints; Convergence threshold; maximum number of iterations
3 1 2 2 120.0 1. 0 !IC; Particle #; coordinate#=2 for Y;element #; wanted value; weigth ; 0
3 1 3 2 0.0 1. 0 !IC; Particle #; coordinate#=3 for T;element #; wanted value; weigth ; 0
'END'
......
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !60 keV U1+ 544.1192/544.1328
3.01 !read from particles from
1 3000 1 !IP1; IP2; IPStep (Read particles numbered IT1 to IT2, step ITStep)
1 1 1 !not used here
1. 1. 2. 0.25 1. 1. 1. * !scaling factors on Y,T,Z,P,S,DP
0 0 0 0 0 0 0 ! not used here
0 ! not sure what this one does
noEspread.beam ! name of the input file
'FAISCNL'
init.fai
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 223. 0.0 0.0 223. !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -223. 0. 0. -223. 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 223. 0.0 0.0 223. !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -223. 0. 0. -223. 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
curv_mass1.fai
'END'
'MATRIX'
2 0 !order of the map; 0:means calculate the matrix here
'FIT' 13
7 !Number of physical parameters to be varied
3 24 3.021 0.1
3 37 -3.021 0.1
3 40 -3.021 0.1
6 21 3.021 0.1
6 24 3.021 0.1
6 37 -3.021 0.1
6 40 -3.021 0.1
1 1e-15 100
2 1 22 8 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'
'MATRIX'
1 0
'END'
'FIT' 13
3 !Number of physical parameters to be varied
3 36 -3.020 0.1 !element#4 (DIPOLE); parameter #36 (exit edge angle); inverse (minus sign) correlation with element#4 parameter#20 (entrance edge angle); relative range
6 20 3.020 0.1
6 36 -3.020 0.1
2 1e-11 100 !Number of constraints; Convergence threshold; maximum number of iterations
1 1 2 8 0.0 1. 0 !transfer matrix; matrix element 1 1; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
1 2 1 8 0.0 1. 0 !transfer matrix; matrix element 2 2; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
'END'
'OBJET'
544.12 !60 keV U1+
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.
'OBJET'
544.12 !BORO[kg.cm]: Brho of 60 keV 238U+ in kG.cm
2 !KOBJ=2: 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'. Note: Brho=BORO*D
1 !1 or -9 (-9 disables the tracking of this particle)
......@@ -2,7 +2,7 @@ zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !60 keV U1+
543.687 !60 keV U1+
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.
......@@ -14,7 +14,7 @@ zgoubi HRS tuto: DIPOLE
'DIPOLE'
2 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
65 4.5307259 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.5 1e99 0.0 0.0 1e99 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
......
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !544.1192 or 544.1328
543.701 !543.687 or 543.701 and not 543.714, since we are aiming for a mass resolution of 20000 and not momentum resolution of 20000
3.01 !read from particles from
1 3000 1 !IP1; IP2; IPStep (Read particles numbered IT1 to IT2, step ITStep)
1 1 1 !not used here
......@@ -11,6 +11,7 @@ zgoubi HRS tuto: DIPOLE
noEspread.beam ! name of the input file
'DRIFT'
36.32357
......@@ -18,55 +19,49 @@ noEspread.beam ! name of the input file
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
65 4.5315270 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.474 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
45. 26.486450 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.474 -227.0 0. 0. -227.0 0. !W-; edge angle; R1; U1; U2; R2
-45. -26.486450 -227.0 0. 0. -227.0 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 127.96721 -0.349067 127.96721 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
0.3 !Integration step size[cm]
2 127.97493 -0.349067 127.97493 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'DRIFT'
36.32357
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
65 4.5315270 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.474 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
45. 26.486450 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.474 -227.0 0. 0. -227.0 0. !W-; edge angle; R1; U1; U2; R2
-45. -26.486450 -227.0 0. 0. -227.0 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 127.96721 -0.349067 127.96721 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
0.3 !Integration step size[cm]
2 127.97493 -0.349067 127.97493 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357 ! drift length=80-120*tan(20/deg.)
'DRIFT'
36.32357
'FAISCNL'
mass1.fai
'REBELOTE'
1 0.1 0 1
1
1 1 544.1328
mass2.fai
'END'
......
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1328 !544.1192 or 544.1328 or 544.1464?? mass resolution vs momentum resolution
3.01 !read from particles from
1 3000 1 !IP1; IP2; IPStep (Read particles numbered IT1 to IT2, step ITStep)
1 1 1 !not used here
1. 1. 1. 0.5 1. 1. 1. * !scaling factors on Y,T,Z,P,S,DP
0 0 0 0 0 0 0 ! not used here
0 ! not sure what this one does
noEspread.beam ! name of the input file
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 1e99 0.0 0.0 1e99 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -1e99 0. 0. -1e99 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 1e99 0.0 0.0 1e99 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -1e99 0. 0. -1e99 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
mass2.fai
'END'
This source diff could not be displayed because it is too large. You can view the blob instead.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "twomasses_edgeCurvature.eps"
set xlabel 'X/mm'
set ylabel 'X`/mrad'
set grid
titletext="Horizontal and vertical emittance: 3{/Symbol m}m, no DE/E, m/Dm=20000 \n"
set title titletext
set key left cent
cm2mm=10.0
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 1 title "^{238}U mass", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 3 title"20001/20000*^{238}U mass"
......@@ -2,9 +2,9 @@ zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !60 keV U1+
543.687 !60 keV U1+
5 !5: generate 11 particles used to calculate transfer matrix ('MATRIX')
0.0001 0.1 0.0001 0.1 .1 .001 !step size in Y; T; Z; P; S; D
0.00001 0.001 0.00001 0.001 .0 .0 !step size in Y; T; Z; P; S; D
0.0 0. 0. 0. 0. 1.
'DRIFT'
......@@ -14,18 +14,18 @@ zgoubi HRS tuto: DIPOLE
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
65 4.5307259 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.575 220.0 0.0 0.0 220.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
45. 26.5 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.575 -220.0 0. 0. -220.0 0. !W-; edge angle; R1; U1; U2; R2
-45. -26.5 -227.0 0. 0. -227.0 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.1 !Integration step size[cm]
0.3 !Integration step size[cm]
2 127.96721 -0.349067 127.96721 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
......@@ -34,15 +34,28 @@ zgoubi HRS tuto: DIPOLE
'MATRIX'
1 0
'FIT'
1 !Number of physical parameters to be varied
3 36 -3.020 0.1 !element#4 (DIPOLE); parameter #36 (exit edge angle); inverse (minus sign) correlation with element#4 parameter#20 (entrance edge angle); relative range
2 1e-11 100 !Number of constraints; Convergence threshold; maximum number of iterations
1 1 1 5 0.0 1. 0 !transfer matrix; matrix element 1 1; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
1 2 2 5 0.0 1. 0 !transfer matrix; matrix element 2 2; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
'FIT'
3 !Number of physical parameters to be varied
3 5 0. 0.2 !element #3 (DIPOLE); vary parameter #5 (BO); coupling switch; relative range
3 64 0. 0.2 !element #3 (DIPOLE); vary parameter #64 (RE); coupling switch; relative range
3 36 -3.020 0.1 !element #3 (DIPOLE); parameter #36 (exit edge angle); inverse (minus sign) correlation with element#4 parameter#20 (entrance edge angle); relative range
4 1e-8 100 !Number of constraints; Convergence threshold; maximum number of iterations
1 1 1 5 0.0 1. 0 !transfer matrix; matrix element 1 1; FIT after element#5 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
1 2 2 5 0.0 1. 0 !transfer matrix; matrix element 2 2; FIT after element#5 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
3 1 2 5 0.0 10. 0 !IC; Particle #; coordinate#=2 for Y;element #; wanted value; weigth ; 0
3 1 3 5 0.0 10. 0 !IC; Particle #; coordinate#=3 for T;element #; wanted value; weigth ; 0
'END'
'FIT'
3 !Number of physical parameters to be varied
3 36 -3.020 0.1 !element #3 (DIPOLE); parameter #36 (exit edge angle); inverse (minus sign) correlation with element#4 parameter#20 (entrance edge angle); relative range
3 5 0. 0.2 !element #3 (DIPOLE); vary parameter #5 (BO); coupling switch; relative range
3 64 0. 0.2 !element #3 (DIPOLE); vary parameter #64 (RE); coupling switch; relative range
4 1e-8 100 !Number of constraints; Convergence threshold; maximum number of iterations
3 1 2 5 0.0 1000. 0 !IC; Particle #; coordinate#=2 for Y;element #; wanted value; weigth ; 0
3 1 3 5 0.0 1000. 0 !IC; Particle #; coordinate#=3 for T;element #; wanted value; weigth ; 0
1 1 1 5 0.0 1. 0 !transfer matrix; matrix element 1 1; FIT after element#5 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
1 2 2 5 0.0 1. 0 !transfer matrix; matrix element 2 2; FIT after element#5 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
This diff is collapsed.
......@@ -2,14 +2,20 @@
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "twomasses.eps"
set xlabel 'X/mm'
set ylabel 'X`/mrad'
set grid
set ylabel "X'/mrad"
#set grid
titletext="Horizontal and vertical emittance: 3{/Symbol m}m, no DE/E, m/Dm=20000 \n"
set title titletext
set key left cent
cm2mm=10.0
set key at -0.5, 5
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 1 title "^{238}U mass", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 3 title"20001/20000*^{238}U mass"
set xrange [-1.5:0.5]
set yrange [-60:60]
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 1 notitle, \
"" w points pt 7 ps 0.7 lc 1 title "^{238}U mass ", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 2 notitle ,\
"" w points pt 7 ps 0.7 lc 2 title "^{238}U mass*20001/20000"
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !60 keV U1+
543.687 !BRHO of 60 keV 238U1+ in kG.cm
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.
......@@ -13,13 +13,13 @@ zgoubi HRS tuto: DIPOLE
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
65 4.5307259 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 200. 0.0 0.0 200. !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
45. 26.543694 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -200. 0. 0. -200. 0. !W-; edge angle; R1; U1; U2; R2
-45. -26.543694 -227.0 0. 0. -227.0 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.
......@@ -61,13 +61,13 @@ zgoubi HRS tuto: DIPOLE
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
65 4.5307259 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 200. 0.0 0.0 200. !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
45. 26.543694 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -200. 0. 0. -200. 0. !W-; edge angle; R1; U1; U2; R2
-45. -26.543694 -227.0 0. 0. -227.0 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.
......
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1328 !544.1192 or 544.1328
3.01 !read from particles from
1 3000 1 !IP1; IP2; IPStep (Read particles numbered IT1 to IT2, step ITStep)
1 1 1 !not used here
1. 1. 1. 0.5 1. 1. 1. * !scaling factors on Y,T,Z,P,S,DP
0 0 0 0 0 0 0 ! not used here
0 ! not sure what this one does
noEspread.beam ! name of the input file
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 220.000 0.0 0.0 220.000 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -220.000 0. 0. -220.000 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 220.000 0.0 0.0 220.000 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -220.000 0. 0. -220.000 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
mass2-ec.fai
'END'
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "twomasses_edgeCurvature.eps"
set xlabel 'X/mm'
set ylabel 'X`/mrad'
set grid
titletext="Horizontal and vertical emittance: 3{/Symbol m}m, no DE/E, m/Dm=20000 \n"
set title titletext
set key left cent
cm2mm=10.0
plot "mass1-ec.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 1 title "^{238}U mass", \
"mass2-ec.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 3 title"20001/20000*^{238}U mass"
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !544.1192 or 544.1328 and not 544.1464, since we are aiming for a mass resolution of 20000 and not momentum resolution of 20000
543.687 !543.687 or 543.701 and not 543.714, since we are aiming for a mass resolution of 20000 and not momentum resolution of 20000
3.01 !read from particles from
1 3000 1 !IP1; IP2; IPStep (Read particles numbered IT1 to IT2, step ITStep)
1 1 1 !not used here
......@@ -19,7 +19,7 @@ noEspread.beam ! name of the input file
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX