Commit 3b9d12e4 authored by Thomas Planche's avatar Thomas Planche

before asking help to Francois about REBELOTE

parent ebe8d90c
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !544.1192 or 544.1328
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
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.474 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.474 -227.0 0. 0. -227.0 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 127.96721 -0.349067 127.96721 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.474 227.0 0.0 0.0 227.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.474 -227.0 0. 0. -227.0 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 127.96721 -0.349067 127.96721 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
mass1.fai
'REBELOTE'
1 0.1 0 1
1
1 1 544.1328
'END'
This source diff could not be displayed because it is too large. You can view the blob instead.
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1328 !544.1192 or 544.1328 or 544.1464?? mass resolution vs momentum resolution
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. 1. 0.5 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
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 1e99 0.0 0.0 1e99 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -1e99 0. 0. -1e99 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.543694 1e99 0.0 0.0 1e99 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.543694 -1e99 0. 0. -1e99 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 127.96748 -0.349067 127.96748 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
mass2.fai
'END'
This source diff could not be displayed because it is too large. You can view the blob instead.
set term postscript eps color enhanced "Times-Roman" 18 lw 2
set output "twomasses_edgeCurvature.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
set key left cent
cm2mm=10.0
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 1 title "^{238}U mass", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 3 title"20001/20000*^{238}U mass"
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1192 !60 keV U1+
5 !5: generate 11 particles used to calculate transfer matrix ('MATRIX')
0.0001 0.1 0.0001 0.1 .1 .001 !step size in Y; T; Z; P; S; D
0.0 0. 0. 0. 0. 1.
'DRIFT'
36.32357
'DIPOLE'
0 !output flag: 0: no outpot, 2: output trajectory to zgoubi.plt, etc.
130 120. !AT[deg.], RM[cm]
65 4.5343276 0. 0. 0. !ACENT[deg.];B0[kG]; N; B; GX
7.0 0. !ENTRANCE FIELD BOUNDARY: fringe field extend[cm]; unused
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
45. 26.575 220.0 0.0 0.0 220.0 !W+[deg.]; edge angle[deg.]; R1[cm]; U1; U2; R2[cm]
7.0 0. !EXIT FIELD BOUNDARY: lambda; =-1 for Enge like
4 0. 1.8 0. 0. 0. 0. 0. !unused; C0...C5; shift
-45. -26.575 -220.0 0. 0. -220.0 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.1 !Integration step size[cm]
2 127.96721 -0.349067 127.96721 0.349067 !KPOS RE[cm]; TE[rad]; RS[cm]; TS[rad]
'DRIFT'
36.32357
'MATRIX'
1 0
'FIT'
1 !Number of physical parameters to be varied
3 36 -3.020 0.1 !element#4 (DIPOLE); parameter #36 (exit edge angle); inverse (minus sign) correlation with element#4 parameter#20 (entrance edge angle); relative range
2 1e-11 100 !Number of constraints; Convergence threshold; maximum number of iterations
1 1 1 5 0.0 1. 0 !transfer matrix; matrix element 1 1; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
1 2 2 5 0.0 1. 0 !transfer matrix; matrix element 2 2; FIT after element#12 ('MATRIX');Wanted value(=0.0); Weigth ; 0:no additional parameters
'END'
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
set key left cent
cm2mm=10.0
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 1 title "^{238}U mass", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 3 title"20001/20000*^{238}U mass"
zgoubi HRS tuto: DIPOLE
'OBJET'
544.1328 !544.1192 or 544.1328
544.1192 !544.1192 or 544.1328 and not 544.1464, 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. 1. 0.5 1. 1. 1. * !scaling factors on Y,T,Z,P,S,DP
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
......@@ -61,7 +61,7 @@ noEspread.beam ! name of the input file
36.32357 ! drift length=80-120*tan(20/deg.)
'FAISCNL'
mass2.fai
mass1.fai
'END'
......
This diff is collapsed.
......@@ -10,6 +10,6 @@ set key left cent
cm2mm=10.0
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 1 title "^{238}U mass", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.5 lc 3 title"20001/20000*^{238}U mass"
plot "mass1.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 1 title "^{238}U mass", \
"mass2.fai" u ($10*cm2mm):11 w points pt 7 ps 0.2 lc 3 title"20001/20000*^{238}U mass"
** TEST MULTIPOLE
'OBJET' 1
546.198 !60 keV U1+ = 546.198
1 !1:multi particles on a grid
11 1 1 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
! forFrancois.fai
! 1
'END'
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
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 with fringe field"
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 "Ex_vs_x.eps"
plot "zgoubi.plt" u 22:37 w p lc 1
set key top right
set output "Ey_vs_x.eps"
plot "zgoubi.plt" u 22:38 w p lc 1
set output "Ez_vs_x.eps"
plot "zgoubi.plt" u 22:39 w p lc 1
@article{MARCHETTO2019,
title = "Status of the CANREB high-resolution separator at TRIUMF",
journal = "Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms",
year = "2019",
issn = "0168-583X",
doi = "https://doi.org/10.1016/j.nimb.2019.05.032",
url = "http://www.sciencedirect.com/science/article/pii/S0168583X19303349",
author = "M. Marchetto and F. Ames and R. Baartman and C. Barquest and S. Kiy and T. Planche and S. Saminathan and J.A. Maloney and S. Brown and M. Corwin and D. Sehayek and J. Laroche",
keywords = "High-resolution separator, Mulitpole, CANREB, ARIEL, TRIUMF",
abstract = "A new ISOL rare isotope beam production facility, ARIEL, is under construction at TRIUMF. ARIEL aims at increasing the delivery of radioactive beams three fold with respect to the present capability of the ISAC facility. Part of ARIEL is the new CANREB equipment that can be described by the two main functionalities: a charge breeding system that includes RFQ cooler, EBIS and Nier separator, and a high-resolution mass separator system. The latter is designed to achieve a resolving power of twenty thousand for a transmitted emittance of three micrometer. The separator optics has been designed with symmetry in order to minimize high-order aberrations. The dispersion of the system is created by two identical ninety degree magnetic dipoles with a field flatness of one part in one hundred thousand. The dipoles are tested and the magnetic field characterized before being installed on line for operation. High-order aberrations can also be corrected by an electrostatic multipole; this features a novel design as well as a new tuning technique. In this paper we will present the latest results from the field characterization and discuss the high level application to tune the multipole."
}
@article{de1983first,
title={A First Order Space Charge Option for TRANSOPTR},
author={De Jong, MS and Heighway, EA},
journal={IEEE Transactions on Nuclear Science},
volume={30},
number={4},
pages={2666--2668},
year={1983},
publisher={IEEE}
}
No preview for this file type
......@@ -219,7 +219,7 @@
Jim warned me that in this file the vertical emittance is 6\,$\mu$m (while the horizontal one is 3\,$\mu$m), and that the aspect ration between vertical size and vertical angle was wrong by a factor of 4: to fix the aspect ratio and bring the vertical emittance down to 3\,$\mu$m I use the scaling factors on line 5 of the {\tt 'OBJECT'} input.
{\it Assignment: write an input file to simulate the entire HRS from source slit to image slit. Track the 3\,000 particles in noEspread.beam and same their final coordinates to a file. Now change the mass of the particle by a factor 1/20\,000 (i.e. change its {\tt BRHO} by a factor 1/10\,000), track gaian and save the final coordinates into another file. Plot the horizontal phase space coordinates from both files, you should get a plot like~\cref{fig:twomassesInit}:}
{\it Assignment: write an input file to simulate the entire HRS from source slit to image slit. Track the 3\,000 particles in noEspread.beam and same their final coordinates to a file. Now change the mass of the particle by a factor 1/20\,000 (i.e. change its {\tt BRHO} by a factor 1/40\,000)\footnote{Take two particles with the same energy $qV=\frac{P^2}{2m}$, but with adifferent mass $s$. Logarithmic differentiation leads to: $0=\text{d}\ln\left(\frac{P^2}{2m}\right)\implies 2\frac{\text{d}B\rho}{B\rho}=\frac{\text{d}m}{m}$}, track gaian and save the final coordinates into another file. Plot the horizontal phase space coordinates from both files, you should get a plot like~\cref{fig:twomassesInit}:}
\begin{figure}[htb]
\centering
......
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