From 040f27bd7891467a75f48575ae74af493844d515 Mon Sep 17 00:00:00 2001
From: Eric Teunis de Boone <ericteunis@deboone.nl>
Date: Thu, 2 Feb 2023 14:27:25 +0100
Subject: [PATCH] ZH: Makefile with config.mk to control properties of the
 beacon

---
 .../airshower_beacon_simulation/.gitignore    |  1 +
 .../airshower_beacon_simulation/Makefile      | 32 +++++++++++++++----
 2 files changed, 26 insertions(+), 7 deletions(-)
 create mode 100644 simulations/airshower_beacon_simulation/.gitignore

diff --git a/simulations/airshower_beacon_simulation/.gitignore b/simulations/airshower_beacon_simulation/.gitignore
new file mode 100644
index 0000000..aee2e4c
--- /dev/null
+++ b/simulations/airshower_beacon_simulation/.gitignore
@@ -0,0 +1 @@
+config.mk
diff --git a/simulations/airshower_beacon_simulation/Makefile b/simulations/airshower_beacon_simulation/Makefile
index d7e3fb1..fa18822 100644
--- a/simulations/airshower_beacon_simulation/Makefile
+++ b/simulations/airshower_beacon_simulation/Makefile
@@ -3,17 +3,32 @@
 
 FIG_DIR := ./figures
 
-all: beacon clocks phases findks reconstruct
+-include config.mk
 
+MAX_K ?= 3
+CLK_DEV ?= 20
+TRACE_N ?= 4096
+TRACE_PRE ?= 500
+NOISE_SIGMA ?= 1e4
+
+BEAC_AMP ?= 1e3 0 0
+BEAC_F ?= 51.53e-3
+BEAC_DECAY ?= 0
+
+PB_LOW ?= 0.03
+PB_HIGH ?= 0.08
+
+
+all: beacon clocks phases findks vary-fixes reconstruct
 
 beacon:
-	./aa_generate_beacon.py | tee figures/aa.log
-	./ab_modify_clocks.py 0 | tee figures/ab.log
+	./aa_generate_beacon.py -N ${TRACE_N} -P ${TRACE_PRE} -n ${NOISE_SIGMA} -a ${BEAC_AMP} -f ${BEAC_F} -l ${PB_LOW} -u ${PB_HIGH} -d ${BEAC_DECAY} | tee ${FIG_DIR}/aa.log
+	./ab_modify_clocks.py 0 | tee ${FIG_DIR}/ab.log
 	./ac_show_signal_to_noise.py --no-show-plots --fig-dir=${FIG_DIR}
-	./view_beaconed_antenna.py 72 -p x -p y -p z --no-show-plots --fig-dir=${FIG_DIR}
+	./view_beaconed_antenna.py 72 -p x -p y -p z -p n -p b --ft --no-show-plots --fig-dir=${FIG_DIR}
 
 clocks:
-	./ab_modify_clocks.py 15 --gaussian
+	./ab_modify_clocks.py ${CLK_DEV} --gaussian
 
 phases:
 	./ba_measure_beacon_phase.py --no-show-plots --fig-dir=${FIG_DIR}
@@ -22,13 +37,16 @@ phases:
 	./bd_antenna_phase_deltas.py --no-show-plots --fig-dir=${FIG_DIR}
 
 findks:
-	./ca_period_from_shower.py --no-show-plots --fig-dir=${FIG_DIR}
+	./ca_period_from_shower.py --max-k ${MAX_K} --no-show-plots --fig-dir=${FIG_DIR}
 	./cb_report_measured_antenna_time_offsets.py --no-show-plots --fig-dir=${FIG_DIR}
 
+vary-fixes:
+	./dc_grid_power_time_fixes.py --no-show-plots --fig-dir=${FIG_DIR}
+
 reconstruct:
 	./da_reconstruction.py --no-show-plots --fig-dir=${FIG_DIR}
 	./db_longitudinal_figure.py --no-show-plots --fig-dir=${FIG_DIR}
-	./dc_grid_power_time_fixes.py --no-show-plots --fig-dir=${FIG_DIR}
+
 dist-clean:
 	rm -f ZH_airshower/antennas.hdf5
 	rm -f ZH_airshower/clocks.csv