\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 tsignals Tree
}
~\\
Amy Connolly\\
\today\\
\end{center}
~\\
~\\

\section{The tsignals Tree}

The tree called ``tsignals'' is stored in outputs/signals.root, which
is written by \texttt{icemc}.  Only events that survive any preselection cuts
to reach the Anita trigger are written to this tree.
Remember that as with the \texttt{passing\_events} tree, each event
has a weight given by \texttt{weight\_test/dnutries} to account for earth absorption and
selection of the interaction position within the horizon of the payload.
This is only filled for the antenna that is closest to the interaction, which
is always on the lowest of the three layers of the antennas (there are no
drop-downs implemented yet).  For all waveforms, there are 512 samples separated
by (1/2.6)~ns.

Where I refer to the first polarization, I mean LCP for Anita~I and vertical polarization for
Anita~II.  And for second polarization I mean RCP for Anita~I and horizontal polarization for
Anita~II.

A channel passes the trigger if the output of the tunnel diode\\
(\texttt{timedomain\_output\_1\_inanita[j]} or 
\texttt{timedomain\_output\_1\_inanita[j]})
goes lower than the (negative) 
threshold within the time window that lies between the 
232$^{\mathrm{nd}}$ and 284$^{\mathrm{th}}$ bins, or between about
89~ns and 109~ns.  Notice that this window
is centered about 15~ns later than the center of the 
signal waveforms since under my
tunnel diode model it imposes a delay.
The thresholds for ANITA~I are -1.251E-14 (Lo),-1.13106E-14 (Mid1),
-7.98583E-15 (Mid2),-1.40634E-14 (Hi) and -3.55069E-11 (full when used).
For ANITA~II, they are -1.15638e-14 (Lo), -1.42725e-14 (Mid),
-2.20845e-14 (Hi), -2.71869e-14 (Full).


\begin{center}
  \begin{tabular}{p{2.2in}p{4.0in}}
\texttt{double dangle} & Viewing angle - Cerenkov angle \\
\texttt{double emfrac} & EM fraction of the shower \\
\texttt{double hadfrac} & Hadronic fraction of the shower \\
\texttt{int inu} & \texttt{icemc}'s event number \\
\multicolumn{2}{p{6.0in}}{The event weight is the quotient \texttt{weight\_test/dnutries}}\\
\texttt{double dnutries} &  The number of ``tries'' for the interaction position across Antarctica before it would have landed within the horizon of the balloon.\\
\texttt{double weight\_test} & Weight due to absorption in the earth.  Is it the probability that the neutrino would have survived its trip through the earth.  When we fill this tree, we haven't stepped through the Earth using the full geographical model to calculate the absorption yet, so instead we just assume the 
neutrino only traverses ice on its way to the interaction.\\
\texttt{double peak\_e\_rx} & Peak of the signal waveform at the antenna receiver, first polarization.  No RFCM amplification or any attenuation of any sort imposed.  Just antenna effective heights imposed on the incident signal.\\
\texttt{double peak\_h\_rx} & Same as previous, but second polarization.\\
\texttt{double volts\_rx\_rfcm\_lab\_e[512]} & Time domain {\it signal} waveform read by the surf in the first polarization\\
\texttt{double peak\_e\_rx\_rfcm\_lab} & Peak of the previous waveform \\
\texttt{double volts\_rx\_rfcm\_lab\_h} & Time domain {\it signal} waveform read by the surf in the second polarization\\
\texttt{double peak\_h\_rx\_rfcm\_lab} & Peak of the previous waveform \\
\end{tabular}
\end{center}

\begin{center}
  \begin{tabular}{p{3.0in}p{3.2in}}
\texttt{double signal\_vpol\_inanita[5][512]} & Time domain {\it signal} waveforms (voltage v. time) for each of the five channels and each of the 512 samples in the vertical polarization only, before coverting to LCP, RCP for ANITA~I. \\
\texttt{double noise\_vpol\_inanita[5][512]} & Same as previous, but just the {\it noise} waveform in the vertical polarization only, with no signal.  This is before converting to LCP, RCP for ANITA~I. \\
\texttt{double total\_vpol\_inanita[5][512]} & Sum of the previous two waveforms. \\
\texttt{double total\_diodeinput\_1\_inanita[5][512]} & Waveform that is input to the tunnel diode for the first polarization.\\
\texttt{double total\_diodeinput\_2\_inanita[5][512]} & Waveform that is input to the tunnel diode for the second polarization.\\
\texttt{double peak\_e[5]} & Peak voltage of the signal waveform in the first polarization \\
\texttt{double peak\_h[5]} & Peak voltage of the noise waveform in the second polarization \\
\texttt{double timedomain\_output\_1\_inanita[5][512]} & Output waveforms of the tunnel diodes for each channel (voltage v. time) for the first polarization \\
\texttt{double timedomain\_output\_2\_inanita[5][512]} & Output waveforms of the tunnel diodes for each channel (voltage v. time) for the second polarization \\
\texttt{int flag\_e\_inanita[5][512]} & 1 when the channel has fired and 0 otherwise, for the first polarization.  This flag is held high for a L1 trigger coincidence window of 11.19~ns. \\
\texttt{int flag\_h\_inanita} & same as previous, but for the second polarization. \\
\texttt{double bwslice\_vrms[5]} & RMS noise voltage for each channel \\
\texttt{double ston[5]} & This is the peak signal voltage divided by rms noise voltage for each channel (\texttt{peak\_e[i]/bwslice\_vrms[i]}).  Not that this is not what is input to the discriminator, it is just one measure of signal strength. \\
\texttt{int channels\_passing\_e[5]} & Which channels in the first polarization pass the single-channel trigger at some time.  1 for yes, 0 for no. \\
\texttt{int channels\_passing\_h[5]/I} & Same as previous, but second polarization. \\
\texttt{int l1\_passing} & Whether this antenna passes the L1 trigger. 1 for yes, 0 for no. \\
 \texttt{double integral\_vmmhz} & Integral of the electric field in V/m/MHz incident on the closest antenna.  For the older frequency domain simulation, this was the measure of signal strength. \\
\end{tabular}
\end{center}




\end{document}

\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}