\documentstyle[12pt,epsfig]{article}

\setlength{\textwidth}{6.5in}
\setlength{\textheight}{8.5in}
%\setlength{\footskip}{0.35in}
\setlength{\topmargin}{0.0in}
\setlength{\oddsidemargin}{0in}
\setlength{\parindent}{5ex}
\setlength{\parskip}{.12in}
\newcommand{\detw}{\mbox{$\frac{\Delta E}{2}$}}
\newcommand{\hs}{\hspace*{5ex}}



%\def\decayright#1{\kern#1em\raise1.1ex\hbox{$|$}\kern-.5em\rightarrow}

\begin{document}


\begin{center}
{\Large
Summary of a Variables in passing\_events Tree
}
~\\
Amy Connolly\\
\today\\
\end{center}
~\\
~\\

\section{The Tree}

The tree called ``passing\_events'' is stored in outputs/icefinal.root, which
is written by \texttt{icemc}.  Only events that pass all three levels
of the Anita trigger are written to the tree.  When using this tree for
plotting, it is important to weight the events by \texttt{weight} (described
below).

All distances are in meters unless otherwise stated.

\begin{center}
  \begin{tabular}{p{2.2in}p{4.0in}}
\multicolumn{2}{p{6.0in}}{The next four variables use a 2d coordinate system 
from the view looking down on the South Pole, 
with the +x axis pointing in the direction of 
+90 longitude and +y in the 0$^{\circ}$ longitude direction }\\
\texttt{double horizcoord} & Position of the neutrino interaction in km along the +x axis \\
\texttt{double vertcoord} & Position of the neutrino interaction in km along the +y axis \\
\texttt{double horizcoord\_bn} & Position of the balloon in km along the +x axis \\
\texttt{double vertcoord\_bn} & Position of the balloon in km along the +y axis \\
\multicolumn{2}{p{6.0in}}{}\\
\texttt{double weight1} & Weight due to absorption in the earth.  Is it the probability that the neutrino would have survived its trip through the earth. \\
\texttt{double weight2} & A phase space factor.  For each balloon position, we choose the neutrino interaction position and direction such that the Cerenkov signal is visible at the balloon.  This phase space factor is supposed to remove the bias that this selection introduces. \\
\texttt{double weight} & Total weight for the neutrino interaction, which is \texttt{weight1$\times$ weight2}\\  
\texttt{double logweight} & log$_{10}$ of \texttt{weight} \\
\end{tabular}
\end{center}

\begin{center}
  \begin{tabular}{p{2.2in}p{4.0in}}
