Difference between revisions of "Sadiq Thesis Latex"

From New IAC Wiki
Jump to navigation Jump to search
Line 1,025: Line 1,025:
 
\subsection{The Positron Beam on Virtual Detectors DQ4 and DD1}   
 
\subsection{The Positron Beam on Virtual Detectors DQ4 and DD1}   
 
    
 
    
The positron beam energy distributions at the entrance of Q4 observed by virtual detector DQ4 is shown in Figure~\ref{STSimS1_En_DQ4_DD1}.  The distribution of positrons energies observed leaving the downstream side of the procturion target T1 as observed by virtual detector DD1UP are also shown. The ration of positrons leaving the target T1 and entering the first quadrupole indicate that nearly 90\% of the positrons are lost.  The large positron divergence is responsible for this loss. The median positron energy does not change suggesting that the beamline simulation is tuned to transport the optimal number of positrons emitted by the positron production target, T1.
+
The positron beam energy distributions at the entrance of Q4 observed by virtual detector DQ4 is shown in Figure~\ref{STSimS1_En_DQ4_DD1}.  The distribution of positrons energies observed leaving the downstream side of the production target T1 as observed by virtual detector DD1UP are also shown. The ratio of positrons leaving the target T1 and entering the first quadrupole indicate that nearly 90\% of the positrons are lost.  The large positron divergence is responsible for this loss. The median positron energy does not change suggesting that the beamline simulation is tuned to transport the optimal number of positrons emitted by the positron production target, T1.
  
 
\begin{figure}[htbp]  
 
\begin{figure}[htbp]  

Revision as of 22:08, 5 March 2014

Chapter 1: Introduction File:Sadiq thesis chapt 1.txt

\chapter{Introduction}

The positron production efficiency using a proposed linac based quadrupole triplet collection system is measured in partnership with the Continuous Electron Beam Accelerator Facility (CEBAF) of Thomas Jefferson National Accelerator Facility (commonly referred to as Jefferson Lab or JLab) and the Idaho Accelerator Center (IAC) of Idaho State University (ISU). There has been substantial interest in developing a polarized positron source for the nuclear physics community~\cite{PolPos}. The main challenge has been increasing the intensity and polarization of the beam. This work investigates a proposed method for collecting positrons using a quadrupole triplet system.


\section{Different Approaches to Produce Positrons} Positrons are used in several disciplines of science, such as chemistry, physics, material science, surface science, biology and nanoscience~\cite{Chemerisov:2009zz}. Positrons may be obtained either from radioactive sources or from pair production of high energy photons. Radioactive sources like Co-58, Na-22, and Cu-64 emit positive beta particles $\beta^{+}$ when they decay. As given in Table~\ref{tab:e+_source}, the positron rate of radioactive sources is about $10^{7} \sim 10^{8}$ e$^{+}$/s because of the source activity limits of 100~mCi~\cite{PosSource1} that are imposed for radiation safety.

Positrons are available for use in experiments at the nuclear reactors~\cite{PosSource2} shown in Table~\ref{tab:e+_source}. Neutrons are used to produce positrons from the photons emitted when neutrons interact with cadmium. Photons pair produce positrons within the tungsten foils positioned in the photon flux. The positron beam created in the nuclear reactor has the intensity shown in Table~\ref{tab:e+_source}.

Another common method used to generate positrons is by using an electron linear accelerator (linac). In this method, electrons, from the linac, incident on a high Z material like tungsten produce positrons by bremsstrahlung and pair production. One of the advantages of the linac based positron beam is its variable energy, intensity, and the ability to pulse the source on and off at specific time intervals as given in Table~\ref{tab:e+_source}.

\begin{table} \centering \caption{Different Positron Sources~\cite{PosSource1,PosSource2}.} \begin{tabular}{llll} \toprule {} & {} & {19:09, 24 December 2013 (MST)$\beta^{+}$ sources} & {} \\ \midrule {Source} & {Activity} & {Rate (e$^{+}$/s)} & {Facility/Source Type} \\ \midrule %{Source} & {Energy Range} & {Intensity} & {Production Method} \\ {Co-58} & {100~mCi} & {$3\times10^{5}$ } & {W-moderator} \\ {Na-22} & {100~mCi} & {$2\times10^{6}$ } & {W-moderator} \\ {Cu-64} & {80~Ci/cm$^2$} &{$10^{7}\sim10^{8}$} & {Brookhaven, reactor beam} \\

\toprule {} & {} & {~~Pair-production} & {} \\ \midrule {Energy } & {} & {Rate (e$^{+}$/s)} & {Facility/Source Type} \\ \midrule {36~MeV} & {} & {$10^{8}$} & {Giessen, LINAC} \\ {100~MeV} & {} & {$10^{10}$} & {Livermore, LINAC} \\ {150~MeV} & {} & {$10^{8}$} & {Oak Ridge, LINAC} \\ {5~keV} & {} & {$5\times10^{6}$} & {NC State University, reactor} \\ \bottomrule \end{tabular} \label{tab:e+_source} \end{table}

\section{Positron Beam Generation from \\ Bremsstrahlung} When a moving charged particle interacts with the electric field of another charged particle, it can be deflected and lose energy in the form of photons, as shown in Figure~\ref{fig:Theo-Brem}. This interaction is known as the bremsstrahlung process. The probability of this interaction increases with the square of the atomic number of the material traversed by the incident charged particle. Figure~\ref{fig:Brems_photon_Ene} shows the photon energy distribution produced when the 12~MeV electron energy distribution from Figure~\ref{fig:Theo-Brems_ele_Ene} interacts with a 1~mm thick tungsten target. As shown in Figure~\ref{fig:Brems_photon_Ene}, the distribution peaks at 0.3~MeV.

\begin{figure}[htbp] \centering \includegraphics[scale=0.40]{1-Introduction/Figures/bremsstrahlung/brems2.eps} \caption{Photon emission from the bremsstrahlung processes.} \label{fig:Theo-Brem} \end{figure}

The bremsstrahlung cross section for the energy range of this experiment given is by~\cite{brms-cors} %\begin{equation} %d \sigma = 4 Z^2r_e^2 \alpha \frac{d \nu}{\nu} \left \{ \left (1 + \left( \frac{E}{E_0} \right )^2 \right ) \left [ \frac{\phi_1(\gamma)}{4} - \frac{1}{3} \ln Z -f(Z)\right ] - \frac{2E}{3E_0} \left [ \frac{\phi_2(\gamma)}{4} - \frac{1}{3} \ln Z -f(Z)\right ] \right \}, %\label{eq:Brem-cross} %\end{equation} \begin{equation} \begin{array}{cl} d\sigma= & 4Z^{2}r_{e}^{2}\alpha\frac{d\nu}{\nu}\{\left(1+\left(\frac{E}{E_{0}}\right)^{2}\right)\left[\frac{\phi_{1}(\gamma)}{4}-\frac{1}{3}\ln Z-f(Z)\right]\\

& -\frac{2E}{3E_{0}}\left[\frac{\phi_{2}(\gamma)}{4}-\frac{1}{3}\ln Z-f(Z)\right]\},

\end{array}

\label{eq:Brem}

\end{equation} \noindent where $E_0$ is initial total energy of the electron, $E$ is the final total energy of the electron, $\nu = \frac{E_0-E}{h}$ is frequency of the emitted photon, and $Z$ is atomic number of the target. $\gamma = \frac{100 m_ec^2 h \nu}{E_0 E Z^{1/3}}$ is the charge screening parameter, and $f(Z)$ is given by \begin{equation} f(Z) = (Z \alpha)^2 \sum_1^{\infty} \frac{1}{ n [ n^2 + (Z \alpha)^2]}, \end{equation} \noindent where $\alpha = \frac{1}{137}$ is the fine-structure constant, $\phi_1$ and $\phi_2$ are screening functions that depend on Z. \begin{figure}[htbp] \centering \includegraphics[scale=0.75]{1-Introduction/Figures/En_photon_dnT1_logY_3.eps} \caption{Bremsstrahlung photon energy distribution produced when the 12~MeV electron energy distribution from Figure~\ref{fig:Theo-Brems_ele_Ene} interacts with the simulation's 1~mm thick tungsten target.} \label{fig:Brems_photon_Ene} \end{figure} \begin{figure}[htbp] \centering \includegraphics[scale=0.75]{1-Introduction/Figures/En_e_upT1.eps} \caption{The electron energy distribution incident on a tungsten foil in simulation.} \label{fig:Theo-Brems_ele_Ene} \end{figure}

There are three competing processes that a photon can undergo when interacting with matter. Figure~\ref{fig:Theo-3pro-in-W} illustrates the cross-sections for the different interactions that can occur when a photon traverses tungsten as a function of photon energy. At electron volt (eV) energies, which are comparable to the electron atomic binding energy, the dominant photon interaction is the photoelectric effect. As the photon energy increases up to the kilo electron volt (keV) range, the Compton scattering process starts to be more dominant. Although the photon is totally absorbed during the photoelectric effect, photons merely lose energy when undergoing Compton scattering. As the photon energy reaches twice the rest mass energy of the electron, $i.e.$ 2 \begin{math} \times \end{math} 511~keV, pair production begins to occur. Pair production becomes the dominant interaction process when photon energies are beyond 5~MeV~\cite{Krane}. In this process, a photon interacts with the electric field of the nucleus or the bound electrons and is converted into an electron and positron pair.

\begin{figure}[htbp] \centering \includegraphics[scale=0.5]{1-Introduction/Figures/xcom/10.eps} \caption{The cross-sections for different types of photon interactions with tungsten as a function of photon energy~\cite{nistxcom}.} \label{fig:Theo-3pro-in-W} \end{figure}

Using natural units, where \begin{math}c \equiv 1\end{math}, the differential cross-section for pair production can be expressed as \begin{equation} \begin{array}{l} \frac{d \sigma}{d \epsilon_1 d \theta_1 d \theta_2} = 8 \left ( \frac{\pi a}{\sinh (\pi a)} \right )^2 \frac{a^2}{2 \pi} \frac{e^2}{\hbar c} \left ( \frac{\hbar}{m_e c }\right )^2 \frac{\epsilon_1 \epsilon_2}{k^3} \theta_1 \theta_2 \\ \\ \times \{ \frac{V^2(x)}{q^4} \left [ k^2 (u^2 + v^2) \xi \eta - 2 \epsilon_1 \epsilon_2 (u^2 \xi^2 + v^2 \eta^2 ) + 2 (\epsilon_1^2 + \epsilon_2^2)uv \xi \eta cos(\phi) \right ] \\ \\ \left . + a^2W^2(x) \xi^2 \eta^2 \left [ k^2(1 - (u^2+v^2)\xi \eta - 2 \epsilon_1 \epsilon_2 (u^2 \xi^2 + v^2 \eta^2) \right . \right .\\ \\ \left . \left .- 2 (\epsilon_1^2 + \epsilon_2^2) u v \xi \eta \cos(\phi)\right ]\right \}, \\ \end{array} \end{equation}

\noindent where $k$ is photon energy, $\theta_{1}$ and $\theta_2$ are the scattering angle of $e^+$ and $e^-$ respectively, $ \phi = \phi_1 - \phi_2$ is the angle between the $e^+$ and $e^-$ pair, $\epsilon_1$ and $\epsilon_2$ are the energy of the positron and electron respectively. Other constants are $u = \epsilon_1 \theta_1$, $v=\epsilon_2 \theta_2$, $\xi = \frac{1}{1+u^2}$, $\eta= \frac{1}{1+v^2}$, $q^2 = u^2 + v^2 + 2 u v \cos(\phi)$, $x= 1-q^2 \xi \eta$, $a = \frac{Ze^2}{\hbar c}$, $V(x) = 1 + \frac{a^2}{(1!)^2} + \frac{a^2 (1+a^2) x^2}{(2!)^2} + \frac{a^2 (1+a^2)(2^2+a^2)x^4 x^2}{(3!)^2} + \cdots$, and $W(x) = \frac{1}{a^2} \frac{d V(x)}{d x}$.

In pair production, positron and electron pairs are created back to back in the center of mass frame as given by \begin{equation} \gamma \rightarrow e^{+}+e^{-}. \label{eq:pair} \end{equation} In the lab frame, electrons and positrons tend to move in the direction of the photon, as shown in Figure~\ref{fig:Theo-pair-pro}. The positron and electron carry away the energy from the photon that is in excess of 1.022~MeV. In the center of mass frame, the kinetic energy is equally shared. Photons with an energy above 1.022~MeV in the bremsstrahlung spectrum of Figure~\ref{fig:Brems_photon_Ene} have the potential to create electron and positron pairs. Figure~\ref{fig:Theo-brem} is the simulation of 10 million 12~MeV mono energetic electrons impinging on a 1.016~mm thick tungsten target. Turning on the annihilation process resulted in a 511~keV peak on top of the bremsstrahlung spectrum. This 511~keV peak represents photons produced when the created positrons from the pair production annihilate with atomic electrons inside the tungsten target.

\begin{figure}[htbp] \centering \includegraphics[scale=0.40]{1-Introduction/Figures/Pair_Production/Pair_Production.eps} \caption{Pair production.} \label{fig:Theo-pair-pro} \end{figure} \begin{figure}[htbp] \centering \includegraphics[scale=0.7]{1-Introduction/Figures/pair/on/Eng3.eps} \\ \caption{Photon spectrum created by one million 12~MeV mono energy electrons impinging on a 1.016~mm thick tungsten target.} \label{fig:Theo-brem} \end{figure}

\section{Generation of Electron Beam Using a \\ Linac} Electrons need to be accelerated, using a linac (linear particle accelerator), up to at lest few mega electron volts (MeV) in order to generate bremsstrahlung photons with enough energy to produce electron and positron pairs. A linac is a device that accelerates charged particles to nearly the speed of light using an electromagnetic wave. The HRRL (High Repetition Rate Linac), located in the Beam Lab of ISU's Physics Department, can accelerate electrons up to 16~MeV peak energy, which is sufficient energy for positron production.

To optimize positron production, one needs to know the beam size and divergence of the electrons impinging the production target. The emittance and the Twiss parameters ($\alpha$, $\beta$, and $\gamma$) quantify these beam properties and are used as the input parameters for accelerator simulation tools. These tools can study beam transport and predict beam properties (beam size and divergence) along the beam line.

\subsection{Emittance and Twiss Parameters} In accelerator physics, a Cartesian coordinate system is used to describe the motion of the accelerated particles. As shown in Figure~\ref{fig:coordinates}, the $s$-axis of the Cartesian coordinate system is defined as the natural coordinate that is oriented along the same direction as the beam momentum. The $x$-axis and $y$-axis are horizontal and vertical coordinates which constitute the transverse beam profile. The transverse beam profiles are described as a function of the longitudinal coordinates, $x(s)$ and $y(s)$.

\begin{figure}[htbp] \centering \includegraphics[scale=0.7]{1-Introduction/Figures/coordinates.eps} \caption{Coordinated system and reference orbit (dashed line)~\cite{Conte}.} \label{fig:coordinates} \end{figure}

The horizontal phase space $x'$ vs. $x$ (similar for $y$ projection), shown in Figure~\ref{fig:phase-space}, of the beam is an ellipse with invariant area along the beamline (under conditions that space charge effects, coherent synchrotron radiation, and wakefield are ignorable)~\cite{Conte}. Here $x'$ is defined as \begin{equation} x'=\frac{dx}{ds}. \label{eq:divergence} \end{equation} \noindent The area of the ellipse is \begin{equation} A = \pi w = \pi (\gamma z^2 + 2\alpha z z^{\prime} + \beta z^{\prime 2}). \label{eq:el-area} \end{equation} \noindent Here $w$ is called the Courant-Snyder invariant~\cite{Conte} and $\alpha$, $\beta$, and $\gamma$ are called Twiss parameters. The transverse emittance $\epsilon$ of the beam is defined to be the area of the ellipse that contains fraction of the particles. The units for the emittance are m$\cdot$rad and mm$\cdot$mrad. Conventionally used unit for emittance is $\mu$m, which is similar as mm$\cdot$mrad (1 $\mu$m = 1 mm$\cdot$mrad = $10^{-6}$ m$\cdot$rad). Twiss parameters are useful because they are related to the beam size and divergence by \begin{equation} \sigma_{x}(s)=\sqrt{\epsilon _x (s) \beta _x (s)},~ \sigma_{x'}(s)=\sqrt{\epsilon _x (s) \gamma _x (s)}, \label{eq:twiss-emit} \end{equation} \noindent where $\epsilon_{x}$ is beam's horizontal emittance, $\sigma_{x}$ is the horizontal rms beamsize, $\sigma_{x'}$ is the horizontal rms beam divergence, and $\beta _x $ and $\gamma _x$ are two of the three Twiss parameters (same for vertical projection). The Twiss parameters are related by \begin{equation} \gamma = \frac{1 + \alpha^2}{\beta}. %\text{ or } \beta _x (s) \gamma _x (s) - \alpha _x (s)^2 = 1. \label{eq:twiss} \end{equation} \begin{figure}[htbp] \centering \includegraphics[scale=0.5]{1-Introduction/Figures/Emttance_Ellipse.eps} \caption{Phase space ellipse and its relation to the emittance and the Twiss parameters $\alpha$,$\beta$, and $\gamma$~\cite{Conte}.} \label{fig:phase-space} \end{figure} The beam emittance is inversely proportional to the beam energy. As the electron beam gains energy during acceleration, the divergence of the beam decreases because the momentum increases in forward direction. This emittance decrease due to the energy gain is called adiabatic damping. Thus, it is more practical to use normalized emittance which is defined as \begin{equation} \epsilon_{\text{n}} = \beta\gamma\epsilon, \label{eq:nor_emi} \end{equation} where $\beta$ is $c/v$, $\gamma$ is relativistic (Lorentz) factor, and $\epsilon$ unnormalized emittance.

