SEXEC=singularity exec -B $TMPDIR gstlal-dev/ #CIT example


rankstatpdf=/home/gstlalcbc/observing/4/dev/runs/trigs.edward_O3_replay_MDC05/dist_stat_pdfs/H1L1V1-GSTLAL_DIST_STAT_PDFS-0-0.xml.gz
tmpweights=/home/anarya.ray/gstlal-pastro/manifold_template_weights
V='{"BBH": 3.720777380922114, "NSBH": 0.13036984478340222, "BNS": 0.01447408360378373}' # O3 sentivity volume
farth=730.5 # FAR threshold of 2 per day in one julian year
rates='{"BBH": 18.0730098234, "NSBH": 52.9087878433, "BNS": 662.448825556}' # O3 merge rates
ratesinj='{"BBH": 4045.100324417181, "BNS": 1086074.923130263, "NSBH": 269801.54621791927}' # MDC injection rates
injdb=/home/anarya.ray/gstlal-pastro/mario_injections/injections-minSNR-4.xml.gz

all: pastro_inj.txt

pastro_model_fgmc.h5:
	$(SEXEC) sgnl-pastro-initialize-fgmc-model --rankstat-filename ${rankstatpdf} --weights-dir ${tmpweights} --V-new ${V} --far-threshold  ${farth} --output pastro_model_fgmc.h5 --rates ${rates} --rates-inj ${ratesinj}

pastro_inj.txt: pastro_fgmc_model.h5
	$(SEXEC) -B $TMPDIR gstlal-dev sgnl-pastro-calculate-fgmc-test-inj --output-file pastro_inj.txt --fgmc-model pastro_model_fgmc.h5 --injection-file ${injdb}

clean:
	rm -rf fgmc_test/*.h5 && rm -rf fgmc_test/*.txt
