“on the applicability of phase space tomography to ......joshua auriemma ([email protected])...
TRANSCRIPT
“On the Applicability of Phase SpaceTomography to Reconstruct thePhase Space of a Space-Charge-
Dominated Electron Beam”
By Joshua Auriemma
Advisor: Philippe Piot
August 4, 2004
Basic Underlying Theory (BUT)
• With knowledge of what the electronbunch looks like at some time (t ) and theenvironmental changes endured by thebunch since the initial time (t0 ) – which inthis case include drift lengths, quadstrength, and quad length – it is possibleto obtain a very good representation ofthe bunch at time t0 .
Experimental Setup
Transfer Matrices
Drifts: The drift is an extremely simple beam line elementto model. In the drift, we do not have to take intoaccount any forces on the particles. The only things thatwe must take into account are x’, y’ and the length (L )of the drift.
Quadrupoles: Quads are used to focus the beam on oneplane and defocus it on the orthogonal plane. Because ofthis, the transformation is more difficult than the driftmatrix, and is dependant on the angle (which isdependant on the quad strength, k )
Why The Matrices?
• As the bunch travels through the quadrupolesand the drift, the beam distribution becomesrotated and skewed.
Enter The Matrix
In our program, the drifts are
given the following matrix:
While quadrupole transfer
matrices are solved via the
matrix to the left.
Two Programs(Both Alike in Dignity)
The Transfer Matrix Program: A program developed by me at Fermilab.
The purpose of the program is to solve the transfer matrices of thebeam as it progresses through space. It then applies the final
transfer matrix to a 1x4 vector of x, x’, y, and y’. The output of this
is later applied to the Max Entropy program.
The Max Entropy Program: The brilliant part of the code was done at
DESY in Germany. My contribution to the code was getting itexecute its algorithms on external (real) data as opposed to the
internal test which it came with. This program will take thehistograms and their associated angles as input, and attempt to
recreate the original bunch distribution.
Transfer Matrix Program
• The Transfer Matrix Program asks how many quads the
beam went through, the strength of the quads, their
lengths, and the length of the drifts. It will then calculate
all of the transfer matrices and apply the final matrix to
the data set. From this final matrix, we can then calculate
the angle of polarization:
_ = arctan(x21/x11)
Sample TMP Input
Transverse Beam Distribution
-6.00E-04
-5.00E-04
-4.00E-04
-3.00E-04
-2.00E-04
-1.00E-04
0.00E+00
1.00E-04
2.00E-04
3.00E-04
-6.00E-04 -4.00E-04 -2.00E-04 0.00E+00 2.00E-04 4.00E-04 6.00E-04 8.00E-04
Sample TMP Output1D Histogram at 1.45466 rad
0
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 50 100 150 200 250 300 350 400 450 500
bins
pa
rtic
le e
ve
nts
1D Histogram at 1.39095 rad
0
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 50 100 150 200 250 300 350 400 450 500
bins
pa
rtic
le e
ve
nts
1D Histogram at 1.39783 rad
0
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 50 100 150 200 250 300 350 400 450 500
bins
pa
rtic
le e
ve
nts
1D Histogram at 1.39316 rad
0
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
0 50 100 150 200 250 300 350 400 450 500
bins
pa
rtic
le e
ve
nts
Max Entropy Program
Recreated Distributions
Conclusions
• Even Gaussian distributions can be veryaccurately recreated by using themaximum entropy algorithm.
• It is likely that an optimum number ofhistograms needed to accurately recreatethe bunch will be obtained with furtherresearch.
Contact, Thanks, Etc.
Joshua Auriemma ([email protected])
University of Massachusetts - Lowell
35 Standish St. MS# 4566
Lowell, MA 01854
Philippe Piot ([email protected])
Thanks to Roger, Erik, and Max for paying us to learn physics at the
best particle accelerator in the world.
And most importantly, thanks to the interns whose interesting
mesh of personalities and corny physics humor will leave
me with fond memories for a long time to come.