In this chapter, the possible ways to generate positrons were described. Positrons from radioactive sources have lower intensities due to safety limits placed on the source activity. Positrons from reactors have intensities comparable to radioactive sources but are low in energy. Generating positrons using a linac can produce higher energy positrons at higher fluxes than reactors.

Chapter 2 discusses the hardware used in this experiment and the tools used to measure emittance, Twiss parameter, and energy of the HRRL. The measured emittance, Twiss parameter and energy, as well as the positron production efficiency are discussed in the Chapter 3. The simulation of the positron production efficiency and transportation process are compared with measurement in Chapter 4. The conclusion from this work is presented in Chapter 5.

Chapter 2: Apparatus File:Sadiq thesis chapt 2.txt

\chapter{Apparatus} This chapter describes the apparatus and associated hardware used in the experiment to produce positrons as described in Chapter 1. The new HRRL beamline was constructed using two dipoles as well as ten quadrupoles to optimize the beam transportation. Beamline elements were aligned using a laser alignment system. Additional beamline elements such as an energy slit, beam viewers, and Faraday cups were added. Two NaI detectors were installed at the end of the beamline and were used to detect 511~keV photons emitted when positrons annihilate in a second tungsten target.

\section{HRRL Beamline} A 16~MeV S-band (2856~MHz RF frequency) standing-wave High Repetition Rate Linac (HRRL) located in the Department of Physics Beam Lab at Idaho State University was used to impinge a 12~MeV electron beam onto a tungsten foil. The energy of the HRRL is tunable between 3 to 16~MeV and its repetition rate is variable from 1 to 300~Hz. The operating parameters of the HRRL is given in Table~\ref{tab:hrrl-par}. As shown in Figure~\ref{fig:app-hrrl-cavity}, the HRRL has a thermionic gun, vertical and horizontal steering magnet sets on two ends, and two solenoid magnets. The startup, shutdown, and beam optimization procedure of the HRRL is given in the Appendix C.

\begin{table} \centering \caption{The Basic Parameters of the HRRL.} \begin{tabular}{lcc} \toprule {Parameter} & {Unit} & {Value} \\ \midrule maximum energy & MeV & 16 \\ peak current & mA & 100 \\ repetition rate & Hz & 300 \\ absolute energy spread & MeV & 25\% \\ macro pulse length & ns & $>$50 \\ RF Frequency & MHz & 2856 \\ \bottomrule \end{tabular} \label{tab:hrrl-par} \end{table}

\begin{figure}[htbp] \centering \includegraphics[scale=0.74]{2-Apparatus/Figures/HRRL_Cavity4.png} \caption{The configuration of the HRRL cavity.} \label{fig:app-hrrl-cavity} \end{figure}

The accelerator's cavity was relocated to the position shown in Figure~\ref{fig:app-hrrl-line} to provide enough space for a beam line that can transport either positrons or electrons. The beam elements are described in Table~\ref{tab:app-hrrl-coordinates}. Quadrupole and dipole magnets, described in Appendix D, were added to the new beam line as well as an OTR and YAG view screen to measure and monitor beam profile and position. Faraday cups and toroids were installed to measure the electron beam current. Energy slits were installed to control the energy/momentum spread of the beam after the first dipole (D1).

A 1.016~mm thick retractable tungsten (99.95\%) foil target (T1) was placed between the 1st and 2nd quadrupole triplets and used to produce positrons when the electron beam interacts with it. The room where the HRRL is located is divided by a wall into two parts; the accelerator side and the experimental cell. A beam pipe at the end of the 90 degree beamline goes through a hole in the wall and delivers the beam from the accelerator side to the experimental cell. The positron detection system consisting of two NaI detectors was placed at the end of the beamline in the experimental cell side as shown in Figure~\ref{fig:app-hrrl-line}.

\begin{sidewaysfigure} \centering \includegraphics[scale=0.265]{2-Apparatus/Figures/HRRL_line2.eps} \caption{The HRRL beamline layout and parts.} \label{fig:app-hrrl-line} \end{sidewaysfigure}

\begin{table} \centering \caption{The HRRL Beamline Parts and Coordinates.} \begin{tabular}{llll} \toprule {Label} & {Beamline Element} & {Distance from} & {} \\ {} & {} & {Linac Exit (mm)} & \\ \midrule Q1 & quadrupole & 335 & \\ Q2 & quadrupole & 575 & \\ Q3 & quadrupole & 813 & \\ T1 & e$^+$ production target & 1204 & \\ Q4 & quadrupole & 1763 & \\ Q5 & quadrupole & 2013 & \\ Q6 & quadrupole & 2250 & \\ D1 & dipole & 2680 & \\ S1 & OTR screen & 3570 & \\ FC1 & Faraday cup & 3740 & \\ EnS & energy slit & 3050 & \\ S2 & YAG screen & 3410 & \\ Q7 & quadrupole & 3275 & \\ D2 & dipole & 3842 & \\ FC2 & Faraday cup & 4142 & \\ Q8 & quadrupole & 4044 & \\ Q9 & quadrupole & 4281 & \\ Q10 & quadrupole & 4571 & \\ T2 & annihilation target & 7381 & \\ \bottomrule \end{tabular} \label{tab:app-hrrl-coordinates} \end{table} The positron production target, T1 was shielded with 8 inches of Fe bricks and 4 inches of Pb bricks to lower photon and electron background. The two dipoles were also shielded with Pb bricks because the beam scrapes the dipole vaccuum chamber when it is transported. An eight inch thick Pb wall was installed both sides of the wall dividing the accelerator and the experimental cell to shield the cell from background generated on the accelerator side. The NaI detectors were shielded with Pb brick to lower background. Positioning the NaI detector in a separate room from the positron production target T1 was done to decrease the high background that is historically present in the accelerator room. While this choice reduces the background there is also substantial signal loss, as will be shown in subsequent chapters, transporting a chromatic positron sample to the experimental cell.

\section{HRRL Beamline Alignment Using Laser} The HRRL beamline was aligned using a laser beam as shown in Figure~\ref{fig:alignment}. The gun of the HRRL was removed so that the laser beam from a laser placed on a table in the experimental side of the HRRL cell would be directed through the cavity. Mirrors were mounted onto holders with horizontal and vertical adjustments. The laser was first adjust with two mirrors on the laser table and focused using two focusing lenses. Two mirrors, one in the experimental side and one at on the linac side, reflected the laser beam through the center of the linac.

\begin{figure}[htbp] \centering \includegraphics[scale=0.24]{2-Apparatus/Figures/HRRL_Alignment.png} \caption{HRRL beamline alignment using laser.} \label{fig:alignment} \end{figure}

The laser beam was shot through the center of the HRRL cavity and the geometrical center of the 0 degree beamline magnets (quadrupoles Q1$\sim$Q6, first dipole D1) downstream were aligned according to the laser beam. The laser beam was reflected by a mirror mounted on a rotator that reflected the beam by 90$^{\circ}$ to the 90 degree beamline and the quadrupole magnets in the 90 degree beamline were aligned relative to the laser beam reference point.

Two irises were placed at the end of the 90 degree beamline and centered on the reflected 90 degree laser beam. A second laser was mounted on the wall in the experimental cell side. The laser on the wall was positioned to pass through the center of the two irises placed at the end of the 90 degree beamline. Thus, the laser on the wall was aligned to the center of the 90 degree beamline and can be used as a reference as shown in Figure~\ref{fig:laser}. \begin{figure}[htbp] \centering \includegraphics[scale=0.23]{2-Apparatus/Figures/BeamlineParts/laser.jpg} \caption{The laser mounted on the wall in experimental cell side. This laser was aligned to the center of the 90 degree beamline.} \label{fig:laser} \end{figure}

At last, the second dipole, D2, was connected to the 90 degree beampipe and beamline elements between D1 and D2 (energy slits, Q7, and 5-way cross holds YAG screen) were placed one after another.

\section{Energy Slit and Flag Controller} A control box was built to open and close the energy slit as well as control the three flags shown in Figure~\ref{fig:controller}. The slit control is on the right side of the box and the maximum width of the slit is 3.47~cm as shown in Figure~\ref{fig:controller}. On the left of the controller are the switches for three flags. Power supplies were installed inside the box and controlled by these switches to remove the targets from the beamline, turn on/off the cameras of the flags, and turn on/off the lights within each flag.

\subsection{Energy Slit Controller} The controller was built to open or close the energy slit (Danfysik water cooled slit model 563 system 5000) based on the design from Danfysik~\cite{Danfysik} as shown in Figure~\ref{fig:controller} (on right side). The wiring diagram of the controller is shown in Figure~\ref{fig:controllerConnection}. When the energy slit is fully open, the width between slits is 3.47~cm as indicated by the LED number display. When the slit is fully open/closed, one of the two LED lights will light up and the motor will stop.

The power source of the energy slit provides 12~VDC, 1.8~A max current, 20~W power, and takes 85$\sim$264~VAC input. The motor and the relay switch both takes 12~VDC. The LED indicating lights take 12~V/50~mA current and the 57~$\Omega$ resistors take 0.5~W power.

A potentiometer is placed inside the energy slit and the resistance of the potentiometer changes with the width of the slit. A 370~mV voltage is applied to the potentiometer. An LED number display is connected to the potentiometer so that the voltage change in the potentiometer is numerically indicated by the LED display. The slit width is indicated by the potentiometer resistance in the circuit.The potentiometer resistance in the circuit is indicated by the voltage across the resistor in the circuit which is displayed on the LED number display. For example, when the width of the energy slit is 3.24~cm, the voltage on the potentiometer read by the LED display is 324~mV, and the LED number display indicates 3.24. The cm unit is labeled on the right of LED display. \begin{figure}[htbp] \centering \includegraphics[scale=0.29]{2-Apparatus/Figures/BeamlineParts/controller.jpg} \caption{Front panel of the energy slit controller (on the right) and flag controllers (on the left). } \label{fig:controller} \end{figure} \begin{figure}[htbp] \centering \includegraphics[scale=0.43]{2-Apparatus/Figures/controllerConnection.png} \caption{Controller wiring diagram of energy slit controller (modified from the design given in the Danfysik water cooled slit model 563 system 5000 manual)~\cite{Danfysik}.} \label{fig:controllerConnection} \end{figure}

\subsection{Flag Controller} On the left of the controller box are the switches for three flags (0 degree/OTR flag, 45 degree/YAG flag, and 90 degree/annihilation target flag) as shown in Figure~\ref{fig:controller}. Four power supplies (one back up) were installed inside the box and controlled by these switches. The switches in the top row insert or remove the targets from the beamline and the ones in the middle turn on/off the cameras of the flags. The bottom row switches turn on/off the lights of the flag which lights up the target.

\section{The OTR Imaging System} The OTR target is a 10~$\mu$m thick aluminum foil with a 1.25 inch diameter. A moving charged particle produces electromagnetic fields that are dependent on the dielectric constant of the media the particle traverses. When a moving charge particle cross the boundary between two media (vacuum and aluminum in this case), electromagnetic radiation is emitted to reorganize the fields in the new media. The radiation is emitted in both forward (beam direction) and backward (image charge direction) as shown in Figure~\ref{fig:otr}. \begin{figure}[htbp] \centering \includegraphics[scale=0.6]{2-Apparatus/Figures/OTR.eps} \caption{OTR pattern when the incident beam is at $45^{\circ}$ angle with respect to the foil~\cite{OTR-Gitter}.} \label{fig:otr} \end{figure}

Three two inch diameter lenses were used for the imaging system to avoid optical distortion at lower electron energies. The focal lengths and positions of the lenses, shown in Figure~\ref{image_sys}, were calculated with the thin lens equation to produce the desired image magnification. The lenses, assembly rods, and lens holder plates for the camera cage system were purchased from Thorlabs, Inc.~\cite{thorlabs}.

The camera used is a JAI CV-A10GE digital 1/2" progressive scan camera with a 767 (horizontal) by 576 (vertical) pixel area and 6.49 (horizontal) by 4.83 (vertical) mm sensing area. It has high speed shutter from 1/60 to 1/300,000 second, edge pre-select, pulse width trigger modes, auto shutter, and smear-less mode. The images were taken by triggering the camera (in edge pre-select mode) synchronously with the electron gun. \begin{figure}[htbp] \centering {\scalebox{0.46} [0.46]{\includegraphics{2-Apparatus/image_sys2.eps}}} {\scalebox{0.5} [0.5]{\includegraphics{2-Apparatus/MOPPR087f3}}} \caption{The OTR imaging system.} \label{image_sys} \end{figure}

\section{Positron Detection} When the electron beam is incident on T1, photons and secondary electrons are created along with positrons. These particles are the main source of noise in the experiment. The positrons were transported to the second tungsten target (T2) which was shielded from this background by the concrete wall and Pb bricks. The setup is shown in Figure~\ref{fig:HRRL-pos-det-setup}. A 6-way cross was placed at the end of the beamline to hold T2. The 6-way cross has three 1~mil (0.0254 millimeters) thick stainless steel windows. The two horizontal windows perpendicular to the beamline allowed the 511~keV photons created from the positron annihilation to escape the beamline with a limited attenuation. A third window at the end of the 90 degree beamline was used as the beam exit. Two NaI detectors were placed facing the two exit windows to detect the photons produced in T2. A 2 inch thick lead brick collimator with a 2 inch diameter hole was placed between the exit window and NaI detector. A scintillator (Scint) and a Faraday cup (FC3) were placed at the end of the beamline and were used to tune the electron and positron beam. When positrons reach T2, they can thermalize and annihilate inside T2. During thermalization, a positron looses its kinetic energy. When it annihilates with an electron, two 511~keV photons are emitted back to back. A triplet coincidence is required between the accelerator RF pulse and the detection of a photon in each NaI detector.

\begin{figure}[htbp] \centering \includegraphics[scale=0.50]{2-Apparatus/HRRL_Pos_detection2.eps} \caption{The positron detection system: T2 was placed with 45$^\circ$ angle to the horizontal plane first, then rotated 45$^\circ$ along the vertical axis.} \label{fig:HRRL-pos-det-setup} \end{figure}

\subsection{NaI Detectors}

NaI crystals, shown in Figure~\ref{fig:PMT}, were used to detect 511~keV photons from positron annihilation. Originally, the detectors had pulse signal lengths around 400~$\mu$s. New PMT base were built to use the HV divider shown in Figure~\ref{fig:PMT_base}. A picture of the constructed bases is shown in Figure~\ref{fig:new_base_made}. The pulse length of the new PMT base is about 1~$\mu$s. The NaI crystal is from Saint-Gobain Crystal \& Detectors (Mod. 3M3/3) with a dimension of $3 \times 3$. Operating high voltage of -1150~V would position the 511~keV photons within the range of the charge sensing ADC (CAEN Mod. V792).

\begin{table} \centering \caption{The Radioactive Sources and Corresponding Photon Peaks.} \begin{tabular}{lccc} \toprule {Radioactive Sources} & Unit & First Peak & Second Peak \\ \midrule Co-60 & keV & 1173 & 1332 \\ Na-22 & keV & 511 & 1275 \\ \bottomrule \end{tabular} \label{tab:Na22_Co60} \end{table}

\begin{figure}[htbp] 
\centering
\includegraphics[scale=0.52]{2-Apparatus/SAINT-GOBAIN_3M33.png}
\caption{The NaI crystal dimension.}
\label{fig:PMT}
\end{figure}

\begin{figure}[htbp] \centering \includegraphics[scale=0.75]{2-Apparatus/Modified_PMT.png} \caption{The modified PMT base design.} \label{fig:PMT_base} \end{figure}

\begin{figure}[htbp] \centering \includegraphics[scale=0.13]{2-Apparatus/IAC_NaI.png} \caption{The NaI crystals and new bases.} \label{fig:new_base_made} \end{figure}

The NaI detectors were calibrated using a Na-22 and a Co-60 source with the photon peaks indicated in Table~\ref{tab:Na22_Co60}. Figure~\ref{fig:NaI_Co60_Scope} is the oscilloscope image of several Co-60 photon pulses observed by the detector with the new PMT. The calibrated NaI detector spectrum from the Na-22 and Co-60 sources is shown in Figure~\ref{fig:NaI-Calb}. The rms values of the fits on the four peaks shown in Figure~\ref{fig:NaI-Calb} are $\sigma_{Na, 511}=18.28\pm0.04$~keV, $\sigma_{Na, 1275}=44.51\pm0.27$~keV, $\sigma_{Co, 1173}=42.49\pm0.24$~keV, and $\sigma_{Co, 1332}=50.30\pm0.39$~keV.

\begin{figure}[htbp] \centering \includegraphics[scale=0.5]{2-Apparatus/NaI_Co60_Scope.png} \caption{Detector output pulses using the Co-60 source and new PMT. The amplitude of the pulse is about 60~mV. The rise time of the pulse is larger than 50~ns, and the fall time is larger than 700~ns.} \label{fig:NaI_Co60_Scope} \end{figure}

