Just some quick directions on how to run a script to get the full 524 kHz DCPD data, and then save it in a useful file format.

The script to save this is in /ligo/home/elenna.capote/OMC_DCPD/H1OMCDPCD_HF.py

Right now it is set to collect data from these two channels (line 62 of the code):

channels = ["H1:OMC-DCPD_A0_OUT", "H1:OMC-DCPD_B0_OUT"]

The amount of time to collect data is set on line 66:

duration = 600  # seconds

To run the code:

> conda activate labutils

> python H1OMCDPCD_HF.py

The files will be saved in /ligo/home/elenna.capote/OMC_DCPD in a new folder labeled with a datetime and "xcorr_data" since I have traditionally used this script to collect data for cross correlation measurements. However, you can change the label on line 71, label="xcorr_data"

The output is a set of 1 second HDF5 frames for each second of data collection (so 600 seconds = 600 files). To put this data into a single .gwf format of the full data length, which can easily be read in with GWPY, use the combine.py script. Copy combine.py into the new directory with the frames and run it. The code generates a file for each OMC DCPD channel at the full rate, and at a decimated 64 kHz rate.

This is a bit odd the way it is set up, but the way this works ensures that even if you have to quit the script early, you don't lose all the data, since it saves each file by the second.