all: bayestar.out

inj.xml:
	lalapps_inspinj -o inj.xml \
	--m-distr fixMasses --fixed-mass1 1.4 --fixed-mass2 1.4 \
	--t-distr uniform --time-step 7200 \
	--gps-start-time 1000000000 --gps-end-time 1000086400 \
	--d-distr volume --min-distance 1 --max-distance 600e3 \
	--l-distr random --i-distr uniform \
	--f-lower 30 --waveform TaylorF2threePointFivePN --disable-spin

psd.xml:
	bayestar_sample_model_psd -o psd.xml \
	--H1=aLIGOZeroDetHighPower \
	--L1=aLIGOZeroDetHighPower \
	--I1=aLIGOZeroDetHighPower \
	--V1=AdvVirgo \
	--K1=KAGRA

coinc.xml: psd.xml inj.xml
	bayestar_realize_coincs --reference-psd psd.xml inj.xml -o coinc.xml \
	--detector H1 L1 V1 I1 K1 --keep-subthreshold

coinc_inj_found.sqlite: coinc.xml inj.xml
	ligolw_add coinc.xml inj.xml | \
	lalapps_inspinjfind | \
	ligolw_sqlite  --preserve-ids --replace \
	--database coinc_inj_found.sqlite /dev/stdin

0.toa_phoa_snr.fits: coinc.xml psd.xml
	bayestar_localize_coincs coinc.xml

bayestar.out: coinc_inj_found.sqlite 0.toa_phoa_snr.fits
	bayestar_aggregate_found_injections \
	coinc_inj_found.sqlite \*.fits -o bayestar.out \
	--contour 50 90

clean:
	rm -f inj.xml psd.xml coinc.xml coinc_inj_found.sqlite *.fits bayestar.out