\multicolumn{2}{p{6.0in}}{The following arrays are defined by a coordinate system whose origin is at the center of the earth, +z is pointed towards the South Pole, +x is parallel to 90$^{\circ}$E at the South Pole and +y is parallel to 0$^{\circ}$ longitude at the South Pole.}\\
\texttt{double posnu[3]} & Position of the neutino interaction point \\
\texttt{double r\_bn[3]} & Position of the balloon \\
\texttt{double n\_bn[3]} & Normalized array pointing from center of the earth to the balloon \\
\texttt{double costheta\_nutraject} & $\cos{\theta_{\nu}}$ where $\theta_{\nu}$ is the zenith angle of the neutrino's incident trajectory \\
\texttt{double phi\_nutraject} & $\phi$ of the neutrino's incident trajectory\\
\texttt{double nnu[3]} & Normalized array pointing in the direction of the neutrino's incident trajectory\\
\texttt{double n\_exit2bn[5][3]} & Five 3-element arrays, each a unit vector pointing from the exit point to the balloon for the ray seen by the balloon.  We find the exit point in three iterations, and each array is the result of a different iteration.  So what you want is \texttt{n\_exit2bn[2]} because that is the third iteration.  The last two are zeros.\\
\texttt{double rfexit[5][3]} & Five 3-element arrays, each a unit vector pointing from the center of the earth to the point where the RF seen at the balloon exits the ice.  \\
\texttt{double n\_exit\_phi} & $\phi$ of \texttt{n\_exit2bn[2]} \\
\multicolumn{2}{p{6.0in}}{}\\
\texttt{double pnu} & Energy of the neutrino in eV\\
\texttt{double altitude\_int} & ``altitude'' of the interaction point relative to the ice surface (it is a negative number because the interaction is always below the surface)\\
\texttt{double elast\_y} & Inelasticity of the interaction (fraction of the neutrino's energy that goes into the hadronic shower) \\
\texttt{double emfrac} & Fraction of neutrino energy that goes into the EM component of the shower \\
\texttt{double hadfrac} & Fraction of the neutrino energy that goes into the hadronic component of the shower \\
\texttt{double sumfrac} & \texttt{emfrac+hadfrac}.  This can be less than unity since some of the energy can be carried away by a lepton.\\
\texttt{int nuflavor} & Neutrino flavor. 1=electron, 2=muon, 3=tau \\
\texttt{int current} & Type of interaction. 1=charged current, 2=neutral current. \\
\texttt{double viewangle} & Angle in the ice between the ray seen by the balloon and the neutrino direction in radians (actually in the firn just below the surface) \\
\texttt{double changle} & Cerenkov angle in the ice or firn, depending on where the interaction occured in radians\\
\texttt{double offaxis} & $|\texttt{viewangle-changle}|$ \\
\end{tabular}
\end{center}

\begin{center}
  \begin{tabular}{p{2.2in}p{4.0in}}
\texttt{int l1trig[3]} & For each trigger layer, which antennas pass L1.  16 bit, 16 bit and 8 bit for trigger layers 1 (top), 2 (bottom) and nadirs\\
\texttt{int l2trig[3][400]} & For each trigger layer, which set of three neighboring antennas pass L2.  16 bit, 16 bit and 8 bit for layers 1 (top) \& 2 (bottom) and nadirs.  Indexed according to the number of the centre antenna.  It is 400 long in the 2nd dimension so that the same array can be used for EeVA.  Only the relevant elements are filled. \\
\texttt{int l3trig[400]} & 16 bit number which says which phi sectors pass L3  \\
\texttt{double arrivaltimes[40]} & Time of arrival of the signal at each antenna, relative to the time the signal hits the first antenna \\
\texttt{double e\_component,h\_component} & Projection of the polarization of the signal along the horizontal (vertical) polarization of the antenna \\
\multicolumn{2}{p{6.0in}}{Each of the three arrays below are filled for the
antenna that is closest to the interaction (with the minimum arrival time).  The first index is the polarization (0=LCP, 1=RCP)}\\
\texttt{double rx0\_signal\_eachband[2][5]} & The signal strength in each band, including noise fluctutions.  Here the signal strength is $\int{V(f)} df$ where 
V(f) is the signal spectrum in Volts/Hz read by each channel 
and df is the frequency bin width. \\
\texttt{double rx0\_noise\_eachband[2][5]} & The average noise voltage that we 
use for each channel.\\
\texttt{double rx0\_threshold\_eachband[2][5]} & The threshold a signal
must exceed for a channel to pass,
in terms of a multiple of the average noise.\\
\texttt{double vmmhz\_max} & Electric field per frequency in V/m/MHz at FREQ\_HIGH=1200~MHz of the signal as it is incident on the payload\\
\texttt{double vmmhz\_min} & Electric field per frequency in V/m/MHz at FREQ\_LOW=200~MHz of the signal as it is incident on the payload\\
\texttt{volts\_rx\_rfcm\_lab\_e\_all[48][512]} & \\
 & Filled voltages in Volts for vertical polarization, 48 antennas and 512 samples each, with the
samples being every 1/2.6 ns so that the waveforms are 196.9 ns long\\
\texttt{volts\_rx\_rfcm\_lab\_h\_all[48][512]} & \\
 & Filled voltages in Volts for horizontal polarization, 48 antennas and 512 samples each, with the
samples being every 1/2.6 ns so that the waveforms are 196.9 ns long\\


\end{tabular}
\end{center}



\end{document}










