From: Nathan Rodning <rodning@relay.phys.ualberta.ca>
Date: Sun, 28 Feb 1999 14:15:58 -0700 (MST)
To: e614tn@relay.phys.ualberta.ca
Subject: TN-28 Injection of Muon Beam into Solenoid

The optimal focus and phase space of the beam as it crosses
the solenoid fringe field is studied, and consequences for
the effective polarization and count rate are shown.


*************************************************
Description: full .ps file , Filename: sol981.ps

Description: figure 8 , Filename: fig8.ps

Description: figure 7 , Filename: fig7.ps

Description: figure 6 , Filename: fig6.ps

Description: figure 5 , Filename: fig5.ps

Description: figure 4 , Filename: fig4.ps

Description: figure 3 , Filename: fig3.ps

Description: figure 2 , Filename: fig2.ps

Description: figure 1 , Filename: fig1.ps

\documentstyle[11pt,rotate]{article}
\renewcommand{\topfraction}{1.0}
\renewcommand{\bottomfraction}{1.0}
\renewcommand{\floatpagefraction}{1.0}

\input epsf
\setlength{\parindent}{0.5cm}
\begin{document}

\title{Injection of a 30 MeV/c surface muon beam into the 2.3 Tesla solenoid 
for TRIUMF experiment E614}
\author{Jaap Doornbos, TRIUMF}
\date{February 27, 1998}
\maketitle


\section{Introduction and summary}


TRIUMF experiment E614 aims to measure the Michel parameters of muon
decay. A surface muon beam is stopped in the center of a 2.3 T solenoid.
 Surface muons have the spin in the direction of
motion because they are produced by the decay of pions at rest in the
surface of the production target. Since G=2 for muons, the spin follows the
direction of motion during transport through magnetic elements.

An essential requirement of the experiment is that the beam is 100 percent
polarized in the longitudinal Z direction. It is inevitable that the muons
have a transverse momentum. Their direction makes an angle $\alpha$ with
the axis of the solenoid and the component of the spin in the z direction
is smaller than unity. If the $\alpha$ is 10 mr, than the depolarization is
$1-cos(0.010)=0.5*10^{-4}$. If $\alpha$ is 15 mr, the depolarization is
$1.1*10^{-4}$,while it is $2.0*10^{-4}$ for $\alpha$ is 20 mr.
The experiment requires a beam which has a small depolarization. $2.0
*10^{-4}$ is too big.

A beam that is injected into a solenoid passes through a fringe field
region, where the beam is being focussed by the radial component of the
magnetic field. The angle in the beam becomes large, giving a large
depolarzation. The injection of the muon beam from TRIUMF beam line M13
needs to be optimized to keep the increase in angle as small as
possible. This note gives some very preliminary results from a study to
understand the behaviour of the muon beam in the solenoid.  

\section{Method of calculation}

Figure 1a shows the Z component of the magnetic field along the axis of the
solenoid as obtained from calculations with the program OPERA2D. Z=0
corresponds to the center of the solenoid. This has been obtained from the
full two dimensional map which extends to a radius of 30 cm.
 The field on the axis is
sufficient to determine the field anywhere in the solenoid using the
Maxwell equations. This is done in the raytracing program ZGOUBI, which
traces particles stepwise through electric and magnetic fields in a
calculation which is accurate to high order. Notice that the field rises
steeply from Z=-160 cm on. Here is where the yoke of the solenoid ends. The
radial component $B_{r}$ of the field off the axis is to good approximation
proportional to $(r/2)dB/dz$, where r is the distance from the
axis. This was checked for the two dimensional field map for the region
where the beam is. Figure 1b
shows that $B_{r}$ has a double maximum.






\section{The behaviour of some individual rays in the solenoid}

Figures 2a and 2c are for for four rays which start out parallel to the axis at
1, 3, 5 and 7 mm from the axis.  Figure 2a shows how the radial position
changes. The fringe field focusses the beam to the axis at a position which
is independent of the inital radius.The radial
focussing field $B_{r}$ is proportional to the distance from the axis.
Figure 2c gives the development of the angle $\alpha$ as the particle enters
the solenoid.
Figures 2b and 2d are for rays converging into the solenoid with initial
angles of 5, 10, 15 and 20 mr. When the solenoid field is turned off they
converge to a point 152 cm from the center. This is the optimum situation,
as will be shown in the next section.

Figure 3 shows X-Y projections of how the indicated rays develop transversely. 

Figure 4 shows the dependence of $\alpha$ in the center of the solenoid as
function of the initial radius and of the initial angle. Even for a radius
of 3 mm the final angle is already 18 mr. An initial angle of 10 mr gives a
final angle of 17.5 mr.


 