\begin{figure}[htbp] \centering \includegraphics[scale=0.77]{2-Apparatus/Figures/NaI_Calbration/NaI_Calb8.eps} \caption{The calibrated NaI spectrum of Na-22 and Co-60 sources.} \label{fig:NaI-Calb} \end{figure}

\subsection{The DAQ Setup}

The data acquisition (DAQ) setup and timing diagram is shown in Figure~\ref{fig:daq-setup}. The last dynode signals from left and right NaI detectors were inverted using a ORTEC 474 inverting amplifier and sent to a CAEN Mod. N842 constant fraction discriminator (CFD). The electron gun pulse from the linac was used to generate a VETO for the CFD to prevent the RF noise from triggering the CFD, otherwise the CFD would generate multiple digital pulses. The PMT output itself could produce multiple CFD output pulses due to the voltage fluctuations within a single pulse. To prevent multiple CFD outputs, a GG~8000-01 octal gate generator was used to create a single 1~$\mu$s wide pulse from the CFD first pulse in order to ignore the multiple CFD pulses produced by a single analog output pulse from the detector. A triple coincidence was formed between the gun pulse and the 1~$\mu$s wide pulse from each detector using a LeCroy model 622 logic module. %\begin{equation} %\text{NaI~Left~\&\&~NaI~Right~\&\&~Gun~Trigger}. %\end{equation}

The ADC requires 5.7~$\mu$s to convert the analog signal to a digital signal. The logic module output was delayed 6~$\mu$s by a dual timer (CAEN Mod. N93B) to accommodate the ADC's conversion time and trigger the DAQ. The ADC (CAEN Mod. V792) converted the NaI detector's analog signals to digital when a 1~$\mu$s gate, created by the gun pulse using a dual timer, was present as shown in the lower part of Figure~\ref{fig:daq-setup}. The ADC was fast cleared unless it received a veto signal from the inverted output of logic module created using a dual timer.

\begin{sidewaysfigure}[htbp] \centering \includegraphics[scale=0.8]{2-Apparatus/Figures/DAQ_Logic_all.eps} \caption{The DAQ setup and timing diagram.} \label{fig:daq-setup} \end{sidewaysfigure}

Chapter 3: Data Analysis File:Sadiq thesis chapt 3.txt

\chapter{Data Analysis} In the first section of this chapter, the measurement of the HRRL emittance and Twiss parameters using the quadrupole scanning method is described. The current in the first quadrupole, Q1, was changed incrementally and the electron beam shapes were observed using the OTR screen located at the end of the 0 degree beamline. The emittance and Twiss parameters are obtained by fitting a parabolic function to the plot $\sigma_{x,y}^2$ vs K$_1$L (K$_1$ is quadrupole strength and L is the quadrupole effective length).

The second section describes the energy distribution measurement of the HRRL when it was tuned to generate an electron beam with its peak at 12~MeV. A Faraday cup, FC2, was placed at the end of the 45 degree beamline to measure the electron beam current when D1 was on and D2 was off. The energy scan data and the MATLAB script used for data analysis are given in the appendix A.

Last section is about the positron produced using the HRRL located in the Beam Lab of Physics Department at Idaho State University. The electron beam from the HRRL with 12~MeV peak incident on the tungsten production target (T1) produced positrons downstream and positrons were transported to annihilation target (T2), where they annihilate and crated back to back scattered 511~keV photon pairs. Two NaI detectors, their shielding, placement, and operation modes are described. The different methods to count positron rates are discussed. The electron beam current of the HRRL during the positron production was monitored by a scintillator. The calibration of the scintillator using a Faraday cup and oscilloscope is described. At last, the ratios of positrons detected using NaI detectors in coincidence mode to the electrons impinging on T1 are given for 1$\sim$5~MeV positron beam.

\section{Emittance Measurement} The HRRL beam emittance was measured by using an optical transition radiation (OTR) screen. This transition radiation was theoretically predicted by Ginzburg and Frank~\cite{Ginzburg-Frank} in 1946 to occur when a charged particle passes the boundary of two media. An OTR based viewer (31.75~mm in diameter and 10~$\mu$m thick aluminum) was installed to observe the electron beam size at 37.2~mA (higher currents might saturate the camera) electron peak currents, available using the HRRL at $15\pm1.5$~MeV with 200~ns macro pulse width. The visible light is produced when a relativistic electron beam crosses the boundary of two media (aluminum and air) with different dielectric constants.

Light is emitted in a conical shape at backward angles with the peak intensity at an apex angle of $\theta = 1/\gamma$ ($\gamma$ is relativistic factor) with respect to the incident electron's angle of reflection. A 15 MeV electron accelerated by the HRRL would emit light at $\theta = 2^\circ$. Orienting the OTR target at 45${^\circ}$ with respect to the incident electron beam will result in the high intensity photons being observed at an angle of 90${^\circ}$ with respect to the incident beam, see Figure~\ref{q-scan-layout}. These backward-emitted photons are observed using a digital camera and can be used to measure the shape and the intensity of the electron beam. Although an emittance measurement can be performed in several ways~\cite{emit-ways, sole-scan-Kim}, the quadrupole scanning method~\cite{quad-scan} was used to measure the emittance and Twiss parameters in this work.

