ericteunis/m-thesis-documentation
1
0
Fork 0
mirror of https://gitlab.science.ru.nl/mthesis-edeboone/m.internship-documentation.git synced 2024-11-22 23:13:35 +01:00
Code Issues Projects Releases Packages Wiki Activity

Thesis: Single Sine (WIP)

Browse source
This commit is contained in:
Eric Teunis de Boone 2023-10-12 12:23:35 +02:00
parent 9ff75c5629
commit 45d8a12031
3 changed files with 206 additions and 137 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

1
documents/thesis/chapters/beacon_discipline.tex
View file

@ -202,6 +202,7 @@ Chapter~\ref{sec:single_sine_sync} shows a special case of this last scenario wh
\\%>>> \\%>>>
\subsection{Array synchronisation}% <<< \subsection{Array synchronisation}% <<<
\label{sec:beacon:array}
% extending to array % extending to array
The idea of a beacon is to synchronise an array of antennas. The idea of a beacon is to synchronise an array of antennas.
As \eqref{eq:synchro_mismatch_clocks} applies for each pair of antennas in the array, all the antennas that record the beacon signal can determine the synchronisation mismatches simultaneously.% As \eqref{eq:synchro_mismatch_clocks} applies for each pair of antennas in the array, all the antennas that record the beacon signal can determine the synchronisation mismatches simultaneously.%

318
documents/thesis/chapters/single_sine_interferometry.tex
View file

