From: Robert MacDonald <rmacdon@relay.phys.ualberta.ca>
Date: Thu, 8 Mar 2001 13:55:27 -0700 (MST)
To: E614TN@relay.phys.ualberta.ca
Subject: TN-52: Summary of M13 Momentum Scans of Surface Muon Edge

  This document describes four momentum scans performed on the M13
  channel---two with a beryllium 1AT1 target, and two with graphite. The
  scans were performed by Glen Marshall and Robert MacDonald.

  The muon and positron rates were measured for a series of magnetic
  field settings of the first dipole magnet (B1) of the M13 channel.
  The B1 magnet selects the channel momentum, and the purpose of these
  scans was to calibrate the B1 magnetic field (according to an NMR
  probe) to the well-known surface-muon edge momentum (29.79 MeV/c).
  Since this edge should be very sharp, the measured width of the edge
  also gives the momentum bite of the channel for the chosen F1 slit
  setting.

Filename: momentumscan_technote.ps

\documentclass[11pt,letterpaper]{article}
\usepackage{graphicx}  % Use this for inserting graphics into figures.
%  Graphics can be inserted with the \includegraphics command.

% Widen the margins (measured from 1in down and 1in from left edge of paper).
\setlength{\topmargin}{-0.5in} \setlength{\headsep}{0.5in} 
\setlength{\footskip}{0.5in}
\setlength{\oddsidemargin}{0in} 
\setlength{\textheight}{9in} \setlength{\textwidth}{6.5in}
% Allow figures to take up to 80% of the page.
\renewcommand{\textfraction}{0.2}
% 1.5 line spacing.
%\renewcommand{\baselinestretch}{1.5}
% Use unindented block paragraphs.
%\setlength{\parskip}{2ex plus0.5ex minus0.5ex}
%\setlength{\parindent}{0em}
\pagestyle{headings}

