Cannot run 4 detectors test analysis because some program do not support K1
Specifically, https://git.ligo.org/lscsoft/gstlal/-/blob/master/gstlal-inspiral/bin/gstlal_inspiral_injection_snr and https://git.ligo.org/lscsoft/gstlal/-/blob/master/gstlal-inspiral/python/stats/inspiral_lr.py#L555 uses hard-coded information for doing calculation and they cause some jobs fails for 4 detectors (H1, L1, V1, K1) analysis with injection jobs.
Tracing back the failed jobs, it appears that gstlal_inspiral_injection_snr computes the injection SNRs using PSDs and hard coding them into alpha4
<-- H1, alpha5
<-- L1, and alpha6
<-- V1 of the SimInspiralTable, and then these injection SNRs are being used in this line for simulating the numerator of the LR. However, I believe alphaX
is not supposed to be used for injection SNRs, so it suggests that it's not just hard coding problem but a design problem.
To solve this problem properly, kipp suggested two solutions:
- don't pre-compute these numbers in SimInprialTable, move this calculation to where it is needed and compute the numbers on-the-fly
- correct the table structure (SimInprialTable) so that it can accommodate the desired instruments (preferred)
It would be great for people to know this problem and discuss about potential solutions.
Related information: https://git.ligo.org/waveforms/o4-injection-file/-/tree/main