@ -13,18 +13,89 @@
% <<< % <<<
As shown in Chapter~\ref{sec:disciplining}, both impulsive and sine beacon signals can synchronise air shower radio detectors to enable the interferometric reconstruction of extensive air showers. As shown in Chapter~\ref{sec:disciplining}, both impulsive and sine beacon signals can synchronise air shower radio detectors to enable the interferometric reconstruction of extensive air showers.
% period multiplicity/degeneracy This chapter will focus on using a single sine beacon to synchronise an array due to the simple setup and analysis required for such a beacon.
Although a sine beacon is the least intrusive, due to its periodicity, it can only synchronise two detectors up to its period Additionally, at \gls{Auger}, a public TV-transmitter is broadcasting at $67.25\MHz$.
(expressed as the $\Delta k_{ij}$ term in \eqref{eq:synchro_mismatch_clocks_periodic} This poses an opportunity to use a ``free'' beacon to synchronise the radio antennas of \gls{AERA} and \gls{AugerPrime}.
\footnote{%<<< \\
Note that $\Delta k_{ij}$ will be referenced in this chapter as $k_j$ since we can take station $i$ as the reference ($k_i = 0$).
}%>>> Due to the periodicity of sine beacons, the ability to synchronise an array is limited up to the beacon period $T$.
). As previously mentioned, the correct periods can be ascertained by choosing a beacon period much longer than the estimated accuracy of another timing mechanism.\footnote{For reference, \gls{GNSS} timing is expected to be below $30\ns$}
As previously mentioned, by choosing a beacon period much longer than the estimated accuracy of another timing mechanism, the correct periods can be ascertained. Likewise, this can be achieved using the beating of multiple frequencies such as the four frequency setup in \gls{AERA}, amounting to a total period of $>1\ns$..
\\ \\
In this chapter, a different method of resolving these period mismatches is investigated by recording an impulsive signal in combination with the sine beacon. In this chapter, a different method of resolving these period mismatches is investigated by recording an impulsive signal in combination with the sine beacon.
This extra signal declares a shared time $\tTrueEmit$ that is common to the stations, after which the periods can be counted (see Figure~\ref{fig:beacon_sync:sine}). Figure~\ref{fig:beacon_sync:sine} shows the steps of synchronisation using this combination.
The extra signal declares a shared time $\tTrueEmit$ that is common to the stations, after which the periods can be counted.
Note that the period mismatch term $\Delta k_{ij}$ in \eqref{eq:synchro_mismatch_clocks_periodic} will be referenced throughout this Chapter as $k$ since we can take station $i$ as reference ($k_i =0$).
\\
\begin{figure}%<<<
\centering
\begin{subfigure}[t]{0.47\textwidth}
\includegraphics[width=\textwidth]{beacon/beacon_sync.pdf}
\caption{
Phase alignment
}
\label{fig:beacon_sync:syntonised}
\end{subfigure}
\hfill
\begin{subfigure}[t]{0.47\textwidth}
\includegraphics[width=\textwidth]{beacon/beacon_sync_period.pdf}
\caption{
Period alignment
}
\label{fig:beacon_sync:period_alignment}
\end{subfigure}
\caption{
Synchronisation scheme for two antennas using a single sine wave beacon (orange) and an impulsive signal (blue).
Vertical dashed lines indicate periods of the sine wave.
\subref{fig:beacon_sync:syntonised} A small time delay $t_\varphi$ is derived from the phase difference of the beacon as measured by the antennas.
\subref{fig:beacon_sync:period_alignment} The period mismatch $k$ is determined from the overlap between the impulsive signals.
%Expecting the impulsive signals to come from the same source, the overlap between the impulsive signals is used to determine the period mismatch $k$.
Note that the impulsive signals do not coincide perfectly due to different propagation delays from the source to the antennas.
}
\label{fig:beacon_sync:sine}
% \begin{subfigure}{\textwidth}
% \centering
% \includegraphics[width=\textwidth]{beacon/08_beacon_sync_timing_outline.pdf}
% \caption{
% Measure two waveforms at different antennas at approximately the same local time (clocks are not synchronised).
% }
% \label{fig:beacon_sync:timing_outline}
% \end{subfigure}
% \begin{subfigure}{\textwidth}
% \centering
% \includegraphics[width=\textwidth]{beacon/08_beacon_sync_synchronised_outline.pdf}
% \caption{
% The beacon signal is used to remove time differences smaller than the beacon's period.
% The detector clocks are now an unknown amount of periods out of sync.
% }
% \label{fig:beacon_sync:syntonised}
% \end{subfigure}
% \begin{subfigure}{\textwidth}
% \centering
% \includegraphics[width=\textwidth]{beacon/08_beacon_sync_synchronised_period_alignment.pdf}
% \caption{
% Lifting period degeneracy ($k=n-m=7$ periods) using the optimal overlap between impulsive signals.
% \protect\Todo{note misaligned overlap due to different locations}
% }
% \label{fig:beacon_sync:period_alignment}
% \end{subfigure}
% \caption{
% Synchronisation scheme for two antennas using a continuous beacon and an impulsive signal, each emitted from a separate transmitter.
% Vertical dashed lines indicate periods of the beacon (orange),
% solid lines indicate the time of the impulsive signal (blue).
% \\
% \subref{fig:beacon_sync:syntonised}: The beacon allows to resolve a small timing delay ($\Delta \tClockPhase$).
% \\
% \subref{fig:beacon_sync:period_alignment}: Expecting the impulsive signals to come from the same source, the overlap between the two impulsive signals is used to lift the period degeneracy ($k=n-m$).
% }
% \label{fig:beacon_sync:sine}
% \protect\Todo{
% Redo figure without xticks and spines,
% rename $\Delta \tClockPhase$
% }
\end{figure}%>>>
% Same transmitter / Static setup % Same transmitter / Static setup
When the beacon transmitter is also used to emit the signal defining $\tTrueEmit$, the number of periods $k$ can be obtained directly from the signal. When the beacon transmitter is also used to emit the signal defining $\tTrueEmit$, the number of periods $k$ can be obtained directly from the signal.
@ -36,57 +107,15 @@ In dynamic setups, such as for transient signals, the time delays change per eve
The time delays must therefore be resolved from the information of a single event. The time delays must therefore be resolved from the information of a single event.
\\ \\
\Todo{text continuity}
% Beacon + Impulsive -> discrete % Beacon + Impulsive -> discrete
In a sine beacon synchronised array, finding this minimum deviation must control for the period defects. As shown in Chapter~\ref{sec:beacon:array}, an impulsive signal allows to reconstruct the direction of origin in a single event depending on the timing resolution of the array.
In general, these can be constrained using estimates of the accuracy of other timing mechanisms (see below). Synchronising the array with a sine beacon, any clock mismatch is discretized into a number of periods $k$.
This allows to improve the reconstruction by iterating the discrete clock mismatches during reconstruction.
\\ \\
With a restricted set of allowed period defects, we can alternate optimising the calibration signal's origin and optimising the set of period time delays of the array. Of course, a limit on the number of periods is required to prevent over-optimisation.
In general, they can be constrained using estimates of the accuracy of other timing mechanisms (see below).
\begin{figure}%<<<
\centering
\begin{subfigure}{\textwidth}
\centering
\includegraphics[width=\textwidth]{beacon/08_beacon_sync_timing_outline.pdf}
\caption{
Measure two waveforms at different antennas at approximately the same local time (clocks are not synchronised).
}
\label{fig:beacon_sync:timing_outline}
\end{subfigure}
\begin{subfigure}{\textwidth}
\centering
\includegraphics[width=\textwidth]{beacon/08_beacon_sync_synchronised_outline.pdf}
\caption{
The beacon signal is used to remove time differences smaller than the beacon's period.
The detector clocks are now an unknown amount of periods out of sync.
}
\label{fig:beacon_sync:syntonised}
\end{subfigure}
\begin{subfigure}{\textwidth}
\centering
\includegraphics[width=\textwidth]{beacon/08_beacon_sync_synchronised_period_alignment.pdf}
\caption{
Lifting period degeneracy ($k=n-m=7$ periods) using the optimal overlap between impulsive signals.
\Todo{note misaligned overlap due to different locations}
}
\label{fig:beacon_sync:period_alignment}
\end{subfigure}
\caption{
Synchronisation scheme for two antennas using a continuous beacon and an impulsive signal, each emitted from a separate transmitter.
Vertical dashed lines indicate periods of the beacon (orange),
solid lines indicate the time of the impulsive signal (blue).
\\ \\
\textit{Middle panel}: The beacon allows to resolve a small timing delay ($\Delta \tClockPhase$). With a restricted set of allowed period shifts, we can alternate optimising the calibration signal's origin and optimising the set of period time delays of the array.
\\
\textit{Lower panel}: Expecting the impulsive signals to come from the same source, the overlap between the two impulsive signals is used to lift the period degeneracy ($k=n-m$).
}
\label{fig:beacon_sync:sine}
\Todo{
Redo figure without xticks and spines,
rename $\Delta \tClockPhase$
}
\end{figure}%>>>
% >>> % >>>
\section{Lifting the Period Degeneracy with an Air Shower}% <<< \section{Lifting the Period Degeneracy with an Air Shower}% <<<
@ -96,65 +125,73 @@ In the case of radio detection of air showers, the very signal of the air shower
This falls into the dynamic setup previously mentioned. This falls into the dynamic setup previously mentioned.
The best period defects must thus be recovered from a single event. The best period defects must thus be recovered from a single event.
\\ \\
When doing the interferometric analysis for a sine beacon synchronised array, waveforms can only be delayed by an integer amount of periods, thereby giving discrete solutions to maximizing the interferometric signal. When doing the interferometric analysis for a sine beacon synchronised array, waveforms can only be delayed by an integer amount of periods, thereby giving discrete solutions to maximising the interferometric signal.
\Todo{add size of shower at plane vs period defects in meters} \Todo{add size of shower at plane vs period defects in meters}
\subsection{Air Shower simulation} \subsection{Air Shower simulation}
% simulation of proton E15 on 10x10 antenna % simulation of proton E15 on 10x10 antenna
To test the idea of combining a single sine beacon with an air shower, we simulated a set of recordings of a single air shower that also contains a beacon signal. To test the idea of combining a single sine beacon with an air shower, we simulated a set of recordings of a single air shower that also contains a beacon signal.
\footnote{\Todo{Url to repository}} \footnote{\url{https://gitlab.science.ru.nl/mthesis-edeboone/m-thesis-introduction/-/tree/main/airshower_beacon_simulation}}
\\ \\
The air shower signal (here a $10^{16}\eV$ proton) is simulated by \acrlong{ZHAires} on a grid of 10x10 antennas with a spacing of $50$\,meters.\Todo{cite ZHAires?} The air shower signal (here a $10^{16}\eV$ proton) was simulated by \acrlong{ZHAires} on a grid of 10x10 antennas with a spacing of $50$\,meters.\Todo{cite ZHAires?}
Each antenna recorded a waveform of $500$ samples with a sample rate of $1\GHz$. Each antenna recorded a waveform of $500$ samples with a samplerate of $1\GHz$ for each of the X,Y and Z polarisations.\Todo{polarisation names}
Figure~\ref{fig:single:proton_waveform} shows the earliest and latest waveforms recorded by the array with their true time. %Figure~\ref{fig:single:proton_waveform} shows the earliest and latest waveforms recorded by the array with their true time.
\Todo{verify numbers in paragraph} \Todo{verify numbers in paragraph}
\\ \\
Figure~\ref{fig:single:proton_grid} shows the maximum electric field measured at each of the antennas.
The ring of antennas with maximum electric fields in the order of $30\mathrm{\;\mu V/m}$ at the center of the array is the Cherenkov ring.
The Cherenkov ring forms due to the forward beaming of the radio emissions of the airshower.\Todo{expand/refer}
Outside this ring, the maximum electric field quickly falls with increasing distance to the array core.
\Todo{geometry, vertical shower}
\\
%% add beacon %% add beacon
A sine beacon ($\fbeacon = 51.53\MHz$) is introduced at a distance of approximately $75\mathrm{\,km}$ northwest of the array. A sine beacon ($\fbeacon = 51.53\MHz$) was introduced at a distance of approximately $75\mathrm{\,km}$ northwest of the array, primarily received in the X polarisation.
The distance between the antenna and the transmitter results in a phase offset with which the beacon is received at each antenna.% The distance between the antenna and the transmitter results in a phase offset with which the beacon is received at each antenna.%
\footnote{%<<< \footnote{%<<<
The beacon's amplitude is also dependent on the distance. Although simulated, this effect has not been incorporated in the analysis as it is neglible for the considered grid and distance to the transmitter. The beacon's amplitude is also dependent on the distance. Although simulated, this effect has not been incorporated in the analysis as it is negligible for the considered grid and distance to the transmitter.
} %>>> } %>>>
The beacon signal is recorded over a longer time, to be able to distinguish the beacon and air shower later in the analysis. The beacon signal was recorded over a longer time ($10240$\;samples), to be able to distinguish the beacon and air shower later in the analysis.
\\ \\
The final waveform of an antenna (see Figure~\ref{fig:single:annotated_full_waveform}) is then constructed by adding its beacon and air shower waveforms and bandpassing with relevant frequencies (here $30$ and $80\MHz$ are taken by default). The final waveform of an antenna (see Figure~\ref{fig:single:annotated_full_waveform}) was then constructed by adding its beacon and air shower waveforms and band-passing with relevant frequencies (here $30$ and $80\MHz$ are taken by default).
Of course, a gaussian white noise component can be introduced to the waveform as a simple noise model (see Figure~\ref{fig:sine:time_accuracy} for a treatise on the timing accuracy of a sine beacon). Of course, a gaussian white noise component is introduced to the waveform as a simple noise model (see Figure~\ref{fig:sine:time_accuracy} for a treatise on the timing accuracy of a sine beacon).
\\ \\
\begin{figure}% <<< \begin{figure}% <<<
\centering \centering
%\begin{subfigure}{0.47\textwidth} \includegraphics[width=0.5\textwidth]{ZH_simulation/array_geometry_shower_amplitude.pdf}
% \includegraphics[width=\textwidth]{ZH_simulation/array_geometry_shower_amplitude.pdf} \caption{
% \caption{} \protect\Todo{text}
% \label{fig:single:proton_grid} The 10x10 antenna grid used for recording the air shower.
%\end{subfigure} Colours indicate the maximum electric field recorded at the antenna.
%\hfill }
\begin{subfigure}[t]{0.47\textwidth} \label{fig:single:proton_grid}
\includegraphics[width=\textwidth]{ZH_simulation/first_and_last_simulated_traces.pdf} \end{figure}% >>>
\caption{} \begin{figure}% <<<
\label{fig:single:proton_waveform} \begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{ZH_simulation/ba_measure_beacon_phase.py.A74.no_mask.zoomed.pdf}
\label{fig:single:annotated_full_waveform}
\end{subfigure} \end{subfigure}
\hfill \hfill
\begin{subfigure}[t]{0.47\textwidth} \begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{ZH_simulation/ba_measure_beacon_phase.py.A74.no_mask.zoomed.pdf} \includegraphics[width=\textwidth]{ZH_simulation/ba_measure_beacon_phase.py.A74.fourier.pdf}
\caption{} \label{fig:single:fourier}
\label{fig:single:annotated_full_waveform}
\end{subfigure} \end{subfigure}
\caption{ \caption{
\textit{Left:} \textit{Left:}
%The 10x10 antenna grid used for recording the air shower.
%Colours indicate the maximum electric field recorded at the antenna.
%\textit{Right:} %\textit{Right:}
Example of the earliest and latest recorded air shower waveforms in the array as simulated by ZHAires. %Example of the earliest and latest recorded air shower waveforms in the array as simulated by ZHAires.
Excerpt of a fully simulated waveform ($N=10240\mathrm{\;samples}$) (blue) containing the air shower (a $10^{16}\eV$~proton), the beacon (green, $\fbeacon = 51.53\MHz$) and noise.
The part of the waveform between the vertical dashed lines is considered airshower signal and masked before measuring the beacon parameters.
\textit{Right:} \textit{Right:}
Excerpt of a fully simulated waveform (blue) containing the air shower (a $10^{16}\eV$~proton), the beacon (orange, $\fbeacon = 51.53\MHz$) and noise. Fourier spectra of the waveforms.
\Todo{indicate cuts?} The green dashed lines indicate the measured beacon parameters.
} }
\label{fig:single:proton} \label{fig:single:proton}
\end{figure}% >>> \end{figure}% >>>
% randomise clocks % randomise clocks
After the creation of the antenna waveforms, the clocks are randomised by sampling a gaussian distribution with a standard deviation of $30\ns$ (see Figure~\ref{fig:simu:sine:periods:repair_none}). After the creation of the antenna waveforms, the clocks are randomised by sampling a gaussian distribution with a standard deviation of $30\ns$. % (see Figure~\ref{fig:simu:sine:periods:repair_none}).
At a beacon period of $\sim 20\ns$, this ensures that multiple antennas have clock defects of at least one beacon period. At a beacon period of $\sim 20\ns$, this ensures that multiple antennas have clock defects of at least one beacon period.
This in turn allows for synchronisation mismatches of more than one beacon period. This in turn allows for synchronisation mismatches of more than one beacon period.
Moreover, it falls in the order of magnitude of clock defects that were found in \gls{AERA}\cite{PierreAuger:2015aqe}. Moreover, it falls in the order of magnitude of clock defects that were found in \gls{AERA}\cite{PierreAuger:2015aqe}.
@ -162,28 +199,29 @@ Moreover, it falls in the order of magnitude of clock defects that were found in
% separate air shower from beacon % separate air shower from beacon
To correctly recover the beacon from the waveform, it must be separated from the air shower. To correctly recover the beacon from the waveform, it must be separated from the air shower.
In Figure~\ref{fig:single:annotated_full_waveform} it is readily identified at $t=500\ns$. Typically, a trigger sets the location of the airshower signal in the waveform.
In our case, the airshower signal is located at $t=500\ns$ (see Figure~\ref{fig:single:annotated_full_waveform}).
Since the beacon can be recorded for much longer than the air shower signal, we mask a window of $500$ samples around the maximum of the trace as the air shower's signal. Since the beacon can be recorded for much longer than the air shower signal, we mask a window of $500$ samples around the maximum of the trace as the air shower's signal.
% measure beacon phase, remove distance phase % measure beacon phase, remove distance phase
The remaining waveform is fed into a \gls{DTFT} \eqref{eq:fourier:dtft} to measure the beacon's phase $\pMeas$ and amplitude. The remaining waveform is fed into a \gls{DTFT} \eqref{eq:fourier:dtft} to measure the beacon's phase $\pMeas$ and amplitude.
Note that due to explicitly including a time axis in a \gls{DTFT}, a number of samples can be omitted without introducing artifacts.
\\ \\
With the obtained beacon parameters, the air shower signal is in turn reconstructed by subtracting the beacon from the full waveform in the time domain. With the obtained beacon parameters, the air shower signal is in turn reconstructed by subtracting the beacon from the full waveform in the time domain.\Todo{cite GAP note?}
\\ \\
The small clock defect $\tSmallClock$ is then finally calculated from the beacon's phase $\pMeas$ by subtracting the phase $\pProp$ introduced by the propagation from the beacon transmitter. The small clock defect $\tClockPhase$ is then finally calculated from the beacon's phase $\pMeas$ by subtracting the phase introduced by the propagation from the beacon transmitter.
\\ \\
% introduce air shower % introduce air shower
From the above, we now have a set of reconstructed air shower waveforms with corresponding clock defects smaller than one beacon period $T$. From the above, we now have a set of reconstructed air shower waveforms with corresponding clock defects smaller than one beacon period $T$.
Shifting the waveforms to remove these small clocks defects, we are left with resolving the correct number of periods $k$ per waveform (see Figure~\ref{fig:simu:sine:periods:repair_phases}). Shifting the waveforms to remove these small clocks defects, we are left with resolving the correct number of periods $k$ per waveform (see Figure~\ref{fig:grid_power:repair_phases}).
\\ \\
\begin{figure}%<<< \begin{figure}%<<<
\centering \centering
\includegraphics[width=0.8\textwidth]{ZH_simulation/findks/ca_period_from_shower.py.run0.i1.kfind.zoomed.peak.pdf} \includegraphics[width=0.8\textwidth]{ZH_simulation/findks/ca_period_from_shower.py.run0.i1.kfind.pdf}
\caption{ \caption{
Finding the maximum correlation for integer period shifts (up to $k=3$) between two waveforms recording the same (simulated) air shower. Finding the maximum correlation for integer period shifts (best $k=-1$) between two waveforms recording the same (simulated) air shower.
Randomising the antenna clocks up to $30\ns$ and $\fbeacon = 51.53\MHz$ corresponds to at most $3$ periods of time difference between two waveforms. Randomising the antenna clocks up to $30\ns$ and $\fbeacon = 51.53\MHz$ corresponds to at most $3$ periods of time difference between two waveforms.
\Todo{location origin}
} }
\label{fig:single:k-correlation} \label{fig:single:k-correlation}
\end{figure}%>>> \end{figure}%>>>
@ -192,23 +230,24 @@ Shifting the waveforms to remove these small clocks defects, we are left with re
\subsection{\textit{k}-finding} % <<< \subsection{\textit{k}-finding} % <<<
% unknown origin of air shower signal % unknown origin of air shower signal
Up until now, the shower axis and thus the origin of the air shower signal here have not been resolved. Up until now, the shower axis and thus the origin of the air shower signal have not been resolved.
This means that the unknown propagation time delays for the air shower ($\tProp$) affect the alignment of the signals in Figure~\ref{fig:beacon_sync:period_alignment} in addition to the unknown clock period defects ($k_j T$). This means that the unknown propagation time delays for the air shower ($\tProp$) affect the alignment of the signals in Figure~\ref{fig:beacon_sync:period_alignment} in addition to the unknown clock period defects ($k_j T$).
As such, both this origin and the clock defects have to be determined simultaneously. As such, both this origin and the clock defects have to be determined simultaneously.
\\ \\
% radio interferometry % radio interferometry
If the antennas had been fully synchronised, radio interferometry as introduced in Section~\ref{sec:interferometry} can be applied to find the origin of the air shower signal, thus resolving the shower axis. If the antennas had been fully synchronised, radio interferometry as introduced in Section~\ref{sec:interferometry} can be applied to find the origin of the air shower signal, thus resolving the shower axis.
Still, a rough first estimate of the shower axis might be made using this or other techniques (see Figure~\ref{fig:dynamic-resolve}). Still, a rough first estimate of the shower axis might be made using this or other techniques.% (see Figure~\ref{fig:dynamic-resolve}).
\\ \\
Starting with an initial grid around this estimated axis, a two-step process zooms in on the shower axis while optimising the interferometric signal. Starting with an initial grid around this estimated axis, a two-step process zooms in on the shower axis while optimising the interferometric signal.
In this process each waveform of the array is allowed to shift by a restricted amount of periods with respect to a reference waveform (taken to be the waveform with the highest maximum). In this process each waveform of the array is allowed to shift by a restricted amount of periods with respect to a reference waveform (taken to be the waveform with the highest maximum).
\\ \\
Note that the grids in Figure~\ref{fig:findks} are defined in shower plane coordinates with $\vec{v}$ the shower axis and \vec{B} the local magnetic field. Note that these grids are defined in shower plane coordinates with $\vec{v}$ the true shower axis and $\vec{B}$ the local magnetic field.
Searching a grid that is slightly misaligned with the true shower axis is expected to give comparable results.
\\ \\
The first step consists of finding the set of period shifts $k_j$ that maximizes the overall maximum amplitude on the grid. The first step consists of finding the set of period shifts $k_j$ that maximises the overall maximum amplitude on the grid.
At each location, after removing propagation delays, each waveform and the reference waveform are summed while varying the time delay $kT$ ($\left| k \right| \leq 3$ in Figure~\ref{fig:single:k-correlation}) to find the maximum attainable amplitude of this combined trace. At each location, after removing propagation delays, each waveform and the reference waveform are summed with different time delays $kT$ to find the maximum attainable amplitude of the combined trace.%
\footnote{%<<< \footnote{%<<<
Note that one could use a correlation method instead of a maximum to select the best time delay. Note that one could use a correlation method instead of a maximum to select the best time delay.
However, for simplicity and ease of computation, this has not been implemented. However, for simplicity and ease of computation, this has not been implemented.
@ -217,12 +256,15 @@ At each location, after removing propagation delays, each waveform and the refer
%Since it is this peak that is of interest, it would have been possible to cut the waveforms such to only correlate the peaks. %Since it is this peak that is of interest, it would have been possible to cut the waveforms such to only correlate the peaks.
} %>>> } %>>>
\\ \\
As shown in Figure~\ref{fig:single:k-correlation}, the maximum possible period shift has been limited to $\pm 3\mathrm{\; periods}$.
This corresponds to the maximum expected time delay between two antennas with a clock randomisation up to $30\ns$ for the considered beacon frequency.
\\
% %
This amplitude optimisation is iterated over the grid (see Figure~\ref{fig:findks:maxima}) resulting in a grid measurement of the maximum attainable amplitude and the corresponding set of period defects $k$. This amplitude optimisation is iterated over the grid (see Figure~\ref{fig:findks:maxima}) resulting in a grid measurement of the maximum attainable amplitude and the corresponding set of period defects $k$.
\\ \\
The second step then uses the period defects belonging to the highest maximum amplitude to measure the interferometric power on the same grid (see Figure~\ref{fig:findks:reconstruction}). The second step applies the obtained period defects and measures the interferometric power on the same grid (see Figure~\ref{fig:findks:reconstruction}).
Afterwards, a new grid zooms in on the power maximum and the process is repeated (Figures~\ref{fig:findks:maxima:zoomed} and \ref{fig:findks:reconstruction:zoomed}) until the set of period defects does not change between grids. Afterwards, a new grid zooms in on the power maximum and the process is repeated (Figures~\ref{fig:findks:maxima:zoomed},~\ref{fig:findks:reconstruction:zoomed}) until the set of period defects is constant between zooming grids.
\\ \\
\begin{figure}%<<< findks \begin{figure}%<<< findks
@ -257,42 +299,63 @@ Afterwards, a new grid zooms in on the power maximum and the process is repeated
} }
\label{fig:findks:reconstruction:zoomed} \label{fig:findks:reconstruction:zoomed}
\end{subfigure} \end{subfigure}
\\
\begin{subfigure}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{ZH_simulation/findks/ca_period_from_shower.py.maxima.run2.pdf}
\caption{
Final test grid obtaining the same $k$ as Figure~\ref{fig:findks:maxima:zoomed}.
}
\label{fig:findks:maxima:zoomed2}
\end{subfigure}
\hfill
\begin{subfigure}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{ZH_simulation/findks/ca_period_from_shower.py.reconstruction.run2.power.pdf}
\caption{
Final power measurement.
}
\label{fig:findks:reconstruction:zoomed2}
\end{subfigure}
\caption{ \caption{
Iterative $k$-finding algorithm: Iterative $k$-finding algorithm:
First (\textit{upper left}), find the set of period shifts $k$ per point on a grid that returns the highest maximum amplitude (blue cross). First \subref{fig:findks:maxima}, find the set of period shifts $k$ per point on a grid that returns the highest maximum amplitude (blue cross).
The grid starts as a $8^\circ$ wide shower plane slice at $X=400\mathrm{\,g/cm}$, centered at the true shower axis (red cross). The grid starts as a $8^\circ$ wide shower plane slice at $X=400\mathrm{\,g/cm}$, centred at the true shower axis (red cross).
Second (\textit{upper right}), perform the interferometric reconstruction with this set of period shifts. Second \subref{fig:findks:reconstruction}, perform the interferometric reconstruction with this set of period shifts.
Finally (\textit{lower panes}), zooming in on the maximum power of the reconstruction, repeat the steps until the set of period shifts does not change. Zooming on the maximum power \subref{fig:findks:maxima:zoomed},\subref{fig:findks:reconstruction:zoomed} repeat the steps until the $k$'s are equal between the zoomed grids \subref{fig:findks:maxima:zoomed2},\subref{fig:findks:reconstruction:zoomed2}.
} }
\label{fig:findks} \label{fig:findks}
\end{figure}%>>> \end{figure}%>>>
% >>> % >>>
\clearpage
%\phantomsection %\phantomsection
\subsection{Strategy / Result} %<<< \subsection{Strategy / Result} %<<<
The initial grid plays an important role in finding the correct axis. Figure~\ref{fig:grid_power_time_fixes} shows the effect of the various synchronisation stages on both the alignment of the air shower waveforms, and the interferometric power measurement near the true shower axis.
\\
The initial grid plays an important role here in finding the correct axis.
Due to selecting the highest maximum amplitude, and the process above zooming in aggressively, wrong candidate axes are selected when there is no grid-location sufficiently close to the true axis. Due to selecting the highest maximum amplitude, and the process above zooming in aggressively, wrong candidate axes are selected when there is no grid-location sufficiently close to the true axis.
Figure~\ref{fig:findks:reconstruction} shows such a potential point near $(-1, 0.5)$ with a maximum that is comparable to the selected maximum near the true axis.
\\ \\
% premature optimisation / degeneracy % premature optimisation / degeneracy
Such locations are subject to differences in propagation delays that are in the order of a few beacon periods. Such locations are subject to differences in propagation delays that are in the order of a few beacon periods.
The restriction of the possible delays is therefore important to limit the number of potential axes. The restriction of the possible delays is therefore important to limit the number of potential axis locations.
\\ \\
% fall in local extremum, maximum % fall in local extremum, maximum
As the number of computations scales linearly with the number of grid points ($N = N_x N_y$), it is favourable to minimise the number of grid locations.
Unfortunately, the above process has been observed to fall into local maxima when a too coarse initial grid ($N_x < 13$ points over $8^\circ$ around the true axis\Todo{words}) was used while restricting the time delays to $\left| k \right| \leq 3$.
\\
In this analysis, the initial grid is defined as $8^\circ$ wide around the true axis.\Todo{why?} In this analysis, the initial grid is defined as $8^\circ$ wide around the true axis.\Todo{why?}
Figure~\ref{fig:findks:reconstruction} shows such a potential point near $(-1, 0.5)$ with a maximum that is comparable to the selected maximum near the true axis. As the number of computations scales linearly with the number of grid points ($N = N_x N_y$), it is favourable to minimise the number of grid locations.
Unfortunately, the above process has been observed to fall into local maxima when a too coarse initial grid ($N_x < 13$ at $X=400\mathrm{\,g/cm}$) was used while restricting the time delays to $\left| k \right| \leq 3$.
\bigskip
Figures~\ref{fig:simu:sine:periods, fig:grid_power_time_fixes} show the effect of the various synchronisation stages on both the alignment of the air shower waveforms, and the interferometric power measurement near the true shower axis.
\\ \\
Of course, this algorithm must be evaluated at proper\Todo{word} atmospheric depths where the interferometric technique can resolve the shower.
Additionally, since the true period shifts are static per event, the evaluation at multiple atmospheric depths allows to compare the obtained sets thereof.
\\
\Todo{directional reconstruction}
\Todo{missing power}
\Todo{res k Figure}
%\begin{figure}%<<< fig:simu:sine:periods %\begin{figure}%<<< fig:simu:sine:periods
% \centering % \centering
% \begin{subfigure}[t]{0.45\textwidth} % \begin{subfigure}[t]{0.45\textwidth}
@ -339,7 +402,7 @@ Figures~\ref{fig:simu:sine:periods, fig:grid_power_time_fixes} show the effect o
\begin{subfigure}[t]{0.9\textwidth} \begin{subfigure}[t]{0.9\textwidth}
\includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.repair_none.scale4d.pdf} \includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.repair_none.scale4d.pdf}
\hfill \hfill
\includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.repair_none.axis.trace_overlap.repair_none.pdf} \includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.repair_none.axis.X400.trace_overlap.repair_none.pdf}
\caption{ \caption{
Randomised clocks Randomised clocks
} }
@ -350,9 +413,9 @@ Figures~\ref{fig:simu:sine:periods, fig:grid_power_time_fixes} show the effect o
\begin{subfigure}[t]{0.9\textwidth} \begin{subfigure}[t]{0.9\textwidth}
\includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.repair_phases.scale4d.pdf} \includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.repair_phases.scale4d.pdf}
\hfill \hfill
\includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.repair_phases.axis.trace_overlap.repair_phases.pdf} \includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.repair_phases.axis.X400.trace_overlap.repair_phases.pdf}
\caption{ \caption{
Clock syntonisation Phase synchronisation
} }
\label{fig:grid_power:repair_phases} \label{fig:grid_power:repair_phases}
\end{subfigure} \end{subfigure}
@ -360,7 +423,7 @@ Figures~\ref{fig:simu:sine:periods, fig:grid_power_time_fixes} show the effect o
\begin{subfigure}[t]{0.9\textwidth} \begin{subfigure}[t]{0.9\textwidth}
\includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.repair_all.scale4d.pdf} \includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.repair_all.scale4d.pdf}
\hfill \hfill
\includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.repair_full.axis.trace_overlap.repair_full.pdf} \includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.repair_full.axis.X400.trace_overlap.repair_full.pdf}
\caption{ \caption{
Fully resolved clocks Fully resolved clocks
} }
@ -370,20 +433,19 @@ Figures~\ref{fig:simu:sine:periods, fig:grid_power_time_fixes} show the effect o
\begin{subfigure}[t]{0.9\textwidth} \begin{subfigure}[t]{0.9\textwidth}
\includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.no_offset.scale4d.pdf} \includegraphics[width=0.45\textwidth]{radio_interferometry/dc_grid_power_time_fixes.py.X400.no_offset.scale4d.pdf}
\hfill \hfill
\includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.no_offset.axis.trace_overlap.no_offset.pdf} \includegraphics[width=0.46\textwidth]{radio_interferometry/trace_overlap/on-axis/dc_grid_power_time_fixes.py.no_offset.axis.X400.trace_overlap.no_offset.pdf}
\caption{ \caption{
True clocks True clocks
} }
\label{fig:grid_power:no_offset} \label{fig:grid_power:no_offset}
\end{subfigure} \end{subfigure}
\caption{ \caption{
Different stages of array synchronisation (unsynchronised, beacon synchronised, $k$-resolved and true clocks). Different stages of array synchronisation (unsynchronised, beacon synchronised, $k$-resolved and true clocks)
\allowbreak and
\textit{Right:} The overlap between the array's waveforms at the true axis. their effect on (\textit{right}) the alignment of the waveforms at the true axis
\allowbreak and (\textit{left}) the interferometric power near the simulation axis (red plus).
\textit{Left:} Power measurements near the simulation axis (red plus).
The maximum power is indicated by the blue cross. The maximum power is indicated by the blue cross.
% \Todo{x-axis relative to reference waveform, remove titles, no SNR} % \protect\Todo{x-axis relative to reference waveform, remove titles, no SNR}
} }
\label{fig:grid_power_time_fixes} \label{fig:grid_power_time_fixes}
\end{figure}%>>> \end{figure}%>>>

