NUM_INJECTIONS ?= 100
SEED ?= 42

# Mass model parameters (override on the command line, e.g. `make MODEL=salpeter`)
MODEL ?= power-law
INDEX ?= -2.0
SNR ?= 10

.PHONY: all clean mass-model-plots

all : mc-test_mismatch_histogram.png mass-model-plots

mc-only-bank.h5 : mc-only-bank.yaml O4_projected_psds.xml.gz
	sgn-manifold-cbc-bank --yaml $< --output-h5 $@ --verbose

mc-test.h5 : mc-only-bank.h5
	sgn-manifold-cbc-bank-test \
		--bank $< \
		--output-h5 $@ \
		--num-injections $(NUM_INJECTIONS) \
		--seed $(SEED) \
		--verbose

mc-test_mismatch_histogram.png : mc-test.h5
	sgn-manifold-cbc-bank-test-plot \
		--output-prefix mc-test \
		--test $< \
		--verbose

mc-only-mass-model.h5 : mc-only-bank.h5
	sgn-manifold-cbc-bank-mass-model \
		--bank $< \
		--output-h5 $@ \
		--model $(MODEL) \
		--index $(INDEX) \
		--verbose

# Produces mc-only-mass-model_snr-$(SNR).0_weight_mchirp.png.  For a 1-D
# log10mc bank only the chirp-mass marginal is non-degenerate, so that is the
# single plot emitted; richer banks additionally get m1-m2 and mc-chi views.
mass-model-plots : mc-only-mass-model.h5
	sgn-manifold-cbc-bank-plot-mass-model \
		--snr $(SNR) \
		--title mc-only-mass-model \
		$<

clean:
	rm -f mc-only-bank.h5 mc-test.h5 mc-test_*.png \
		mc-only-mass-model.h5 mc-only-mass-model*.png
