Commit f2f64fc7 authored by Thomas Planche's avatar Thomas Planche

TOSCA+MAP2D-E is now working

parent 7393f29b
This diff is collapsed.
../TOSCA/TOSCAMOD22.map
\ No newline at end of file
** HSR using TOSCA map
'OBJET'
543.687 !543.687 or 543.701 and not 543.714, since we are aiming for a mass resolution of 20000 and not momentum resolution of 20000
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
'PARTICUL'
22172.3 1.602176487E-19 0.0 0.0 0.0
'DRIFT'
21.516804
'TOSCA'
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2 (no print = 0)
1.0318 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 93 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
TOSCAMOD22.map !File name
0 0 0 0. ! ID (see manual);
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.2 !Integration step size[cm]
2 !KPOS, normally=2
180.12483 0. 180.12483 0. !RE fitted to be equal to RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'MAP2D-E'
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
0.00000001 0.1 0.1!Normalization factor for E (to be put in MV/cm), X (cm), and Y(cm)
HEADER_0 !this actually only works if you have n=0 header lines. It crashes if n>0 and I could not figure out why
41 27 !Number of longitudinal and horizontal-transverse nodes of the mesh (the Z elevation is arbitrary)
Multipole2D.map !file name
0 0. 0. 0.
4
1.
1 0. 0. 0.
'TOSCA'
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2 (no print = 0)
1.0318 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 93 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
TOSCAMOD22.map !File name
0 0 0 0. ! ID (see manual);
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.2 !Integration step size[cm]
2 !KPOS, normally=2
180.12483 0. 180.12483 0. !RE fitted to be equal to RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'DRIFT'
21.516804
'FAISCNL'
twomasses.fai
'REBELOTE'
1 0.1 0 1
1
1 1 543.701
'END'
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
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="HRS 'pure separetor'"
plot "zgoubi.fai" u 10:11 with points lc 0 pt 5 ps 0.5 title "Multipole OFF"
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
cm2mm=10.0
set key at -0.15, 5
set xrange [-1.:1.0]
set yrange [-60:60]
plot "twomasses.fai" u ($10*cm2mm):($37==543.687 ? $11:1/0) w points pt 7 ps 0.2 lc 1 notitle, \
"" w points pt 7 ps 0.7 lc 1 title "^{238}U mass ", \
"" u ($10*cm2mm):($37==543.701 ? $11:1/0) w points pt 7 ps 0.2 lc 2 notitle ,\
"" w points pt 7 ps 0.7 lc 2 title "^{238}U mass*20001/20000"
../TOSCA-Vshaped-map/voltages.pdf
\ No newline at end of file
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