ericteunis/m-thesis-introduction
1
0
Fork 0
mirror of https://gitlab.science.ru.nl/mthesis-edeboone/m-thesis-introduction.git synced 2024-12-22 11:33:32 +01:00
Code Issues Projects Releases Packages Wiki Activity

ZH: allow to modify tracelengths in beacon_phases

Browse source
This commit is contained in:
Eric Teunis de Boone 2022-11-25 11:31:19 +01:00
parent dd00aebe3b
commit 22ccf7e3e3
1 changed files with 15 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

20
simulations/airshower_beacon_simulation/ba_beacon_phases.py
View file

@ -99,6 +99,14 @@ if __name__ == "__main__":
test_traces = traces[1:] test_traces = traces[1:]
orients = ['Ex', 'Ey', 'Ez', 'B'] orients = ['Ex', 'Ey', 'Ez', 'B']
# modify the length of the traces
if False:
t_trace = t_trace[:len(t_trace)//2]
half_traces = []
for trace in test_traces:
half_traces.append( trace[:len(trace)//2])
test_traces = half_traces
# Do Fourier Transforms # Do Fourier Transforms
# to find phases and amplitudes # to find phases and amplitudes
if True: if True:
@ -132,11 +140,13 @@ if __name__ == "__main__":
if show_plots and (i == 60 or i == 72): if show_plots and (i == 60 or i == 72):
fig, ax = plt.subplots() fig, ax = plt.subplots()
trace_amp = max(traces[-1]) - min(traces[-1]) ax.set_xlabel("t [ns]")
ax.set_ylabel("Amplitude")
myt = np.linspace(min(traces[0]), max(traces[0]), 10*len(traces[0])) myt = np.linspace(min(t_trace), max(t_trace), 10*len(t_trace))
ax.plot(t_trace, traces[-1], marker='.', label='trace') for i, trace in enumerate(test_traces):
ax.plot(myt, lib.sine_beacon(frequency, myt, amplitude=amplitude, phase=phase), ls='dashed', label='simulated') ax.plot(t_trace, test_traces[i], marker='.', label='trace '+orients[i])
ax.plot(myt, lib.sine_beacon(frequency, myt, amplitude=amplitude, phase=phase), ls='dashed', label='simulated beacon')
ax.set_title(f"Beacon at antenna {h5ant.attrs['name']}\nF:{frequency:.2e}, P:{phase:.4f}, A:{amplitude:.1e}") ax.set_title(f"Beacon at antenna {h5ant.attrs['name']}\nF:{frequency:.2e}, P:{phase:.4f}, A:{amplitude:.1e}")
ax.legend() ax.legend()
@ -178,4 +188,4 @@ if __name__ == "__main__":
ax.set_ylabel("Counts") ax.set_ylabel("Counts")
ax.hist(found_data[:,2], bins='sqrt', density=False) ax.hist(found_data[:,2], bins='sqrt', density=False)
plt.show() plt.show()