\subsection{Emittance Measurement Using Quadrupole Scanning Method} Fig.~\ref{q-scan-layout} illustrates the beamline components used to measure the emittance for the quadrupole scanning method. A quadrupole is positioned at the exit of the linac to focus or de-focus the beam as observed on the OTR view screen. The 3.1~m distance between the quadrupole and the screen was chosen in order to minimize chromatic effects and to satisfy the thin lens approximation. %The quadrupole and the screen are located far away to minimize chromatic effects and to increase the veracity of the thin lens approximation used to calculate beam optics. \begin{figure}[htbp] \centering \includegraphics[scale=0.65]{1-Introduction/Figures/quad_scan_setup2.eps} \caption{Apparatus used to measure the beam emittance.} \label{q-scan-layout} \end{figure} Assuming the thin lens approximation, $\sqrt{k_1}L << 1$, is satisfied, the horizontal transfer matrix of a quadrupole magnet may be expressed as % thin lens approximation (sqrt{k1}*L << 1). In our case sqrt{k1}*L =0.07 \begin{equation} \label{quad-trans-matrix} \mathrm{\mathbf{Q}}=\Bigl(\begin{array}{cc} 1 & 0\\ -k_{1}L & 1 \end{array}\Bigr)=\Bigl(\begin{array}{cc} 1 & 0\\ -\frac{1}{f} & 1 \end{array}\Bigr), \end{equation} where $L$ is the length of quadrupole and $f$ is the focal length. $k_{1}$ is the quadrupole strength given by \begin{equation} \label{eq:k1} k_1 = 0.2998 \frac{g\text{(T/m)}}{p\text{(GeV/c)}} = 0.2998 \frac{ \frac{ B(I(\text{A})) }{ R_{\text{q}} } }{p\text{(GeV/c)}}, \end{equation} where $g$ is the gradient of the quadrupole with the Bore aperture radius $R_q$ at a given coil current $I$ and $p$ is the momentum of the electron beam. A matrix representing the drift space between the quadrupole and screen is given by \begin{equation} \label{drift-trans-matrix} \mathbf{\mathbf{S}}=\Bigl(\begin{array}{cc} 1 & l\\ 0 & 1 \end{array}\Bigr), \end{equation} where $l$ is the distance between the scanning quadrupole and the screen. The transfer matrix $\mathbf{M}$ of the scanning region is given by the matrix product $\mathbf{SQ}$. In the horizontal plane, the beam matrix at the screen ($\mathbf{\sigma_{s}}$) is related to the beam matrix of the quadrupole ($\mathbf{\sigma_{q}}$) using the similarity transformation \begin{equation} \mathbf{\mathbf{\sigma_{s}=M\mathrm{\mathbf{\mathbf{\sigma_{q}}}}}M}^{\mathrm{T}}. \end{equation} where the $\mathbf{\sigma_{s}}$ and $\mathbf{\sigma_{q}}$ are defined as~\cite{SYLee} \begin{equation} \mathbf{\mathbf{\sigma_{s,\mathnormal{x}}=}}\Bigl(\begin{array}{cc} \sigma_{\textnormal{s},x}^{2} & \sigma_{\textnormal{s},xx'} \\ \sigma_{\textnormal{s},xx'} & \sigma_{\textnormal{s},x'}^{2} \end{array}\Bigr) ,\; \mathbf{\mathbf{\sigma_{q,\mathnormal{x}}}}=\Bigl(\begin{array}{cc} \sigma_{\textnormal{q},x}^{2} & \sigma_{\textnormal{q},xx'}\\ \sigma_{\textnormal{q},xx'} & \sigma_{\textnormal{q},x'}^{2} \end{array}\Bigr). \end{equation} \noindent %By defining the new parameters~\cite{quad-scan}, $A \equiv \sigma_{11},~B \equiv \frac{\sigma_{12}}{\sigma_{11}},~C \equiv\frac{\epsilon_{x}^{2}}{\sigma_{11}}$ By defining the new parameters~\cite{quad-scan} \begin{equation} A \equiv l^2\sigma_{\textnormal{q},x}^{2},~B \equiv \frac{1}{l} + \frac{\sigma_{\textnormal{q},xx'}}{\sigma_{\textnormal{q},x}^{2}},~\text{and}~C \equiv l^2\frac{\epsilon_{x}^{2}}{\sigma_{\textnormal{q},x}^{2}}, \end{equation} the matrix element $\sigma_{\textnormal{s},x}^{2}$, the square of the beam size's rms at the screen, may be expressed as a parabolic function of the product of $k_1$ and $L$ \begin{equation} \sigma_{\textnormal{s},x}^{2}=A(k_{1}L)^{2}-2AB(k_{1}L)+(C+AB^{2}). \label{par_fit} \end{equation}

The emittance measurement was performed by changing the quadrupole current, which changes $k_{1}L$, and measuring the corresponding beam image size on the view screen. The measured two-dimensional beam image was projected along the image's abscissa and ordinate axes. A super Gaussian fitting function is used on each projection to determine the rms value, $\sigma_{\textnormal{s}}$ in Eq.~(\ref{par_fit}). Measurements of $\sigma_{\textnormal{s}}$ for several quadrupole currents ($k_{1}L$) are then fit using the parabolic function in Eq.~(\ref{par_fit}) to determine the constants $A$, $B$, and $C$. The unnormalized projected rms emittance ($\epsilon$) and the Twiss parameters ($\alpha$ and $\beta$) can be found using Eq.~(\ref{emit-relation}) \begin{equation} \epsilon=\frac{\sqrt{AC}}{l^2},~\beta=\sqrt{\frac{A}{C}},~\alpha=\sqrt{\frac{A}{C}}(B+\frac{1}{l}). \label{emit-relation} \end{equation}

\subsection{HRRL Emittance Measurement Experiment} A quadrupole scanning method was used to measure the accelerator's emittance. The quadrupole current is changed to alter the strength and direction of the quadrupole magnetic field such that a measurable change in the beam shape is seen by the OTR system. Initially, the beam was steered by the quadrupole indicating that the beam was not entering along the quadrupole's central axis. Several magnetic elements upstream of this quadrupole were adjusted to align the incident electron beam with the quadrupole's central axis. First, the beam current observed by a Faraday cup located at the end of beam line was maximized using upstream steering coils within the linac nearest the gun. Second, the first solenoid nearest the linac gun was used to focus the electron beam on the OTR screen. Steering coils were adjusted to maximize the beam current to the Faraday cup and minimize the deflection of the beam by the quadrupole. A second solenoid and the last steering magnet shown in Figure~\ref{fig:app-hrrl-cavity}, both near the exit of the linac, were used in the final step to optimize the beam spot size on the OTR target and maximize the Faraday cup current. A configuration was found that minimized the electron beam deflection when the quadrupole current was altered during the emittance measurements.

The emittance measurement was performed using an electron beam energy of 15~MeV and a 200~ns long macro pulse of 40~mA electron current. The current in the first quadrupole after the exit of the linac was changed from $-$~5~A to $+$~5~A with an increment of 0.2~A. Seven measurements were taken at each current step in order to determine the average beam width and the variance. Background measurements were taken by turning the linac's electron gun off while keep the RF on. OTR images of the beam is shown in Figure~\ref{fig:eBeam}. The digitized OTR images before and after background subtraction are shown in Figure~\ref{bg}. A small dark current is visible in Figure~\ref{bg} (b) that is known to be generated when electrons are pulled off the cavity wall and accelerated.

\begin{figure}[htbp] \centering \includegraphics[scale=0.6]{4-Experiment/Figures/ElectronBeam26} \caption{The 15~MeV electron beam observed using the OTR screen when dipole coil current was at 0. The macro pulse was 200~ns and the electron peak current was 40~mA .} \label{fig:eBeam} \end{figure}

\begin{figure}[htbp] \begin{tabular}{ccc} \centerline{\scalebox{0.42} [0.42]{\includegraphics{4-Experiment/Figures/MOPPR087f4.eps}}} \\ (a)\\ \centerline{\scalebox{0.42} [0.42]{\includegraphics{4-Experiment/Figures/MOPPR087f5.eps}}}\\ (b)\\ \centerline{\scalebox{0.42} [0.42]{\includegraphics{4-Experiment/Figures/MOPPR087f6.eps}}}\\ (c) \end{tabular} \caption{Digital image from the OTR screen; (a) an image taken with the beam on, (b) a background image taken with the RF on but the electron gun off, (c) The background subtracted beam image ((a)-(b)).} \label{bg} \end{figure}

The electron beam energy was measured using a dipole magnet downstream of the quadrupole used for the emittance measurements. Prior to energizing the dipole, the electron micro-pulse bunch charge passing through the dipole was measured using a Faraday cup located approximately 50~cm downstream of the OTR screen. The dipole current was adjusted until a maximum beam current was observed on another Faraday cup located just after the 45 degree exit port of the dipole. A magnetic field map of the dipole indicates that the electron beam energy was 15~$\pm$~1.6~MeV.

%\subsection{Data Analysis and Results} Images from the JAI camera were calibrated using the OTR target frame. An LED was used to illuminate the OTR aluminum frame that has a known inner diameter of 31.75~mm. Image processing software was used to inscribe a circle on the image to measure the circular OTR inner frame in units of pixels. The scaling factor can be obtained by dividing this length with the number of pixels observed. The result is a horizontal scaling factor of 0.04327~$\pm$~0.00016 mm/pixel and vertical scaling factor of 0.04204~$\pm$~0.00018 mm/pixel. Digital images from the JAI camera were extracted in a matrix format in order to take projections on both axes and perform a Gaussian fit. The observed image profiles were not well described by a single Gaussian distribution.

The profiles may be described using a Lorentzian distribution, however, the rms of the Lorentzian function is not defined. A super Gaussian distribution was used~\cite{sup-Gau}, because it has a sharper distribution than Gaussian and, unlike Lorentzian, the rms values may be directly extracted. The super Gaussian reduced the Chi-square per degree of freedom by a factor of ten compared to a typical Gaussian fit. The beam spot, beam projections, and fits are shown in Figure~\ref{Gau-SupGaus-fits}. In a typical Gaussian distribution, $\phi(x)=\frac{1}{\sqrt{2\pi}}e^{-\frac{1}{2}x^2}$, the independent variable $x$ is raised to the 2nd power while in the super Gaussian it is raised to a power smaller than 2. For example, the beam projections shown in Figure~\ref{Gau-SupGaus-fits} were fitted with the super Gaussian distributions. The variable $x$ on the exponent was raised to $0.9053$ ($x$-projection) and $1.0427$ ($y$-projection).

\begin{figure}[htbp] \begin{tabular}{lr} {\scalebox{0.38} [0.38]{\includegraphics{4-Experiment/Figures/Gau_SupGau/Gau_ChiSqaure.eps}}} {\scalebox{0.38} [0.38]{\includegraphics{4-Experiment/Figures/Gau_SupGau/SupGau_ChiSqaure.eps}}} \end{tabular} \caption{Gaussian and super Gaussian fits for beam projections. The beam images is background subtracted image and taken when quadrupole magnets are turned off. Left image is Gaussian fit and right image is super Gaussian fit.} \label{Gau-SupGaus-fits} \end{figure}

Figure~\ref{fig:par-fit} shows the square of the rms ($\sigma^2_{\textnormal{s}}$) $vs$ $k_1L$ for $x$ (horizontal) and $y$ (vertical) beam projections along with the parabolic fits using Eq.~\ref{fig:par-fit}. The emittance and Twiss parameters from these fits are summarized in Table~\ref{tab:results}. The MATLAB~\cite{MATLAB} scripts used to calculate emittance and Twiss parameters are given in appendix B. \begin{figure}[htbp] \begin{tabular}{lr} {\scalebox{0.27} [0.27]{\includegraphics{4-Experiment/Figures/par_fit_x.eps}}} {\scalebox{0.27} [0.27]{\includegraphics{4-Experiment/Figures/par_fit_y.eps}}} \end{tabular} \caption{Square of rms values and parabolic fittings. As the quadrupole current changes, so does quadrupole strength times quadrupole legnth, k$_1$L, and the square of the beam rms changes accordingly.} \label{fig:par-fit} \end{figure}

\begin{table} \centering \caption{Emittance Measurement Results} \begin{tabular}{lcc} \toprule {Parameter} & {Unit} & {Value} \\ \midrule unpolarized projected emittance $\epsilon_x$ & $\mu$m & $0.37 \pm 0.02$ \\ unpolarized projected emittance $\epsilon_y$ & $\mu$m & $0.30 \pm 0.04$ \\ % normalized \footnote{normalization procedure assumes appropriate beam chromaticity.} emittance $\epsilon_{n,x}$ & $\mu$m & $10.10 \pm 0.51$ \\ %normalized emittance $\epsilon_{n,y}$ & $\mu$m & $8.06 \pm 1.1$ \\ $\beta_x$-function & m & $1.40 \pm 0.06$ \\ $\beta_y$-function & m & $1.17 \pm 0.13$ \\ $\alpha_x$-function & rad & $0.97 \pm 0.06$ \\ $\alpha_y$-function & rad & $0.24 \pm 0.07$ \\ micro-pulse charge & pC & 11 \\ micro-pulse length & ps & 35 \\ energy of the beam $E$ & MeV & 15 $\pm$ 1.6 \\ relative energy spread $\Delta E/E$ & \% & 10.4 \\ \bottomrule \end{tabular} \label{tab:results} \end{table}

\section{Measurement of HRRL Electron Beam \\Energy Spread at 12~MeV} The HRRL energy profile was measured when it was tuned to accelerate electrons to 12~MeV peak energy. A Faraday cup, FC2, was placed at the end of the 45 degree beamline to measure the electron beam current when D1 was on and D2 was off. The dipole coil current for D1 was changed in 1~A increments. As the dipole coil current was changed, the energy of the electrons transported to the Faraday cup would change as described in Appendix A. Figure~\ref{fig:En-Scan} illustrates a measurement of the 12~MeV peak illustrating the observed low energy tail. The HRRL energy profile can be described by overlapping two skewed Gaussian fits~\cite{sup-Gau}. The fit function is given by \begin{equation} G(En)=A_{1}e^{\frac{-\left(En-\mu_{1}\right)^{2}}{2\left\{ \sigma_{1}\left[1+sgn\left(En-\mu_{1}\right)\right]\right\} ^{2}}}+A_{2}e^{\frac{-\left(En-\mu_{2}\right)^{2}}{2\left\{ \sigma_{2}\left[1+sgn\left(En-\mu_{2}\right)\right]\right\} ^{2}}}, \label{eq:skew-Gau} \end{equation} where $sgn$ is the sign function, that is defined as $$ sgn=\left\{ \begin{array}{ll}

  \mbox -1\ & {(x<0)} \\
  \mbox 0 \ & {(x=0)} \\
  \mbox 1 \ & {(x>0)}.\\
                  \end{array}\right.

$$ The other variables are defined as $\sigma_{1}=\frac{\sigma_{r,1}+\sigma_{l,1}}{2}$, $\sigma_{2}=\frac{\sigma_{r,2}+\sigma_{l,2}}{2}$, $E_{1}=\frac{\sigma_{r,1}-\sigma_{l,1}}{\sigma_{r,1}+\sigma_{l,1}}$, and $E_{2}=\frac{\sigma_{r,2}-\sigma_{l,2}}{\sigma_{r,2}+\sigma_{l,2}}$. The measurement results and fits are shown in Figure~\ref{fig:En-Scan} and in Table~\ref{tab:En-Scan_resluts}. \begin{figure}[htbp] \centering \includegraphics[scale=0.5]{4-Experiment/Figures/En_Fit_Assym_Gau.eps} \caption{HRRL energy scan (dots) and fit (line) with two skewed Gaussian distribution.} \label{fig:En-Scan} \end{figure} \begin{table} \centering \caption{Fit Parameters for Two Skewed Gaussian.} \begin{tabular}{lclcc} \toprule {Parameter} & Notation & Unit & {First Gaussian} & {Second Gaussian} \\ \midrule amplitude & A & mA & ~2.14 & 10.88 \\ mean & $\mu$ & MeV & 12.07 & 12.32 \\ sigma left & $\sigma_L$ & MeV & ~4.47 & ~0.70 \\ sigma right & $\sigma_R$ & MeV & ~1.20 & ~0.45 \\ \bottomrule \end{tabular} \label{tab:En-Scan_resluts} \end{table}

% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % \section{Positron Production Experiment} Positron produced using the HRRL located in the Beam Lab of Physics Department at Idaho State University. The electron beam of the HRRL was incident on T1, positron production target, and the positrons produced downstream of T1 was transported to T2, annihilation target, where they annihilate and crated back to back scattered 511~keV photon pairs. Two NaI detectors, shielded with Pb bricks, were placed horizontally and operated in coincidence with electron electron gun pulse to detect 511~keV photon pairs came out of T2. The electron beam current of the HRRL was monitored by a scintillator was placed between Q9 and Q10. The ROOT script used to calculate number/charge of positrons/electrons and the ratio of the two is given in the appendix.

%Positrons transported to the scintillator located at the end of 90 degree the beamline. \subsection{The Electron Beam Current Measurement} A scintillator was placed between Q9 and Q10, as shown in Figure~\ref{fig:Scint_e-}, to monitor the electron beam current. The electron beam current was changed incrementally to measure the correlation between the scintillator and the accelerated electron beam. The beam current was measured using FC1 and the output was integrated using an oscilloscope. The scintillator output was integrated using an ADC (CAEN Mod. V792). As the electron beam current was decreased, the signal observed on the oscilloscope decreased and the ADC measured less charge from the scintillator as shown in Figure~\ref{fig:ADC-CH9}. A linear fit to data resulted a linear relation \begin{equation} Q_{\text{e}^-}(i) = (0.0186 \pm 0.0028)i + 2.79 \pm 0.08~\text{C}, \end{equation} where $i$ is ADC channel number and $Q$ is the accelerated electron beam charge. The fit is shown in Figure~\ref{fig:calb-fit}.

\begin{figure}[htbp] \centering \includegraphics[scale=0.4]{4-Experiment/Figures/HRRL_line.eps} \caption{The electron beam monitor.} \label{fig:Scint_e-} \end{figure}

\begin{figure}[htbp] \centering \includegraphics[scale=0.76]{4-Experiment/Figures/ScintilatorCalibration/ADC_CHAN4.eps} \caption{The photon flux detected using scintillator. The mean of the ADC channel decreased linearly as the electron beam current was decreased. The electron beam current was measured using the Faraday cup at the end of the 0 degree beamline, FC1, and integrated using oscilloscope.} \label{fig:ADC-CH9} \end{figure}

\begin{figure}[htbp] \centering \includegraphics[scale=0.70]{4-Experiment/Figures/ScintilatorCalibration/calb.eps} \caption{Fit for accelerator beam current $v.s.$ the ADC channel.} \label{fig:calb-fit} \end{figure}

\begin{table} \centering \caption{Scintillator Calibration Data.} \begin{tabular}{lll} \toprule {Run Number} & {Faraday Cup} & {Mean of ADC} \\ { } & {Charge Area (nVs)} & {Channel 9} \\ \midrule 3703 & $1201 \pm 10$ & $1126 \pm 0.8$ \\ 3705 & $777 \pm 110$ & $791.8 \pm 0.6$ \\ 3706 & $367.7 \pm 2.3$ & $242.1 \pm 0.3$ \\ \bottomrule \end{tabular} \label{tab:scint_calb} \end{table}

When running in coincidence mode, the electron beam current is sampled by the scintillator only when a coincidence event occurs causing a trigger that gates the ADC and measures the scintillator output for that positron coincidence event. The charge measured by the ADC is the total charge of the electron pulses that created the positron events. The number of beam pulses are counted using a scaler. The total charge on T1 for the entire run is estimated using \begin{equation} Q_{\text{C}} = \left ( \overset{N}{\underset{i}{\sum}} 0.0186i\times(\text{bin~content}[i]) + 2.79 \right ) \times \frac{\text{(\# of beam pulses)}}{\text{(\# of events)}}. \label{eq:q_calc1} \end{equation} \noindent

\subsection{NaI Detectors} The NaI detectors were calibrated using radioactive sources as mention in the Chapter 2. However, in the positron production experiment conducted using HRRL, the photon peak observed by the NaI detector was not at the channels corresponding the 511~keV energy region (according to the calibration with sources). To determine the peak observed is the 511~keV peak, the Na-22 source was placed between two NaI detectors to measure the spectrum when the RF was on, represented by dotted line (no electrons were fired from accelerator gun), and off, represented by dashed line in Fig~\ref{fig:NaI-peak-shift}. The solid line in Fig~\ref{fig:NaI-peak-shift} represents the photon spectrum created by 3~MeV positrons impinging on T2.

When the RF was on, no difference was observed for photon peaks created by the 3~MeV positrons and the 511~keV photon peak from Na-22 source. However, when the RF was off, the 511~keV photon peak from Na-22 source occupied at different channels. It seems the 511~keV peak would shift when the RF is turned on. The 511 keV peak in the right NaI detector was shifted to the right side of the spectrum by 17 channels when the RF was on as shown in Figure~\ref{fig:NaI-peak-shift} (a) while in the left NaI detector the peak shifted to the left by 28 channels as shown if Figure~\ref{fig:NaI-peak-shift} (b).

\begin{figure}[htbp] \begin{tabular}{ccc} \centerline{\scalebox{0.6} [0.6]{\includegraphics{4-Experiment/Figures/RunDec2012/PeakShiftNaIL9.eps}}}\\ (a) Left NaI detector. \\ \\ \\ \centerline{\scalebox{0.6} [0.6]{\includegraphics{4-Experiment/Figures/RunDec2012/PeakShiftNaIR9.eps}}} \\ (b) Right NaI detector.\\ \end{tabular} \caption{The 511~keV peak observed using NaI detectors shifted when accelerator RF was on. The spectrum were taken with RF on (dotted line) and with RF off (dashed line). The solid line represents the photon spectrum created by 3~MeV positrons impinging on T2.} \label{fig:NaI-peak-shift} \end{figure}

\subsection{Positron Rate Estimation in AND Mode}

The background subtracted and normalized photon energy spectra observed using two NaI detectors for $3.00 \pm 0.07$~MeV positrons are shown in Figure~\ref{fig:pos_NaILR}. Figure~\ref{fig:pos_NaILR} (a) and (b) are the background subtracted spectrum with no coincidence or energy cut (OR mode). In OR mode, no coincidence between two detectors are required. Figure~\ref{fig:pos_NaILR} (c) and (d) illustrate events observed in coincidence and within a energy window of $511\pm75$~keV for both detectors (AND or coincidence mode). The measured positron rate using NaI detectors in AND mode was 0.25~Hz for $3.00 \pm 0.07$~MeV positrons when HRRL was operated at 300~Hz repetition rate, 100~mA peak current, and 300~ns (FWHM) macro pulse length. %The uncertainty of the power supplies used for dipoles is 0.1~A which creates 0.06~MeV uncertainty in energy of the beam bent by dipoles. \begin{figure}[htbp] \centering \begin{tabular}{cc} {\scalebox{0.74} [0.74]{\includegraphics[scale=0.40]{4-Experiment/Figures/subtract/whole/NaI_L/r3735_sub_r3736.png}}} & {\scalebox{0.74} [0.74]{\includegraphics[scale=0.40]{4-Experiment/Figures/subtract/whole/NaI_R/r3735_sub_r3736.png}}} \\ (a) & (b) \\ & \\ & \\ {\scalebox{0.74} [0.74]{\includegraphics[scale=0.40]{4-Experiment/Figures/subtract/511_peak/NaI_L/r3735_sub_r3736.png}}} & {\scalebox{0.74} [0.74]{\includegraphics[scale=0.40]{4-Experiment/Figures/subtract/511_peak/NaI_R/r3735_sub_r3736.png}}} \\ (c) & (d) \\

& \\

\end{tabular} \caption{Photon spectrum of NaI detectors after background subtraction created by 3~MeV positrons incident on T2. (a) and (b) are spectrum after background subtractions. (c) and (d) are the spectrum of events coincident in both detectors in 511~keV peaks.} \label{fig:pos_NaILR} \end{figure}

%The electron beam was transported to a phosphorous screen at the end of the 90 degree beamline to find the errors on the positron beam energy. The positron beam current was too low to be observed on phosphorous screen. %The beam centered on the phosphorous screen and then steered to the edge by changing current of D2 by 0.1~A. %The phosphorous screen is twice as large as T2 in horizantal direction.

\subsection{The Positron Production Runs} The annihilation target T2 is insertable into the beamline and placed inside a 6-way cross that has horizontal sides vacuum sealed with 1~mil stainless steel windows. Positrons intercepting T2 thermalise, annihilate, and produce 511~keV photon pairs back-to-back. These photons are detected by two NaI detectors facing T2 and shielded with Pb bricks as shown in Figure~\ref{fig:HRRL-En-Scan}. The background was measured by retracting T2 thereby allowing positrons to exit the beamline and be absorbed by the beam dump.

\begin{figure}[htbp] \centering \includegraphics[scale=0.58]{4-Experiment/Figures/PositronDetection/NaI_Setup4.png} \caption{Positron detection using T2 and NaI detectors.} \label{fig:HRRL-En-Scan} \end{figure}

As shown in Figure~\ref{fig:Dipole-in}, a 511~keV peak was observed (solid line) when T2 was in. A permanent dipole magnet was placed on the beamline after Q10 to deflect charged particles on the accelerator side preventing them from entering the shielded cell. The peak was not observed (dashed line) when T2 was in and the permanent magnet was used. The peak was also not observed (doted dashed line) when T2 and the permanent magnet were removed. Thus, one may argue that the observed peak is due to positrons annihilating in T2.

\begin{figure}[htbp] \centering \begin{tabular}{cc} {\scalebox{0.319} [0.319]{\includegraphics{4-Experiment/Figures/SweepingDipole/LNaI3.png}}} & {\scalebox{0.319} [0.319]{\includegraphics{4-Experiment/Figures/SweepingDipole/RNaI3.png}}} \\ (a) Original spectrum on left NaI. & (b) Original spectrum on right NaI.\\ \end{tabular} \caption{Photon spectrum when a permanent dipole magnet is inserted along with T2 (dashed line), dipole out and T2 in (solid line), and dipole removed and T2 out (dotted dashed line). The positron energy incident on the T2 was $2.15\pm0.06$~MeV.} \label{fig:Dipole-in} \end{figure}

The normalized photon energy spectra observed by two NaI detectors are shown in Figure~\ref{fig:in-out-runs} when T2 was both in (signal) and out (background) of the beamline. Figure~\ref{fig:in-out-runs} (c) and (d) illustrate the events observed when 511~keV photons were required in both detectors. The signal is indicated by the solid line and background by dashed line.

\begin{figure}[htbp] \centering \begin{tabular}{ll} {\scalebox{0.29} [0.29]{\includegraphics{4-Experiment/Figures/NaI_L1/r3735_sub_r3736_2.png}}} & {\scalebox{0.29} [0.29]{\includegraphics{4-Experiment/Figures/NaI_R1/r3735_sub_r3736_2.png}}} \\ (a) Original spectrum on left NaI. & (b) Original spectrum on right NaI.\\ & \\ & \\ {\scalebox{0.28} [0.28]{\includegraphics{4-Experiment/Figures/NaI_L2/r3735_sub_r3736_2.png}}} & {\scalebox{0.28} [0.28]{\includegraphics{4-Experiment/Figures/NaI_R2/r3735_sub_r3736_2.png}}} \\

(c) Spectrum with cut on left NaI. & (d) Spectrum with cut on right NaI.\\ & \\ \end{tabular} \caption{The time normalized spectra of photons created by 3~MeV positrons incident on T2. In the top row are original spectrum and in the bottom row are spectrum of incidents happened in the 511~keV peak coincidently in both detectors. The positron beam energy incident on the T2 was $3.00\pm0.06$~MeV.} \label{fig:in-out-runs} \end{figure}

%\begin{table} %\centering %\caption{Run Parameters of The Run No. 3735.} %\begin{tabular}{lll} %\toprule %{Parameter} & {Unit} & {Value} \\ %\midrule %run number & & 3735 \\ %repetition rate & Hz & 300 \\ %run time & s & 1002 \\ %pulses & & 301462 \\ %events & & 9045 \\ %e$^+$ Counts on NaI Detectors & & 256 $\pm$ 16\\ %\bottomrule %\end{tabular} %\label{tab:run3735} %\end{table}

\subsection{Positron to Electron Ratio in AND Mode} Gaussian distributions were fit to 511~keV photon peaks observed on NaI detectors as shown in Figure~\ref{fig:511Fit} and the fits result in $\sigma_{NaI, 511} = 37.5 \pm 3.0$~keV. Events in the experiment was observed within a 2$\sigma_{NaI, 511}$ energy window (436$-$586~keV) for both detectors. The ratios of positrons, detected using NaI detectors in coincidence mode, to the electrons impinging on T1 at different energies are given in Table~\ref{tab:e+2e-} and the errors on the ratio are statistical. The systematic errors are discussed in the nest section and results are plotted in Figure~\ref{fig:e+2e-}. %error on total counts = sqt(total counts). Rate = Counts/time. Error on rates = sqrt(d(rate)/d(counts)^2*(error on the counts)^2)=sqrt((error on the counts)^2/time^2)=sqrt(counts/time^2)=sqrt(rate/time). \begin{figure}[htbp] \centering \begin{tabular}{cc} {\scalebox{0.76} [0.76]{\includegraphics[scale=0.48]{4-Experiment/Figures/systimatics/NaI_L/3MeV_511Fit_r3735_sub_r3736L.eps}}} & {\scalebox{0.76} [0.76]{\includegraphics[scale=0.48]{4-Experiment/Figures/systimatics/NaI_R/3MeV_511Fit_r3737_sub_r3736R.eps}}} \\ (a) Left NaI detector.& (b) Right NaI detector.\\ & \\ \end{tabular} \caption{} \label{fig:511Fit} \end{figure}

\begin{table} \centering \caption{Positron to Electron Rate Ratio.} \begin{tabular}{cc} \toprule {Energy} & {Positron to Electron Ratio} \\ \midrule $1.02 \pm 0.03$ & $(1.7 \pm 0.6)\times10^{-16}$ \\ $2.15 \pm 0.06$ & $(8.6 \pm 1.5) \times10^{-16}$ \\ $3.00 \pm 0.07$ & $(8.52 \pm 0.54)\times10^{-15}$ \\ $4.02 \pm 0.07$ & $(3.11 \pm 0.31)\times10^{-15}$ \\ $5.00 \pm 0.06$ & $(3.32 \pm 0.89)\times10^{-15}$ \\ \bottomrule \end{tabular} \label{tab:e+2e-} \end{table}

\subsection{Positron Rate Estimation in Software OR Mode}

The ratio of 511~keV photons counted in software OR mode to the electrons incident on T1 are given in Table~\ref{tab:OrMode_e+/e-}. In these runs listed in the table, NaI detectors were operated in triple coincidence between NaI detectors and gun pulse. In the software of OR mode, there is no software coincidence are required between two NaI detectors. In OR mode, left and right NaI detectors observe different 511~keV photon rates (see section 4.4 of Chapter 4 for the detailed analysis on the count asymmetry in NaI detectors). In the experiment, the upstream side of the T2 was facing right NaI detector and right NaI detected more 511~keV photons. The asymmetry is defined as \begin{equation} \text{Asymmetry} = \frac{N_{r}-N_{l}}{N_{r}+N_{l}}\times 100\%, \end{equation} where $N_{r}$ and $N_{l}$ are the number of 511~keV photons detected by right and left NaI detectors respectively. The 511~keV photons are counted in 1, 2, and 3 $\sigma_{NaI, 511}$ ($\sigma_{NaI, 511} = 37.5 \pm 3.0$~keV) energy windows. \begin{table} \centering \caption{The Ratio of 511~keV Photons Measured in OR Mode to the Electrons Incident on T1 Measured in Different $\sigma_{NaI, 511}$ Energy Windows.} \begin{tabular}{cccc} \toprule 1 $\sigma_{NaI, 511}$ Cut \\ \midrule {Energy (MeV)} & {e$^{+}$/e$^{-}$ on NaI Right} & {e$^{+}$/e$^{-}$ on NaI Left} & {Asymmetry} \\ \midrule $1.02 \pm 0.03$ & $( 2.5 \pm 0.3 ) \times10^{-15}$ & $( 0.0 \pm 0.0 ) \times10^{-15}$ & $ 96.0 \pm 2.8 \%$ \\ $2.15 \pm 0.06$ & $( 8.5 \pm 0.6 ) \times10^{-15}$ & $( 1.0 \pm 0.2 ) \times10^{-15}$ & $ 79.2 \pm 4.0 \%$ \\ $3.00 \pm 0.07$ & $( 75.1\pm 1.5 ) \times10^{-15}$ & $( 14.2\pm 0.6 ) \times10^{-15}$ & $ 68.2 \pm 1.3 \%$ \\ $4.02 \pm 0.07$ & $( 42.0\pm 1.2 ) \times10^{-15}$ & $( 6.0 \pm 0.4 ) \times10^{-15}$ & $ 75.1 \pm 1.7 \%$ \\ $5.00 \pm 0.06$ & $( 21.5\pm 2.3 ) \times10^{-15}$ & $( 4.0 \pm 1.0 ) \times10^{-15}$ & $ 68.5 \pm 7.0 \%$ \\ \midrule

2 $\sigma_{NaI, 511}$ Cut  \\

\midrule {Energy (MeV)} & {e$^{+}$/e$^{-}$ on NaI Right} & {e$^{+}$/e$^{-}$ on NaI Left} & {Asymmetry} \\ \midrule $1.02 \pm 0.03$ & $( 2.9 \pm 0.3 ) \times10^{-15}$ & $( 0.3 \pm 0.1 ) \times10^{-15}$ & $ 81.5 \pm 5.2 \%$ \\ $2.15 \pm 0.06$ & $( 11.8 \pm 0.7 ) \times10^{-15}$ & $( 1.4 \pm 0.2 ) \times10^{-15}$ & $ 79.0 \pm 3.4 \%$ \\ $3.00 \pm 0.07$ & $( 98.0 \pm 1.7 ) \times10^{-15}$ & $( 19.0 \pm 0.7 ) \times10^{-15}$ & $ 67.6 \pm 1.2 \%$ \\ $4.02 \pm 0.07$ & $( 53.9 \pm 1.3 ) \times10^{-15}$ & $( 8.1 \pm 0.5 ) \times10^{-15}$ & $ 73.9 \pm 1.5 \%$ \\ $5.00 \pm 0.06$ & $( 28.6 \pm 2.6 ) \times10^{-15}$ & $( 5.5 \pm 1.1 ) \times10^{-15}$ & $ 67.6 \pm 6.1 \%$ \\ \midrule

3 $\sigma_{NaI, 511}$ Cut  \\

\midrule {Energy (MeV)} & {e$^{+}$/e$^{-}$ on NaI Right} & {e$^{+}$/e$^{-}$ on NaI Left} & {Asymmetry} \\ \midrule $1.02 \pm 0.03$ & $( 3.4 \pm 0.3 ) \times10^{-15}$ & $( 0.4 \pm 0.1 ) \times10^{-15}$ & $ 78.4 \pm 5.1 \%$ \\ $2.15 \pm 0.06$ & $( 13.4\pm 0.7 ) \times10^{-15}$ & $( 1.5 \pm 0.2 ) \times10^{-15}$ & $ 79.3 \pm 3.1 \%$ \\ $3.00 \pm 0.07$ & $( 106.3\pm1.8 ) \times10^{-15}$ & $( 21.4\pm 0.8 ) \times10^{-15}$ & $ 66.5 \pm 1.1 \%$ \\ $4.02 \pm 0.07$ & $( 58.9\pm 1.4 ) \times10^{-15}$ & $( 9.4 \pm 0.5 ) \times10^{-15}$ & $ 72.5 \pm 1.5 \%$ \\ $5.00 \pm 0.06$ & $( 16.4\pm 1.4 ) \times10^{-15}$ & $( 2.9 \pm 0.6 ) \times10^{-15}$ & $ 70.0 \pm 5.6 \%$ \\ \bottomrule \end{tabular} \label{tab:OrMode_e+/e-} \end{table}

The errors given in Table~\ref{tab:OrMode_e+/e-} are statistical. Counting positrons in AND mode underestimates the positron rates by factor of 10 approximately compared to the OR mode. The asymmetry tends to decrease as the positron energy increases. The average ratios are given in Table~\ref{tab:Sys_OrMode_e+/e-}. In the table, the first error on ratios are systematic and second ones are statistical.

\begin{sidewaystable} \centering \caption{The Ratio of 511~keV Photons Measured in OR Mode to the Electrons Incident on T1.} \begin{tabular}{cccc} \toprule {Energy (MeV)} & {e$^{+}$/e$^{-}$ on NaI Right} & {e$^{+}$/e$^{-}$ on NaI Left} & {Asymmetry} \\ \midrule

$1.02 \pm 0.03$ & $(2.9\pm 0.5 \pm0.3) \times10^{-15}$ & $(0.3 \pm 0.2 \pm 0.1) \times10^{-15}$ & $85.3\pm 9.4\pm5.2\%$ \\ $2.15 \pm 0.06$ & $(11.2\pm2.5 \pm0.7) \times10^{-15}$ & $(1.3 \pm 0.3 \pm 0.2) \times10^{-15}$ & $79.2\pm 0.2\pm3.4\%$ \\ $3.00 \pm 0.07$ & $(93.1\pm16\pm1.7) \times10^{-15}$ & $(18.2\pm 3.7 \pm 0.7) \times10^{-15}$ & $67.4\pm 0.9\pm1.2\%$ \\ $4.02 \pm 0.07$ & $(51.6\pm8.7 \pm1.3) \times10^{-15}$ & $(7.8 \pm 1.7 \pm 0.5) \times10^{-15}$ & $73.8\pm 1.3\pm1.5\%$ \\ $5.00 \pm 0.06$ & $(22.2\pm6.1 \pm2.6) \times10^{-15}$ & $(4.1 \pm 1.3 \pm 1.1) \times10^{-15}$ & $68.7\pm 1.2\pm6.1\%$ \\

\bottomrule \end{tabular} \label{tab:SysOrMode_e+/e-} \end{sidewaystable}

\subsection{Positron Rate Estimation in Hardware OR Mode} The positron rate was also measured by running NaI detectors in OR mode and the photon spectrum are shown in Figure~\ref{fig:Or-mode} for $2.15 \pm 0.06$~MeV positrons incident on T2. Two NaI detectors were operated in hardware OR mode in these runs. The rate was $0.21$~Hz for left and $0.35$~Hz right NaI detectors in OR mode. However, when the software coincidence are required between two detectors by only drawing the events happened in 511~keV peaks of both detectors, the rate dropped $0.028$~Hz.

\begin{figure}[htbp] \centering \begin{tabular}{cc} {\scalebox{0.76} [0.76]{\includegraphics[scale=0.40]{4-Experiment/Figures/NaIOrMode/r3679_sub_r3680_NaI_L.png}}} & {\scalebox{0.76} [0.76]{\includegraphics[scale=0.40]{4-Experiment/Figures/NaIOrMode/r3679_sub_r3680_NaI_R.png}}} \\ (a) Left NaI detector.& (b) Right NaI detector.\\ & \\ \end{tabular} \caption{The photon rates for $2.15 \pm 0.06$~MeV incident positrons measured by running NaI detectors in OR mode. The rate was $0.21$~Hz for left and $0.35$~Hz right NaI detectors in OR mode while the rate in coincidence (AND) mode was $0.028$~Hz.} \label{fig:Or-mode} \end{figure}

\section{Systematic Error Analysis for AND Mode Rate} The systematic errors in the positron production experiment could be introduced by the uncertainty in the magnetic fields and misalignment of the quadrupoles. These systematic errors are discussed in section 6 of Chapter 4.%The uncertainty of the emittances and Twiss parameters of the electron beam would also create systematic error as well. The positron and electron rates are normalized by time which can introduce systematic error as well. However, the uncertainty in run duration time is very small (0.1$-$0.2$\%$) and can be ignored (the duration of positron production runs are around 500$\sim$1000 seconds and the uncertainty is about 1 second). The uncertainty in the 511~keV peaks detected by the NaI detectors would also contribute to the systematic error.

%\subsection{Systematic Error Introduced by the Uncertainty of The 511~keV Peak} Gaussian distributions were fit to the 511~keV photon peaks observed on NaI detectors as shown in Figure~\ref{fig:511Fit} and the fits result in $\sigma_{NaI, 511} = 37.5 \pm 3.0$~keV. Events in the experiment was observed in coincidence and within a 2$\sigma_{NaI, 511}$ energy window (436$-$586~keV) for both detectors. The systematic error of the 511~keV photon counts introduced by the uncertainty in the 511~keV peak is studied by counting the photons with 1$\sigma_{NaI, 511}$ (473.5$-$548.5~keV) and 3$\sigma_{NaI, 511}$ (398.5$-$623.5~keV) energy windows. The positron to electron rate ratio for different energy windows are given in Table~\ref{tab:sys2}. The results are plotted in Figure~\ref{fig:e+2e-}.

\begin{sidewaystable} \centering \caption{Positron to Electron Rate Ratio: Systematic Error Introduced by The Uncertainty of The 511~keV Peak.} \begin{tabular}{llll} \toprule {Energy} & {Positron to Electron Ratio} & {Positron to Electron Ratio} & {Positron to Electron Ratio} \\ {(MeV)} & {Energy Window: 2$\sigma_{NaI, 511}$} & {Energy Window: 3$\sigma_{NaI, 511}$} & {Energy Window: 1$\sigma_{NaI, 511}$} \\ \midrule $1.02 \pm 0.03$ & $(1.7 \pm 0.6) \times10^{-16}$ & $(1.7 \pm 0.6) \times10^{-16}$ & $(0.5 \pm 0.4) \times10^{-16}$ \\ $2.15 \pm 0.06$ & $(8.6 \pm 1.5) \times10^{-16}$ & $(9.2 \pm 2.0) \times10^{-16}$ & $(3.0 \pm 0.9) \times10^{-16}$ \\ $3.00 \pm 0.07$ & $(8.52 \pm 0.54)\times10^{-15}$ & $(9.90 \pm 0.58)\times10^{-15}$ & $(4.20 \pm 0.35)\times10^{-15}$ \\ $4.02 \pm 0.07$ & $(3.11 \pm 0.31)\times10^{-15}$ & $(3.66 \pm 0.34)\times10^{-15}$ & $(2.07 \pm 0.26)\times10^{-15}$ \\ $5.00 \pm 0.06$ & $(3.32 \pm 0.89)\times10^{-15}$ & $(3.55 \pm 0.92)\times10^{-15}$ & $(1.66 \pm 0.63)\times10^{-15}$ \\

\bottomrule \end{tabular} \label{tab:sys2} \end{sidewaystable}

\begin{figure}[htbp] \centering \includegraphics[scale=0.84]{4-Experiment/Figures/Ratio/R.eps} \caption{The ratios of positrons detected by NaI detectors in coincidence mode to the electrons impinging T1. The solid error bars statistical and the dashed ones are systematic.} \label{fig:e+2e-} \end{figure}

Chapter 4: Simulation File:Sadiq thesis chapt 4.txt

\chapter{Simulation} The ratios of positrons to the electrons impinging the tungsten production target, T1, was substantially smaller after they traversed the beam line than what was produced at the target (10$^{-15}$ instead of 10$^{-3}$). Simulations were performed using G4beamline to better understand the losses from the transportation of positrons to the experimental cell. ``G4beamline is a particle tracking and simulation program based on the GEANT4~\cite{geant4} toolkit that is specifically designed to easily simulate beamlines and other systems using single-particle tracking~\cite{muonsinc}. A sample G4beamline script for positron generation using the new HRRL beamline is given in the appendix F.

The simulation predicts that at least one positron per 1000 incident 10~MeV electrons is produced using a 2mm thick tungsten target. The simulation revealed that the number of electrons decreased by orders of magnitude as they were transported through the beamline magnets. As a result of this beam loss, the simulation was divided into three steps to increase the beam line simulation efficiency. Each step generates particles at different locations along the beamline where the beam loss was found to be substantial. While the first simulation step used the measured electron energy profile, subsequent steps would generate particles based on the particle phase space observed at the end of previous step. The method decreased the simulation time so a sample of more than one million events could be procuded within a single day.

The first step in the simulation generated an electron beam with the energy distribution observed in the experiment, see Figure~\ref{fig:En-Scan}. The electrons were focused by three quadrupoles onto the positron production target T1 (see Figure~\ref{fig:app-hrrl-line}). Electrons traversing T1 produced bremsstrahlung photons of sufficient energy to produce $e^+e^-$ pairs that would escape the downstream side of the target and be collected by a second quadrupole triplet. The second step simulated the collection and transportation of positrons exiting T1 to the entrance of the first dipole D1. The last step transported positrons from the entrance of D1 all the way to the annihilation target T2, the interactions of positrons with T2, and the detection of the resulting 511~keV photon pairs.

The two targets (T1 and T2) were positioned at different angles with respect to the incident beam momentum vector. The beamline coordinate system aligns the "z" axis to point along the incident electrons momentum vector and the "x" axis along the horizontal plane. The positron production target, T1, was placed such that the upstream side of T1 was facing vertically down at a 45$^{\circ}$ angle with respect to the vertical(rotated 45$^{\circ}$ counter clockwise about x-axis). The positron conversion target, T2, was rotated about two different axes. The first rotation positioned the upstream side of the target so it was facing down by 45$^{\circ}$ (rotated 45$^{\circ}$ clockwise about x-axis like T1). The second rotation was 45$^{\circ}$ clockwise about y-axis. As a result, target T2's upstream face was directed towards the beam left NaI detector.

\section{Step 1 - The Electron Beam Generation and Transportation to T1} In the first simulation step, an electron beam was generated with an energy distribution that was observed in the experiment. The emittance, the Twiss parameters, and the energy distribution of the electron beam were measured experimentally. The energy distribution of the electron beam is shown in Figure~\ref{fig:En-Scan}. The distribution was fit using two skewed Gaussian distributions. The fit parameters given in Table~\ref{tab:En-Scan_resluts} were used by the simulation to generated electrons.

A series of virtual detectors were placed along the beamline to sample the beam. As an example, three virtual circular detectors and T1 are shown in Figure~\ref{fig:T1_UpD_DwD2}. The electron beam was detected by a virtual detector DUPT1 (Detector UPstream of T1) placed 25.52~mm upstream of T1. Positrons, electrons, and photons generated during the interaction of the electron beam with T1 were observed by virtual detectors DT1 (Detector of T1) and DDNT1 (Detector DowNstream of T1) placed 25.52~mm downstream of T1. %In Figure~\ref{fig:SimS1_pos_En_DDNT1}, $13.8 \times 10^{10}$ electrons shot at T1 and generated positrons positrons shown in blue.

\begin{figure}[htbp] \centering \includegraphics[scale=0.55]{3-Simulation/Figures/sim_setup_T1_UpD_DwD2.png} \caption{T1 is the positron production target. DUPT1 is a virtual detector located upstream of T1 to detect the incoming electron beam. DDNT1 is a virtual detector downstream of T1. DT1 is a virtual detector that is placed right after T1 and parallel to it.} \label{fig:T1_UpD_DwD2} \end{figure}


\subsection{The Positron Beam on DDNT1}

In the first step, $1.38 \times 10^{10}$ electrons were generated with the energy distribution shown by the dotted-dashed line in Figure~\ref{fig:SimS1_T1UPDN}. These electrons were transported to T1 where they produced photons that produced the positron distribution shown by the solid the line in Figure~\ref{fig:SimS1_T1UPDN} by pair-production. The dashed line is the electron energy distribution observed by DDNT1. The simulation result using $1.38 \times 10^{7}$ electrons incident on T1 is drawn in Figure~\ref{fig:SimS1_T1UPDN}. The incident electrons were detected by the virtual detector DUPT1 and downstream positrons and electrons were detected by DDNT1. As shown in Figure~\ref{fig:SimS1_T1UPDN}, the electrons pass through T1 loosing approximately 4 MeV while the positrons escape the downstream side of T1 with a mean energy of about 3 MeV. The beam line was set to transport positrons with this mean energy as a result of this prediction.

\begin{figure}[htbp] \centering \includegraphics[scale=0.75]{3-Simulation/Figures/s/s1/overlay8.eps} \caption{The incident electron energy distribution (dotted dashed line), the distribution of electrons after T1 (dashed line), and the distribution of positrons produced(solid line). The incident electron distribution counts were weighted by 0.001.} \label{fig:SimS1_T1UPDN} \end{figure}

\begin{figure}[htbp] \begin{tabular}{cc} {\scalebox{0.364} [0.364]{\includegraphics{3-Simulation/Figures/X_e+_DDNT1.eps}}} & {\scalebox{0.364} [0.364]{\includegraphics{3-Simulation/Figures/Y_e+_DDNT1.eps}}} \\ (a) $x$ projection. & (b) $y$ projection. \\ {\scalebox{0.364} [0.364]{\includegraphics{3-Simulation/Figures/XP_e+_DDNT1.eps}}} & {\scalebox{0.364} [0.364]{\includegraphics{3-Simulation/Figures/YP_e+_DDNT1.eps}}} \\ (c) $x$ projection of the divergence. & (d) $y$ projection of the divergence.\\ {\scalebox{0.364} [0.364]{\includegraphics{3-Simulation/Figures/XY_e+_DDNT1.png}}} & {\scalebox{0.364} [0.364]{\includegraphics{3-Simulation/Figures/XY_e+_DDNT1_zoom.png}}} \\ (e) The transverse beam profile. & (f) Zoom in of (e). \\ \end{tabular} \caption{The transverse beam projections and angular distributions of positrons detected. The virtual detector diameter is 30~mm while the beam pipe is only 24 mm.} \label{fig:DDNT1_results} \end{figure}

The positron spatial and angular distribution detected by DDNT1 is shown in Figure~\ref{fig:DDNT1_results}. The $y$ $vs.$ $x$ spatial distribution of the beam is shown in Figure~\ref{fig:DDNT1_results}~(e) and Figure~\ref{fig:DDNT1_results} (f). As can be seen from Figure~\ref{fig:DDNT1_results} (b) and (d), the $y$ spatial distribution and divergence, defined in equation~\ref{eq:divergence}, of the positron beam have a sharp drop in counts in the region between $-25.8$~mm and $-27.2$~mm from the beam center. Figure~\ref{fig:sim-DDNT1-T1-geo} shows the geometry and location of T1 and DDNT1. If the size of T1 were to be increased, it would eventually intersect with DDNT1 at a distance between $25.8$~mm and $27.2$~mm from the beam center, $i.e.$ the edge of the T1 is facing this $1.4$~mm wide low count area. %This is the result of the target's thickness of 1.016 mm and the 45$^{\circ}$ angle of intersection ($1.016\sqrt{2}=1.44$). The edge of the target does not produce many positrons compared to the face of the target. \begin{figure}[htbp] \centering \includegraphics[scale=1]{3-Simulation/Figures/sharp_drop2.eps} \caption{The geometry of the target T1 and the virtual detector DDNT1.} \label{fig:sim-DDNT1-T1-geo} \end{figure}

As shown in Figure~\ref{fig:DDNT1_YpY}, the $y$ distribution count decreases at $\theta$ = 45$^{\circ}$. Positrons were emitted from both the downstream and upstream side of T1. Positrons from the downstream side of T1 intersected the detector at angles below 45$^{\circ}$ while positrons from the upstream side of T1 begin to hit the detector at angles beyond 45$^{\circ}$. Neither positrons upstream nor downstream of T1 traveled to the 1.4~mm wide low count area. Only positrons created on the edge of T1 reached the low count area between 25.8~mm~$<x<$~27.2~mm. As a result, the counts in this area are comparatively lower.

\begin{figure}[htbp] \begin{tabular}{cc} {\scalebox{0.39} [0.39]{\includegraphics{3-Simulation/Figures/STSimS1_YYP_DDNT1.png}}} & {\scalebox{0.39} [0.39]{\includegraphics{3-Simulation/Figures/STSimS1_YYP_DDNT1_zoom.png}}} \\ (a) y' $vs.$ y. & (b) y' $vs.$ y zoom. \\ \end{tabular} \caption{The positron beam $y'$ $vs$. $y$ detected by DDNT1.} \label{fig:DDNT1_YpY} \end{figure}

\subsection{The Positron Beam on Virtual Detectors DQ4 and DD1}

The positron beam energy distributions at the entrance of Q4 observed by virtual detector DQ4 is shown in Figure~\ref{STSimS1_En_DQ4_DD1}. The distribution of positrons energies observed leaving the downstream side of the production target T1 as observed by virtual detector DD1UP are also shown. The ratio of positrons leaving the target T1 and entering the first quadrupole indicate that nearly 90\% of the positrons are lost. The large positron divergence is responsible for this loss. The median positron energy does not change suggesting that the beamline simulation is tuned to transport the optimal number of positrons emitted by the positron production target, T1.

\begin{figure}[htbp] \begin{tabular}{ccc} \centerline{\scalebox{0.8} [0.8]{\includegraphics{3-Simulation/Figures/En_e+_DQ4.eps}}} \\ (a) The positron energy distribution on DQ4. \\ \centerline{\scalebox{0.8} [0.8]{\includegraphics{3-Simulation/Figures/En_e+_DD1.eps}}}\\ (b) The positron energy distribution on DD1. \\ \end{tabular} \caption{The positron beam energy distribution detected downstream of T1 just before entering the first quadrupole (a) and after leaving the production target upstream of this quadrupole (b).} \label{STSimS1_En_DQ4_DD1} \end{figure}

\section{Step 2 - The Transportation of The Positron Beam from DDNT1 to The Entrance of The First Dipole}

The simlation's second step was designed to predict the amount of beam loss that would result when the postiron beam is transported from the entrace of Q4 to the entrance of the first energy selecting dipole D1. The beam observed on detector DDNT1 in step 1 was used to generate positrons directed towards the entrance of the quadrupole Q4. At detector DDNT1, the higher energy positrons tend to have smaller polar angles and are closer to the beam center. Positrons were generated in 1~keV/c momentum bins with different weights, spatial, and angular distributions to reproduce the distributions observed in step 1. Figure~\ref{fig:STSimSetupS2} illustrates the quadurpole triplet system and the dipole magnet D1 used in the beamline. Positrons generated at DDNT1 are transported to the entrance of D1 through this quadrupole triplet system. The virtual detectors were placed at the entrance of Q4 (DQ4) and D1 (DD1UP) to track the positrons. The number of positrons is predicted to decrease by a factor of ten as they travel from through the quad triple system to the entrance of the the first Dipole D1, see Table 4.2.

\begin{figure}[htbp] \centering \includegraphics[scale=0.75]{3-Simulation/Figures/STSimSetupS2_3.png} \caption{The generation and transportation of the positron beam in step 2. The virtual detectors were used to track the positrons.} \label{fig:STSimSetupS2} \end{figure}

\section{Step 3 - The Transportation of Positrons from the Entrance of The First Dipole to T2 and The Detection of 511~keV Photons}

In this step, the positron beam was generated at the entrance of the first dipole D1 and observed using virtual detector DD1UP. The beam was deflected $45^{\circ}$ by D1 and passed through the energy slit. The postrons diverged along the horizon according to their energy when traversing D1's magnetic field. The energy slit constrained the spatial distribution of the beam by blocking the beam using a 34.8~mm wide gap oriented vertically. A quadrupole Q7, located after the energy slit, was set to 3.5~A (similar to the experiment) in order to focus the beam before it enters a second dipole D2. The beam was deflected another $45^{\circ}$ by D2 and sent through three quadrupoles towards the annihilation target T2 located at the end of the beamline as shown in Figure~\ref{fig:T2}.

\begin{figure}[htbp] \centering \includegraphics[scale=0.4]{3-Simulation/Figures/HRRL_T2_2.png} \caption{T2 and virtual detectors located upstream (DT2UP) and downstream (DT2DN) of T2 are shown at the center of the figure. NaI detectors and Pb shielding are located horizontally at two sides.} \label{fig:T2} \end{figure}

As shown in Figure~\ref{fig:T2}, two virtual circular detectors DT2UP and DT2DN with a 48~mm diameter (48~mm is the inner diameter of the beam pipe) were placed upstream and downstream of T2 to detect positrons. Two other virtual detectors, DT2L and DT2R (not shown) with the same diameters as the annihilation target T2 were placed on the left and right side of the beam and parallel to T2 to detect positrons. Two additional virtual detectors were placed horizontally at the locations of NaI detectors, which were 170~mm away from the beamline center, to detect photons. Two-inch-thick Pb bricks with 2-inch diameter circular openings were positioned between T2 and the NaI virtual detectors. Each detector was surrounded by 2$$ of Pb. When a positron annihilates inside T2, two back-to-back scattered 511~keV photons are generated. An event is registered in the simulation when both NaI detectors observe a 511~keV photon.


\subsection{Positrons Detected by The Detection System} The photons observed by each NaI detector are subject to the detector efficiency. The detector efficiency chart (shown in Figure~\ref{fig:NaI_Ef}), obtained from Saint-Gobain Crystals~\cite{NaI-Eff}, indicates that the NaI crysta,l used in this experiment, has an efficiency of 68\% for 511~keV photons. If two detectors are operated in coincidence mode, the detection efficiency of the system is $68\% \times 68\%~=~46.24\%$.

Figure~\ref{fig:e+_Generated_and_Detected} shows the number of 511~keV photon pairs detected in coincidence mode (multiplied by 46.24\%) and overlaid with the positrons detected on DDNT1. In step 3 of the simulation, the quadrupole current for Q7 was simulated for 0, 3.5, and 10~A to study the effect of Q7 on the positron energy distribution. As shown in Figure~\ref{fig:e+_Generated_and_Detected}, the shape of the predicted positron energy distribution more closely resembling the oen observed in the experiment when Q7 was set to the value used in the actual experiment, 3.5~A. Fewer positrons are predicted to reach T2 when the Q7 currents are 0 and 10~A. The simulation's prediction for the experiment based on the beamline parameters used during the experiment are given in Table~\ref{tab:Sim_e+/e-} as the ratio of the number of e$^+$ detected in coincidence mode by NaI detectors to the number of e$^-$ incident on the T1 for the measured positron energies. The values in Table~\ref{tab:Sim_e+/e-} are compared with the experiment in the next chapter.

\begin{figure}[htbp] \centering \includegraphics[scale=0.53]{3-Simulation/Figures/NaI_Ef_3.png} \caption{NaI detector efficiency obtained from Saint-Gobain Crystals~\cite{NaI-Eff}. The lines corresponding to the different crystal sizes (in inches) are shown on the right side of the figure.} \label{fig:NaI_Ef} \end{figure}

\begin{figure}[htbp] \centering \includegraphics[scale=0.77]{3-Simulation/Figures/overlay/e+_Generated_and_Detected.eps} \caption{Positrons detected on virtual detector DDNT1 and 511~keV photon pairs detected by the NaI detectors in coincidence mode when Q7 set as at 0~A, 3.5~A and 10~A.} \label{fig:e+_Generated_and_Detected} \end{figure}

\begin{table} \centering \caption{Positron to Electron Ratio Estimated by the Simulation.} \begin{tabular}{ll} \toprule {Energy (MeV)} & {Positron to Electron Ratio} \\ \midrule $1.02 \pm 0.03$ & $(8.90 \pm 5.1) \times10^{-17}$ \\ $2.15 \pm 0.06$ & $(10.10 \pm 0.06) \times10^{-15}$ \\ $3.00 \pm 0.07$ & $(9.96 \pm 0.06) \times10^{-15}$ \\ $4.02 \pm 0.07$ & $(9.67 \pm 0.03) \times10^{-15}$ \\ $5.00 \pm 0.06$ & $(7.80 \pm 0.05) \times10^{-15}$ \\ \bottomrule \end{tabular} \label{tab:Sim_e+/e-} \end{table}


The remaining sections of this chapter use the above simulation to investigate beam loss, an observed asymmetry in the photon detectors, and systematic errors. Understanding the source of the beam line through the beam line was critical to determining the veracity of the data. The large discrepancy of nine order of magnitude between the observed positron production rate and the expected rate was difficult to believe. As shown in the sections below, the poor beam quality from a medical linac combined with a poor positron collection system had a large impact on this efficiency. There is also an indication that the beam line alignment also plays an important role than one might naively expect. The follwing sections attempt to explains these effects as well as investigate their impact on the systematic error.

\subsection{Beam Loss Study} The three step simulation method described in the previous section was used to predict the amount of beam loss at several locations along the beam line. Table~\ref{tab:sim-S2E} quantifies these predictions. In particular, columns 2 and 3 in Table~\ref{tab:sim-S2E} represent the number of positrons that would need to be generated in order to observe the number of positrons in the remaining columns. The number of 511~keV photon pairs detected by NaI detectors in coincidence mode is shown in the last column of Table~\ref{tab:sim-S2E}. The number of positrons observed by virtual detectors which were placed along the beamline and 511~keV photon pairs detected are shown in Figure~\ref{fig:TransEff}.

\begin{sidewaystable} \centering \caption{Simulation of $7.253 \times 10^{16}$ Electrons Incident on the T1: The Number of Positrons Transported and the Number of 511~keV Photons Detected in Coincidence Mode.} \begin{tabular}{lccccccccc} %\begin{tabular}{lllllllllll} \toprule Energy & e$^{+}$ on & e$^{+}$ enter & e$^{+}$ enter & e$^{+}$ exit & e$^{+}$ enter & e$^{+}$ enter & e$^{+}$ exit & e$^{+}$ reach & 511~keV $\gamma$ \\ (MeV) & DDNT1 & Q4 & D1 & D1 & Q7 & D2 & D2 & T2 & by NaI \\

&  &  &  &  &  &   &   &   &  detectors \\

\midrule $1.02 \pm 0.25$ & $2.5 \times 10^{12} $ & $1.4 \times 10^{10} $ & $1.4 \times 10^{9} $ & $1.2 \times 10^{8}$ & $2.6 \times 10^{7} $ & $2.7 \times 10^{6} $ & $1.2 \times 10^{6} $ & $4.3 \times 10^{3} $ & $7$\\ $1.50 \pm 0.25$ & $4.9 \times 10^{12} $ & $2.7 \times 10^{10} $ & $2.8 \times 10^{9} $ & $1.8 \times 10^{8}$ & $9.1 \times 10^{7} $ & $2.7 \times 10^{7} $ & $1.2 \times 10^{7} $ & $5.1 \times 10^{4} $ & $97 $\\ $2.15 \pm 0.25$ & $6.5 \times 10^{12} $ & $3.8 \times 10^{10} $ & $3.8 \times 10^{9} $ & $3.8 \times 10^{8}$ & $1.7 \times 10^{8} $ & $7.8 \times 10^{7} $ & $3.6 \times 10^{7} $ & $4.4 \times 10^{5} $ & $734 $\\ $2.50 \pm 0.25$ & $6.8 \times 10^{12} $ & $4.1 \times 10^{10} $ & $4.1 \times 10^{9} $ & $4.4 \times 10^{8}$ & $2.3 \times 10^{8} $ & $1.1 \times 10^{8} $ & $5.0 \times 10^{7} $ & $4.6 \times 10^{5} $ & $794 $\\ $3.00 \pm 0.25$ & $6.6 \times 10^{12} $ & $4.1\times 10^{10} $ & $4.1 \times 10^{9} $ & $4.6 \times 10^{8}$ & $2.7 \times 10^{8} $ & $1.3 \times 10^{8} $ & $6.1 \times 10^{7} $ & $4.5 \times 10^{5} $ & $723 $\\ $3.50 \pm 0.25$ & $6.1 \times 10^{12} $ & $3.9 \times 10^{10} $ & $3.9 \times 10^{9} $ & $4.4 \times 10^{8}$ & $2.8 \times 10^{8} $ & $1.4 \times 10^{8} $ & $6.7 \times 10^{7} $ & $4.4 \times 10^{5} $ & $716 $\\ $4.02 \pm 0.25$ & $5.3 \times 10^{12} $ & $3.5 \times 10^{10} $ & $3.5 \times 10^{9} $ & $4.1 \times 10^{8}$ & $2.7 \times 10^{8} $ & $1.4 \times 10^{8} $ & $6.9 \times 10^{7} $ & $4.3 \times 10^{5} $ & $702 $\\ $4.50 \pm 0.25$ & $4.6 \times 10^{12} $ & $3.1 \times 10^{10} $ & $3.1 \times 10^{9} $ & $3.6 \times 10^{8}$ & $2.5 \times 10^{8} $ & $1.3 \times 10^{8} $ & $6.7 \times 10^{7} $ & $4.0 \times 10^{5} $ & $646 $\\ $5.00 \pm 0.25$ & $3.8 \times 10^{12} $ & $2.7 \times 10^{10} $ & $2.7 \times 10^{9} $ & $3.1 \times 10^{8}$ & $2.2 \times 10^{8} $ & $1.2 \times 10^{8} $ & $6.2 \times 10^{7} $ & $3.6 \times 10^{5} $ & $566 $\\ $5.50 \pm 0.25$ & $3.0 \times 10^{12} $ & $2.2 \times 10^{10} $ & $2.2 \times 10^{9} $ & $2.6 \times 10^{8}$ & $1.9 \times 10^{8} $ & $1.0 \times 10^{8} $ & $5.6 \times 10^{7} $ & $3.3 \times 10^{5} $ & $533 $\\ \bottomrule \end{tabular} \label{tab:sim-S2E} \end{sidewaystable}

\begin{figure}[htbp] \centering \includegraphics[scale=0.75]{3-Simulation/Figures/Transporation_Efficiency/Ef.eps} \caption{Predicted number of positrons transported. Black cube: positrons incident on DDNT1. Red cube: positrons entered Q4. Blue cube: positrons entered D1. Magenta cube: positrons exited D1. Black circle: positrons entered Q7. Red circle: positrons entered D2. Blue circle: positrons exited D2. Magenta circle: positrons incident on DT2UP. Black triangle: 511~keV photons detected by NaI detectors in coincidence mode.} \label{fig:TransEff} \end{figure}

The simulation's beam loss prediction between the positron production target, T1, and the first collection quadrupole, Q4, may be contrasted with a simple solid angle argument to demonstrate the veracity of the prediction. As shown in Table~\ref{tab:sim-S2E}, the positron energy distribution is divided into 10 bins. Positrons were detected at both detector DDNT1 (28.5~mm downstream T1) and detector DQ4 (484.4 mm downstream T1). The ratio of positrons detected by DDNT1 to the ones detected by DQ4 is about 157. The solid angle that the entrance of Q4 makes with respect to T1, if one assumes that T1 is a point source of positrons is about 0.6 $\pi$~steradians. The distance between T1 and the virtual detector DQ4 is $484.4$~mm and the radius of DQ4 is $24$~mm. The solid angle of the DQ4 is $\Omega_{\text{Q4}}=\frac{\pi r^2}{d^2}=\frac{\pi 24^2}{484.4^2}$~steradian. Positrons make up a cone with a 45$^\circ$ half angle, which is $\Omega_{\text{beam}}=0.6\pi$~steradian in solid angle. The ratio of the two solid angles, $\Omega_{\text{Q4}}/\Omega_{\text{beam}}$, is 1:244, $i.e.$ 1 out of 244 positrons makes it from T1 to DQ4, assuming that the positron beam is isotropic inside the cone. However, the positron beam peaks at a smaller angle as shown in Figure~\ref{fig:DDNT1_results} (c) and (d), and as a result more positrons are transported from DDNT1 to DQ4 which may result in a ratio closer to 1:157.

%In a simulation where the dipoles were set to transport 3~MeV positron, the beam energy distribution at the exit of D1, entrance %of Q7, exit of Q7, entrance of D2, exit of D2, and T2 are shown in Figure~\ref{fig:dipole-trans} and Table~\ref{tab:pos-beam-loss}. %In Table~\ref{tab:pos-beam-loss}, the relative counts reported in the last column are obtained by dividing the second %column over 723, the number of 511~keV photons detected by NaI detectors in coincidence mode.

The first dipole is another region where substantial beam loss is predicted by the simulation. The positron counts dropped one order of magnitude when positrons are transported from the entrance of D1 to the exit of D1. This could be explained by beam scraping on the vacuum chamber. The width of the dipole chamber is 18~mm and the beam pipe diameter is 48~mm. Positrons scraping on the top and bottom of the chamber would be lost. When the dipole was set to transport 3~MeV positrons, positrons having an energy range between 2.8~MeV and 3.3~MeV are transported to the exit D1 as shown in Figure~\ref{fig:dipole-trans}. This small energy range of 0.5 MeV is only one tenth of the full range subtended somewhat uniformly by the positrons and could easily explain the order of magnitude drop in the number of positrons that travers the first dipole.

\begin{table} \centering \caption{The Positron Beam Loss along the Beamline When Dipoles Set to Transport 3~MeV Positrons.} \begin{tabular}{lcc} \toprule {Beam Sample} & {Absolute Counts} & {Relative Counts} \\ {Locations} & {} & {} \\ \midrule e$^+$ on DDNT1 & $6.6\times10^{12}$ & $9.4\times10^{10}$ \\ e$^+$ enter Q4 & $4.1\times10^{10}$ & $5.8\times10^{7}$ \\ e$^+$ enter D1 & $4.1\times10^{9}$ & $5.8\times10^{6}$ \\ e$^+$ exit D1 & $4.6\times10^{8}$ & $6.6\times10^{5}$ \\ e$^+$ enter Q7 & $2.7\times10^{8}$ & $3.8\times10^{5}$ \\ e$^+$ enter D2 & $1.3\times10^{8}$ & $1.9\times10^{5}$ \\ e$^+$ exit D2 & $6.1\times10^{7}$ & $8.7\times10^{4}$ \\ e$^+$ on T2 & $4.4\times10^{5}$ & $6.3\times10^{2}$ \\ e$^+$ annihilated in T2 & $2.5\times10^{5}$& $3.6\times10^{2}$ \\ 511~keV photons on NaI & 723 & 1 \\

\bottomrule \end{tabular} \label{tab:pos-beam-loss} \end{table}

%The dipoles in the simulation were set to transport 3~MeV positrons to the annihilation target T2, where 58\% of the positrons %annihilated and created 511~keV photon pairs. 0.7\% of the photon pairs are lost when passing through the vacuum windows. The left %NaI detector observed 2897 511~keV photons and the right one detected 4100. The ratio of 511~keV photon pairs created in T2 to %those that reached the right and left NaI are 88:1 and 62:1 respectively. This is comparable to the ratio made by out-going %photons to the solid angle made by a NaI detector.

%The distance between the NaI detectors and the beamline center is 170~cm. The Pb shielding has a 2-inch-diameter hole facing T2. %Assuming that positrons annihilated at the center of T2, the solid angle made by a NaI detector is $\Omega=\frac{\pi %r^2}{d^2}=\frac{\pi 25.4^2}{170^2}$~steradian. The 511~keV photons created during the annihilation are opposite in direction and %emitted from the two surfaces of T2 (ignoring the ones escaping from the edge). Photons emitted in each side make a solid angle of %a hemisphere, $2\pi$~steradian. The ratio of $2\pi$~steradian to solid angle of the NaI detector is about 90:1.

%1518 511~keV photon pairs are detected by NaI detectors in coincidence mode. The ratio of 511~keV photon pairs created to the ones %detected in coincidence mode is 168:1. If one counts photons in coincidence mode with the 46.24\% detection efficiency of the %system, it cuts the rate in half.

\begin{figure}[htbp] \centering \includegraphics[scale=0.74]{3-Simulation/Figures/Transporation_Efficiency/En2.eps} \caption{Energy distribution of positrons transported when dipoles were set to bend 3~MeV positrons along the beamline.} \label{fig:dipole-trans} \end{figure}


\section{Analysis of the Photon Count Asymmetry in NaI Detectors}

The beam right NaI detector observed higher rates that the beam left detector in both the experiment and the simulation. Eight virtual NaI detectors were placed as shown in Figure~\ref{fig:NaI-detectors} to study this count asymmetry. The orientation of T2 with respect to the incident positron beam was simulated for three cases. In the first case, the photon count rate was predicted by the simulation when the T2 is perpendicular to the beam (its area vector is parallel to the beam). In the second case, T2 was rotated clockwise (starting with the T2 location in the first case) about the x-axis (the axis pointing beam left) by 45$^\circ$ making the upstream side of T2 face downward as shown in Figure~\ref{fig:NaI-detectors} (a). In the third case, T2 was rotated clockwise about the y-axis (axis points beam up) by 45$^\circ$ so it is facing the beam right NaI detector as shown in Figure~\ref{fig:NaI-detectors} (b). The setup in the third case is similar to the setup in the experiment and the simulation described in the beginning of this chapter.

\begin{figure}[htbp] \begin{tabular}{ccc} \centerline{\scalebox{0.45} [0.45]{\includegraphics{3-Simulation/Figures/Photon_Counts_Asymmetry/T2FaceDown.png}}}\\ (a) \\ \centerline{\scalebox{0.45} [0.45]{\includegraphics{3-Simulation/Figures/Photon_Counts_Asymmetry/T2Exp.png}}} \\ (b)\\ \end{tabular} \caption{NaI detector locations around T2. The positron beam (blue line) is traveling along the z-axis (into the paper in the right figures). (a) T2 was rotated counter-clockwise about the x-axis by 45$^\circ$ positioning the upstream side of the T2 such that it faces the bottom NaI detector. (b) T2 was positioned as in the experiment. It was first rotated to the position as in (a), then it was rotated clockwise about the y-axis by 45$^\circ$, positioning the upstream side of the T2 such that it faces the beam right NaI detector.} \label{fig:NaI-detectors} \end{figure}

\begin{table} \centering \caption{Number of 511~keV photons observed by the NaI detectors.} \begin{tabular}{lcccccc} \toprule {T2 Placement} & Exp. & Perp. & Face & Face & Face & Face \\ & & to Beam& Down & Down & Down & Down \\ \midrule {Energy} & 3~MeV & 3~MeV & 3~MeV & 1~MeV & 6~MeV & 10~MeV \\ \midrule {NaI Right} & 18085 & 7610 & 7160 & 10315 & 6209 & 4436 \\ {NaI Left} & 12798 & 7651 & 7114 & 10254 & 6111 & 4487 \\ {NaI Top Right} & 7050 & 7580 & 12964 & 12371 & 14698 & 10636 \\ {NaI Bottom Left} & 7084 & 7609 & 18563 & 20238 & 15989 & 10239 \\ {NaI Top} & 12687 & 7599 & 14810 & 14332 & 16131 & 11479 \\ {NaI Bottom} & 18008 & 7609 & 18874 & 20193 & 16950 & 11181 \\ {NaI Top Left} & 14632 & 7656 & 12818 & 12268 & 14812 & 10735 \\ {NaI Bottom Right}& 18764 & 7623 & 18415 & 20197 & 16004 & 10317 \\ \bottomrule \end{tabular} \label{tab:photon-counts-asym} \end{table}

The 511~keV photon counts observed by eight virtual NaI detectors for one million positrons are given in Table~\ref{tab:photon-counts-asym} for each configuration. When T2 was placed perpendicular to the incoming positron beam, the 511~keV photons created inside T2 have the same probability to escape from T2 and reach any one of the detectors. As shown in the third column of Table~\ref{tab:photon-counts-asym}, all eight detectors observed a similar number of photons. In a separate simulation, T2 was positioned as in the first case and impinged by 1, 3, 6, 10~MeV positrons. The average distance traveled by positrons inside T2 before annihilation was $0.0847 \pm 0.0001$, $0.3006 \pm 0.0004$, $0.511 \pm 0.001$, $0.564 \pm 0.008$~mm for the four energies respectively. In another GEANT4 simulation, T2 was positioned according to the second case and two virtual NaI detectors were placed on both the top and the bottom of T2. The bottom detector observed more photons than the top one when 511~keV photons were generated 0.3006~mm inside T2 isotropically. Photons are more likely to reach the bottom detector because they would travel through a thinner layer of tungsten to reach it.

In the second case, shown in Figure~\ref{fig:NaI-detectors} (a), the left and the right NaI detectors had the lowest counts for 3~MeV positrons as shown in the fourth column of Table~\ref{tab:photon-counts-asym}. There is less detection probability for a photon traversing T2 in the radial direction towards the left/right detectors than the top/bottom surfaces of T2 due to the amount of material. The average distance traveled by 3~MeV positrons inside T2 before annihilation was $0.3006 \pm 0.0004$~mm. In this case, positrons annihilated near the upstream face of T2. For this reason, the bottom, bottom right, and bottom left NaI detectors (facing the upstream side of T2) observed more photons than the top, top right, and top left as shown in the fourth column of Table~\ref{tab:photon-counts-asym}.

The lower the positron beam energy, the shorter the annihilation depth, and the bigger asymmetry in the counts. As shown in the fourth and fifth columns of Table~\ref{tab:photon-counts-asym}, more/less 511~keV photons were observed on the bottom/top detectors with the 1~MeV positron beam than with the 3~MeV. As the positron beam energy increases, as shown in the sixth and seventh columns of Table~\ref{tab:photon-counts-asym}, the top and bottom detectors observed a similar number of photons, because positrons annihilate more uniformly inside T2 and the asymmetry in the counts decreases. With the increasing positron beam energy, fewer positrons were annihilated inside T2 and more penetrated through which was shown in the top two rows of Table~\ref{tab:photon-counts-asym}.

For the third case, shown in Figure~\ref{fig:NaI-detectors} (b), the top right and left bottom detectors observed the lowest counts, because a photon would need to travel in the radial direction to reach these two detectors as shown in the second column of Table~\ref{tab:photon-counts-asym}. The right, bottom and bottom right (facing upstream face of T2) observed more photons than the left, top, top left.

According to the simulation, the asymmetry in the photon counts was due to the average positron annihilation depth and the photon attenuation inside T2. Low energy positrons tend to annihilate and produce photon pairs near the incident surface. The created photons are more likely to be detected from the incident surface. In the experiment and the simulation, the right NaI detector was facing the upstream side of T2 and observed more 511~keV photons than the left.

\section{Quadrupole Triplet Collection Efficiency \\Study}

A G4beamline simulation was carried out to study the collection efficiency of the second quadrupole triplet (Q4, Q5, and Q6). %There is no magnet scanning option in G4beamline as in the ELEGANT code. The second quadrapole triplet magnets were set to the similar setting as in the experiment. In this simulation, 5,475,869,400 positrons were generated at DDNT1 and transported to DD1 to study the quadrupole triplet positron collection and transportation efficiency. Six quadrupole current settings of the triplet system were simulated as shown in Table~\ref{tab:triplet-eff}. For different quadrupole current settings, no significant differences were observed in the number of positrons, transverse beam profiles, and momentum distributions. The ratio of positrons generated at DDNT1 to the ones that enter D1 is 1525:1.

\begin{sidewaystable} \centering \caption{Quadrupole Triplet System Collection and Transportation Efficiency Data.} \begin{tabular}{cccccccccccccc} \toprule Q4 & Q5 & Q6 & Entries & $x$ & $\sigma_{x} $ & y & $\sigma_{y}$ & $P_{x}$ & $\sigma_{P_{x}}$ & $P_{y}$ & $\sigma_{P_{y}}$ & $P_{z}$ & $\sigma_{P_{z}}$ \\ \midrule A & A & A & & mm & mm & mm & mm & MeV & MeV & MeV & MeV & MeV & MeV \\ \midrule -1 & 2 & -1 & 3587220 & -0.005 & 12 & 0.030 & 12 & $~~3.0 \times 10^{-5}$ & 0.0461 & -0.00216 & 0.04553 & 3.848 & 1.875 \\ -2 & 4 & -2 & 3591423 & -0.012 & 12 & 0.049 & 12 & $~~2.2 \times 10^{-5}$ & 0.0461 & -0.00211 & 0.04554 & 3.848 & 1.875 \\ 1 & -2 & 1 & 3591509 & -0.009 & 12 & 0.040 & 12 & $-1.5 \times 10^{-5}$ & 0.0462 & -0.00216 & 0.04557 & 3.849 & 1.876 \\ 1 & 2 & 1 & 3589854 & -0.005 & 12 & 0.034 & 12 & $-1.8 \times 10^{-5}$ & 0.0462 & -0.00216 & 0.04556 & 3.849 & 1.876 \\ 2 & -4 & 2 & 3592977 & -0.007 & 12 & 0.032 & 12 & $-3.3 \times 10^{-6}$ & 0.0462 & -0.00217 & 0.04549 & 3.849 & 1.876 \\ 2 & 4 & 2 & 3589495 & -0.004 & 12 & 0.033 & 12 & $~~3.0 \times 10^{-6}$ & 0.0462 & -0.00218 & 0.04554 & 3.849 & 1.875 \\ % & & & & & & & & & & & & & \\ \bottomrule \end{tabular} \label{tab:triplet-eff} \end{sidewaystable}


\section{Systematic Errors Study using Simulation} Different sources of systematic errors in the positron production experiment was studied using G4beamline simulation. The power supply of the magnets might fluctuate $\pm$0.1~A which will change the magnetic field strength of the magnets. The misalignment of the magnets would also contribute to the systematic error of the experiment.

\subsection{Systematic Error Created by The Uncertainty in The Magnetic Fields Strength of The Magnets} The field strength of the magnets are dependent on the current provided by the power supply. The uncertainty of the magnet power supply is 0.1~A. Systematic errors for the positron counts were estimated by carrying out simulations with different magnetic field settings as shown in Table~\ref{tab:sim-error}.

The magnet settings are indicated in the top three rows of Table~\ref{tab:sim-error}. The transported positron energy is given by the first column. The 511~keV photon pairs counted in coincidence mode ($i.e.$ the original count multiplied by 46.42\%) for different magnet settings are given in the corresponding columns below. The ``Max/``Min in the table refers to the Maximum/Minimum magnetic field strength when the magnet coil current is at $I_{\text{max}}$/$I_{\text{min}}$, where $I_{\text{max}}$ = $I_{\text{def}} + 0.1$~A and $I_{\text{min}}$ = $I_{\text{def}} - 0.1$~A. ``def refers to the default magnetic field strength of the magnet.

The average counts and fractional errors are shown in the last two columns of Table~\ref{tab:sim-error}. The fractional error in the 511~keV photon pair counts is calculated by dividing the standard deviation of counts by the counts in default magnet setting, $i.e.$ $\frac{\text{standard deviation of counts}}{\text{counts in default setting}}$. %The fractional error in the 511~keV photon pair counts is calculated by dividing the standard deviation of counts at different settings by the counts in default magnet setting, $i.e.$ $\frac{\text{standard deviation of counts}}{\text{counts in default setting}}$.

\begin{sidewaystable} \centering \caption{Systematic Error Study: Counts of 511~keV Photon Pairs for Different Magnet Settings.} %\begin{tabular}{cccccccccccc} \begin{tabular}{lllllllllllll} \toprule D1 & Max & Min & Max & Min & Max & Max & Min & Min & Def & & \\ Q7 & Def & Def & Def & Def & Max & Min & Min & Def & Def & & \\ D2 & Max & Min & Def & Def & Def & Def & Def & Def & Def & & \\ \midrule Energy (MeV)& & & & & & & & & & Average & Fractional Error \\ \midrule 0.765~1.265 & 12 & 2 & 3 & 1 & 1 & 3 & 3 & 0 & 6 & 4 & 57.3 \% \\ 1.25~-~1.75 & 95 & 75 & 70 & 86 & 85 & 92 & 129 & 86 & 97 & 90 & 17.6 \% \\ 1.85~-~2.35 & 674 & 724 & 737 & 769 & 686 & 747 & 803 & 681 & 734 & 728 & 5.8 \% \\ 2.25~-~2.75 & 781 & 755 & 753 & 759 & 772 & 963 & 794 & 763 & 794 & 793 & 8.3 \% \\ 2.75~-~3.25 & 739 & 737 & 765 & 752 & 738 & 698 & 738 & 757 & 723 & 739 & 2.7 \% \\ 3.25~-~3.75 & 713 & 699 & 747 & 712 & 707 & 751 & 718 & 705 & 716 & 719 & 2.5 \% \\ 3.77~-~4.27 & 690 & 708 & 715 & 676 & 704 & 658 & 706 & 715 & 701 & 697 & 2.7 \% \\ 4.25~-~4.75 & 627 & 666 & 634 & 646 & 635 & 675 & 637 & 653 & 646 & 646 & 2.5 \% \\ 4.75~-~5.25 & 558 & 595 & 594 & 603 & 606 & 688 & 582 & 582 & 566 & 597 & 6.7 \% \\ 5.25~-~5.75 & 513 & 536 & 535 & 522 & 525 & 535 & 518 & 544 & 533 & 529 & 1.9 \% \\

\bottomrule \end{tabular} \label{tab:sim-error} \end{sidewaystable}

\subsection{Systematic Error Introduced by The Mis-alignment of The Magnets and The Uncertainty in The Electron Beam} The misalignment of the magnets are one of the main sources of systematic error. The 0 and 90 degree beamlines were aligned with a laser beam, while the 45 degree beamline was placed without any reference laser. Therefore, a misalignment was most likely to occur at the 45 degree beamline. The misalignment was estimated to be 5~mm. To study the effect of misalignments from different sources, beam line elements were shifted as shown in Table~\ref{tab:sys-sim-single}. The measured positron to electron ratios were given as well. Among the different sources, the systematic error produced when the first dipole, D1, is misaligned appeard to have the largest impact.

\begin{longtable}{|c|c|c|}

\caption{Positron to Electron Ratio When a Single Beamline Element Is Mis-Aligned and the Percent Error Compared to the Ratio When No Mis-Alignment.} \label{tab:sys-sim-single}\\ \hline

{Energy (MeV)} & {45 degree line moved } & {percent error} \\ {} & {to beam right by 3~mm} & { } \\

\hline \endfirsthead \multicolumn{3}{c}% {\tablename\ \thetable\ -- \textit{Continued from previous page}} \\ %\hline %{Energy (MeV)} & { } & {percent error} \\ %{} & {} & { } \\ %\hline \endhead \hline \multicolumn{3}{r}{\textit{Continued on next page}} \\ \endfoot \hline \endlastfoot $1.02 \pm 0.03$ & $( 8.3 \pm 3.4) \times10^{-17}$ & -6.7 \% \\ $2.15 \pm 0.06$ & $( 7.4 \pm 0.3) \times10^{-15}$ & -26.7 \% \\ $3.00 \pm 0.07$ & $( 8.1 \pm 0.3) \times10^{-15}$ & -18.7 \% \\ $4.02 \pm 0.07$ & $( 8.2 \pm 0.3) \times10^{-15}$ & -15.2 \% \\ $5.00 \pm 0.06$ & $( 7.2 \pm 0.3) \times10^{-15}$ & -7.7 \% \\

\midrule {Energy (MeV)} & {D1 raised up by 5~mm } & {percent error} \\ \midrule $1.02 \pm 0.03$ & $( 3.2 \pm 5.7) \times10^{-17}$ & -64.0 \% \\ $2.15 \pm 0.06$ & $( 6.1 \pm 0.9) \times10^{-15}$ & -39.6 \% \\ $3.00 \pm 0.07$ & $( 7.2 \pm 1.0) \times10^{-15}$ & -27.7 \% \\ $4.02 \pm 0.07$ & $( 5.5 \pm 0.9) \times10^{-15}$ & -43.1 \% \\ $5.00 \pm 0.06$ & $( 4.8 \pm 0.8) \times10^{-15}$ & -38.5 \% \\


\midrule {Energy (MeV)} & {D2 lowered down by 5~mm } & {percent error} \\ \midrule $1.02 \pm 0.03$ & $( 1.6 \pm 2.3) \times10^{-17}$ & -82.0 \% \\ $2.15 \pm 0.06$ & $( 6.9 \pm 1.0) \times10^{-15}$ & -31.7 \% \\ $3.00 \pm 0.07$ & $( 8.6 \pm 1.1) \times10^{-15}$ & -13.7 \% \\ $4.02 \pm 0.07$ & $( 6 \pm 0.9) \times10^{-15}$ & -38.0 \% \\ $5.00 \pm 0.06$ & $( 7.7 \pm 1.0) \times10^{-15}$ & -1.3 \% \\


\midrule {Energy (MeV)} & {detectors moved up by 3~mm, downstream } & {percent error} \\ {} & {by 3~mm, and right by 3~mm} & { } \\ \midrule $1.02 \pm 0.03$ & $( 9.6 \pm 8.1) \times10^{-17}$ & 7.9 \% \\ $2.15 \pm 0.06$ & $( 7.4 \pm 1.0) \times10^{-15}$ & -26.7 \% \\ $3.00 \pm 0.07$ & $( 8.2 \pm 1.1) \times10^{-15}$ & -17.7 \% \\ $4.02 \pm 0.07$ & $( 9.6 \pm 1.2) \times10^{-15}$ & -0.7 \% \\ $5.00 \pm 0.06$ & $( 7.7 \pm 1.0) \times10^{-15}$ & -1.3 \% \\


\midrule {Energy (MeV)} & {T2 moved up by 3~mm, downstream } & {percent error} \\ {} & {by 3~mm, and right by 3~mm} & { } \\ \midrule $1.02 \pm 0.03$ & $( 6.4 \pm 6.6) \times10^{-17}$ & -28.1 \% \\ $2.15 \pm 0.06$ & $( 7.7 \pm 1.0) \times10^{-15}$ & -23.8 \% \\ $3.00 \pm 0.07$ & $( 7.9 \pm 1.1) \times10^{-15}$ & -20.7 \% \\ $4.02 \pm 0.07$ & $( 7.9 \pm 1.0) \times10^{-15}$ & -18.3 \% \\ $5.00 \pm 0.06$ & $( 7.7 \pm 1.0) \times10^{-15}$ & -1.3 \% \\


\midrule {Energy (MeV)} & {T2 rotated c.w. about x-axis by 5$^\circ$ } & {percent error} \\ {} & {and c.c.w. about y-axis by 5$^\circ$} & { } \\ \midrule $1.02 \pm 0.03$ & $( 6.4 \pm 6.6) \times10^{-17}$ & -28.1 \% \\ $2.15 \pm 0.06$ & $( 7.2 \pm 1.0) \times10^{-15}$ & -28.7 \% \\ $3.00 \pm 0.07$ & $( 6.4 \pm 0.9) \times10^{-15}$ & -35.7 \% \\ $4.02 \pm 0.07$ & $( 7.2 \pm 1.0) \times10^{-15}$ & -25.5 \% \\ $5.00 \pm 0.06$ & $( 6.1 \pm 0.9) \times10^{-15}$ & -21.8 \% \\ \end{longtable}


There are many combinations of magnet misalignments (Q1$-$Q10, D1, and D1), the ``worst case scenario" is considered here to estimate the largest source of systematic error. A ``worst case scenario" was created by mis-aligning the beamline as given in Table~\ref{tab:mis-sim}. The results is given in Table~\ref{tab:sys-sim} in terms of the ratio of the number of e$^+$ detected in coincidence mode by NaI detectors to the number of e$^-$ incident on the T1. Beam line misalignments far beyond what can be reasonably argued indicate drops in positron rate dropped by 60\% to 80\%. A more realistic mis-alignment of beamline elements dropped this rate around 20\%. Since the chance of the beamline misalignment is large, the systematic error of the simulation is argued to be at least around 20\%.

\begin{table} \centering \caption{Mis-alignment of The Beam in The Worst Case Scenario.} \begin{tabular}{ll} \toprule {Beam Element} & {Placement} \\ \midrule D1 & raised up by 5~mm \\ D2 & lowered down by 5~mm \\ 45 degree line & shifted to beam right by 3~mm \\ NaI detectors & raised up by 3~mm, moved 3~mm beam downstream and 3~mm \\

 &  to beam right     \\ 

T2& raised up by 3~mm, moved 3~mm beam downstream and 3~mm \\

 &  to beam right, rotated clockwise 5$^\circ$ about x-axis and counter \\
 &  clockwise 5$^\circ$ about y-axis    \\ 

\bottomrule \end{tabular} \label{tab:mis-sim} \end{table}

\begin{table} \centering \caption{Positron to Electron Ratio Predicted When Beamline is Mis-aligned as Given in Worst Case Scenario and the Percent Error Compared to the Ratio When No Misalignment..} \begin{tabular}{lll} \toprule {Energy (MeV)} & {Positron to Electron Ratio} & percent error \\ \midrule $1.02 \pm 0.03$ & $( 3.4 \pm 3.5) \times10^{-18}$ & -61.8 \% \\ $2.15 \pm 0.06$ & $( 1.79 \pm 0.50) \times10^{-15}$ & -83.2 \% \\ $3.00 \pm 0.07$ & $( 3.39 \pm 0.68) \times10^{-15}$& -66.9 \% \\ $4.02 \pm 0.07$ & $( 3.14 \pm 0.66) \times10^{-15}$& -67.9 \% \\ $5.00 \pm 0.06$ & $( 2.82 \pm 0.62) \times10^{-15}$& -64.1 \% \\ \bottomrule \end{tabular} \label{tab:sys-sim} \end{table}

Chapter 5: Conclusion File:Sadiq thesis chapt 5.txt

\chapter{Conclusions and Suggestions} A new High Repetition Rate Linac (HRRL) beamline located in the ISU's Physics Department beam lab has been successfully reconstructed to produce and transport positrons to the experimental cell. The electron beam energy profile and emittance of the HRRL were measured using a Faraday cup and an OTR based diagnostic system. The positron production rate was measured for positron energies between 1 and 5 MeV. The results are shown in Figure~\ref{fig:e+2e-exp-sim} along with the prediction made by a GEANT4 simulation of the beamline.

The production of positrons using an electron linac was done in several steps. First, positrons are emitted from the downstream side of a tungsten target (T1) when electrons impinge on the upstream side an produce photons of sufficient energy to pair produce within the tungsten target. Postrons escaping the downstream side of the target were collected by the quadrupole triplet. The positrons are then deflected by two dipoles in order to measure the positron rate as a function of the positron energy. Positrons that traversed the two dipoles would annihilate in a second tungsten target (T2) producing back-to-back 511 keV photons that were measured using two NaI detectors. The positron rate was mesured by requiring a coincidence between both NaI detectors and the electron beam pulse.

The positron beam creation, beam loss in the transportation, and detection process were studied using the simulation package G4beamline and compared to this experiment. The simulated e$^+$/e$^-$ ratios are shown in Figure~\ref{fig:e+2e-exp-sim} for the energies measured in this experiment. The simulation includes electron beam generation with the measured electron energy profile, beam losses during transportation, positron annihilation in the tungsten target (T2), and the detection of 511~keV photons in coincidence by the two NaI detectors. While the simulation result agrees with the experiment in that the peak energy distribution is near 3~MeV, it predicts a higher positron to electron ratio as shown in Figure~\ref{fig:e+2e-exp-sim}.

The simulation was used to study the systematic errors in the experiment. The simulation predicts that a realistic misalignment of the beamline can reduce the e$^+$/e$^-$ ratios by 20\% to 30\%. In the worst case scenario, the ratios dropped by 60\% to 80\%. The systematic errors in the experiment bring it into agreement with the simulation.


The ratio of the positrons contained within the 90 degree beampipe to the 511~keV photons detected in coincidence mode when the dipoles were set to bend 3~MeV positrons is predicted to be 1655:1 by the simulation. The ratio of the positrons on T2 to the 511~keV photons is 620:1 under the same conditions as above. The 3~MeV positron rate measured in the experiment was $0.25\pm0.2$~Hz when the HRRL was operated at a 300~Hz repetition rate, 100~mA peak current, and 300~ns (FWHM) RF macro pulse length. Based on this simulation, a measured $0.25 \pm 0.02$~Hz coincidence rate by the NaI detectors would correspond to a $155 \pm 12$~Hz positron rate incident on T2.

In the simulation, the number of positrons collected was insensitive to the quadrupole triplet collection field setting (see section 4.5). The ratio of solid angles subtended by the quadrupole (Q4) and dipole (D1) entrance windows approximated the ratios of positron transported. Dipoles defocus in one plane and defocus in the other. Thus, one can only collect positrons in one plane while loses occur in the other. Solenoids, on other hand, focus in both planes. A solenoid may be a better option to improve the collection efficiency. Positioning the target T1 at the entrance of the solenoid may be the optimal choice for capturing positrons.

%7. Experimental results show quadrupole magnets are not efficient in collecting positrons, since positrons have large angular distribution. Solenoid might be able to improve the collection efficiency of positrons~\cite{kim-bindu-solenoid} and should be placed as close the production target as possible for better efficiency.

\begin{figure}[htbp] \includegraphics[scale=0.79]{5-Conclusion/Figures/Overlay_Exp-Sim-Ratio/R.eps} \caption{Ratio of positrons detected to electrons measured in the experiment (hollow diamond) and simulation (full circle) in coincidence mode. The black solid error bars are statistical and dashed ones are systematic. The experimental systematic errors (red dashed lines) are discussed in section 3.4 of Chapter 3 and the systematic errors in simulation (red dashed lines) estimation is described section 4.6 of Chapter 4.} \label{fig:e+2e-exp-sim} \end{figure} %An OTR based diagnostic tool was designed, constructed, and used to measure the beam emittance of the HRRL. The electron spatial profile measured using the OTR system was not described by a Gaussian distribution but by a super Gaussian or Lorentzian distribution. The unnormalized projected emittances of the HRRL were measured to be less than 0.4~$\mu$m by the OTR based tool using the quadrupole scanning method when accelerating electrons to an energy of 15~MeV. %OTR used, Not Guassian, Changed magnet, measured emttiance

pdf file: File:Sadiq hesis Latex.pdf