Skip to content

add spectrogram during data generation

Avi Vajpeyi requested to merge avi.vajpeyi/bilby_pipe:add_gwpy_spectogram into master
NOTE: Requires !147 (merged) to be completed.

Summary

Add GWpy spectograms during data_generation if create_plots=True.

  • Creates plots for H1 L1 PSDs and data segment.
  • Labels files based on the duration of data plotted.
  • Allows GWpy to determine xlim and ylim.
  • If injection inserted, plots made before and after injection.

Examples of plots

outdir_timeshift_GW150914_data_H1_timeshift_GW150914_data0_1126259462-4_generation_D4_data

outdir_timeshift_GW150914_data_L1_timeshift_GW150914_data0_1126259462-4_generation_D4_data

outdir_timeshift_GW150914_data_L1_timeshift_GW150914_data0_1126259462-4_generation_D128_data

outdir_timeshift_GW150914_data_H1_timeshift_GW150914_data0_1126259462-4_generation_D128_data

Current issues with plots

Strain data needs to be bandpassed

Not sure of bandpassing that can generally apply to all strain data.

q-transform plot needs to have lower y-lim

Currently cant reduce the freq ylim for the 4-second duration. If ylim manually reduced, the plot is just blank below the default ylim. The freq ylim can be reduced if more data (start time - 10, end time +10) is acquired. Then during plotting the plot will not be blank if the freq ylim is reduced.

However, this will require an extra gwpy DataFetch call (which is slow).

Edited by Avi Vajpeyi

Merge request reports