6
documents/thesis/preamble.tex
View file

@ -101,6 +101,10 @@
\newcommand{\GeV}{\text{\,G\eV}} \newcommand{\GeV}{\text{\,G\eV}}
\newcommand{\TeV}{\text{\,T\kern-0.1ex\eV}} \newcommand{\TeV}{\text{\,T\kern-0.1ex\eV}}
\newcommand{\metre}{\text{\,m}}
\newcommand{\mV}{\text{\,mV}}
\newcommand{\ns}{\text{\,ns}} \newcommand{\ns}{\text{\,ns}}
\newcommand{\us}{\text{\,\textmu s}} \newcommand{\us}{\text{\,\textmu s}}
@ -158,6 +162,8 @@
\newacronym{GNSS}{GNSS}{Global Navigation Satellite System} \newacronym{GNSS}{GNSS}{Global Navigation Satellite System}
\newacronym{EAS}{EAS}{Extensive Air Shower} \newacronym{EAS}{EAS}{Extensive Air Shower}
\newacronym{UHECR}{UHECR}{Ultra High Energy Cosmic Ray} \newacronym{UHECR}{UHECR}{Ultra High Energy Cosmic Ray}
\newacronym{DU}{DU}{Detector Unit}
\newacronym{PPS}{PPS}{Pulse Per Second}
%% <<<< Names %% <<<< Names
\newacronym{GRAND}{GRAND}{Giant Radio Array for Neutrino Detection} \newacronym{GRAND}{GRAND}{Giant Radio Array for Neutrino Detection}