% Formatting:
\newcommand{\textcal}[1]{\ensuremath{\mathcal{#1}}}
% Typing shortcuts:
\newcommand{\e}[1]{\ensuremath{\times 10^{#1}}}
\newcommand{\musr}{$\mu$SR }
\newcommand{\twist}{\textcal{TWIST} }

%%%%%%%%%%%%%%%%
\begin{document}


\title{Summary of M13 Momentum Scans of Surface Muon Edge\\
  (\twist Technote 52)}
\author{Robert MacDonald\\\twist\\University of Alberta}
\date{8 March, 2001}
\maketitle

\pagenumbering{roman}

\begin{abstract}

  This document describes four momentum scans performed on the M13
  channel---two with a beryllium 1AT1 target, and two with graphite.
  The scans were performed by Glen Marshall and Robert MacDonald.

  The muon and positron rates were measured for a series of magnetic
  field settings of the first dipole magnet (B1) of the M13 channel.
  The B1 magnet selects the channel momentum, and the purpose of these
  scans was to calibrate the B1 magnetic field (according to an NMR
  probe) to the well-known surface-muon edge momentum (29.79 MeV/c).
  Since this edge should be very sharp, the measured width of the edge
  also gives the momentum bite of the channel for the chosen F1 slit
  setting.

  Summary of the momentum scan results:

  \begin{center}
    \begin{tabular}{|l|r|r|}
      \hline
      \textbf{Beryllium target} & \textbf{30 Oct.} & \textbf{5 Dec.} \\
      \hline
      B1 Field at Edge (Gauss) & $879.6 \pm 0.4$ & $876.5 \pm 0.4$\\
      Momentum bite & 2.1\% & 1.3\%\\
      \hline
      \hline
      \textbf{Graphite target} & \textbf{13 Nov.} & \textbf{20 Nov.}\\
      \hline
      B1 Field at Edge (Gauss) & $878.7 \pm 0.4$ & $877.8 \pm 0.4$ \\
      Momentum bite & 1.6\% & 1.7\% \\
      \hline
    \end{tabular}
  \end{center}

  The momentum bite results are in agreement with the (very rough)
  estimate of 2\% momentum bite for the F1 horizontal slit setting of 25
  mm.

\end{abstract}

%\tableofcontents
%\listoffigures
%\listoftables
%\clearpage

\pagenumbering{arabic}
\setcounter{page}{1}

\section{About the Momentum Scans}
\label{sec:pscan}

The momentum scans were taken using two different targets at 1AT1.
The beryllium target was in place for the 30 October and 5 December
scans.  The graphite target was in place for the 13 and 20 November
scans.  (All scans took place in 2000.)  The scans were performed by
Glen Marshal and/or Robert MacDonald.

An initial tune, at a momentum slightly below the surface muon edge,
was chosen to be the ``reference momentum.''  The channel momentum was
scaled by percentages relative to this reference.  The B1 and B2
bender magnets were scaled by their NMR readings; the quadrupoles were
scaled by their DAC values.

About the November 13 scan: A coarse scan was performed, starting at
the reference momentum and going up to 20\% higher, in various step
sizes.  Additional points were then filled in on the surface muon
edge, as well as points down to 10\% below the reference momentum.

The same scale percentages were used for the later two scans (20 Nov.
and 5 Dec.)  described here, but scaling monotonically from -10\% to
+20\%.  It was hoped this would reduce any hysteresis in the
quadrupole magnets.

The first scan (30 Oct.) used different (slightly coarser)
percentages, but was also done monotonically in an attempt to avoid
hysteresis effects.

For all but the 30 Oct. scan, the detector at the end of the M13
channel consisted primarily of a set of scintillators and an aluminum
target, arranged as shown in figure \ref{fig:apparatus}.  $M1$ is a
thin scintillator used for counting muons, and is relatively
insensitive to positrons.  Surface muons stop in the aluminum.  For
these reasons, $(M1\cdot\bar{P})$ was used as a muon ID.  Positrons
pass through the entire setup undeterred, so a coincidence of $E
\equiv (P \cdot E1 \cdot E2)$ is a good positron ID.  

\begin{figure}
  \begin{center}
    \includegraphics{pscan_apparatus.eps}
    \caption[Schematic diagram of apparatus.]{Schematic diagram of
      apparatus.  M1, P, E1, and E2 are scintillators.  WC1 and WC2 are wire
      chambers (WC1 is pushed against the end of the beam pipe).  Al is a
      set of aluminum sheets.  (Not to scale.)}
    \label{fig:apparatus}
  \end{center}
\end{figure}

Two wire chambers were also part of the apparatus, but were used for
other beam studies.  They are included in the diagram for
completeness.

The 30 Oct. scan used a similar apparatus, but without the wire
chambers (the apparatus used was for the \musr studies).

The data was taken using visual scalers, set to accumulate over 10
second intervals.  Scalers counted $M1$, ``clean'' muons
$(M1\cdot\bar{P})$, positrons, and the T1 ion chamber (which puts out
pulses at a rate proportional to the proton current in the BL1A
beamline).  The rates were then normalized to proton current by
dividing by the T1Ion scaler.


\section{Muon Flux, Surface Muon Edge, and Momentum Bite}
\label{sec:muflux}

%\subsection{Beryllium 1AT1 Target}
%\label{sec:Be}

%\subsection{Graphite 1AT1 Target}
%\label{sec:C}

Figures \ref{fig:pscan_muons_Be} and \ref{fig:pscan_muons_C} show the
normalized muon flux as a function of the B1 NMR reading (in Gauss),
for both scans on each target.  The shape is approximately as
predicted: the muon rate rises with momentum until it reaches the
surface muon edge, at which point it falls sharply.  It resumes
rising, but more slowly, after the edge---at these momenta, the only
muons are ``cloud'' muons.
\begin{figure}
  \begin{center}
    \includegraphics[angle=90, width=3in, height=3in]{pscan_muons_001030.ps}
    \includegraphics[angle=90, width=3in, height=3in]{pscan_muons_001205.ps}
  \end{center}
  \caption[Normalized muon flux vs. B1 NMR reading, Be
  target.]{Normalized muon flux vs. B1 NMR reading, \textbf{Be
      target}.  Left: October~30.  Right: December~5.  Muon flux is
    normalized against the T1 Ion Chamber rate, which is proportional
    to the BL1A proton rate.}
  \label{fig:pscan_muons_Be}
\end{figure}
\begin{figure}
  \begin{center}
    \includegraphics[angle=90, width=3in, height=3in]{pscan_muons_001113.ps}
    \includegraphics[angle=90, width=3in, height=3in]{pscan_muons_001120.ps}
  \end{center}
  \caption[Normalized muon flux vs. B1 NMR reading, C
  target.]{Normalized muon flux vs. B1 NMR reading, \textbf{C target}.
    Left: November~13.  Right: November~20.  Muon flux is normalized
    against the T1 Ion Chamber rate, which is proportional to the BL1A
    proton rate.}
  \label{fig:pscan_muons_C}
\end{figure}

The minima and maxima of these graphs are summarized in table
\ref{tab:pscan_muons_minmax}.  Of particular interest is the graphite
(C) target data.  The minima (mostly cloud muons) are roughly the same
between the two scans, but the maxima (mostly surface muons) differ by
a factor of almost two.  The graphite target split in two during the
November beam studies, and it is not known when this happened; the
difference in surface muon flux may be due to this break.
\begin{table}[tbp]
  \begin{center}
    \begin{tabular}{|l|l|l|l|}
      \hline
      Date & Target & Min. Flux & Max. Flux \\
      \hline
      Oct. 30 & Be & 0.066 & 0.73 \\
      Dec. 5  & Be & 0.043 & 0.74 \\
      \hline
      Nov. 13 & C  & 0.14  & 6.3  \\
      Nov. 20 & C  & 0.13  & 2.7  \\
      \hline
    \end{tabular}
    \caption{Minima and maxima of normalized muon flux graphs from
      figures \ref{fig:pscan_muons_Be} and \ref{fig:pscan_muons_C}.}
    \label{tab:pscan_muons_minmax}
  \end{center}
\end{table}

The curves are not quite as smooth as expected.  Some of the jitter in
the November 13 data could be attributed to hysteresis (the data was
not taken in order of momentum), but not so for the other scans.
Moreover, the shape of these jitter ``features'' is reproduced in
three of the four graphs (and the fourth---from 30 Oct.---was taken
using different equipment).

Figures \ref{fig:pscan_muons_zoom_Be} and
\ref{fig:pscan_muons_zoom_C} show an expanded view of the surface
muon edge for each scan.  The ``10\%'' and ``90\%'' levels correspond
to 10\% and 90\% of the range between the minimum and maximum muon
fluxes.
\begin{figure}
  \begin{center}
    \includegraphics[angle=90, width=3in,
    height=3in]{surfacemu_closeup_001030data.ps}
    \includegraphics[angle=90, width=3in,
    height=3in]{surfacemu_closeup_001205data.ps}
  \end{center}
  \caption[Closeup of surface muon edge, Be target.]{Closeup of surface
    muon edge, \textbf{Be target}.  Normalized as in figure
    \ref{fig:pscan_muons_Be}.  The horizontal dashed lines show the
    10\% and 90\% levels, and the vertical dashed lines show where the
    momentum scan crosses those levels.  The boxed point shows the
    halfway-point between those levels, ie. the surface muon edge (at
    $879.6 \pm 0.4$ Gauss (left graph) and $876.5 \pm 0.4$ Gauss
    (right graph)).}
  \label{fig:pscan_muons_zoom_Be}
\end{figure}
\begin{figure}
  \begin{center}
    \includegraphics[angle=90, width=3in,
    height=3in]{surfacemu_closeup_001113data.ps}
    \includegraphics[angle=90, width=3in,
    height=3in]{surfacemu_closeup_001120data.ps}
  \end{center}
  \caption[Closeup of surface muon edge, C target.]{Closeup of surface
    muon edge, \textbf{C target}.  Normalized as in figure
    \ref{fig:pscan_muons_C}.  The horizontal dashed lines show the
    10\% and 90\% levels, and the vertical dashed lines show where the
    momentum scan crosses those levels.  The boxed point shows the
    halfway-point between those levels, ie. the surface muon edge (at
    $878.7 \pm 0.4$ Gauss (left graph) and $877.8 \pm 0.4$ Gauss
    (right graph)).}
  \label{fig:pscan_muons_zoom_C}
\end{figure}

Let $B$ represent the magnetic field in B1 (that is, the B1 NMR).
Then $B_{10}$, $B_{50}$, and $B_{90}$ are the $B$ values for the 10\%,
50\%, and 90\% levels.

The 10\% and 90\% levels were chosen because the derivative of the
edge shape should approximate a Gaussian curve.  In this case, the B1
NMR values for the 10\% and 90\% levels should roughly correspond to
the half-maximum points of the Gaussian's peak.  Therefore, the
difference $B_{90} - B_{10}$ should give the ``width'' of the edge,
and $B_{50}$ should occur at exactly the surface muon edge.

The 10\% and 90\% levels are shown as dashed horizontal lines in
figures \ref{fig:pscan_muons_zoom_Be} and \ref{fig:pscan_muons_zoom_C}; the
corresponding B1 NMR values are shown as dashed vertical lines.  The
boxed point shows the position of $B_{50}$---the surface muon edge of
29.79 MeV/c.  This point occurs at the field values shown in table
\ref{tab:pscan_results}.  (The uncertainty in the field values
represents the precision to which the B1 magnet can be set, as limited
by the coarseness of the power supply control (DAC).  A detailed error
analysis has not yet been done, so this may not be the only source of
error on these points.)  The fields for the graphite target momentum
scans agree within error.  The beryllium target results, however, vary
significantly.  
\begin{table}[tbp]
  \begin{center}
    \begin{tabular}{|l|r|r|}
      \hline
      \textbf{Beryllium target} & \textbf{30 Oct.} & \textbf{5 Dec.} \\
      \hline
      B1 Field at Edge (Gauss) & $879.6 \pm 0.4$ & $876.5 \pm 0.4$\\
      Momentum bite & 2.1\% & 1.3\%\\
      \hline
      \hline
      \textbf{Graphite target} & \textbf{13 Nov.} & \textbf{20 Nov.}\\
      \hline
      B1 Field at Edge (Gauss) & $878.7 \pm 0.4$ & $877.8 \pm 0.4$ \\
      Momentum bite & 1.6\% & 1.7\% \\
      \hline
    \end{tabular}
    \caption{Summary of momentum scan results, both targets.}
    \label{tab:pscan_results}
  \end{center}
\end{table}

One possible explanation for the difference in beryllium target
results is beam tuning: the 5~December scan was taken during beam tune
studies, and the tune used was probably not exactly the same as that
used for the other scans.  (In fact, it was probably a better tune
(closer to what \twist will likely use).)

Figures \ref{fig:deriv_Be} and \ref{fig:deriv_C} show plots of the
``derivative'' of the surface muon edge.  The derivative was
calculated as the difference in y positions (normalized flux) of
adjacent points, divided by the difference in x positions (B1 NMR).
They are plotted at the centres of their B1 NMR intervals.  Vertical
lines show the positions of $B_{10}$, $B_{50}$, and $B_{90}$.  All of
these plots are vaguely Gaussian, some moreso than others; with all
curves, with both curves $B_{10}$ and $B_{90}$ occur near the half-max
points of the curves, and $B_{50}$ occurs roughly in the middle of the
peak.  Note that uncertainties in these ``derivatives'' were not
calculated.
\begin{figure}
  \begin{center}
    \includegraphics[angle=90, width=3in, height=3in]{naive_deriv_001030.ps}
    \includegraphics[angle=90, width=3in, height=3in]{naive_deriv_001205.ps}
  \end{center}
  \caption[``Derivative'' of surface muon edge, Be
  target.]{``Derivative'' (slope of segments between data points) of
    surface muon edge, \textbf{Be target}.  Vertical lines show
    $B_{10}$, $B_{50}$, and $B_{90}$ (see text).  These plots do not
    include uncertainties.}
  \label{fig:deriv_Be}
\end{figure}
\begin{figure}
  \begin{center}
    \includegraphics[angle=90, width=3in, height=3in]{naive_deriv_001113.ps}
    \includegraphics[angle=90, width=3in, height=3in]{naive_deriv_001120.ps}
  \end{center}
  \caption[``Derivative'' of surface muon edge, both
  scans, C target.]{``Derivative'' (slope of segments between data
    points) of surface muon edge, \textbf{C target}.  Vertical lines
    show $B_{10}$, $B_{50}$, and $B_{90}$ (see text).  These plots do
    not include uncertainties.}
  \label{fig:deriv_C}
\end{figure}

The magnetic field $B$ should be directly proportional to the channel
momentum $p$.  Therefore 
\begin{equation}
  \label{eq:bite}
  \frac{dB}{B} = \frac{dp}{p} .
\end{equation}
Thus we can calculate the momentum bite $dp/p$ from the width of the
surface muon edge: 
\begin{equation}
  \label{eq:dB}
  \frac{dp}{p} \approx \frac{|B_{10} - B_{90}|}{B_{50}}
\end{equation}
The momentum bite results are listed in table \ref{tab:pscan_results}.
No error analysis has yet been done on momentum bite, but the results
are in good agreement with the 2\% momentum bite estimated by Jaap
Doornbos.  (The width of the F1 horizontal slit defines the momentum
bite.  Jaap Doornbos calculates a momentum bite of (very roughly) 1\%
per 1.25 cm slit width; the slit was set to 2.5 cm wide for these
scans.)


\section{Positron Flux}
\label{sec:eflux}

The surface muon edge should have little to no effect on the M13
positron flux---the flux should increase approximately linearly with
momentum---so this flux can be used as a check.  Figures
\ref{fig:pscan_e_Be} and \ref{fig:pscan_e_C} show the normalized
positron flux as a function of B1 NMR.  The flux was normalized
against the T1 ion chamber as with the muon flux.
\begin{figure}[p]
  \begin{center}
    \includegraphics[angle=90, width=3in,
    height=3in]{pscan_positrons_001030.ps}
    \includegraphics[angle=90, width=3in,
    height=3in]{pscan_positrons_001205.ps}
    \caption[Normalized positron flux vs. B1 NMR reading, Be
    target..]{Normalized positron flux vs. B1 NMR reading, \textbf{Be
        target}.  Positron flux is normalized against the T1 Ion
      Chamber rate, which is proportional to the BL1A proton rate.}
    \label{fig:pscan_e_Be}
  \end{center}
\end{figure}
\begin{figure}[p]
  \begin{center}
    \includegraphics[angle=90, width=3in,
    height=3in]{pscan_positrons_001113.ps}
    \includegraphics[angle=90, width=3in,
    height=3in]{pscan_positrons_001120.ps}
    \caption[Normalized positron flux vs. B1 NMR reading, C
    target.]{Normalized positron flux vs. B1 NMR reading, \textbf{C
        target}.  Positron flux is normalized against the T1 Ion
      Chamber rate, which is proportional to the BL1A proton rate.}
    \label{fig:pscan_e_C}
  \end{center}
\end{figure}

As the graphs show, the positron flux does seem to scale linearly with
channel momentum.  There are kinks in the curves, however.  For the
November 13 data, these may be due to some hysteresis in the
quadrupoles, since the data points were not taken in a single series.
(The dipole magnets B1 and B2 should not be subject to hysteresis,
since the magnetic fields for those were set by NMR.)  

The kinks in the other scans are probably not due to hysteresis, since
the data was taken in order of momentum---any hysteresis should affect
all points in roughly the same way.  They could be due to some
non-linearity or scaling problems in the quadrupoles.

It is worth noting that most of the features of these curves are in
the same places on both graphs for a given target (though different
for the different targets).  This suggests that hysteresis is probably
not the cause.

\end{document}
Filename: pscan_apparatus.eps

Filename: pscan_muons_001030.ps

Filename: pscan_muons_001113.ps

Filename: pscan_muons_001120.ps

Filename: pscan_muons_001205.ps

Filename: surfacemu_closeup_001030data.ps

Filename: surfacemu_closeup_001113data.ps

Filename: surfacemu_closeup_001120data.ps

Filename: surfacemu_closeup_001205data.ps

Filename: naive_deriv_001030.ps

Filename: naive_deriv_001113.ps

Filename: naive_deriv_001120.ps

Filename: naive_deriv_001205.ps

Filename: pscan_positrons_001030.ps

Filename: pscan_positrons_001113.ps

Filename: pscan_positrons_001120.ps

Filename: pscan_positrons_001205.ps


TN-52: Summary of M13 Momentum Scans of Surface Muon Edge / Robert MacDonald

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