Gaussian MDC: All IFOs were using same random seed for frame generation
Hanford, Livingston, and Virgo data were all using same randomseed in pycbc_condition_strain
script. Thus the noise was identical in the detectors. This was resulting in large number of triggers being generated due to correlated noise. This is fixed by changing the randomseed variable for each detector:
Diff of the wrapper script that generates the frames:
shaon@alpha:~/projects/analysis/mock-data-challenge$ git diff tools_to_inject/Gaussian_MDC/wrapper_script.py
diff --git a/tools_to_inject/Gaussian_MDC/wrapper_script.py b/tools_to_inject/Gaussian_MDC/wrapper_script.py
index 69d24da..8827133 100644
--- a/tools_to_inject/Gaussian_MDC/wrapper_script.py
+++ b/tools_to_inject/Gaussian_MDC/wrapper_script.py
@@ -35,6 +35,8 @@ parser.add_argument("-i", "--index", type=str, action='store',
help="Key of the JSON file")
parser.add_argument("-o", "--outdir", type=str, action='store', default="frames",
help="Name of the output directory")
+parser.add_argument("-s", "--seed", type=int, action="store", default=100,
+ help="Random seed for pycbc_condition_strain")
args = parser.parse_args()
@@ -60,9 +62,9 @@ for t in tstarts:
thisiter0 = time.time()
cmd = "pycbc_condition_strain --injection-file {} --channel-name {}1:MDC-STRAIN \
--output-strain-file {}_STRAIN-{}-{}.gwf --fake-strain-from-file {} \
- --fake-strain-seed 100 --low-frequency-cutoff 11 --fake-strain-flow 11 \
+ --fake-strain-seed {} --low-frequency-cutoff 11 --fake-strain-flow 11 \
--frame-duration {} --gps-start-time {} --gps-end-time {} \
- --sample-rate 16384".format(args.injection_file, args.det, det, x, y, args.psdfile, seg_len, t, t + seg_len)
+ --sample-rate 16384".format(args.injection_file, args.det, det, x, y, args.psdfile, args.seed, seg_len, t, t + seg_len)
print('Creating frame file for segment {} ---> {}'.format(t, t + seg_len))
os.system(cmd)
The new function call looks like:
python wrapper_script_git.py -d 4096 -f $timefile -i $i -I $injfile -p $psdfile -s $seed -D $detector -o $output