m-thesis-introduction/simulations/airshower_beacon_simulation
2022-12-22 16:51:39 +01:00
..
figures ZH: save figures to dir 2022-12-05 17:49:47 +01:00
lib ZH: upgrade lib test scripts 2022-12-19 19:02:35 +01:00
ZH_airshower Got ZHaires simulated signal 2022-09-22 09:36:37 +02:00
aa_generate_beacon.py ZH: move c_light value to lib 2022-12-15 14:40:20 +01:00
ab_modify_clocks.py ZH: set ab_modify_clock back to default situation 2022-12-15 15:10:42 +01:00
ba_measure_beacon_phase.py ZH: show phase due to t_trace[0] from measurement 2022-12-12 20:32:04 +01:00
bb_measure_true_phase.py ZH: plot showing (measured-actual) residuals for true_phases 2022-12-20 15:26:10 +01:00
bc_baseline_phase_deltas.py ZH: small tweak to figure for baseline phases 2022-12-22 16:51:39 +01:00
bd_antenna_phase_deltas.py ZH: show sigma_phase matrix as a plot 2022-12-22 16:33:19 +01:00
ca_beacon_periods.py ZH: rename scripts 2022-12-02 16:37:57 +01:00
ca_period_from_shower.py ZH: period_from_shower sort by maxima 2022-12-08 15:29:05 +01:00
Makefile ZH: rename b* scripts 2022-12-08 15:05:24 +01:00
README.md ZH: rename b* scripts 2022-12-08 15:05:24 +01:00
reconstruction.py ZH: rename scripts 2022-12-02 16:37:57 +01:00
show_beacon_amplitude_antennas.py ZH: improve some plots 2022-11-22 15:37:55 +01:00
view_beaconed_antenna.py ZH: improve some plots 2022-11-22 15:37:55 +01:00
view_orig_ant0.py ZH: tweaks to view scripts 2022-11-18 19:36:48 +01:00

Airshower + Beacon simulation

Simulate receiving an airshower and beacon in an array. Using the beacon, the clocks can be calibrated.

The ZHaires simulated airshower is stored in ./ZH_airshower. The produced files can be read using ./earsim.

All quantities are stored in HDF5 files with the airshower data.

Steps:

  1. Setup

    1. Beacon ([./aa_generate_beacon.py])

      1. Define tx position
      2. Read in antennas
      3. Sample beacon at each antenna
      4. Add to relevant polarisation
      5. Save antenna traces
    2. Timeoffset ([./ab_modify_clocks.py])

      1. Generate timeoffsets for each antenna
      2. Modify time samples
  2. Beacon analysis

    1. Find beacon frequency and phase in antenna traces ([./ba_measure_beacon_phase.py])
    2. Remove phase due to distances
    3. Analyse the sigma between antenna pairs
    4. Remove 0j part and take the mean of (:,j)
    5. Assign mean as clock shift phase of antenna
  3. Find k*2\pi phase offsets (periods) ([./ca_periods_from_showers.py])

  4. Rewrite clocks