Commit 2b8174a2 authored by Thomas Planche's avatar Thomas Planche

cleaning up...

parent f5297376
Pipeline #6629 canceled with stages
**HRS FIT DIPOLE
'OBJET'
544.12 !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)
'DRIFT'
80.0 !80 cm long drift
'CHANGREF'
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)
'FAISCNL'
checkY.fai
'END'
'DIPOLE'
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'
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.1 !in element #5 (i.e. DIPOLE); vary parameter #6 (i.e. B0); coupling switch OFF=0; relative range +/- 20%
2 1e-10 1 !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'
!!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.
This diff is collapsed.
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "BzThplot.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 notitle
**HRS FIT TOSCA
'OBJET'
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'
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
200.0 0. 200.0 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'FIT'
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'
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "RThplot.eps"
set xlabel 'Theta [deg.]'
set ylabel 'R [m]'
#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
set xrange[-50:50]
set yrange[1.8:2.2]
set key bot cent
plot "Bscale1_zgoubi.plt" using ($22/pi*180.):($10/100.) with p pt 7 ps 0.3 lc 3 title "B scale = 1",\
"zgoubi.plt" using ($22/pi*180.):($10/100.) with p pt 7 ps 0.3 lc 1 title "B scale = 1.0365",\
1.2+0.6*sqrt(2.) w l lw 0.5 lt 2 lc 0 title "R=1.2+0.6*sqrt(2.)"
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "thinBeam.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 right
cm2mm=10.0
plot "zgoubi_R1450.fai" u ($10*cm2mm):11 with points lc 1 title "EdgeR=1.450 m",\
'zgoubi_R2238.fai' u ($10*cm2mm):11 with points lc 3 title "EdgeR=2.2379 m",\
'zgoubi.fai' u ($10*cm2mm):11 with points lc 0 pt 5 ps 0.5 title "2.2379 m, finer mesh"
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "xyplot.eps"
set xlabel 'y [m]'
set ylabel 'x [m]'
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 bot rig
cm2mm=10.0
set xrange[-2:0]
set yrange[0:2]
set sample 2000
set size ratio -1
plot "zgoubi.plt" using ($10/100.*sin($22-pi/4.)):($10/100.*cos($22-pi/4.)) with p pt 7 ps 0.3 title "reference trajectory",\
sqrt(-(x+0.6)**2+1.2**2)+0.6 w l lw 0.5 lt 2 lc 0 title "designed bending radius",\
-sqrt(-(x+0.6)**2+1.2**2)+0.6 w l lw 0.5 lt 2 lc 0 notitle
This diff is collapsed.
**HRS FIT TOSCA
'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' 2
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
200.0 0. 200.0 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'FIT' 3
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' 4
************************************************************************************************************************************
1 Keyword, label(s) : OBJET
MAGNETIC RIGIDITY = 546.198 kG*cm
TRAJECTOIRY SETTING UP
OBJET (2) BUILT UP FROM 1 POINTS
************************************************************************************************************************************
2 Keyword, label(s) : TOSCA
NDIM = 2 ; Number of data file sets used is 1 ; Stored in field array # IMAP = 1 ;
Value of MOD.MOD2 is 22.1
3-D map. MOD=22 or 23. Single field map, with field coefficient value :
1.000000E+00
New field map(s) now used, polar mesh (MOD .ge. 20) ; name(s) of map data file(s) are :
../TOSCAMOD22.dat
----
Map file number 1 ( of 1) :
../TOSCAMOD22.dat map, FORMAT type : regular. Field multiplication factor : 1.00000000E+00
HEADER (2 lines) :
165.852814 1.500000 0.50000 0.000000
Y[cm] Z[cm] X[cm] Br[kG] Bth[kG] Bz[kG]
R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.000000E-01 0.000000E+00
Field map limits, angle : min, max, max-min (rad) : -0.78539816E+00 0.78539816E+00 0.15707963E+01
Field map limits, radius : min, max, max-min (cm) : 0.16585281E+03 0.23485281E+03 0.69000000E+02
Min / max fields drawn from map data : -1.8960 / 4.2372
@ X-node, Y-node, z-node : 0.314 232. 0.00 / 0.777 185. 0.00
Normalisation coeff. BNORM : 1.03650
Field min/max normalised : -1.9652 / 4.3919
Nbre of nodes in X/Y/Z : 181/ 47/ 1
Node distance in X/Y/Z : 8.726646E-03/ 1.50000 / 0.00000
Option for interpolation : 4
Smoothing using 25 point
Integration step : 0.1000 cm (i.e., 4.9912E-04 rad at mean radius RM = 200.4 )
KPOS = 2. Position of reference orbit on mechanical faces
at entrance RE = 200.000 cm TE = 0.00000 rad
at exit RS = 200.000 cm TS = 0.00000 rad
A 1 1.0000 0.000 0.000 0.000 0.000 0.785 182.857 -0.118 0.000 0.000 1
Cumulative length of optical axis = 0.00000000 m ; Time (for ref. rigidity & particle) = 0.00000 s
************************************************************************************************************************************
3 Keyword, label(s) : FIT
FIT procedure launched. Method is 1
variable # 1 IR = 2 , ok.
variable # 1 IP = 70 , ok.
variable element 2, prmtr # 70 : coupled with element 2, prmtr # 72
variable # 1 XC.= 2 , ok.
variable # 1 .XC= 72 , ok.
constraint # 1 IR = 2 , ok.
constraint # 1 I = 1 , ok.
FIT variables and constraints in good order, FIT will proceed.
Final FIT status will NOT be saved. For so, use the 'save [FileName]' command
STATUS OF VARIABLES (Iteration # 0 / 999 max.)
LMNT VAR PARAM MINIMUM INITIAL FINAL MAXIMUM STEP NAME LBL1 LBL2
2 1 70 100. 180. 180.13968 300. 1.129E-05 TOSCA - -
2 1 72 100. 180. 180.13968 300. 1.129E-05
STATUS OF CONSTRAINTS (Target penalty = 1.0000E-10)
TYPE I J LMNT# DESIRED WEIGHT REACHED KI2 NAME LBL1 LBL2 Nb param. [value]
3 1 2 2 0.000000E+00 1.000E+00 -9.469218E-07 1.00E+00 TOSCA - - 0
Fit reached penalty value 8.9666E-13
************************************************************************************************************************************
MAIN PROGRAM : FIT completed. Now doing a last run using variable values from FIT.
************************************************************************************************************************************
1 Keyword, label(s) : OBJET
MAGNETIC RIGIDITY = 546.198 kG*cm
TRAJECTOIRY SETTING UP
OBJET (2) BUILT UP FROM 1 POINTS
************************************************************************************************************************************
2 Keyword, label(s) : TOSCA
NDIM = 2 ; Number of data file sets used is 1 ; Stored in field array # IMAP = 1 ;
Value of MOD.MOD2 is 22.1
3-D map. MOD=22 or 23. Single field map, with field coefficient value :
1.000000E+00
No new map file to be opened. Already stored.