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: mask the beacon where we have signal (optional)

Browse source
This commit is contained in:
Eric Teunis de Boone 2022-12-08 12:02:48 +01:00
parent b11e6f3756
commit 888d81c1f6
1 changed files with 14 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
simulations/airshower_beacon_simulation/ba_beacon_phases.py
View file

@ -87,6 +87,20 @@ if __name__ == "__main__":
test_traces = [ traces[1] ] test_traces = [ traces[1] ]
orients = ['E_AxB'] orients = ['E_AxB']
# TODO: refine masking
# use beacon but remove where E_AxB-Beacon != 0
if True:
if True:
t_mask = np.isclose(h5ant['E_AxB'][1], h5ant['traces'][4], rtol=1e-3, atol=1e-3)
else:
t_mask = np.ones(len(t_trace), dtype=bool)
t_mask[1500:3000] = False # magic numbers from aa_generate_beacon
t_trace = t_trace[t_mask]
for j, t in enumerate(test_traces):
test_traces[j] = t[t_mask]
orients[j] = orients[j] + ' masked'
# use separate polarisations # use separate polarisations
else: else:
if 'traces' not in h5ant.keys(): if 'traces' not in h5ant.keys():