Burst style MDCs
Instructions for running burst-style MDCs in which a signal is injected from a .gwf frame file.
Making an MDC cache
BayesWave requires a path to the cache file of the MDC frames.
If you have specific MDC frames you want to run on, you can make a cache file via
$ lal_path2cache [directory_with_frames]/*gwf > MDC.cache
Injecting into real noise
In this example we will inject a white noise burst into real H1 and L1 data from O2. Here we already have an MDC cache made on CIT.
Todo
Put parameters of WNB here
Config file
[input]
dataseed=1234
seglen=4.0
window=1.0
flow=32
srate=2048
ifo-list=['H1','L1']
[engine]
install=/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py39
bayeswave=%(install)s/bin/BayesWave
bayeswave_post=%(install)s/bin/BayesWavePost
megaplot=%(install)s/bin/megaplot.py
[datafind]
channel-list={'H1':'H1:DCH-CLEAN_STRAIN_C02','L1':'L1:DCH-CLEAN_STRAIN_C02'}
frtype-list={'H1':'H1_CLEANED_HOFT_C02','L1':'L1_CLEANED_HOFT_C02'}
url-type=file
veto-categories=[1]
[injections]
mdc-cache=/home/meg.millhouse/Development_BW/bayeswave-docs-and-examples/MDC/MDC_caches/MDC.cache
mdc-channels={'H1':'H1:SCIENCE','L1':'L1:SCIENCE'}
mdc-prefactor=1
[bayeswave_options]
; command line options for BayesWave. See BayesWave --help
bayesLine=
updateGeocenterPSD=
waveletPrior=
Dmax=100
Niter=5000
[bayeswave_post_options]
; command line options for BayesWavePost. See BayesWavePost --help
bayesLine=
0noise=
[condor]
; see e.g., https://ldas-gridmon.ligo.caltech.edu/ldg_accounting/user
universe=vanilla
checkpoint=
bayeswave-request-memory=4000MB
bayeswave_post-request-memory=4000MB
bayeswave-request-disk=100MB
bayeswave_post-request-disk=100MB
datafind=/usr/bin/gw_data_find
ligolw_print=/usr/bin/ligolw_print
segfind=/usr/bin/ligolw_segment_query_dqsegdb
accounting-group = ligo.dev.o4.burst.paramest.bayeswave
Call to bayeswave_pipe
#!/bin/bash
workdir=WNB_MDC_simNoise
trigtime=1174112418.84
bayeswave_pipe \
--workdir ${workdir} \
-t ${trigtime} \
--sim-data \
MDC_simNoise.ini
Injecting into simulated data
In this example we will inject a white noise burst simulated H1 and L1 data.
Config file
[input]
dataseed=1234
seglen=4.0
window=1.0
flow=32
srate=2048
ifo-list=['H1','L1']
[engine]
install=/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py39
bayeswave=%(install)s/bin/BayesWave
bayeswave_post=%(install)s/bin/BayesWavePost
megaplot=%(install)s/bin/megaplot.py
[datafind]
channel-list={'H1':'H1:LALSimAdLIGO','L1':'L1:LALSimAdLIGO'}
frtype-list={'H1':'LALSimAdLIGO','L1':'LALSimAdLIGO'}
url-type=file
veto-categories=[1]
[injections]
mdc-cache=/home/meg.millhouse/Development_BW/bayeswave-docs-and-examples/MDC/MDC_caches/MDC.cache
mdc-channels={'H1':'H1:SCIENCE','L1':'L1:SCIENCE'}
mdc-prefactor=1
[bayeswave_options]
; command line options for BayesWave. See BayesWave --help
updateGeocenterPSD=
waveletPrior=
Dmax=100
; Note this is a short number of iterations for testing purposes only
Niter=5000
; To use multithreading, un-comment the lines below
; Nchain=20
; Nthreads=4
[bayeswave_post_options]
; command line options for BayesWavePost. See BayesWavePost --help
0noise=
[condor]
; see e.g., https://ldas-gridmon.ligo.caltech.edu/ldg_accounting/user
universe=vanilla
checkpoint=
bayeswave-request-memory=4000MB
bayeswave_post-request-memory=4000MB
bayeswave-request-disk=100MB
bayeswave_post-request-disk=100MB
datafind=/usr/bin/gw_data_find
ligolw_print=/usr/bin/ligolw_print
segfind=/usr/bin/ligolw_segment_query_dqsegdb
accounting-group = ligo.dev.o4.burst.paramest.bayeswave
Call to bayeswave_pipe
#!/bin/bash
workdir=WNB_MDC_simNoise
trigtime=1174112418.84
bayeswave_pipe \
--workdir ${workdir} \
-t ${trigtime} \
--sim-data \
MDC_simNoise.ini