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: fixup ab1550c8: somehow removed beacon from trace generation

Browse source
This commit is contained in:
Eric Teunis de Boone 2022-11-22 15:57:00 +01:00
parent f6b474373d
commit 39cff267ac
1 changed files with 15 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

21
simulations/airshower_beacon_simulation/aa_generate_beacon.py
View file

@ -159,11 +159,20 @@ if __name__ == "__main__":
fname = "ZH_airshower/mysim.sry" fname = "ZH_airshower/mysim.sry"
tx = Antenna(x=-500,y=0,z=0,name='tx') # m tx = Antenna(x=-500,y=0,z=0,name='tx') # m
f_beacon = 51.53e-3 # GHz if not True: # slowest beacon to be found:
f_beacon = 10e-3 # GHz
low_bp = 5e-3 # GHz
high_bp = 80e-3 # GHz
else: # original wanted beacon
f_beacon = 51.53e-3 # GHz
# Bandpass for E field blockfilter # Bandpass for E field blockfilter
low_bp = 0.03 # GHz low_bp = 30e-3 # GHz
high_bp = 0.08 # GHz high_bp = 80e-3 # GHz
# Disable block_filter
if False:
block_filter = lambda x, dt, low, high: x
beacon_amplitudes = 1e-6*np.array([1e2, 0, 0]) # mu V/m beacon_amplitudes = 1e-6*np.array([1e2, 0, 0]) # mu V/m
beacon_radiate_rsq = True beacon_radiate_rsq = True
@ -205,8 +214,8 @@ if __name__ == "__main__":
# add to relevant polarisation # add to relevant polarisation
# and apply block filter # and apply block filter
dt = antenna.t[1] - antenna.t[0] dt = antenna.t[1] - antenna.t[0]
for j, _ in enumerate(beacon_amplitudes): for j, amp in enumerate(beacon_amplitudes):
traces[j] = block_filter(traces[j], dt, low_bp, high_bp) traces[j] = block_filter(traces[j] + amp*beacon, dt, low_bp, high_bp)
append_antenna_hdf5( antennas_fname, antenna, traces, name='traces', prepend_time=True, attrs_dict=dict(t0=t0)) append_antenna_hdf5( antennas_fname, antenna, traces, name='traces', prepend_time=True, attrs_dict=dict(t0=t0))