Commit 11af7bb1 authored by Thomas Planche's avatar Thomas Planche

layout plot works!

parent 94fd4554
......@@ -6,6 +6,7 @@
*.fai
*.plt
*.ps
*.res
# logs and backup files #
#########################
......@@ -14,3 +15,4 @@
# Other stuffs#
#########################
tempRes
# COORDINATES - STORAGE FILE, 09-08-2017 08:17:14. ** HSR using TOSCA, FIT drift length
# 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 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 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.63746041E+02 1.63746041E+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.46198000E+02 1 1 'OBJET ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.2896661909804306E+00 6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 2.1507338596905175E+01 0.0000000000000000E+00 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.63746041E+02 1.63746041E+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.46198000E+02 1 3 'DRIFT ' ' ' ' ' ''
# COORDINATES - STORAGE FILE, 09-08-2017 08:17:14. ** HSR using TOSCA, FIT drift length
# 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 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.3404629803896739E+00 -5.9984720697026646E+01 0.0000000000000000E+00 0.0000000000000000E+00 7.0902380160285850E+02 2.1598824295367084E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 11 'TOSCA ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 5.1125998649097060E-02 -5.9984720697026646E+01 0.0000000000000000E+00 0.0000000000000000E+00 7.3053112046156923E+02 2.1598824295367084E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 13 'DRIFT ' ' ' ' ' ''
This diff is collapsed.
This diff is collapsed.
# COORDINATES - STORAGE FILE, 09-08-2017 08:17:14. ** HSR using TOSCA, FIT drift length
# 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 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 9.8426931356122509E+00 1.5233413466594317E-01 0.0000000000000000E+00 0.0000000000000000E+00 3.4522695176766820E+02 1.0798123986520483E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 5 'TOSCA ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 9.8457398183291360E+00 1.5233413466594317E-01 0.0000000000000000E+00 0.0000000000000000E+00 3.6522695199972509E+02 1.0798123986520483E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 7 'DRIFT ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 9.8487865010460212E+00 1.5233413466594317E-01 0.0000000000000000E+00 0.0000000000000000E+00 3.8522695223178198E+02 1.0798123986520483E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 9 'DRIFT ' ' ' ' ' 'â'
# COORDINATES - STORAGE FILE, 09-08-2017 09:19:26. ** HSR using TOSCA, FIT drift length
# 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 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 -6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 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.63746041E+02 1.63746041E+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.46198000E+02 1 1 'OBJET ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -1.2896661909804306E+00 -6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 2.1507338596905175E+01 0.0000000000000000E+00 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.63746041E+02 1.63746041E+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.46198000E+02 1 3 'DRIFT ' ' ' ' ' ''
# COORDINATES - STORAGE FILE, 09-08-2017 09:19:27. ** HSR using TOSCA, FIT drift length
# 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 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -1.2867311074642487E+00 5.9999799092314603E+01 0.0000000000000000E+00 0.0000000000000000E+00 6.5066552426924420E+02 1.9652201713180419E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 11 'TOSCA ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 2.9307547371528742E-03 5.9999799092314603E+01 0.0000000000000000E+00 0.0000000000000000E+00 6.7217286260657886E+02 1.9652201713180419E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 13 'DRIFT ' ' ' ' ' ''
This diff is collapsed.
This diff is collapsed.
# COORDINATES - STORAGE FILE, 09-08-2017 09:19:27. ** HSR using TOSCA, FIT drift length
# 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 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -9.5300459081509814E+00 9.2447027313427299E-03 0.0000000000000000E+00 0.0000000000000000E+00 3.1608430372047633E+02 9.8260298837665929E-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.63746041E+02 1.63746041E+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.46198000E+02 1 5 'TOSCA ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -9.5298610140963493E+00 9.2447027313427299E-03 0.0000000000000000E+00 0.0000000000000000E+00 3.3608430372133097E+02 9.8260298837665929E-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.63746041E+02 1.63746041E+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.46198000E+02 1 7 'DRIFT ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -9.5296761200417173E+00 9.2447027313427299E-03 0.0000000000000000E+00 0.0000000000000000E+00 3.5608430372218561E+02 9.8260298837665929E-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.63746041E+02 1.63746041E+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.46198000E+02 1 9 'DRIFT ' ' ' ' ' 'â'
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "layout.eps"
set xlabel 'X [m]'
set ylabel 'Y [m]'
set grid
#titletext="HRS TOSCA -- particle trajectories"
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
initer(x) = ((!exists("first")) ? first = x : first, first=NaN) ; #just a hack to get the first value of a column. Not very elegent but will do. I added the first=NaN to reset "first" so that initer() can be used multiple times
set size ratio -1
dl=21.468637/100.0 #length of the first/last drift
dl2=20/100. #drift length between the two TOSCA maps
RE=180.13967/100 #fitted RE (=RS) value in 'TOSCA'
plot "ref_entrance.fai" u (-$10/100.):($14/100.) w l lt 1 lc 1 notitle,\
"ref_mag1.plt" using ($10/100.*sin($22-pi/4.)+RE):($10/100.*cos($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
"ref_mid.fai" u (($14-initer($14))/100.+RE):($10/100.+RE+dl) w l lt 1 lc 1 notitle,\
"ref_mag2.plt" using ($10/100.*cos($22-pi/4.)+RE):(-$10/100.*sin($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
"a+60mrad_entrance.fai" u (-$10/100.):($14/100.) w l lt 1 lc 1 notitle,\
"a+60mrad_mag1.plt" using ($10/100.*sin($22-pi/4.)+RE):($10/100.*cos($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
"a+60mrad_mid.fai" u (($14-initer($14))/100.+RE):($10/100.+RE+dl) w l lt 1 lc 1 notitle,\
"a+60mrad_mag2.plt" using ($10/100.*cos($22-pi/4.)+RE):(-$10/100.*sin($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
# r0=1.2013 #reference bending radius
# dle=0.255 #drift length before entering into the TOSCA map
# dlmap=r0*tan(25*pi/180.) #length of the drift inside the map before/after the effective edge
# dlc=0.24043 #drift length between the exit of the TOSCA map and the center of the multipole
# Sdr1=3.26067388585 #value of S at the entrance of the drift after the first dipole (copied from center_ref.fai)
# Sdr2=6.74719321425 #value of S at the entrance of the drift after the second dipole (copied from exit_ref.fai)
#set yrange [-0.05:0.05]
#set xrange [-0.05:0.05]
# plot "entrance_ref.fai" u (-$10/100.):($14/100.) w l lt 1 lc 1 notitle,\
# "mag1_ref.plt" using ($10/100.*sin($22-pi/4.)+r0):($10/100.*cos($22-pi/4.)+dlmap+dle) with l lt 1 lc 1 notitle,\
# "mag2_ref.plt" using ($10/100.*sin($22+pi/4.)+2*dlmap+r0+2*dlc):($10/100.*cos($22+pi/4.)+dlmap+dle) with l lt 1 lc 1 notitle,\
# "entrance+60.fai" u (-$10/100.):($14/100.) w l lt 1 lc 3 notitle,\
# "mag1+60.plt" using ($10/100.*sin($22-pi/4.)+r0):($10/100.*cos($22-pi/4.)+dlmap+dle) with l lt 1 lc 3 notitle,\
# "mag2+60.plt" using ($10/100.*sin($22+pi/4.)+2*dlmap+r0+2*dlc):($10/100.*cos($22+pi/4.)+dlmap+dle) with l lt 1 lc 3 notitle,\
#"exit_ref.fai" u ($10/100.):($14/100.-Sdr2) w l lt 1 lc 3 notitle,\
This diff is collapsed.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "layout.eps"
set xlabel '[m]'
set ylabel '[m]'
set grid
titletext="HRS layout (using TOSCA map)"
labeltext="(c) `whoami`, "
today="`date +%Y/%b/%d`"
set label labeltext.today at screen .01, screen .02
set title titletext
set key cent cent
cm2mm=10.0
set yrange [0:2.5]
findFirst(x) = ($0 > 0 ? first : first = x) #just a hack to get the first value of a column.
set size ratio -1
dl=21.468637/100.0 #length of the first/last drift
mcl=40./100. #multipole corrector length
RE=180.13967/100 #fitted RE (=RS) value in 'TOSCA'
plot "ref_entrance.fai" u (-$10/100.):($14/100.) w l lt 2 lc 1 title "reference trajectory",\
"ref_mag1.plt" using ($10/100.*sin($22-pi/4.)+RE):($10/100.*cos($22-pi/4.)+dl) with l lt 2 lc 1 notitle,\
"ref_mid.fai" u (findFirst($14),($14-first)/100.+RE):($10/100.+RE+dl) w l lt 2 lc 1 notitle,\
"ref_mag2.plt" using ($10/100.*cos($22-pi/4.)+RE+mcl):(-$10/100.*sin($22-pi/4.)+dl) with l lt 2 lc 1 notitle,\
"ref_exit.fai" u ($10/100.+2*RE+mcl):(findFirst($14),-($14-first)/100.+dl) w l lt 2 lc 1 notitle,\
"a+60mrad_entrance.fai" u (-$10/100.):($14/100.) w l lt 1 lc 1 title "+/- 60 mrad",\
"a+60mrad_mag1.plt" using ($10/100.*sin($22-pi/4.)+RE):($10/100.*cos($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
"a+60mrad_mid.fai" u (findFirst($14),($14-first)/100.+RE):($10/100.+RE+dl) w l lt 1 lc 1 notitle,\
"a+60mrad_mag2.plt" using ($10/100.*cos($22-pi/4.)+RE+mcl):(-$10/100.*sin($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
"a+60mrad_exit.fai" u ($10/100.+2*RE+mcl):(findFirst($14),-($14-first)/100.+dl) w l lt 1 lc 1 notitle,\
"a-60mrad_entrance.fai" u (-$10/100.):($14/100.) w l lt 1 lc 1 notitle,\
"a-60mrad_mag1.plt" using ($10/100.*sin($22-pi/4.)+RE):($10/100.*cos($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
"a-60mrad_mid.fai" u (findFirst($14),($14-first)/100.+RE):($10/100.+RE+dl) w l lt 1 lc 1 notitle,\
"a-60mrad_mag2.plt" using ($10/100.*cos($22-pi/4.)+RE+mcl):(-$10/100.*sin($22-pi/4.)+dl) with l lt 1 lc 1 notitle,\
"a-60mrad_exit.fai" u ($10/100.+2*RE+mcl):(findFirst($14),-($14-first)/100.+dl) w l lt 1 lc 1 notitle,\
"mdm100+_entrance.fai" u (-$10/100.):($14/100.) w l lt 1 lc 3 title "particle mass +1%",\
"mdm100+_mag1.plt" using ($10/100.*sin($22-pi/4.)+RE):($10/100.*cos($22-pi/4.)+dl) with l lt 1 lc 3 notitle,\
"mdm100+_mid.fai" u (findFirst($14),($14-first)/100.+RE):($10/100.+RE+dl) w l lt 1 lc 3 notitle,\
"mdm100+_mag2.plt" using ($10/100.*cos($22-pi/4.)+RE+mcl):(-$10/100.*sin($22-pi/4.)+dl) with l lt 1 lc 3 notitle,\
"mdm100+_exit.fai" u ($10/100.+2*RE+mcl):(findFirst($14),-($14-first)/100.+dl) w l lt 1 lc 3 notitle,\
"mdm100-_entrance.fai" u (-$10/100.):($14/100.) w l lt 1 lc 3 notitle,\
"mdm100-_mag1.plt" using ($10/100.*sin($22-pi/4.)+RE):($10/100.*cos($22-pi/4.)+dl) with l lt 1 lc 3 notitle,\
"mdm100-_mid.fai" u (findFirst($14),($14-first)/100.+RE):($10/100.+RE+dl) w l lt 1 lc 3 notitle,\
"mdm100-_mag2.plt" using ($10/100.*cos($22-pi/4.)+RE+mcl):(-$10/100.*sin($22-pi/4.)+dl) with l lt 1 lc 3 notitle,\
"mdm100-_exit.fai" u ($10/100.+2*RE+mcl):(findFirst($14),-($14-first)/100.+dl) w l lt 1 lc 3 notitle,\
** HSR using TOSCA, FIT drift length
'OBJET'
551.660 !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. 60. 0. 0. 0. 1. ’o’ !Y; T; Z; P; S; D; 'marker'
1 !1 or -9 (-9 disables the tracking of this particle)
'FAISTORE'
mdm100+_entrance.fai
1
'DRIFT'
21.468637
'FAISTORE'
mdm100+_entrance.fai
1
'TOSCA'
0 2 !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
'FAISTORE'
mdm100+_mid.fai
1
'DRIFT'
20.
'FAISTORE'
mdm100+_mid.fai
1
'DRIFT'
20.
'FAISTORE'
mdm100+_mid.fai
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
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
'FAISTORE'
mdm100+_exit.fai
1
'DRIFT'
21.468637
'FAISTORE'
mdm100+_exit.fai
1
'END'
# COORDINATES - STORAGE FILE, 09-08-2017 09:30:10. ** HSR using TOSCA, FIT drift length
# 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 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 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.65383507E+02 1.65383507E+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.51660000E+02 1 1 'OBJET ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.2896661909804306E+00 6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 2.1507338596905175E+01 0.0000000000000000E+00 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.65383507E+02 1.65383507E+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.51660000E+02 1 3 'DRIFT ' ' ' ' ' ''
# COORDINATES - STORAGE FILE, 09-08-2017 09:30:11. ** HSR using TOSCA, FIT drift length
# 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 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 6.7909863082628590E+00 -5.7507889339043111E+01 0.0000000000000000E+00 0.0000000000000000E+00 7.1875264052571049E+02 2.1923173342056904E-02 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.65383507E+02 1.65383507E+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.51660000E+02 1 11 'TOSCA ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 5.5550074814248456E+00 -5.7507889339043111E+01 0.0000000000000000E+00 0.0000000000000000E+00 7.4025682659024596E+02 2.1923173342056904E-02 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.65383507E+02 1.65383507E+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.51660000E+02 1 13 'DRIFT ' ' ' ' ' ''
This diff is collapsed.
This diff is collapsed.
# COORDINATES - STORAGE FILE, 09-08-2017 09:30:11. ** HSR using TOSCA, FIT drift length
# 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 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.2141391481069121E+01 1.5967231328486704E+01 0.0000000000000000E+00 0.0000000000000000E+00 3.4599855278287527E+02 1.0823861825969634E-02 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.65383507E+02 1.65383507E+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.51660000E+02 1 5 'TOSCA ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.2460763249641273E+01 1.5967231328486704E+01 0.0000000000000000E+00 0.0000000000000000E+00 3.6600110257850287E+02 1.0823861825969634E-02 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.65383507E+02 1.65383507E+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.51660000E+02 1 7 'DRIFT ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 1.2780135018213425E+01 1.5967231328486704E+01 0.0000000000000000E+00 0.0000000000000000E+00 3.8600365237413047E+02 1.0823861825969634E-02 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.65383507E+02 1.65383507E+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.51660000E+02 1 9 'DRIFT ' ' ' ' ' 'â'
# COORDINATES - STORAGE FILE, 09-08-2017 09:28:28. ** HSR using TOSCA, FIT drift length
# 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 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 -6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 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.65383507E+02 1.65383507E+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.51660000E+02 1 1 'OBJET ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -1.2896661909804306E+00 -6.0000000000000000E+01 0.0000000000000000E+00 0.0000000000000000E+00 2.1507338596905175E+01 0.0000000000000000E+00 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.65383507E+02 1.65383507E+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.51660000E+02 1 3 'DRIFT ' ' ' ' ' ''
# COORDINATES - STORAGE FILE, 09-08-2017 09:28:29. ** HSR using TOSCA, FIT drift length
# 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 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 3.2272391164710825E+00 6.0022233761048845E+01 0.0000000000000000E+00 0.0000000000000000E+00 6.5843975812041799E+02 1.9911391041960862E-02 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.65383507E+02 1.65383507E+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.51660000E+02 1 11 'TOSCA ' ' ' ' ' ''
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 4.5173843591518397E+00 6.0022233761048845E+01 0.0000000000000000E+00 0.0000000000000000E+00 6.7994712544849835E+02 1.9911391041960862E-02 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.65383507E+02 1.65383507E+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.51660000E+02 1 13 'DRIFT ' ' ' ' ' ''
This diff is collapsed.
This diff is collapsed.
# COORDINATES - STORAGE FILE, 09-08-2017 09:28:28. ** HSR using TOSCA, FIT drift length
# 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 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -7.6307655092656717E+00 1.4023078052438061E+01 0.0000000000000000E+00 0.0000000000000000E+00 3.1670804348909496E+02 9.8468356029224290E-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.65383507E+02 1.65383507E+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.51660000E+02 1 5 'TOSCA ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -7.3502855628222301E+00 1.4023078052438061E+01 0.0000000000000000E+00 0.0000000000000000E+00 3.3671001011741322E+02 9.8468356029224290E-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.65383507E+02 1.65383507E+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.51660000E+02 1 7 'DRIFT ' ' ' ' ' 'â'
1 0.00000000E+00 0.00000000E+00 -6.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.0000000000000000E+00 -7.0698056163787886E+00 1.4023078052438061E+01 0.0000000000000000E+00 0.0000000000000000E+00 3.5671197674573148E+02 9.8468356029224290E-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.65383507E+02 1.65383507E+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.51660000E+02 1 9 'DRIFT ' ' ' ' ' 'â'
# COORDINATES - STORAGE FILE, 08-08-2017 16:35:09. ** HSR using TOSCA, FIT drift length
# 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 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 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.63746041E+02 1.63746041E+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.46198000E+02 1 1 'OBJET ' ' ' ' ' ''
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 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 0.0000000000000000E+00 2.1468637000000001E+01 0.0000000000000000E+00 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.63746041E+02 1.63746041E+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.46198000E+02 1 3 'DRIFT ' ' ' ' ' ''
# COORDINATES - STORAGE FILE, 08-08-2017 16:35:10. ** HSR using TOSCA, FIT drift length
# 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 -3.3846372929247082E-05 -2.8472542925152625E-04 0.0000000000000000E+00 0.0000000000000000E+00 6.7852159174405188E+02 2.0582671053854040E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 11 'TOSCA ' ' ' ' ' ''
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 -3.9959039814517448E-05 -2.8472542925152625E-04 0.0000000000000000E+00 0.0000000000000000E+00 6.9999022874405273E+02 2.0582671053854040E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 13 'DRIFT ' ' ' ' ' ''
This diff is collapsed.
This diff is collapsed.
# COORDINATES - STORAGE FILE, 08-08-2017 16:35:09. ** HSR using TOSCA, FIT drift length
# 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 8.2426950882563688E-06 -2.6382734930597532E-04 0.0000000000000000E+00 0.0000000000000000E+00 3.2999514381570276E+02 1.0291336509063031E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 5 'TOSCA ' ' ' ' ' 'â'
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 2.9661481021367396E-06 -2.6382734930597532E-04 0.0000000000000000E+00 0.0000000000000000E+00 3.4999514381570344E+02 1.0291336509063031E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 7 'DRIFT ' ' ' ' ' 'â'
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 -2.3103988839828896E-06 -2.6382734930597532E-04 0.0000000000000000E+00 0.0000000000000000E+00 3.6999514381570413E+02 1.0291336509063031E-02 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.63746041E+02 1.63746041E+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.46198000E+02 1 9 'DRIFT ' ' ' ' ' 'â'
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.
** HSR using TOSCA, FIT drift length
'OBJET' 1
551.660 !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. 60. 0. 0. 0. 1. o !Y; T; Z; P; S; D; 'marker'
1 !1 or -9 (-9 disables the tracking of this particle)
'FAISTORE' 2
mdm100+_entrance.fai
1
'DRIFT' 3
21.468637
'FAISTORE' 4
mdm100+_entrance.fai
1
'TOSCA' 5
0 2 !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
'FAISTORE' 6
mdm100+_mid.fai
1
'DRIFT' 7
20.
'FAISTORE' 8
mdm100+_mid.fai
1
'DRIFT' 9
20.
'FAISTORE' 10
mdm100+_mid.fai
1
'TOSCA' 11
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
'FAISTORE' 12
mdm100+_exit.fai
1
'DRIFT' 13
21.468637
'FAISTORE' 14
mdm100+_exit.fai
1
'END' 15
************************************************************************************************************************************
1 Keyword, label(s) : OBJET
MAGNETIC RIGIDITY = 551.660 kG*cm
TRAJECTOIRY SETTING UP
OBJET (2) BUILT UP FROM 1 POINTS
************************************************************************************************************************************
2 Keyword, label(s) : FAISTORE
OPEN FILE mdm100+_entrance.fai
FOR PRINTING COORDINATES
Print will occur at element[s] labeled :
************************************************************************************************************************************
3 Keyword, label(s) : DRIFT
Drift, length = 21.46864 cm
TRAJ #1 IEX,D,Y,T,Z,P,S,time : 1 0.000000E+00 1.289666E+00 6.000000E+01 0.000000E+00 0.000000E+00 2.1507339E+01 0.00000E+00
Cumulative length of optical axis = 0.214686370 m ; Time (for reference rigidity & particle) = 4.125382E-09 s
************************************************************************************************************************************
4 Keyword, label(s) : FAISTORE
mdm100+_entrance.fai
already open...
Print will occur at element[s] labeled :
************************************************************************************************************************************
5 Keyword, label(s) : TOSCA
OPEN FILE zgoubi.plt
FOR PRINTING TRAJECTORIES
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 = 180.140 cm TE = 0.00000 rad
at exit RS = 180.140 cm TS = 0.00000 rad
A 1 1.0000 0.000 60.000 0.000 0.000 0.785 192.281 0.016 0.000 0.000 1
CONDITIONS DE MAXWELL ( 3245. PAS ) :
DIV(B) LAPLACIEN(B) ROTATIONNEL(B)
0.000 0.000 0.000
0.000 0.000
0.000 0.000
LAPLACIEN SCALAIRE = 0.000
Cumulative length of optical axis = 0.214686370 m ; Time (for ref. rigidity & particle) = 4.125382E-09 s
************************************************************************************************************************************
6 Keyword, label(s) : FAISTORE
OPEN FILE mdm100+_mid.fai
FOR PRINTING COORDINATES
Print will occur at element[s] labeled :
************************************************************************************************************************************
7 Keyword, label(s) : DRIFT
Drift, length = 20.00000 cm
TR