\section{Monte Carlo calculations for a collection of rays}

The angle $\alpha$ in the center of the solenoid depends on the way the beam
enters the solenoid. It is important to have the beam come to a waist
somewhere in the region where the focussing by the fringe field is
strongest. In order to study this we take a beam starting from a waist with
a radius of 5 mm and a maximum angle with the axis of 20 mr. The projected emittance
in both the vertical and the horizontal plane is then $\pi *5*20=100 \pi$
mm.mr. This is the phase space area of the beam projected on either the
vertical or the horizontal plane.
 This beam is traced back and then injected
into the solenoid.

\subsection{Where to focus the beam?}
 In figure 5 $S$ is a measure of where the beam would come
to a waist if the solenoid is turned off. If $S=0$, the beam comes to a
waist 100 cm from the center of the solenoid. If $S=100$ the beam comes to
a waist at 200 cm from the center. The total number of traced particles is
617. Figures 5a and 5b show the number of particles which arrive in the
center with angles less than 10 or 20 mr, while figure 5c shows the maximum
angle in the beam. The optimum situation occurs is when the beam is aimed
at a point 152 cm from the center. As figure 1b indicates, this is just
behind the first bump in the radial field. It turns out that the optimum situation does not
depend very much on the choice of the emittance.

\subsection{Emittance dependence}
How does the intensity within a chosen angle range in the center depend on
the projected emittance?

This is shown in figure 6. To vary the emittance we change the radius from
3 to 10 mm in steps of 1 mm and, simultaneously, the angle from 12 to 40
mr in steps of 4 mr. The projected emittance varies from $36 \pi$ mm.mr to
 $400 \pi$ mm.mr. Figures 6a and 6b show the intensity within 10 and 20 mr in
arbitrary units. For easy comparison, the intensity is normalized to 100
for an emittance of $100 \pi$ mm.mr. The intensity is almost a linear
function of the projected emittance. The total phase space and intensity of
the incoming beam increases, in fact, as the square of the projected emittance.
Figures 6c and 6d show the ratio of the total number of particles and the
particles within 10 or 20 mr. For example for an emittance of $100 \pi$
mm.mr the number of particles within 10 mr is about 8 percent.

\subsection{Influence of the shape of the emittance}
Does the number of particles within the chosen angle range depend on the
shape of the projected emittance?

 Results are given in figure 7.
The emittance is $100 \pi$ mm.mr which can be obtained with a radius of 3 mm
and an angle of 33 mr, a radius of 5 mm and an angle of 20 mr, a radius of
10 mm and an angle of 10 mr, etcetera. In figure 7 the variation of the
shape of the emittance is represented by the radius R. Figures 7a and 7b
show the percentage of the total intensity within 10 or 20 mr, while 7c
shows the maximum angle in the center of the solenoid. The conclusion is
that 5 mm radius and 20 mr angle is an optimum choice.
\subsection{Angle distribution for the optimum situation}
Figure 8 shows for the choice of 5 mm and 20 mr the percentage of particles
within the indicated range of radius or angle in the center. The bottom two
figures give the accumulated percentages. 




\section{Consequences for the M13 beam}

\subsection{Intensity of M13 beam in 1 percent $\Delta$P/P momentum bite}

M13 takes of from proton beamline BL1A at an angle of 135 degrees. The
production target is 1.0 cm long graphite. The proton spot on the target
is of the order of 2 mm in diameter. 
This presents a source to the M13 channel which is horizontally $(10.0 +2)sin(45)=8.4$
mm. Vertically we take a size of 5
mm. This is bigger than 2 mm, because what counts is not the size of the
proton beam but the effective size of the surface area of the target which
contributes to the surface muon beam. M13 accepts vertical angles of up to
160 mr and horizontal angles up to 60 mr. The horizontal emittance is 
$60*8.4/2*\pi=250\pi$ mm.mr. The vertical emittance is
$160*5.0/2*\pi=400\pi$ mm.mr. The momentum acceptance is 8.0
percent$\Delta$P/P. The measured total intensity of surface muons in this phase
space is 1,600 K per second per 150 $\mu$A proton beam. The experiment needs
a small momentum bite, not larger than 1.0 percent. This gives an intensity
of 200 K per second.

\subsection{Intensity of M13 in $100 \pi$ mm.mr emittance
and 1 percent $\Delta$P/P momentum bite}

 It is not advisable to use the full emittance of M13,
