Commit 6759d7e4 authored by Thomas Planche's avatar Thomas Planche
Browse files

added examples/TRIUMF-HRS

parent b846c20c
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !BORO: Brho of 60 keV 238U+ = 544.12 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)
'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
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]
0. 0. !EXIT FIELD BOUNDARY: (not used for now)
0. 0. 0. 0. 0. 0. 0. 0.
0. 0. -1e99 0.0 0.0 -1e99 0.
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 0. 0.0 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'FAISCNL'
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
'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)
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !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.
'DRIFT'
36.32357
'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
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]
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.5 -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.96721 -0.349067 127.96721 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'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
'END'
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 200. 0.0 0.0 200. !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
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
'MATRIX'
2 0 !order of the map; 0:means calculate the matrix here
'FIT'
4 !Number of physical parameters to be varied
3 5 0. 0.01
3 64 3.66 0.01
2 1 0.0 0.05
3 21 -3.037 1.0
4 1e-10 100
3 1 2 4 0.0 10. 0
3 1 3 4 0.0 10. 0
1 2 2 5 0.0 3. 0
2 2 22 5 0.0 1. 0
'END'
'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 200. 0.0 0.0 200. !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
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
# COORDINATES - STORAGE FILE, 21-09-2017 09:11:12. zgoubi HRS tuto: DIPOLE
# 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , 10, 11, 12, 13, 14, 15 , 16 , 17 , 18 , 19 , 20 , 21 , 22 , 23 , 24 , 25 , 26, 27 , 28 , 29 ,30, 31 , 32 , unused, 34 , 35 , 36, 37 , 38 , 39 , 40 , 41/lbl1, 42/lbl2, 43
# KEX, Do-1, Yo, To, Zo, Po, So, to, D-1, Y, T, Z, P, S, time, SXo, SYo, SZo, modSo, SX, SY, SZ, modS, ENEKI, ENERG, IT, IREP, SORT, M, Q, G, tau, unused, RET, DPR, PS, BORO, IPASS, NOEL, KLEY, LABEL1, LABEL2, LET
# int, float, cm, mrd, cm, mrd, cm, mu_s, float, cm, mrd, cm, mrd, cm, mu_s, float,float,float, float, float,float,float,float, MeV, MeV, int, int, cm, MeV/c2, C, float, float, float, float, float, float, kG.cm, int, int, string, string, string, string
1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.1999994623915468E+02 -1.0292160113640624E-03 0.0000000000000000E+00 0.0000000000000000E+00 1.3802143228181038E+02 4.6038994177035227E-03 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 1.63122832E+02 1.63122832E+02 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.44119200E+02 1 2 'DIPOLE ' ' ' ' ' 'o'
1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.1999999997587769E+02 1.0137501949003536E-08 0.0000000000000000E+00 0.0000000000000000E+00 1.3802144849479606E+02 4.6038999585105121E-03 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 1.63122832E+02 1.63122832E+02 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.44119200E+02 1 2 'DIPOLE ' ' ' ' ' 'o'
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
zgoubi HRS tuto: DIPOLE
'OBJET'
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
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'
mass1.fai
'END'
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "traj.eps"
set multiplot layout 1,2
set xlabel "Theta/deg."
set ylabel "R/cm" offset 2,0,0
set xrange[0:65]
set size ratio 1
unset key
plot "zgoubi.plt" using ($22*180./pi):($10) w l
th0=pi*65./180.
set xlabel "R/cm*sin(th)"
set ylabel "R/cm*cos(th)" offset 2,0,0
set xrange[-125:0]
set yrange[0:125]
set size ratio -1
unset key
plot "zgoubi.plt" using ($10*sin($22-th0)):($10*cos($22-th0)) w l
This diff is collapsed.
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
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"
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
zgoubi HRS tuto: DIPOLE
'OBJET' 1
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' 2
36.32357
'DIPOLE' 3
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65.0000000 4.5343313 0.0000000 0.0000000 0.0000000 ! ACNT, HNORM, indices
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]
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
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 128.0223520 -0.3490670 127.9674800 0.3490670 ! KPOS, RE, TE, RS, TS
'DRIFT' 4
36.32357
'MATRIX' 5
2 0 !order of the map; 0:means calculate the matrix here
! 'FIT' 6
! 4 !Number of physical parameters to be varied
! 3 5 0. 0.01
! 3 64 3.66 0.01
! 2 1 0.0 0.05
! 3 21 -3.037 1.0
! 4 1e-10 100
! 3 1 2 4 0.0 10. 0
! 3 1 3 4 0.0 10. 0
! 1 2 2 5 0.0 3. 0
! 2 2 22 5 0.0 1. 0
'END' 7
This diff is collapsed.
This diff is collapsed.
** HSR using TOSCA, FIT drift length
'OBJET'
546.198 !60 keV U1+ = 546.198
5 !5: generate 11 particles used to calculate transfer matrix ('MATRIX')
0.001 0.1 0.001 0.1 0. 0. !step size in Y; T; Z; P; S; D
0.0 0. 0. 0. 0. 1.
'DRIFT'
20.0
'TOSCA'
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling (adjusted so that R at the center=120cm+60cm*sqrt(2), 120cm is reference radius of the HRS, 60cm is the amound by which the map has been enlarge to include fringe fields); X coordinate scaling; Y; Z
HEADER_2 ! Title
181 47 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
TOSCAMOD22.dat !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
180.13967 0. 180.13967 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'DRIFT'
20.
'DRIFT'
20.
'TOSCA'
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling (adjusted so that R at the center=120cm+60cm*sqrt(2), 120cm is reference radius of the HRS, 60cm is the amound by which the map has been enlarge to include fringe fields); X coordinate scaling; Y; Z
HEADER_2 ! Title
181 47 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
TOSCAMOD22.dat !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
180.13967 0. 180.13967 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'DRIFT'
20.0
'MATRIX'
1 0 !order of the matrix/map; 0:means calculate the matrix here
'FIT'
1 !Number of physical parameters to be varied
7 1 2.001 0.5 !element#1 (DRIFT); parameter #1; correlation set so that RE=RS (which is parameters#72 of element #2); allowed variation +/- 50%
1 1e-12 100 !Number of constraints; Convergence threshold; maximum number of iterations
1 1 2 8 0.0 1. 0 !transfer matrix; matrix element 1 2; FIT after element#8 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
! 'FAISTORE'
! zgoubi.fai
! 1
'END'
This diff is collapsed.
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