Skip to content

GitLab

Commit 659db7fe authored by Thomas Planche's avatar Thomas Planche
Browse files

cleaning up for Dan

parent e2de0f76
Hide whitespace changes
Inline Side-by-side
examples/TRIUMF-HRS/TOSCA/FIT_DRIFT/FIT_DRIFT.in
View file @ 659db7fe
......@@ -7,14 +7,14 @@
0.0 0. 0. 0. 0. 1.
'DRIFT'
20.0
20.1
'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
../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]
......@@ -33,7 +33,7 @@ TOSCAMOD22.dat !File name
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
../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]
......@@ -42,7 +42,7 @@ TOSCAMOD22.dat !File name
'DRIFT'
20.0
20.1
'MATRIX'
1 0 !order of the matrix/map; 0:means calculate the matrix here
......@@ -53,8 +53,11 @@ TOSCAMOD22.dat !File name
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
'END'
! 'FAISTORE'
! zgoubi.fai
! 1
'END'
examples/TRIUMF-HRS/TOSCA/FIT_TOSCA/tempRes deleted 100644 → 0
View file @ e2de0f76
**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.
Skip reading field map file : ../TOSCAMOD22.dat
Pgm toscap, restored mesh coordinates for field map # 1, name : ../TOSCAMOD22.dat
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 = 180.140 cm TE = 0.00000 rad
at exit RS = 180.140 cm TS = 0.00000 rad
A 1 1.0000 0.000 0.000 0.000 0.000 0.785 180.140 -0.000 0.000 0.000 1
Cumulative length of optical axis = 0.00000000 m ; Time (for ref. rigidity & particle) = 0.00000 s
************************************************************************************************************************************
3 Keyword FIT[2] is skipped since this is the (end of) last run following the fitting procedure.
Now carrying on beyond FIT keyword.
************************************************************************************************************************************
4 Keyword, label(s) : END
1 particles have been launched
Made it to the end : 1
************************************************************************************************************************************
ZGOUBI RUN COMPLETED.
Zgoubi, author's dvlpmnt version.
Job started on 08-08-2017, at 16:06:17
JOB ENDED ON 08-08-2017, AT 16:06:18
CPU time, total : 0.23859899999999998
examples/TRIUMF-HRS/TOSCA/FIT_TOSCA/zgoubi.res deleted 100644 → 0
View file @ e2de0f76
**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.
Skip reading field map file : ../TOSCAMOD22.dat
Pgm toscap, restored mesh coordinates for field map # 1, name : ../TOSCAMOD22.dat
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 = 180.140 cm TE = 0.00000 rad
at exit RS = 180.140 cm TS = 0.00000 rad
A 1 1.0000 0.000 0.000 0.000 0.000 0.785 180.140 -0.000 0.000 0.000 1
Cumulative length of optical axis = 0.00000000 m ; Time (for ref. rigidity & particle) = 0.00000 s
************************************************************************************************************************************
3 Keyword FIT[2] is skipped since this is the (end of) last run following the fitting procedure.
Now carrying on beyond FIT keyword.
************************************************************************************************************************************
4 Keyword, label(s) : END
1 particles have been launched
Made it to the end : 1
************************************************************************************************************************************
ZGOUBI RUN COMPLETED.
Zgoubi, author's dvlpmnt version.
Job started on 08-08-2017, at 16:06:17
JOB ENDED ON 08-08-2017, AT 16:06:18
CPU time, total : 0.23859899999999998
An updated version of the input data file, with variables in FIT'ed state, has been saved in zgoubi.FIT.out.dat.
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