Commit cce4b0f5 authored by Thomas Planche's avatar Thomas Planche

added details about how I tested MAP2D-E and set up the multipole field map

parent ed3df184
Pipeline #6631 canceled with stages
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set xlabel 'Z [cm]'
set ylabel 'E [MV/cm??]'
set grid
titletext="Testing EMULT vs MAP2D-E"
labeltext="(c) `whoami`, "
today="`date +%Y/%b/%d`"
set label labeltext.today at screen .01, screen .02
set title titletext
#set xrange[-22:22]
set output "plots/Ez_vs_z.eps"
plot "zgoubi.plt" u 22:37 w lp lc 1
set output "plots/Ex_vs_z.eps"
plot "zgoubi.plt" u 22:38 w lp lc 1
set output "plots/Ey_vs_z.eps"
plot "zgoubi.plt" u 22:39 w lp lc 1
** TEST MULTIPOLE
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1 !1:multi particles on a grid
1 1 1 1 1 1
0.3 1.0 0.3 0.0 0.0 1.
0.3 0.0 0.3 0.0 0.0 1.
'PARTICUL' 2
22172.3 1.602176487E-19 0.0 0.0 0.0
! 'OPTIONS'
! 1 1
! CONSTY ON
'MAP2D-E' 5
1 2 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
0.00001 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
45 27 !Number of longitudinal and horizontal-transverse nodes of the mesh (the Z elevation is arbitrary)
test.map !file name
0 0. 0. 0.
4
1.
1 0. 0. 0.
'DRIFT'
20.0
'FAISTORE' 3
zgoubi_map.fai
1
'END'
** TEST MULTIPOLE
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1 !1:multi particles on a grid
11 1 11 1 1 1
0.3 1.0 0.3 0.0 0.0 1.
0.0 0.0 0.0 0.0 0.0 1.
'PARTICUL' 2
22172.3 1.602176487E-19 0.0 0.0 0.0
'ELMULT' 10
2 !IL=7 for output in zgoubi.impdev.out
10. 3. 0. 1000. 20000. 0. 0. 0. 0. 0. 0. 0. ! Length[cm]; tip radius[cm]; dipole tip field[V/m]; quad[V/m]; sextu[V/m]; etc...
20. 10. 1. 1. 1. 1. 1. 1. 1. 1. 1. ! ENTRANCE FACE: Integration zone[cm]; Dipole fringe field extend (FFE)[cm]; quad FFE relative to dipole FFE [unitless]; sextu FFE/dipole FFE [unitless]; etc...
4 0. 1.8 0. 0. 0. 0. !unused; C0...C5
20. 10. 1. 1. 1. 1. 1. 1. 1. 1. 1. ! EXIT FACE: Integration zone[cm]; Dipole fringe field extend (FFE)[cm]; quad FFE relative to dipole FFE [unitless]; sextu FFE/dipole FFE [unitless]; etc...
4 0. 1.8 0. 0. 0. 0. !unused; C0...C5
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. !skew angles[rad] of each component
1.0 !Integration step size[cm]
1 0. 0. 0. !KPOS (1=no misalagnment ); XCE[cm]; YCE[cm]; ALE[rad]
'DRIFT'
40.0
'FAISTORE' 3
zgoubi_emult.fai
1
'END'
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set xlabel 'X [cm]'
set ylabel 'E [MV/cm??]'
set grid
titletext="Testing EMULT vs MAP2D-E"
labeltext="(c) `whoami`, "
today="`date +%Y/%b/%d`"
set label labeltext.today at screen .01, screen .02
set title titletext
set key bot left
set output "plots/Ex_vs_x.eps"
plot "ELMULT.plt" u 22:37 w p lc 1 ,\
"MAP2D.plt" u 22:37 w p lc 3
set key top right
set output "plots/Ey_vs_x.eps"
plot "ELMULT.plt" u 22:38 w p lc 1 ,\
"MAP2D.plt" u 22:38 w p lc 3
set output "plots/Ez_vs_x.eps"
plot "ELMULT.plt" u 22:39 w p lc 1 ,\
"MAP2D.plt" u 22:39 w p lc 3
** TEST MULTIPOLE
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1 !1:multi particles on a grid
11 1 11 1 1 1
0.3 1.0 0.3 0.0 0.0 1.
0.0 0.0 0.0 0.0 0.0 1.
'PARTICUL' 2
22172.3 1.602176487E-19 0.0 0.0 0.0
! 'OPTIONS'
! 1 1
! CONSTY ON
'MAP2D-E' 5
1 2 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1. 1. 1.
HEADER_0
51 15
impdev2EMap.out
0 0. 0. 0.
2
1.
1 0. 0. 0.
'DRIFT'
20.0
'FAISTORE' 3
zgoubi_map.fai
1
'END'
** HSR no quad using DIPOLE
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1 !1:multi particles on a grid
15 1 1 1 1 1
0.3 1.0 0.0 0.0 0.0 1.
0.0 0.0 0.0 0.0 0.0 1.
'PARTICUL' 2 2
22172.3 1.602176487E-19 0.0 0.0 0.0
'OPTIONS'
1 1
CONSTY ON
'ELMULT' 10
7 !IL=7 for output in zgoubi.impdev.out
10. 3. 0. 1000. 20000. 0. 0. 0. 0. 0. 0. 0. ! Length[cm]; tip radius[cm]; dipole tip field[V/m]; quad[V/m]; sextu[V/m]; etc...
20. 10. 1. 1. 1. 1. 1. 1. 1. 1. 1. ! ENTRANCE FACE: Integration zone[cm]; Dipole fringe field extend (FFE)[cm]; quad FFE relative to dipole FFE [unitless]; sextu FFE/dipole FFE [unitless]; etc...
4 0. 1.8 0. 0. 0. 0. !unused; C0...C5
20. 10. 1. 1. 1. 1. 1. 1. 1. 1. 1. ! EXIT FACE: Integration zone[cm]; Dipole fringe field extend (FFE)[cm]; quad FFE relative to dipole FFE [unitless]; sextu FFE/dipole FFE [unitless]; etc...
4 0. 1.8 0. 0. 0. 0. !unused; C0...C5
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. !skew angles[rad] of each component
1.0 !Integration step size[cm]
1 0. 0. 0. !KPOS (1=no misalagnment ); XCE[cm]; YCE[cm]; ALE[rad]
! 'FAISTORE' 3
! zgoubi.fai
! 1
'END'
This diff is collapsed.
This diff is collapsed.
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 "plots/thinBeam.eps"
set xlabel 'Y [cm]'
set ylabel 'Y` [mrad]'
set grid
titletext="Testing EMULT"
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[-80:80]
#set yrange[-60:60]
plot "zgoubi_emult.fai" u 10:11 with points,\
"zgoubi_map.fai" u 10:11 with points lc 3
set output "plots/thinBeam_z.eps"
set xlabel 'Z [cm]'
set ylabel 'Z` [mrad]'
plot "zgoubi_emult.fai" u 12:13 with points,\
"zgoubi_map.fai" u 12:13 with points lc 3
$const #x_step 10.
$const #z_step 10.
$const #x_width 130.*2
$const #z_width 220*2
GRID X0=-#x_width/2. Y0=0 Z0=-#z_width/2. DXG=#x_step DYG=1 DZG=#z_step NXG=INT(#x_width/#x_step)+1 NYG=1 NZG=INT(#z_width/#z_step)+1 FILE='/home/tplanche/text/designs/HRS/opera/multipole/from_Carla/test.table' BINARY=NO FORMAT=2 F1=X UNIT1G=LENGU F2=Y UNIT2G=LENGU F3=Z UNIT3G=LENGU F4=Ex UNIT4G=1 F5=Ey UNIT5G=1 F6=Ez UNIT6G=1 F7=INT((X+#x_width/2.)/#x_step) UNIT7G=1 F8=1 UNIT8G=1 F9=INT((Z+#z_width/2.)/#z_step) UNIT9G=1 UNIT10G=1 UNIT11G=1 UNIT12G=1
y
$const #a INT(#z_width/#z_step)+1
$const #a INT(#x_width/#x_step)+1
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
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 "plots/thinBeam.eps"
set xlabel 'Y [cm]'
set ylabel 'Y` [mrad]'
set grid
titletext="Testing EMULT"
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[-80:80]
#set yrange[-60:60]
plot "zgoubi_emult.fai" u 10:11 with points,\
"zgoubi_map.fai" u 10:11 with points lc 3
set output "plots/thinBeam_z.eps"
set xlabel 'Z [cm]'
set ylabel 'Z` [mrad]'
plot "zgoubi_emult.fai" u 12:13 with points,\
"zgoubi_map.fai" u 12:13 with points lc 3
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