because it would give a large ratio of unwanted to wanted muons, the latter
being muons within a 10 mr angle $\alpha$ in the center of the solenoid. For the
sake of argument, assume that an emittance of $100\pi$ mm.mr in both planes
is a good choice. In that case there will be 12 times as many muons
with $\alpha$ larger than 10 mr, than muons with $\alpha$ smaller than 10 mr. There are
$(100/250)*(100/400)*200 K=20 K$ muons per second in the beam. Of these
muons, there are 8 percent or 1.6 K within a 10 mr angle $\alpha$.
It is important to emphasize that this estimate depends crucially on the
assumption that the vertical source size at the production target is 5
mm. If, in fact, the source size is larger, then the intensity within $100
\pi$ mm.mr is less than 1.6 K per second. It is necessary to find out the
source size.

\subsection{Optimum distance between M13 and the E614 solenoid}
Some calculations have been done to determine whether it is possible to get
an optimum emittance shape for M13. It turns out that a distance of about
1.1 m between the last quadrupole of M13 and the coils of the solenoid ia
optimal. In that case the size of the vertical spot would be 9 mm. The half
size, which is used for calculating the emittance, is 4.5 mm, close to the
optimum 5 mm. Horizontally, the size is 13 mm. The half size is 6.5 mm, not
so far from the optimal 5 mm.
If the distance is increased, the horizontal spot will increase, which is
undesirable. It is not possible to decrease the horizontal spot size,
because of the size of the production target. There is no freedom to vary
the horizontally focussing quadrupoles in the beam line because they are
determined by the requirement of making a dispersed focus at the momentum
slit, and making the beam achromatic after the second bending magnet.



 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig1.ps}}
     \caption{Longitudinal and radial field of the solenoid.}
     \label{fig:1}
   \end{center}
 \end{figure}

 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig2.ps}}
      \caption{The radius and the angle $\alpha$ along the solenoid for selected
               rays. Figures 2a and 2c are for four rays parallel to the
               axis, 1, 3, 5 
               and 7 mm
               from the axis. Figures 2b and 2d are for four rays aimed at a point 152
               cm from the center of the solenoid, starting with angles of 5, 10, 15
               and 20 mr.}
     \label{fig:2}
   \end{center}
 \end{figure}



 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig3.ps}}
      \caption{The transverse X-Y motion of a subset of the rays of figure 2.}
     \label{fig:3}
   \end{center}
 \end{figure}


 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig4.ps}}
       \caption{The dependence of $\alpha$ in the center of the solenoid on the
                initial angle in figure 4b, and on the initial distance to the axis in
                figure 4a}
     \label{fig:4}
   \end{center}
 \end{figure}


 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig5.ps}}
       \caption{Dependence of the intensity within $\alpha$ less than 10 or 20 mr
                in figures 5a and 5b. S indicates where the beam is focussed when the
                solenoid is turned off. $S=0$ corresponds to a point 100 cm upstream of
                the center of the solenoid, $S=100$ corresponds to a point 200 cm upstream
                of the center. This calculation is for a beam with a $100 \pi$ mm.mr
                projected emittance, a radius of 5 mm, and an angle of 20 mr. Figure 5c
                gives the maximum angle in the center}
     \label{fig:5}
   \end{center}
 \end{figure}

 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig6.ps}}
       \caption{The intensity in the center of the solenoid with $\alpha$ less
                than 10 or 20 mr in figures 6a and 6b, as function of the projected
                emittance. The intensities are normalized to 100 for an emittance of
                $100 \pi$ mm.mr. In fact, there are roughly four times as many muons
                in the case of 20 mr, than in the case of 10 mr. Figures 6c and 6d show
                the ratio of the total number of muons to the muons in the indicated
                angle range.}
     \label{fig:6}
   \end{center}
 \end{figure}

 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig7.ps}}
       \caption{The intensity for different shapes of a $100 \pi$ mm.mr
                emittance in figures 7a and 7b. The angle varies inversely
                to the radius to keep the
               emittance constant. If R is 5 mm, the angle is 20 mr. This is the
               optimum choice. The maximum angle in the center of the solenoid is shown
               in figure 7c.}
     \label{fig:7}
   \end{center}
 \end{figure}


 \begin{figure}[htb]
   \begin{center}
     \leavevmode
     \rotate[l]{\epsfysize=5.1in \epsfbox{fig8.ps}}
       \caption{ Figures 8a and 8b give for a $100 \pi$ mm.mr emittance,
                 radius of 5 mm and angle of 20
                 mr, the percentage of beam as function of R or $\alpha$ in the center of
                 the solenoid. Figures 8c and 8d give the accumulated percentages.}
     \label{fig:8}
   \end{center}
 \end{figure}


\end{document}

TN-28 Injection of Muon Beam into Solenoid / Nathan Rodning

Created for the The Center for Subatomic Research E614 Project Projects Page.
Created by The CoCoBoard.