Commit d98bd95d authored by Thomas Planche's avatar Thomas Planche

untracked DIPOLE/*.eps

parent ee9c0b82
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !60 keV U1+
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'
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 242. 0.0 0.0 242. !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 -242. 0. 0. -242. 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 242. 0.0 0.0 242. !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 -242. 0. 0. -242. 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.)
'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'
\ No newline at end of file
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !60 keV U1+ 544.1192/544.1328
544.1192 !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
......@@ -66,48 +66,3 @@ mass1.fai
'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)
zgoubi HRS tuto: DIPOLE
'OBJET'
544.12 !60 keV U1+
3.01 !read from particles from file
1 3000 1
1 1 1
1 1 2.0 0.25 1 1 1 *
0 0 0 0 0 0 0
0
input_noEspread.fai
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'
......@@ -19,13 +21,13 @@ input_noEspread.fai
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343183 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
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]
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 -1e99 0. 0. -1e99 0. !W-; edge angle; R1; U1; U2; R2
-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.
......@@ -42,13 +44,13 @@ input_noEspread.fai
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343183 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
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]
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 -1e99 0. 0. -1e99 0. !W-; edge angle; R1; U1; U2; R2
-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.
......@@ -60,8 +62,38 @@ input_noEspread.fai
'DRIFT'
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
mass1.fai
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'
......
This diff is collapsed.
This diff is collapsed.
......@@ -5,15 +5,16 @@ 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"
labeltext="(c) `whoami`, "
today="`date +%Y/%b/%d`"
#set label labeltext.today at screen .01, screen .02
set title titletext
set key
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"
set key right
set output "test.eps"
plot "init.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 0 title "init",\
"curv_mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 1 title "curved edge",\
"curv_mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 3 title "curved edge"
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment