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: use full band for noise power determination

Browse source
This commit is contained in:
Eric Teunis de Boone 2023-01-31 14:52:27 +01:00
parent 82ad45730e
commit 0641048f67
1 changed files with 6 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
simulations/airshower_beacon_simulation/ac_show_signal_to_noise.py
View file

@ -70,14 +70,14 @@ if __name__ == "__main__":
## Beacon vs Noise SNR ## Beacon vs Noise SNR
## ##
if True: if True:
beacon_snrs = [ lib.signal_to_noise(myfilter(beacon_amp*ant.beacon), myfilter(ant.noise), samplerate=1/dt, signal_band=beacon_pb) for ant in antennas ] beacon_snrs = [ lib.signal_to_noise(myfilter(beacon_amp*ant.beacon), myfilter(ant.noise), samplerate=1/dt, signal_band=beacon_pb, noise_band=pb) for ant in antennas ]
fig, ax = plt.subplots() fig, ax = plt.subplots()
ax.set_title("Maximum Beacon/Noise SNR") ax.set_title("Maximum Beacon/Noise SNR")
ax.set_xlabel("Antenna no.") ax.set_xlabel("Antenna no.")
ax.set_ylabel("SNR") ax.set_ylabel("SNR")
ax.plot([ int(ant.name) for ant in antennas], beacon_snrs, 'o', ls='none') ax.plot([ int(ant.name) for ant in antennas], beacon_snrs, 'o', ls='none')
if fig_dir: if fig_dir:
fig.savefig(path.join(fig_dir, path.basename(__file__) + f".beacon_vs_noise_snr.pdf")) fig.savefig(path.join(fig_dir, path.basename(__file__) + f".beacon_vs_noise_snr.pdf"))
@ -85,7 +85,7 @@ if __name__ == "__main__":
## Beacon vs Total SNR ## Beacon vs Total SNR
## ##
if True: if True:
beacon_snrs = [ lib.signal_to_noise(myfilter(beacon_amp*ant.beacon), ant.E_AxB, samplerate=1/dt, signal_band=beacon_pb) for ant in antennas ] beacon_snrs = [ lib.signal_to_noise(myfilter(beacon_amp*ant.beacon), ant.E_AxB, samplerate=1/dt, signal_band=beacon_pb, noise_band=pb) for ant in antennas ]
fig, ax = plt.subplots() fig, ax = plt.subplots()
ax.set_title("Maximum Beacon/Total SNR") ax.set_title("Maximum Beacon/Total SNR")
@ -99,16 +99,16 @@ if __name__ == "__main__":
## Airshower signal vs Noise SNR ## Airshower signal vs Noise SNR
## ##
if True: if True:
shower_snrs = [ lib.signal_to_noise(ant.E_AxB, myfilter(ant.noise), samplerate=1/dt, signal_band=pb) for ant in antennas ] shower_snrs = [ lib.signal_to_noise(ant.E_AxB, myfilter(ant.noise), samplerate=1/dt, signal_band=pb, noise_band=pb) for ant in antennas ]
fig, ax = plt.subplots() fig, ax = plt.subplots()
ax.set_title("Total (Signal+Beacon+Noise)/Noise SNR") ax.set_title("Total (Signal+Beacon+Noise)/Noise SNR")
ax.set_xlabel("Antenna no.") ax.set_xlabel("Antenna no.")
ax.set_ylabel("SNR") ax.set_ylabel("SNR")
ax.plot([ int(ant.name) for ant in antennas], shower_snrs, 'o', ls='none') ax.plot([ int(ant.name) for ant in antennas], shower_snrs, 'o', ls='none')
if fig_dir: if fig_dir:
fig.savefig(path.join(fig_dir, path.basename(__file__) + f".total_snr.pdf")) fig.savefig(path.join(fig_dir, path.basename(__file__) + f".total_snr.pdf"))
if show_plots: if show_plots:
plt.show() plt.show()